Displaying 1 result from an estimated 1 matches for "uppval".
Did you mean:
uppal
2011 Apr 13
0
ddply and nlminb
...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((Aest-A)^2 )
}
# ------------ Define parameter starting values and limits
startval = c(0.05,-0.01,20, -5,15,0.5)
lowval = c(0.01,-0.05, 5,-15, 7,0.1)
uppval = c(0.2, 0.02,50, 0,25,0.99)
# ------------ Fit using entire data set
R3 <- nlminb(startval, fnonRecHypT, control=list(trace=1),
lower=lowval, upper = uppval)
# ------------ estimate fit parameter standard errors
x = R3$par
D3 = hessian(fnonRecHypT,x)
e3 = sqrt(diag(solve(D3)))
# -----...