Bonjour à tous,
j'ai un problème: je dois simuler n variables aléatoires de loi N(m,s²) où s=0.5 sur le logiciel R pour tester H0: m=m0 contre m>m0, renouveller l'expérience 500 fois, calculer Tobs et noter succès si on rejette H0!!
J'ai fait un petit programme mais je ne sais pas si c'est bon:
>for (i in 1:500) {
x>-rnorm(n,5,0.5)
}
>t.test(x, alternative='greater')
Et je voulais savoir aussi comment on récupère la p-value.
Après il faut que je calcule la puissance du test. J'ai fait:
> power.t.test( n, delta=5-mean(x), sd=0.5, sig.level=0.05)
Et dernière question, il faut que je refasse pareil dans le cas non gaussien (avec la loi exponentielle par exemple) j'ai fait la même début j'ai mi rexp(100,8) au lieu de rnorm(10,5,0.5) mais après je ne sais pas quoi faire comme test (enfin je ne connais pas le nom)
Si quelqu'un pouvait m'aider ça serait vraiment gentil car je suis totalement paumée!
Merci d'avance
-----