在 R 中绘制幂函数(假设检验)

Plotting power function in R (hypothesis testing)

提问人:John hall 提问时间:9/15/2023 最后编辑:SalvatoreJohn hall 更新时间:9/21/2023 访问量:48

问:

我想在 R 中绘制一个幂函数。符号检验通用的幂函数的定义是:如果 pvalue > 0.05,则为 alpha,否则为 1。

我拥有的代码是:

# Power function do sign test
condition <- ifelse(results_df$sign_test_p_value>0.05,0.05,1)
plot(condition,type="l", main="Power Function",lwd=2,xlab= "Mean",ylab="Power of the test", xlim=c(9,15),ylim=c(0,1))
lines(condition_1, col='black',type="l", lty=2)

当我绘图时,x 轴从 12 开始,一直到 16。这里的问题是 x 轴给出了真实均值的值。平均样本值约为 -0.00004。这里的问题是 x 轴上的值不正确。我不知道我的代码是否有问题/我做情节的方式。

如果有人想检查代码,符号测试的 p 值为:

results_df <- data.frame(sign_test_p_value = 
                          c(1.000000e+00, 1.000000e+00, 1.000000e+00, 
                            1.000000e+00, 1.000000e+00, 1.000000e+00, 
                            1.000000e+00, 1.000000e+00, 1.000000e+00, 
                            9.964203e-01, 5.628114e-01, 1.046356e-02, 
                            2.184458e-06, 6.853518e-12, 1.110223e-16, 
                            1.110223e-16, 1.110223e-16, 1.110223e-16, 
                            1.110223e-16, 1.110223e-16, 1.110223e-16))
r 假设检验

评论

0赞 pookpash 9/15/2023
目前,您只是将条件中的所有变量一个接一个地绘制,因此 x 轴是该值的索引,y 绘制 的值。如果按照代码中的 xlim 设置为 9 到 15,则该图将仅显示第 9 到 15 个值。另外,由于您没有指定它,因此我们不知道是什么。但是我真的不明白你想达到什么目的,所以请详细说明。我猜你想在 x 轴上有平均值,但不确定。如果是这样的话,我们还需要平均值,而不仅仅是 p 值。conditioncondition_1

答: 暂无答案