[R-br] predict com minpack.lm
Gilson Sanchez
gilsonsch em gmail.com
Quinta Janeiro 26 11:29:17 BRST 2012
Pessoal mais uma vez na procura de ajuda...
Estou utilizando o pacote minpack.lm para ajustar uma equação muito grande
que com o nls não consegui fazer, agora estou querendo ver os valores
ajustados pelo modelo, alguém pode ajudar...
segue o CMR o qual estou trabalhando
luz=c(0,25,50,75,100,500,1000,1500,2000)
Pn=c(-0.60375,
0.04300,0.27000,0.90175,1.39225,2.70350,4.93100,4.31000, 3.80075)
ps<-list(a=1,b=1)
model1=function(ps,xx){(((0.01940*xx+ps$a+0.60375)-(((0.01940*xx+ps$a+0.60375)^2)-4*0.01940*xx*ps$b*(ps$a+0.60375))^0.5)/2*ps$b)-0.60375}
resid1=function(observed,ps,xx){observed-model1(ps,xx)}
lin1=nls.lm(ps,resid1,observed=Pn, xx=PAR)
summary(lin1)
predict????
Agradeço desde já!
--
Bach. Agr.Gilson Sánchez Chia, MSc.
Laboratório de Fisiologia Vegetal
Embrapa Amazônia Ocidental
Fone: (92) 3303-7841
http://lattes.cnpq.br/3684553005529634
---
Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente.
Please consider the environment before printing this email.
Renew, Reduce and Recycle. The Planet Knows.***
Salva un Árbol.* *
***
*
PIENSA**** EN EL MUNDO PARA TUS HIJOS** *
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120126/964088d6/attachment.html>
Mais detalhes sobre a lista de discussão R-br