আমাকে পরিসংখ্যানগুলির একটি ভূমিকা লিখতে বলা হয়েছে এবং পি-ভ্যালু এবং পাওয়ার সম্পর্কিত হওয়ার উপায়টি কীভাবে গ্রাফিকালি দেখাতে হবে তা নিয়ে আমি লড়াই করছি। আমি এই গ্রাফটি নিয়ে এসেছি:
আমার প্রশ্ন: এটি প্রদর্শনের আরও ভাল উপায় আছে কি?
এই আমার আর কোড
x <- seq(-4, 4, length=1000)
hx <- dnorm(x, mean=0, sd=1)
plot(x, hx, type="n", xlim=c(-4, 8), ylim=c(0, 0.5),
ylab = "",
xlab = "",
main= expression(paste("Type II (", beta, ") error")), axes=FALSE)
axis(1, at = c(-qnorm(.025), 0, -4),
labels = expression("p-value", 0, -infinity ))
shift = qnorm(1-0.025, mean=0, sd=1)*1.7
xfit2 <- x + shift
yfit2 <- dnorm(xfit2, mean=shift, sd=1)
# Print null hypothesis area
col_null = "#DDDDDD"
polygon(c(min(x), x,max(x)), c(0,hx,0), col=col_null)
lines(x, hx, lwd=2)
# The alternative hypothesis area
## The red - underpowered area
lb <- min(xfit2)
ub <- round(qnorm(.975),2)
col1 = "#CC2222"
i <- xfit2 >= lb & xfit2 <= ub
polygon(c(lb,xfit2[i],ub), c(0,yfit2[i],0), col=col1)
## The green area where the power is
col2 = "#22CC22"
i <- xfit2 >= ub
polygon(c(ub,xfit2[i],max(xfit2)), c(0,yfit2[i],0), col=col2)
# Outline the alternative hypothesis
lines(xfit2, yfit2, lwd=2)
axis(1, at = (c(ub, max(xfit2))), labels=c("", expression(infinity)),
col=col2, lwd=1, lwd.tick=FALSE)
legend("topright", inset=.05, title="Color",
c("Null hypoteses","Type II error", "True"), fill=c(col_null, col1, col2), horiz=FALSE)
abline(v=ub, lwd=2, col="#000088", lty="dashed")
arrows(ub, 0.45, ub+1, 0.45, lwd=3, col="#008800")
arrows(ub, 0.45, ub-1, 0.45, lwd=3, col="#880000")
হালনাগাদ
ভয়ঙ্কর উত্তরের জন্য আপনাকে ধন্যবাদ। আমি কিছু কোড পরিবর্তন করেছি:
# Print null hypothesis area
col_null = "#AAAAAA"
polygon(c(min(x), x,max(x)), c(0,hx,0), col=col_null, lwd=2, density=c(10, 40), angle=-45, border=0)
lines(x, hx, lwd=2, lty="dashed", col=col_null)
...
legend("topright", inset=.015, title="Color",
c("Null hypoteses","Type II error", "True"), fill=c(col_null, col1, col2),
angle=-45,
density=c(20, 1000, 1000), horiz=FALSE)
আমি নাল অনুমানের ড্যাশড, কিছুটা অস্পষ্ট ছবি পছন্দ করি কারণ এটি ইঙ্গিত দেয় যে এটি সত্যিকার অর্থে সেখানে নেই। আমি স্বচ্ছতা এবং আলফা যুক্ত করার বিষয়ে চিন্তা করেছি তবে আমি একটি ছবিতে খুব বেশি তথ্য পাওয়ার বিষয়ে চিন্তিত এবং তাই এটি না করার সিদ্ধান্ত নিয়েছি।
মুদ্রিত নিবন্ধগুলির সীমাবদ্ধতা আমাকে পাঠকদের পরীক্ষা করতে দেয় না। আমি টিচারিংডেমোসের সাথে @ গ্রেগ স্নোয়ের জবাবটি আমার উত্তর হিসাবে বেছে নিয়েছি যেহেতু আমি দুটি ত্রুটিটি ওভারল্যাপিং না করে ধারণাটি পছন্দ করি।