[Reading Notes]R in Action

楼主辛苦

1. Listing 8.2. Polynomial regression
   
2. 
   
3. > fit2 <- lm(weight ~ height + I(height^2), data=women)
   
4. > summary(fit2)
   
5. > plot(women$height,women$weight,
   
6.        xlab="Height (in inches)",
   
7.        ylab="Weight (in lbs)")
   
8. > lines(women$height,fitted(fit2))
   
9. library(car)
   
10. scatterplot(weight ~ height, data=women,
    
11.   spread=FALSE, smoother.args=list(lty=2), pch=19,
    
12.   main="Women Age 30-39",
    
13.   xlab="Height (inches)",
    
14.   ylab="Weight (lbs.)")

复制代码