<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=utf-8">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Times;
panose-1:2 2 6 3 5 4 5 2 3 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EstiloDeEmail17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=PT-BR link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal><span style='font-family:"Arial","sans-serif";color:#1F497D'>Eu
tentei rodar e não reconheceu a variável ‘x’, não está no seu data set.</span><span
style='font-family:"Arial","sans-serif"'><o:p></o:p></span></p>
<div>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'><span style='color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Olá,<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>eu gostaria de saber se é possível usar estruturas de
correlação nos resíduos quando ajusto modelos não-lineares em mensurações de
dados de produção de gases <i>in vitro</i> (caracterizado por
mensurações acumuladas ao longo do tempo) e sem repetição.<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'><o:p> </o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Breve descrição do banco de dados de gases <i>in vitro</i>:
o experimento foi conduzido em frascos de (125 mL) que contem fruído
ruminal com feno de alfalfa. Esses frascos são acoplados aos sensores que medem
pressão dos gases no interior dos frascos com o conversor digital que registra
automaticamente a cada 5 minutos, durante o período de 48 horas. As mensurações
dos gases são, portanto acumuladas ao longo do tempo. Tenho quatro compostos
(amostras) diferentes que chamei de (y, y2, y3 e y4), em que cada uma amostra
contem 577 mensurações ao longo do tempo (x) que é dado em horas.<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'><o:p> </o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Eu preciso ajustar modelos não-lineares para cada
amostra (exponencial, logístico e logístico bicompartimental). O modelo
logístico bicompartimental foi o que melhor se ajustou sem usar estrutruras de
correlação, porém, a análise residual possuí uma forte dependência (visto que
os dados são acumulados, ou seja, a mensuração no tempo (t+1) depende da
mensuração anterior do tempo t.<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Então, eu teria que entrar com estruturas de correlações no
modelo logístico bicompartimental para melhorar a análise residual? <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'><o:p> </o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> Eu tentei algumas estruturas usando o argumento
“correlation”, da função gnls. No quadro abaixo segue o seguinte resumo das
estruturas inseridas no modelo, com os respectivos resultados gerados pela FAC
ou os erros retornados pelo software R.<o:p></o:p></p>
<table class=MsoNormalTable border=0 cellspacing=0 cellpadding=0 width=619
style='width:464.5pt;border-collapse:collapse'>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>modelo<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border:solid black 1.0pt;
border-left:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>correlation<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border:solid black 1.0pt;
border-left:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Erros<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border:solid black 1.0pt;
border-left:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>FAC - residual<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M1<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>sem estrutrura de correlação<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Sem erro<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Tendência : possivelmente independência e heterogeneidade<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M2<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>AR1: corAR1()<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Erro: step halving factor reduced below minimum in NLS step<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M3<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>AR1: corAR1(form=~x)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Sem erro<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Tendência : possivelmente independência e heterogeneidade<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M4<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>AR2: corAR1(p=2, q=0)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Coefficient matrix not invertible<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M5<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>ARIMA: corARMA(p=2,q=1)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>step halving factor reduced below minimum in NLS step<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M6<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>ARIMA: corARMA(p=2,q=2)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Coefficient matrix not invertible<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M7<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>ARIMA: corARMA(p=3,q=1)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Coefficient matrix not invertible<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
<tr>
<td width=151 valign=top style='width:113.15pt;border:solid black 1.0pt;
border-top:none;padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>M8<o:p></o:p></p>
</td>
<td width=187 valign=top style='width:140.35pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>AR(): corAR1(form=~1)<o:p></o:p></p>
</td>
<td width=130 valign=top style='width:97.8pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>Erro: step halving factor reduced below minimum in NLS step<o:p></o:p></p>
</td>
<td width=151 valign=top style='width:113.2pt;border-top:none;border-left:
none;border-bottom:solid black 1.0pt;border-right:solid black 1.0pt;
padding:0cm 5.4pt 0cm 5.4pt'>
<p class=MsoNormal align=center style='mso-margin-top-alt:auto;text-align:
center'>-<o:p></o:p></p>
</td>
</tr>
</table>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Observação: Eu fiz a análise usando a metodologia de
séries temporais (modelo linear não estacionário – com uma diferença d=1) e a
FAC e FACP mostraram que os resíduos eram ruído branco, como nos modelos
ARIMA(p=1, d=1, q=0), ARIMA(p=2, d=1, q=0) e ARIMA(p=3, d=1, q=0). Seria
possível ajustar modelos não-lineares (logístico bicompartimental) usando
estruturas de correlações como o ARIMA?<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Os dados e o script seguem em anexo com todos os modelos
que tentei ajustar (somente para uma amostra, ou seja, correspondente as
colunas chamadas de y e x). Agradeço se puderem me dar uma orientação.<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>rm(list=ls(all=TRUE))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#pacotes necessários para a análise.<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>require(manipulate)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>require(nlstools)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>require(car)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>require(nlme)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>require(lmtest)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># modelos estudados<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#a * (1 - exp(-b * (t - c))) # modelo exponencial<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#a /(1+exp(2+4*b*(c-t))) # modelo
logistico<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#a/(1+exp(2+4*b*(c-t)))+d/(1+exp(2+4*e*(c-t))) #
modelo logistico 2 pool<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#-----------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>setwd('E:\\autoregressivo_modelo')<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>gas01<-read.csv("gas.csv",sep = ",",
header = TRUE)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>gas01<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>head(gas01)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>attach(gas01)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>t=Time/60 # para deixar o tempo em horas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>dim(gas01)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>gas=cbind(gas01,t)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>names(gas) <-
c("Time","y","y2","y3", "y4",
"x") # variáveis que quero ajustar y (variavel dependente gases) e x
o tempo em horas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>head(gas)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>dim(gas)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>attach(gas)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#===================================================================<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste usando o gnls<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#Ajuste do modelo exponencial -lote I considerando os erros
independentes<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod1 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas, start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367,
e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod1) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod1), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod1)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod1)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando - AR()<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod2 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corAR1(),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod2) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod2), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod2)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod2)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(x,mod2$residuals)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod2)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#
------------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando - AR()<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod3 <- gnls(y ~ a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corAR1(form=~x),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod3) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod3), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod3)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod3)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(x,mod3$residuals)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot( ACF(mod3, maxLag =24), alpha = 0.05 )<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#
------------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando - AR2()<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod4 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corARMA(p=2,q=0),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=0.78, b=0.06, c=-5.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod4) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod4), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod4)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod4)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot( ACF(mod4, maxLag =24), alpha = 0.05 )<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#
------------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando -
ARIMA(p=2,q=1)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod5 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corARMA(p=2,q=1),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod5) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod5), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod5)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod5)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot( ACF(mod5, maxLag =24), alpha = 0.05 )<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#
------------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando -
ARIMA(p=2,q=2)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod6 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corARMA(p=2,q=2),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod6) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod6), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod6)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod6)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot( ACF(mod6, maxLag =24), alpha = 0.05 )<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>#
------------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'># Ajuste do modelo Logístico 2 pool considerando -
ARIMA(p=3,q=1)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>mod7 <- gnls(y ~
a/(1+exp(2+4*b*(c-x)))+d/(1+exp(2+4*e*(c-x))),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>data=gas,correlation=corARMA(p=3,q=1),<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>start=c(a=1.1179, b=0.0562, c=-9.5923, d=0.367, e=0.0226))<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>summary(mod7) # quadro de estimativas<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>R2(residuals(mod7), gas$y)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>residuosnI<-summary(mod7)$resid<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>shapiro.test(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>ks.test(rnorm(577),residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot(mod6)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>hist(residuosnI)<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>plot( ACF(mod7, maxLag =24), alpha = 0.05 )<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'> <o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Agradeço desde já,<o:p></o:p></p>
<p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;
text-align:justify'>Adriele.<o:p></o:p></p>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<p class=MsoNormal>-- <o:p></o:p></p>
<div>
<div>
<p class=MsoNormal><span style='font-family:"Times","serif"'> Adriele
Giaretta Biase.<o:p></o:p></span></p>
</div>
<div>
<p class=MsoNormal><span style='font-family:"Times","serif"'> Mestre
em Estatística e Experimentação Agropecuária - UFLA. <br>
Doutoranda em Estatística e Experimentação Agronômica - ESALQ/
USP<o:p></o:p></span></p>
</div>
<div>
<p class=MsoNormal><span style='font-family:"Times","serif"'> Contato:
(19) 8861-0619.<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="https://www.avast.com/antivirus">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" alt="Avast logo" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Este email foi escaneado pelo Avast antivírus.
<br><a href="https://www.avast.com/antivirus">www.avast.com</a>
</p>
</td>
</tr>
</table>
<br />
</body>
</html>