<div dir="ltr">Olá pessoal,<div><br></div><div><br></div><div>eu enviei uma dúvida sobre qual melhor estrutura de correlação residual devo considerar num modelo não linear (logístico bicompartimental) para ajuste de dados com produção de gases in vitro, que possui uma alta dependência residual.</div><div><br></div><div>Eu organizei melhor o script e os dados, em anexo (novo anexo), em que os dados que estão na planilha, a primeira coluna é o x variável independente representando o tempo. E o y são as mensurações de gases.</div><div>Então a nova planilha está exatamente as informações no qual preciso obter o ajuste do conjunto.</div><div><br></div><div>Tentei inserir algumas estruturas de correlação, mas não deu certo. Alguém sabe como devo ajustar esse conjunto de dados (qual melhor estrutura de correlação residual que posso considerar) usando modelos não lineares na função gnls?</div><div><br></div><div><br></div><div><br></div><div>Adriele</div><div><br></div><div><br></div><div><div><br></div><div>#-----------------------------------------------------------------------</div><div>setwd('C:\\Users\\Adriele\\Desktop')</div><div>gas<-read.csv("gases_lista_R.csv",sep = ",", header = TRUE)</div><div>head(gas)</div><div>attach(gas)</div><div>#-----------------------------------------------------------------------</div><div><br></div><div># Ajuste usando o gnls</div><div>#Ajuste do modelo exponencial -lote I considerando os erros independentes ( o que não ocorre!!)</div><div>mod1 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas, start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))</div><div>summary(mod1) # quadro de estimativas</div><div><br></div><div># função para o cálculo do R2</div><div>R2 <- function(residuals, observed){</div><div>1-(sum(residuals)^2)/sum((observed-mean(observed))^2)</div><div>}</div><div><br></div><div>R2(residuals(mod1), gas$y)</div><div>residuosnI<-summary(mod1)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod1)</div><div>hist(residuosnI)</div><div>plot( ACF(mod1, maxLag =24), alpha = 0.05 )</div><div><br></div><div># ------------------------------------------------------------------------</div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - AR()</div><div>mod2 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corAR1(),</div><div>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))</div><div>summary(mod2) # quadro de estimativas</div><div>R2(residuals(mod2), gas$y)</div><div>residuosnI<-summary(mod2)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod2)</div><div>hist(residuosnI)</div><div>plot(x,mod2$residuals)</div><div>plot(mod2)</div><div>plot( ACF(mod2, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - AR()</div><div>mod3 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corAR1(form=~x),</div><div>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))</div><div>summary(mod3) # quadro de estimativas</div><div>R2(residuals(mod3), gas$y)</div><div>residuosnI<-summary(mod3)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod3)</div><div>hist(residuosnI)</div><div>plot(x,mod3$residuals)</div><div>plot( ACF(mod3, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - AR2()</div><div>mod4 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corARMA(p=2,q=0),</div><div>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))</div><div>summary(mod4) # quadro de estimativas</div><div>R2(residuals(mod4), gas$y)</div><div>residuosnI<-summary(mod4)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod4)</div><div>hist(residuosnI)</div><div>plot( ACF(mod4, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - ARIMA(p=2,q=1)</div><div>mod5 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corARMA(p=2,q=1),</div><div>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))</div><div>summary(mod5) # quadro de estimativas</div><div>R2(residuals(mod5), gas$y)</div><div>residuosnI<-summary(mod5)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod5)</div><div>hist(residuosnI)</div><div>plot( ACF(mod5, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - ARIMA(p=2,q=2)</div><div>mod6 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corARMA(p=2,q=2),</div><div>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))</div><div>summary(mod6) # quadro de estimativas</div><div>R2(residuals(mod6), gas$y)</div><div>residuosnI<-summary(mod6)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod6)</div><div>hist(residuosnI)</div><div>plot( ACF(mod6, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - ARIMA(p=3,q=1)</div><div>mod7 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corARMA(p=3,q=1),</div><div>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))</div><div>summary(mod7) # quadro de estimativas</div><div>R2(residuals(mod7), gas$y)</div><div>residuosnI<-summary(mod7)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod6)</div><div>hist(residuosnI)</div><div>plot( ACF(mod7, maxLag =24), alpha = 0.05 )</div><div># ------------------------------------------------------------------------</div><div><br></div><div># Ajuste do modelo Logístico 2 pool considerando - AR() a ordem das obsevações</div><div>mod8 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),</div><div>data=gas,correlation=corAR1(form=~1),</div><div>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))</div><div>summary(mod8) # quadro de estimativas</div><div>R2(residuals(mod8), gas$y)</div><div>residuosnI<-summary(mod8)$resid</div><div>shapiro.test(residuosnI)</div><div>ks.test(rnorm(577),residuosnI)</div><div>plot(mod6)</div><div>hist(residuosnI)</div><div><br></div>-- <br><div class="gmail_signature"><div style="font-family:times,serif;margin:0cm 0cm 0.0001pt;font-size:16px"><span style="white-space:pre">    </span>Adriele Giaretta Biase.</div><div style="font-family:times,serif;margin:0cm 0cm 0.0001pt;font-size:16px"><span style="white-space:pre">      </span>Mestre em  Estatística e Experimentação Agropecuária - UFLA. <br><span style="white-space:pre">        </span>Doutoranda em Estatística e Experimentação Agronômica - ESALQ/ USP</div><div style="font-family:times,serif;margin:0cm 0cm 0.0001pt;font-size:16px"><span style="font-size:12pt;white-space:pre">       </span><span style="font-size:12pt">Contato: (19) 8861-0619.</span></div></div>
</div></div>