Displaying 1 result from an estimated 1 matches for "pconfr".
2013 May 17
2
zigzag confidence interval in a plot
...7042
60,11 1535,1
70,04 2244,963303
79,97 1954,507042
100,31 2679,140625
> plot(cd$CHAOsep12RNA,cd$Depth, ylim = rev(range(0:100)), xlab="CHAO", ylab="Depth", pch=15, las=2, main="Sep12-RNA", cex.main=1)
> lmR <- lm(cd$Depth~cd$CHAOsep12RNA)
> abline(lmR)
pconfR <- predict(lmR,interval="confidence")
matlines(cd$CHAOsep12RNA,pconfR[,c("lwr","upr")], col=1, lty=2)
I also tried
> newx <- seq(min(cd$CHAOsep12RNA), max(cd$CHAOsep12RNA), length.out=11)
> a <- predict(lmR, newdata=data.frame(CHAO=newx), interval=c(&q...