[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