我试图在R中创建t临界值表.我使用了以下代码,但是有些错误,因为我无法获得正确的值:
t.power <- function(n) {
alpha=c(0.50,0.75,0.80,0.85,0.90,0.95,
0.975,0.99,0.995,0.999,0.9995,1.0)
power=qt(alpha, df=1:n)
return(power)
}
Run Code Online (Sandbox Code Playgroud)
然后我试着测试它 n=15
t.power(15)
Run Code Online (Sandbox Code Playgroud)
我得到了
[1] 0.0000000 0.8164966 0.9784723 1.1895669 1.4758840 1.9431803 2.3646243
[8] 2.8964594 3.2498355 4.1437005 4.4369793 Inf 0.0000000 0.6924171
[15] 0.8662450
Run Code Online (Sandbox Code Playgroud)
谁能告诉我我做错了什么?