Displaying 1 result from an estimated 1 matches for "waterpot".
Did you mean:
waterport
2010 Mar 29
1
getting CI's for certain y of nls fitted curve
...ata.frame(press =
press_x)),"gradient"),1,
function(x) sum(vcov(mod1)*outer(x,x))))
###plot points and curve with ci's
matplot(press_x, predict(mod1, data.frame(press =
press_x))+outer(se.fit,qnorm(c(.5, .025,.975))),
type="l",lty=c(1,3,3),ylab="transp",xlab="waterpot")
with(por, matpoints(press,tr,pch=1))
###the desired y is the tenth of the y asymptote
asym<-coef(mod1)[1]
tenth<-asym*0.1
abline(h=asym)
abline(h=tenth,lty=5)
text(10,65,"10% of max. transpiration",adj=0)
###finding corresponding x by y with inverse function:
Asym=coef(m...