Question

对于这个项目，我需要使用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
}

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
}

此外，使用“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的恢复。

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

1 个答案: