dvar <- read.csv("http://commres.net/wiki/_media/elemapi2_.csv", fileEncoding="UTF-8-BOM") mod <- lm(api00 ~ ell + acs_k3 + avg_ed + meals, data=dvar) summary(mod) anova(mod) attach(dvar) da1 <- data.frame(api00, ell, acs_k3, avg_ed, meals) da2 <- data.frame(api00, ell, avg_ed, meals) da1 <- na.omit(da1) da2 <- na.omit(da2) spcor(da1) spcor(da2)