使用for循环模拟t检验p值

时间:2010-05-19 11:10:16

标签: r loops statistics

对于这个项目,我需要使用R脚本来模拟t检验的有效性。我必须使用for循环将用于执行以下2000次:

循环看起来像这样吗

i <- 1
for (i <= 2001) { 
    x <-rf(5,df1=5,df2=10)
    b <- df2
    p.value <-t.test(x,mu=(b/(b-2))$p.value
    i <- i+1
}

1 个答案:

答案 0 :(得分:6)

在你编写它的方式中,它将是一个“while”循环。

R中的For循环具有以下语法:

for (i in 1:2000) {
    df1 <- 5
    df2 <- 10
    x <-rf(5, df1=df1, df2=df2)
    b <- df2
    p.value <- t.test(x, mu=(b/(b-2)))$p.value
}

此外,使用“apply”构造可能更有效,例如使用replicate,并将df包含为函数参数:

get.p.value <- function(df1, df2) {
        x <- rf(5, df1=df1, df2=df2)
        p.value <- t.test(x, mu=(df2/(df2-2)))$p.value
    }
replicate (2000, get.p.value(df1 = 5, df2 = 10))

这并非总是如此,但它简化了p.values的恢复。