Bonjour, je suis actuellement en stage.
Je dois simmuler des données qui me serviront pour réaliser un modèle mixte (lme). Je dois simuler des données sur des oiseaux sur le variable quantité lait ingurgité (il s'agit d'oiseaux receuillis qui sont malade) pour des oiseaux seins la quantité normal est 29.65 (en moyenne) ml.
Je part de la formule suivante Y=m+sigma*N(0,1)
m=moyenne des quelques données que j'ai à ma disposition.
sigma=exp^(-alpha*i) avec alpha=1 et i va de 1 à 7
On veut voir l'évolution de la quantité de nouriture prise durant 7 ans avec vérification chaque année pour les oiseaux recevant un traitement pour ouvrir l'appéti.
voici ce que j'ai fait :
n<-60 # population de 60 oiseaux
k<-7 # 7 relevé (relevé annuel pendant 7 ans)
err<-matrix(rnorm(n*k,0,1),nr=n,nc= k) # simulation de la loi normale sur les 60 oiseaux (n) pour les 7 années (k)
err
alpha<-(-1)
for (i in 1:7){ # i = numéro de l'année allant de 1 à 7
y<-22 * exp(alpha*i)+ exp(alpha*i)* err # 22=moyenne de nouriture ingurgité pour les qq oiseaux que j'ai.
}
y
Trait1<-cbind(1,(1:60)) # colonne traitement 1 pr les 60 oiseaux
Trait1
Y<-matrix(cbind(Trait1,y),ncol=9) # jeu de données final
colnames(Y)<-c('Traitement','Oiseau','P1',' P2','P3','P4','P5','P6','P7')# jeu de données final avec le nom des colonnes
Y
le truc c'est qu'il parte d'une quantité de 22 environ pour augmenté chaque année jusqu'à se rapprocher de 29,65 sans jamais l'égaler etant donnée qu'il s'agit d'oiseaux malades.
Je pense que je me suis trompé dans la formul du y car toute les données sont très petites (0,0456...) et qu'il n'y a pas d'évolution.
quelqu'un peut-il m'aider?....
Merci
A.
-----