search for: predict_amount

Displaying 1 result from an estimated 1 matches for "predict_amount".

2009 Mar 15
0
Axes crossing at origin
...predict(fit,newdata) plot.new() plot.window(xlim=c(0,0.6*max(fit$df$x)),ylim=c(15.7,17.5)) lines(newdata$x, newy) axis(1) axis(2) abline(v=0, h=0) mtext("FITC-insulin [mol]",1,2.4) mtext("log intensity [log(LAU)]",2,2.4) myydata=c(16.35,16.5,16.65,16.95,17.1,17.25) source("predict_amount.r") myxdata=predict_amount(fit,myydata,uselog=TRUE) for(i in 1:length(myydata)){ lines(c(myxdata[i],myxdata[i],max(newdata)),c(min(newy),myydata[i],myydata[i]),lty=3) } y1=seq(16.95,17.25,0.003) y2=seq(16.35,16.65,0.003) x1=predict_amount(fit,y1,TRUE) x2=predict_amount(fit,y2,TRUE) norm1=...