<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
<div dir="ltr">
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
<div dir="ltr">
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
<div dir="ltr"><font class="Apple-style-span" face="'Courier New'" size="2">Boa noite a todos<br id="FontBreak"><span class="Apple-style-span" style="background-color: rgb(255, 255, 255); "><div>Estou fazendo doutorado na UEL no curso de Ciência animal. Ainda sou novato no uso do R e também não tenho base estatística aprofundada em modelos não lineares.</div><div>E preciso de um help, estou adaptando uma fórmula de modelo não linear do proc do SAS para o R, mas li o help, busquei dicas no site <a href="http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase30.html#x32-20100030" target="_blank" style="color: rgb(0, 101, 204); ">http://leg.ufpr.br/~<wbr>paulojus/embrapa/Rembrapa/<wbr>Rembrapase30.html#x32-20100030</a><wbr> e fui seguindo o modelo explicado por esse autor, mas tem algumas coisas que não entendi e não estou conseguindo rodar o modelo.</div><div>Primeiramente explicarei resumidamente o que se trata o modelo. Esse modelo é para poder determinar cinco outras variáveis, através da produção de gás pela fermentação de bacteriais do líquido ruminal in vitro, simulando assim o rumen do animal. As variáveis que tenho em mãos são o volume de gás produzido e o tempo da leitura, a partir deles preciso determinar:</div><div><span style="background-color: rgb(255, 255, 255); ">VcnF = volume de </span><span style="background-color: rgb(255, 255, 255); ">gás produzido pelo carboidratos não fibrosos (mL de gás)</span></div><div></div><div><span style="background-color: rgb(255, 255, 255); ">kdcnf = taxa de degradação dos </span><span style="background-color: rgb(255, 255, 255); ">carboidratos<wbr> não fibrosos (%/hora)</span></div><div><span style="background-color: rgb(255, 255, 255); ">VcF = </span><span style="background-color: rgb(255, 255, 255); ">volume de </span><span style="background-color: rgb(255, 255, 255); ">gás produzido pelo carboidratos fibrosos (mL de gás)</span></div><div><span style="background-color: rgb(255, 255, 255); ">kdcf = </span><span style="background-color: rgb(255, 255, 255); ">taxa de degradação dos </span><span style="background-color: rgb(255, 255, 255); ">carboidratos fibrosos (%/hora)</span></div><div><span style="background-color: rgb(255, 255, 255); ">L= tempo de latência (h), este é o tempo que as bactérias gastam para começar a degradar o substrato.</span></div><div>Mas ao rodar o modelo abaixo aparece o seguinte erro:</div><div><span style="background-color: rgb(255, 255, 255); "><div><span class="Apple-tab-span" style="white-space:pre"> </span>fit30=nls(Y~ (VcnF/(1+exp(2-4*kdcnf*I(L-Tempo)))+ VcF/I(1+exp(2-4*kdcf*(L-Tempo)))),</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>+ data=cra30, </div><div><span class="Apple-tab-span" style="white-space:pre"> </span>+ start=list(VcnF = 10,kdcnf = 0.01,L =0.1,VcF = 50,kdcf =0.01, Tempo=1),</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>+ alpha= (0.05))</div><div>Error in qr.qty(QR, resid) : NA/NaN/Inf in foreign function call (arg 5)</div><div>In addition: Warning message:</div><div>In Ops.factor(lhs, rhs) : - not meaningful for factors</div><div><br></div><div>Aproveitando o ensejo gostaria de perguntar como faço para limitar os parâmetros pois no proc do SAS está assim:</div><div><p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">VcnF = </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">10</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> to </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">200</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> by </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">50</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><span> </span>kdcnf
= </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.01</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> to </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.8</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">
by </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.05</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><span> </span>L = </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.1</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> to </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">5</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> by </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.2</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><span> </span>VcF =
</span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">50</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> to </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">150</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">
by </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">30</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><span> </span>kdcf
= </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.01</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> to </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.8</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">
by </span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0.05</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">;<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><span> </span></span><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">bounds</span><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "> L >=</span><b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">0</span></b><span lang="EN-US" style="background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">;<o:p></o:p></span></p><p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="color: black; background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; "><br></span></p><p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="color: black; background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">Caso alguém tenha materiais relativos a modelos não linear para que possa aprofundar melhor, também contribuirá bastante. </span></p><p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><span lang="EN-US" style="color: black; background-image: initial; background-attachment: initial; background-origin: initial; background-clip: initial; background-color: white; background-position: initial initial; background-repeat: initial initial; ">Agradeço a todos, e espero que alguém possa me ajudar! </span></p><p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;mso-layout-grid-align:none;text-autospace:none"><br></p></div></span></div><div><span style="background-color: rgb(255, 255, 255); "><br></span></div></span></font><br><span class="ecxApple-style-span" style="font-size: 10pt; font-family: Tahoma; "> </span><font class="ecxApple-style-span" face="Arial" color="#ff0000" style="font-family: Tahoma; font-size: 10pt; "> </font><font class="ecxApple-style-span" face="Arial" color="#002060" style="font-family: Tahoma; font-size: 10pt; "> Eduardo Lucas Terra Peixoto</font><div style="font-family: Tahoma; font-size: 10pt; "><span class="ecxApple-style-span" style="color:rgb(0, 32, 96);font-family:Arial">Discente de doutorado - Ciência Animal</span></div><div style="font-family: Tahoma; font-size: 10pt; "><font class="ecxApple-style-span" color="#002060" face="Arial">Universidade Estadual de Londrina<br></font><div><font class="ecxApple-style-span" face="Arial" color="#002060">Cel. (43) 9653-5574</font></div></div></div>
</div>
</div>
</div></body>
</html>