> plasma <- read.table("F:\\sta4210\\CH03TA08.txt",header=F,col.names=c("X","Y","log10Y"))
> 
> attach(plasma)
The following object(s) are masked from 'plasma (position 3)':

    log10Y, X, Y
> 
> pdf("C:\\Rmisc\\graphs\\CH03TA08.pdf")
> 
> plot(X,Y)
> 
> plot(X,log10Y)
> 
> plasma.reg1 <- lm(log10Y ~ X)
> summary(plasma.reg1)

Call:
lm(formula = log10Y ~ X)

Residuals:
      Min        1Q    Median        3Q       Max 
-0.116974 -0.035356  0.002626  0.027908  0.168180 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  1.134820   0.021641   52.44  < 2e-16 ***
X           -0.102282   0.008835  -11.58 4.51e-11 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 0.06247 on 23 degrees of freedom
Multiple R-squared: 0.8535,     Adjusted R-squared: 0.8472 
F-statistic:   134 on 1 and 23 DF,  p-value: 4.506e-11 

> 
> Res_log10Y <- 100*resid(plasma.reg1)
> 
> plot(X,Res_log10Y,main="Residuals vs X")
> 
> qqnorm(Res_log10Y,main="Residuals Normal Probability Plot")
> qqline(Res_log10Y)
>