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** *