Dear R-list

 

Suppose I have a data set stored in hmet, for which I did get confidence
limit on a linear regression as shown below.

My question is how I can subset only data points which are within the
confidence limit.

 

Thank you. 

Keun-Hyung 

 

---------------------------------------------------------------

Al=rnorm(100, 3)

Cd=rnorm(100, 0.2)

hmet=data.frame(Al=Al, Cd=Cd)

plot(Al, Cd, xlab="Al", ylab= "Cd", , pch=16)

cd.lm=lm(Cd~ Al, na.action=na.omit, data= hmet)

abline(cd.lm$coeff, lty=1)

xy= data.frame (Al =pretty(hmet$Al, 20))

yhat=predict(cd.lm, newdata=xy, interval="confidence")

ci=data.frame(lower=yhat[,"lwr"], upper=yhat[,"upr"])

lines(xy$Al, ci$lower, lty= 2, lwd=2, col="grey")

lines(xy$Al, ci$upper, lty=2, lwd=2, col="grey")


        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to