在R中模拟t检验的显著性水平
下面的代码在我创建的自定义发行版上生成10000个t测试在R中模拟t检验的显著性水平,r,statistics,simulation,R,Statistics,Simulation,下面的代码在我创建的自定义发行版上生成10000个t测试 > x <- replicate(10000,{ t.test(rcn(20,.25,25), mu=0, alternative="greater") }) >x您可以直接访问replicate中每个t.test的测试统计信息。例如: x <- replicate(10000, { t.test(rcn(20,.25,25), mu=0, alternative="greater")$statistic
> x <- replicate(10000,{
t.test(rcn(20,.25,25), mu=0, alternative="greater")
})
>x您可以直接访问replicate
中每个t.test
的测试统计信息。例如:
x <- replicate(10000, {
t.test(rcn(20,.25,25), mu=0, alternative="greater")$statistic
})
请给出rcn
@Beasterfield的定义,使您的示例具有可复制性。我的问题代表任何分布。如何在数千个测试中选择具有相同特征的t测试?谢谢。我明白了,但是你可以从stats
中选择一个分布函数,让人们有机会复制你的代码。请看,很好!这就是美元符号的作用吗?它细化了函数?@JohnK$
用于访问元素。在本例中,t.test
返回一个列表。$
符号提取名为statistic
的元素。
crit <- 1.729
sum(x > crit)