[R-br] homocedasticidade modelo não linear

[Fernando Antonio de souza]

---------- Mensagem encaminhada ----------
De: Fernando Antonio de souza <nandodesouza em gmail.com>
Data: 14 de março de 2012 13:50
Assunto: duvida sobre presunção de homocedasticidade
Para: r-br em listas.c3sl.ufpr.br


Caros amigos
Estou trabalhando com dados de desenvolcimento de fetos ovinos. Estou
modelando o seu crescimento ao longo da gestação. Para tal tenho utilizado
o modelo de gompertz da forma log(Y)= a-b*exp(-c*x) sendo x a idade
gestacional em dias e Y o peso do feto. Os meus dados atendem o presuposto
de normalidade, porém não  o faz para o pressuposto de igualdade das
variâncias. Tenho lido na literatura a possibilidade de modelar a variância
nestes casos. A variância residual observada graficamente  se eleva e em
seguida rebaixa conforme as médias  do modelo se elevam. Na própria
literatura comenta que este tipo de comportamento pode ser modelado
utilizando o método VarPower() na função glns do pacote nlme. Eu assim o
fiz entretanto o novo modelo não  obteve melhor ajuste que o anterior
quando comparados pela função Anova (p-value e AIC,BIC). Minha dúvida é a
seguinte sem a igualdade das variâncias eu não posso validar o modelo.
Existe alguma forma de corrigir este problema? abaixo disponho do CMR

library(nlme)
dput(gestsimp)
gompertz1<-nls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),data=gestsimp)
gompertz2<-gnls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),weights=varPower(),data=gestsimp)
summary(gompertz1)
plot(log10(PESO.FETO)~Gest,data=gestsimp,main="Peso do feto em função da
idade de cordeiros")
curve(5.697742-8.360311*exp(-0.003531*x),add=T,col=2)

caso o dput()não funcione segue a saída dele

dput(gestsimp)structure(list(Animal = c(3L, 51L, 83L, 89L, 91L, 92L, 93L, 48L,
52L, 55L, 69L, 76L, 78L, 11L, 25L, 41L, 42L, 46L, 50L, 53L, 57L,
10L, 19L, 20L, 30L, 38L, 45L), Manejo = c(1L, 2L, 2L, 2L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 1L, 1L, 1L, 1L), Gest = c(90L, 90L, 90L, 90L, 90L, 90L,
90L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 130L, 130L, 130L,
130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L),
    Fetos = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
    ), PV = c(51.1, 46.4, 53.5, 58.2, 30.7, 43.2, 41.5, 42, 43.4,
    48.1, 42.2, 47.5, 32, 64.1, 49.9, 40.5, 42, 40.4, 50, 49,
    52, 60.1, 55.9, 57, 51.5, 45, 42.8), PCVZ = c(36.93, 36.135,
    46.125, 48.02, 23.095, 33.75, 32.575, 34.99, 34.465, 39.578,
    30.84, 38.045, 27.455, 54.797, 43.595, 32.12, 33.841, 33.16,
    42.06, 40.605, 42.325, 50.436, 46.515, 48.845, 43.075, 36.55,
    35.435), PESO.FETO = c(0.215, 0.555, 0.49, 0.53, 0.225, 0.465,
    0.61, 1.465, 0.585, 1.755, 0.62, 1.615, 1.045, 0.985, 2.79,
    2.42, 2.286, 2.41, 2.3, 3.1, 3.105, 5.38, 4.065, 4.258, 3.878,
    3.49, 3.032), FETOGKG = c(60.5065, 50.8621, 67.9386, 32.7586,
    41.7391, 46.1815, 47.7506, 51.0421, 66.0087, 63.5711, 50.8569,
    60.5784, 49.8817, 62.1197, 60.1357, 0, 0, 62.2886, 0, 55.4768,
    53.5301, 68.0826, 58.413, 59.6224, 61.1382, 50.2531, 62.5518
    )), .Names = c("Animal", "Manejo", "Gest", "Fetos", "PV",
"PCVZ", "PESO.FETO", "FETOGKG"), row.names = c(11L, 13L, 15L,
19L, 20L, 21L, 22L, 30L, 31L, 32L, 35L, 36L, 38L, 42L, 45L, 46L,
47L, 48L, 49L, 51L, 55L, 63L, 66L, 67L, 71L, 76L, 80L), class =
"data.frame", na.action = structure(22L, .Names = "62", class =
"omit"))
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120315/cb7ef675/attachment.html>