Displaying 1 result from an estimated 1 matches for "uselog".
Did you mean:
userlog
2009 Mar 15
0
Axes crossing at origin
...x(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=normpdf(y1,17.1,0.045)
norm1=min(newy)+norm1/max(no...