Displaying 1 result from an estimated 1 matches for "fnonrechypt".
2011 Apr 13
0
ddply and nlminb
...oth the parameter estimates and the parameter errors.
Thank you
Robert
# ------------ END OF CODE
------------------------------------------------------------------------------------------
# ------------ Data loaded in data frame M2 and attached - not shown.
# ------------ Define function
fnonRecHypT <- function(x){
qeo = x[1]
dqe = x[2]
Fmo = x[3]
TL = x[4]
To = x[5]
TH = 45
theta = x[6]
qe = qeo + dqe*T
theta =
phi = (TH-To)/(To-TL)
Fm = Fmo*((T-TL)*(TH-T)^ phi) / ((To-TL)*(TH-To)^ phi)
Aest = (qe*PAR + Fm - ((qe*PAR + Fm)^2 - 4*theta*qe*PAR*Fm)^0.5)/(2.*theta)
result = sum...