[R-br] predict com minpack.lm

[Ivan Bezerra Allaman]

Com ?nls.lm você resolveria seu problema. No entanto, segue um exemplo:

f <-function(luz,a,b){
    (((0.01940*luz+a+0.60375)-(((0.01940*luz+a+0.60375)^2)-4*0.01940*luz*b*(a+0.60375))^0.5)/2*b)-0.60375
}
preditos <- do.call('f',c(list(luz=luz),lin1$par))

plot(luz,Pn,ylim=c(-1,5),type='l')
par(new=TRUE)
plot(luz,preditos,ylim=c(-1,5),col='red',xlab='',ylab='',type='l')


(S,f,P)
Allaman

 

\begin{signature}
<<>>=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5076
E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
@
\end{signature}
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120126/51ce8322/attachment.html>