Simuler test t valeurs p en utilisant une boucle

Question

Pour ce projet, je suis obligé d'utiliser un script R pour simuler l'efficacité du test t. Je dois utiliser une boucle sera utilisée pour effectuer les 2000 heures suivantes:

Est-ce que l'apparence de la boucle quelque chose comme ça

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

Dans la façon dont vous l'avez écrit, ce serait une boucle « while » .

Pour boucles R ont la syntaxe suivante:

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 
} 

De plus, il pourrait être plus efficace d'employer un construire « appliquer », par exemple avec répliquée et comprennent les df comme arguments de la fonction:

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)) 

Ce n'est pas toujours vrai, mais il simplifie la récupération des p.values.