<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">Quero comparar modelos para duas espécies de peixes. Consigo fazer usando nls mas
não gnls. Alguém tem alguma idéia de como fazer? <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">dados<-structure(list(Peixes = structure(c(1L, 1L,
1L, 1L, 1L, 1L, 1L, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L),
.Label = c("1", <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">"2"), class = "factor"), Tempo =
c(0L, 0L, 0L, 30L, 30L, 30L, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L,
270L, 270L, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L,
90L, 240L, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">240L, 240L,
270L, 270L, 270L), Peso = c(6.09, 6.09, 6.09, 65.73, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">52, 53.13,
102.88, 115.8, 129.31, 161.06, 174.94, 213.63, 233.67, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">291.93,
301.33, 258.47, 326.84, 337.04, 8.92, 8.92, 8.92, 78.8, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">72.07, 85.8,
167.82, 146.13, 179.31, 255.75, 230.73, 207.06, <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">320.67, 280.2,
243.93, 376.55, 441.96, 349.84)), .Names = c("Peixes", <o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">"Tempo",
"Peso"), row.names = c(NA, -36L), class = "data.frame")<o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US">library(nlme)<o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-ansi-language:EN-US"><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span class="aesltt-10x-x-120"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";mso-ansi-language:EN-US">m1</span></span><span class="aesltt-10x-x-120"><span lang="EN-US" style="font-size:10.0pt;font-family:
"Courier New";mso-fareast-font-family:Calibri;mso-ansi-language:EN-US"> <-gnls(</span></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:
Calibri;mso-ansi-language:EN-US">Peso~a*(1-b*(exp(-c*Tempo)))</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";mso-ansi-language:
EN-US">^3</span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";
mso-fareast-font-family:Calibri;mso-ansi-language:EN-US"> <span class="aesltt-10x-x-120">, </span></span><span lang="EN-US" style="font-size:
10.0pt;font-family:"Courier New";color:#454545;background:white;mso-ansi-language:
EN-US">weights=varPower(), data=dados,</span><span class="aesltt-10x-x-120"><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:
Calibri;mso-ansi-language:EN-US"> </span></span><span lang="EN-US" style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri;
mso-ansi-language:EN-US">start=c(a=350,b=1, c=0.01)<span class="aesltt-10x-x-120">)<o:p></o:p></span></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span lang="EN-US" style="font-size:10.0pt;
font-family:"Courier New";mso-fareast-font-family:Calibri;mso-ansi-language:
EN-US"><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New"">m2</span></span><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New";
mso-fareast-font-family:Calibri"><- nls(</span></span><span style="font-size:10.0pt;font-family:"Courier New"">Peso</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">~a*(1-b*(exp(-c*Tempo)))</span><span style="font-size:10.0pt;font-family:"Courier New"">^3,</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#454545;
background:white">data=dados,</span><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri"> </span></span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">start=c(a=300,b=10,
c=0.05)<span class="aesltt-10x-x-120">)<o:p></o:p></span></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri"><o:p> </o:p></span></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span style="font-size:10.0pt;font-family:"Courier New"">m3<-gnls(Peso</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))</span><span style="font-size:10.0pt;font-family:"Courier New"">^3,</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">
<span class="aesltt-10x-x-120">, </span></span><span style="font-size:10.0pt;
font-family:"Courier New";color:#454545;background:white">weights=varPower(), data=dados,</span><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New";
mso-fareast-font-family:Calibri"> </span></span><span style="font-size:
10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">start=list(a=c(300,300),b=c(10,10),
c=c(0.05,0.05))<span class="aesltt-10x-x-120">)</span><o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; "><span style="font-size:10.0pt;font-family:"Courier New"">m4<-nls(Peso</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))</span><span style="font-size:10.0pt;font-family:"Courier New"">^3,</span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">
</span><span style="font-size:10.0pt;font-family:"Courier New";color:#454545;
background:white">data=dados,</span><span class="aesltt-10x-x-120"><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri"> </span></span><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-font-family:Calibri">start=list(a=c(300,300),b=c(10,10),
c=c(0.05,0.05))<span class="aesltt-10x-x-120">)</span><o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New""><o:p> </o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">anova(m1,m3)<o:p></o:p></span></div>
<div class="MsoNormal" style="margin-bottom: 0.0001pt; "><span style="font-size:10.0pt;font-family:"Courier New"">anova(m2,m4)<o:p></o:p></span></div></div></div></body></html>