Discussion :

[damien29]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[A. D.]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[damien29]

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

[A. D.]

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.

[damien29]

Merci pour la réponse