<åt sidan class = "s-meddelande s-meddelande__info js-post-notice mb16" role = "status">

Stängt. Denna fråga är utanför ämnet . För närvarande accepteras inte svar.

Kommentarer

  • Lägg till ett reproducerbart exempel för människor att arbeta med.
  • Vid denna punkt har du ' inte beskrivit heteroscedasticitetsstrukturen inom gruppen i din modell än. Kontrollera ?gls och titta på argumentet weights. Om du behöver mer än det, ge ett fungerande exempel som @gung påpekade.
  • @Stefan Tack för ditt förslag. Jag tittade på argumentet weights och fann att det måste anges varFunc och etc. På grund av min begränsade förmåga kan jag dock inte ta reda på hur jag gör det specifikt i mitt fall. Jag vill bara göra min avvikelse mer konstant. Jag har också lagt till min ursprungliga dataset och just nu kan det vara ett lämpligt reproducerbart exempel?

Svar

Vid den här tiden har du inte beskrivit heteroscedasticitetsstrukturen inom gruppen i din modell än. Försök weights=varPower() som visas i exemplet i ?gls . Det blir av med heteroscedasticiteten i ditt fall.

Jämför: 

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

ange bildbeskrivning här 

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

ange bildbeskrivning här

Även om du tittar i kapitel 5.2.1 (sidan 208) i modeller för blandade effekter i S och S -Plus av Pinheiro och Bates 2000, det finns mycket information om variansfunktionerna i nlme. Det här svaret kan också vara till hjälp: Regressionsmodellering med ojämn varians .

Lämna ett svar

Din e-postadress kommer inte publiceras. Obligatoriska fält är märkta *