Uzavřeno. Tato otázka je mimo téma . Momentálně nepřijímá odpovědi.

Komentáře

  • Přidejte reprodukovatelný příklad , se kterým mohou lidé pracovat.
  • V tomto okamžiku jste ještě ' dosud ve svém modelu nepopsali strukturu heteroscedasticity uvnitř skupiny. Zkontrolujte ?gls a podívejte se na argument weights. Pokud potřebujete víc než to, uveďte proveditelný příklad, jak zdůraznil @gung.
  • @Stefan Děkujeme za váš návrh. Podíval jsem se na argument weights a zjistil jsem, že je třeba specifikovat varFunc atd. Vzhledem k mé omezené schopnosti však nemohu přijít na to, jak to udělat konkrétně v mém případě. Jen chci, aby moje variance byla stálejší. Také jsem přidal svůj původní datový soubor a právě teď by to mohl být vhodný reprodukovatelný příklad?

Odpověď

V tomto okamžiku jste ve svém modelu dosud nepopsali strukturu heteroscedasticity v rámci skupiny. Zkuste weights=varPower(), jak je ukázáno v příkladu v ?gls . Tím se ve vašem případě zbavíte heteroscedasticity.

Porovnat: 

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

sem zadejte popis obrázku 

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

zde zadejte popis obrázku

Rovněž pokud se podíváte v kapitole 5.2.1 (strana 208) v Mixed Effects Models v S a S -Plus, Pinheiro a Bates 2000, existuje mnoho informací o variantních funkcích v nlme. Tato odpověď může být také užitečná: Regresní modelování s nerovnoměrnou odchylkou .

Napsat komentář

Vaše e-mailová adresa nebude zveřejněna. Vyžadované informace jsou označeny *