使用for循环模拟t-test p值

Question

对于这个项目，我需要使用一个R脚本来模拟t-test的有效性。我必须使用一个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
}

Answer 1

按照您编写的方式，它将是一个"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
}

此外，使用“应用”结构可能更有效，例如使用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的恢复。