Fechado. Esta pergunta está fora do tópico . Atualmente não está aceitando respostas.
Comentários
Resposta
Neste ponto, você ainda não descreveu a estrutura de heterocedasticidade dentro do grupo em seu modelo. Tente weights=varPower()
conforme mostrado no exemplo em ?gls
. Isso elimina a heterocedasticidade do seu caso.
Compare:
m1 <- gls(salary ~ age*sex) plot(m1)
m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2)
Também se você olhar no Capítulo 5.2.1 (página 208) em Modelos de efeitos mistos em S e S -Além disso, por Pinheiro e Bates 2000, há muitas informações sobre as Funções de Variância em nlme
. Esta resposta também pode ser útil: Modelagem de regressão com variância desigual .
?gls
e dê uma olhada no argumentoweights
. Se você precisar de mais do que isso, forneça um exemplo viável como @gung apontou.weights
e descobri que ele precisa especificarvarFunc
e etc ,. No entanto, devido à minha capacidade limitada, não consigo descobrir como fazer isso especificamente no meu caso. Só quero tornar minha variância mais constante. Além disso, adicionei meu conjunto de dados original e agora ele pode ser um exemplo reproduzível apropriado?