2

Preciso ajustar uma função não linear através das rotinas nls ou nlsLM, entretanto, meu dados geram resíduos heterocedásticos quando calculados com essas rotinas, o que não me permite ajustar diretamente através da soma dos quadrados dos resíduos

Já tenho a equação que torna os resíduos homocedásticos, mas não sei como fazer para que a rotina utilize essa forma para calcular os resíduos no lugar do cálculo direto de (observado-esperado)^2

8
  • Então José, você já tem alguma coisa pronta? Algum trecho do código? Poderia tbm explicar melhor o que quer dizer com resíduos heterocedásticos e resíduos homocedásticos isso me parece ter mais a ver com Regressão Linear -- Matematica. Commented 10/02/2015 às 16:04
  • Olá, @MarceloBonifazio! Eu estou usando diretamente a função nlsLM, apenas insiro meu conjunto de dados e mando rodar. Resíduos heterocedásticos ocorrem quando a diferença entre o observado e o esperado (resíduos) apresentam uma tendência em sua distribuição ao longo de x. As rotinas {nls} e {nlsLM} ajustam o modelo através da minimização do quadrado dos resíduos, preciso alterar essa parte do código Commented 10/02/2015 às 16:13
  • José, o código do nls é aberto e pode ser visto digitando nls no console do R. Caso você queira editar pode recriar a função, fazendo nls <- código mostrado... No entanto, isso não deve ser fácil. Existem outros jeitos de corrigir a heterocedasticidade mudando o modelo e não código, tem certeza que é isso que você precisa fazer? Commented 11/02/2015 às 0:28
  • @DanielFalbel, preciso alterar o código porque minha matriz de dados tem esse problema. Eu cheguei a abrir o código do nls, mas não consegui identificar a parte que calcula os resíduos pra fazer a alteração. Sabe me dizer isso? Commented 11/02/2015 às 11:03
  • Acho que não compensa alterar o código por custo benefício. Por que você não tenta ajustar um modelo com uma função para a variância tbm? A função gnls ajusta modelos com uma função para a variância passada para o argumento weights Commented 11/02/2015 às 16:31

1 Resposta 1

1

Fiz um exemplo que deve te ajudar:

# dados do modelo
DNase1 <- subset(DNase, Run == 1)

# ajuste do mdoelo
fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)

# definicao de a e b
a <- 1
b <- 2

# cálculo do resíduo
residuo <- (DNase1$density - predict(fm1DNase1))/((a*DNase1$density)^b)

Na última parte do código, coloquei a variável resposta no lugar do seu x, mas não sei se é exatamente isso que você queria.

3
  • A ideia é essa mesmo, mas o que preciso é fazer com que a rotina nls ajuste o modelo utilizando essa equação para os resíduos no lugar de simplesmente (observado-esperado)^2 Commented 11/02/2015 às 20:52
  • você poderia postar o código do modelo que você está ajustanto? acho que vai ficar mais fácil de ajudar Commented 12/02/2015 às 0:19
  • Claro. Esse é o código que estou usando para ajustar o modelo: poly.2p.fitted<-gnls(Weight ~ (a1 * Length^b1) * (1/(1 + exp(rate * (Length -scp))))+(a2 * Length^b2) * (1-(1/(1 + exp(rate * (Length-scp))))), data=dados, start=list(a1=0.1,b1=2.5,a2=0.1,b2=3,rate=10,scp=11), control=gnlsControl(maxIter=99999)) Commented 12/02/2015 às 0:22

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .