Question

> prop.test(6,10)

    1-sample proportions test with continuity correction

data:  6 out of 10, null probability 0.5
X-squared = 0.1, df = 1, p-value = 0.7518
alternative hypothesis: true p is not equal to 0.5
95 percent confidence interval:
    0.2736697 0.8630694
sample estimates:
    p 
0.6

请帮帮我，如何手动计算？

这是我尝试过的，

> se =  sqrt(0.6*0.4*10/9)/sqrt(10)
> mean = 6/10
> z =  qnorm(0.975)
> from95 = mean - z*se
> to95 = mean + z*se
> from95
[1] 0.2799392
> to95
[1] 0.9200608

Answer 1

您需要的所有信息都在函数的源代码中。您可以使用as.list(body(prop.test))将其分成多个部分。或者将其发送到sink的文件，以便您可以单独研究它也可能有所帮助。

例如，您可以使用

查看函数体的第31个元素

> as.list(body(prop.test))[[31]]
# if (alternative == "two.sided") PVAL <- pchisq(STATISTIC, PARAMETER, 
#     lower.tail = FALSE) else {
#     if (k == 1) 
#         z <- sign(ESTIMATE - p) * sqrt(STATISTIC)
#     else z <- sign(DELTA) * sqrt(STATISTIC)
#     PVAL <- pnorm(z, lower.tail = (alternative == "less"))
# }

prop.test（x，n）如何在内部工作？

1 个答案: