<span class="Apple-style-span" style="border-collapse:separate;color:rgb(0,0,0);font-family:Arial;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><span class="Apple-style-span" style="font-family:'Lucida Console';font-size:13px;line-height:15px;text-align:-webkit-left;white-space:pre-wrap"><pre tabindex="0" class="GJWPQFQDK4" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2">
segue comando: (pode ser devido aos NA's)<br><br><span class="Apple-style-span" style="border-collapse:separate;color:rgb(0,0,0);font-family:Arial;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><span class="Apple-style-span" style="font-family:'Lucida Console';font-size:13px;line-height:15px;text-align:-webkit-left;white-space:pre-wrap"><pre tabindex="0" class="GJWPQFQDK4" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2">
<span class="GJWPQFQDB4 ace_keyword" style="color:blue">dput(simple)</span><br></pre></span></span>structure(list(Manejo = c(2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 1, 
1, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA), Gest = c(90, 90, 90, 90, 
90, 90, 110, 110, 110, 110, 110, 110, 110, 130, 130, 130, 130, 
130, 130, 130, 130, 140, 140, 140, 140, 140, 140, 140, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), PESO.FETO = c(0.555, 
0.49, 0.53, 0.225, 0.465, 0.61, 0.215, 1.465, 0.585, 1.755, 0.62, 
1.615, 1.045, 2.79, 2.42, 2.286, 2.41, 2.3, 3.1, 3.105, 2.56, 
5.38, 0.985, 4.065, 4.258, 3.878, 3.49, 3.032, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA)), .Names = c("Manejo", "Gest", 
"PESO.FETO"), row.names = c(NA, -57L), class = "data.frame")</pre></span></span><br><br><div class="gmail_quote">Em 10 de março de 2012 21:12, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ainda nao foi dessa vez... mas esta' facil.. Execute o seu codigo,<br>
dai' veja no R o que vai aparecer... vai aparecer algo esquisito como<br>
abaixo:<br>
<br>
> dput(data.frame(x=rnorm(10)))<br>
structure(list(x = c(-0.543332683249016, -1.72683892829404, -1.46539343776498,<br>
-1.39513469551024, 0.809122589671721, -1.58459611384001, -0.818826156554395,<br>
-0.391288255253397, -0.579453920817577, 0.473369713399808)), .Names =<br>
"x", row.names = c(NA,<br>
-10L), class = "data.frame")<br>
<br>
<br>
dai' vc copia pra gente esse resultado estranho (q no meu caso comeca<br>
com "structure")... :)<br>
<br>
b<br>
<br>
2012/3/11 Fernando Antonio de souza <<a href="mailto:nandodesouza@gmail.com">nandodesouza@gmail.com</a>>:<br>
<div class="HOEnZb"><div class="h5">> Veja se agora foi<br>
><br>
><br>
> ###cria a função gompertz e chama a função manipulate para encontrar valores<br>
> do chute<br>
> simple<-as.data.frame(cbind(Manejo=gestsimp$Manejo,Feto=gestsimp$Fetus,Gest=gestsimp$Gest,PESO.FETO=gestsimp$PESO.FETO))<br>
> dput(simple)<br>
> plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=simple)<br>
> attach(simple)<br>
> gompertz<-function(x,a,b,c){a-b*exp(-c*x)}<br>
> gompertz(c(90,110,130,140),1.5,20,0.026)<br>
> library(manipulate)<br>
> start <- list()<br>
> manipulate({plot(log(PESO.FETO)~Gest,data=simple)<br>
><br>
>             curve(gompertz(x,a=a,b=b,c=c),add=TRUE)<br>
>             start <- list(a=a,b=b,c=c)},<br>
>            a=slider(0.00001, 10, initial=0.001),<br>
>            b=slider(-0.0001, 50, initial=0.01),<br>
>            c=slider(0.00007, 0.5, initial=0.0001))<br>
> ##----------------------------------------------------------------------------<br>
> #calcula os parâmetros do modelo de gompertz<br>
> gompertz<-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=2.12,b=14.99,c=0.018),data=simple)<br>
> summary(gompertz)#visualiza os parâmetros estimados<br>
><br>
><br>
><br>
><br>
><br>
> Em 10 de março de 2012 20:39, Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
> escreveu:<br>
><br>
>> Fernando, faltou a saida do dput()...<br>
>><br>
>> 2012/3/10 Fernando Antonio de souza <<a href="mailto:nandodesouza@gmail.com">nandodesouza@gmail.com</a>>:<br>
>> > Segue o CMR, os datos encontra-se no dput(). Gostaria de acrescentar<br>
>> > algo<br>
>> > mais. Meus  dados estão bastante dispersos. Pode ser por isto que não<br>
>> > está<br>
>> > convergindo?Há alguma solução?<br>
>> ><br>
>> > ###cria a função gompertz e chama a função manipulate para encontrar<br>
>> > valores<br>
>> > do chute<br>
>> > plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=gestsimp)<br>
>> > attach(gestsimp)<br>
>> > gompertz<-function(x,a,b,c){a-b*exp(-c*x)}<br>
>> > gompertz(c(90,110,130,140),1,5,0.0006)<br>
>> > library(manipulate)<br>
>> > start <- list()<br>
>> > manipulate({plot(log(PESO.FETO)~Gest,data=gestsimp)<br>
>> >             curve(gompertz(x,a=a,b=b,c=c),add=TRUE)<br>
>> >             start <- list(a=a,b=b,c=c)},<br>
>> >            a=slider(0.00001, 10, initial=0.001),<br>
>> >            b=slider(-0.0001, 50, initial=0.01),<br>
>> >            c=slider(0.00007, 0.5, initial=0.0001))<br>
>> ><br>
>> > ##----------------------------------------------------------------------------<br>
>> > #calcula os parâmetros do modelo de gompertz<br>
>> ><br>
>> > gompertz<-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=1.68,b=24.599,c=0.026),subset=c(Gest>0,Fetos==1,Manejo==1),data=tabela)<br>
>> > summary(gompertz)#visualiza os parâmetros estimados<br>
>> > dput(gestsimp)<br>
>> ><br>
>> > Em 10 de março de 2012 19:20, Ivan Bezerra Allaman<br>
>> > <<a href="mailto:ivanalaman@yahoo.com.br">ivanalaman@yahoo.com.br</a>><br>
>> > escreveu:<br>
>> >><br>
>> >> Provavelmente seus chutes estão errados! Acesse o site<br>
>> >> <a href="http://ridiculas.wordpress.com/" target="_blank">http://ridiculas.wordpress.com/</a> que provavelmente irás resolver o seu<br>
>> >> problema. Do contrário, nos forneça um CMR.<br>
>> >><br>
>> >> (S,f,P)<br>
>> >> Allaman<br>
>> >><br>
>> >><br>
>> >> \begin{signature}<br>
>> >> <<>>=<br>
>> >> Prof. Dr. Ivan Bezerra Allaman<br>
>> >> Universidade Estadual de Santa Cruz<br>
>> >> Departamento de Ciências Exatas e Tecnológicas<br>
>> >> Ilhéus/BA - Brasil<br>
>> >> Fone: <a href="tel:%2B55%2073%203680-5076" value="+557336805076">+55 73 3680-5076</a><br>
>> >> E-mail: <a href="http://ivanalaman@yahoo.com.br/ivanalaman@gmail.com" target="_blank">ivanalaman@yahoo.com.br/ivanalaman@gmail.com</a><br>
>> >> @<br>
>> >> \end{signature}<br>
>> >><br>
>> >> _______________________________________________<br>
>> >> R-br mailing list<br>
>> >> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> >> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>> >> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> >> código mínimo reproduzível.<br>
>> ><br>
>> ><br>
>> ><br>
>> > _______________________________________________<br>
>> > R-br mailing list<br>
>> > <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> > <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>> > Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> > código<br>
>> > mínimo reproduzível.<br>
>> _______________________________________________<br>
>> R-br mailing list<br>
>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</div></div></blockquote></div><br>