Fechado. Esta pergunta está fora do tópico . Atualmente não está aceitando respostas.

Comentários

  • Adicione um exemplo reproduzível para as pessoas trabalharem.
  • Neste ponto, você ainda ' t descreveu a estrutura de heterocedasticidade dentro do grupo em seu modelo. Verifique ?gls e dê uma olhada no argumento weights. Se você precisar de mais do que isso, forneça um exemplo viável como @gung apontou.
  • @Stefan Obrigado por sua sugestão. Eu dei uma olhada no argumento weights e descobri que ele precisa especificar varFunc 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?

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) 

insira a descrição da imagem aqui

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

insira a descrição da imagem aqui

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 .

Deixe uma resposta

O seu endereço de email não será publicado. Campos obrigatórios marcados com *