Création de fonction

Bonjour,
voila mon problème, je dois réaliser trois graphiques qui montrent l'évolution des trois variables que ma fonction doit créer pour les n individus allant de 10 à 50000 de 4 à 4.
Ma fonction fonctionne, mais elle me crée que une valeur pour chaque variable, alors qu'il faudrait plus de 1200.
Et ensuite les représenter dans un graphique.

Voici mon code:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

stat <- function(n){ grillex = rnorm(n, mean = 175, sd = 10) xbar = mean(grillex) varemp = var(grillex) qemp = qnorm(0.975, mean = 175, sd = 10) return(list(xbar = xbar,varemp = varemp,qemp = qemp)) } n = seq(10,5000,4) n stat(n) plot(n, xbar, type = "l", col = 2) plot(n, varemp, type = "l", col = 3) plot(n, qemp, type = "l", col = 4)

---

Merci pour vos réponses

Bonjour,

Effectivement, votre fonction "stat" est faite pour fonctionner avec un entier n.
Pour l'appliquer à la séquence de tous les entiers entre 10 et 50000 espacés de 4, il faut avoir recours aux fonctions de la famille "apply".

Par exemple :

Code:

---

1 2 3

n = seq(10,5000,4)   z<-apply(n,stat)

---

Ici, l'élément z contiendra une liste de listes (car votre fonction "stat" retourne une liste).
Personnellement, je ferais en sorte que la sortie de "stat" soit plutôt un vecteur, car cela facilitera les représentations graphiques je pense.
Par exemple :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

stat_2 <- function(n){ grillex = rnorm(n, mean = 175, sd = 10) xbar = mean(grillex) varemp = var(grillex) qemp = qnorm(0.975, mean = 175, sd = 10) return(v=c(xbar,varemp,qemp)) } n = seq(10,5000,4)   z<-sapply(n,stat_2)   dim(z) [1] 3 1248   plot(n, z[1,], type = "l", col = 2) plot(n, z[2,], type = "l", col = 3) plot(n, z[3,], type = "l", col = 4)

---

Après, il est possible de peaufiner en donnant des noms aux éléments de z, etc...

Bonne continuation :)


Cordialement,

A.D.

Merci pour la réponse

Je voudrais maintenant réaliser des intervalles de confiance (c'est à dire moyenne +ou- le quantile* l'écart type/racine de n.
et enfin de les représenter en graphique.

Merci

De rien, contente d'avoir pu aider :)

Quel problème rencontrez-vous pour la construction de l'intervalle de confiance? Dans la mesure où vous connaissez la formule à utiliser, vous devriez arriver à construire ces intervalles (attention néamoins à la loi suivie par vos variables).

Pour ce qui est des représentations graphiques, pour ajouter une ligne sur un graphique déjà créé (avec la fonction "plot"), vous disposez de la fonction "lines".

N'hésitez pas à préciser le(s) problème(s) rencontré(s) si vous ne parvenez pas à vos fins.
Bonne continuation :)


Cordialement,

A.D.

Merci pour la réponse