Fermé. Cette question est hors sujet . Il naccepte pas les réponses actuellement.

Commentaires

  • Veuillez ajouter un exemple reproductible pour que les gens puissent travailler avec.
  • À ce stade, vous navez ' pas décrit la structure dhétéroscédasticité intra-groupe dans votre modèle. Vérifiez ?gls et regardez largument weights. Si vous avez besoin de plus que cela, veuillez fournir un exemple réalisable comme @gung la souligné.
  • @Stefan Merci pour votre suggestion. Jai jeté un coup dœil à largument weights et jai trouvé quil fallait spécifier varFunc et etc. Cependant, en raison de mes capacités limitées, je ne peux pas comprendre comment faire cela spécifiquement dans mon cas. Je veux juste rendre ma variance plus constante. De plus, jai ajouté mon ensemble de données dorigine et pour le moment, cela pourrait être un exemple reproductible approprié?

Réponse

À ce stade, vous navez pas encore décrit la structure dhétéroscédasticité intra-groupe dans votre modèle. Essayez weights=varPower() comme indiqué dans lexemple de ?gls . Cela élimine lhétéroscédasticité dans votre cas.

Comparez:

m1 <- gls(salary ~ age*sex) plot(m1) 

entrez la description de limage ici

m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2) 

entrez la description de limage ici

Aussi si vous regardez au chapitre 5.2.1 (page 208) dans Modèles deffets mixtes en S et S -Plus par Pinheiro et Bates 2000, il y a beaucoup dinformations sur les fonctions de variance dans nlme. Cette réponse peut également être utile: Modélisation de régression à variance inégale .

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *