dat <- read.table("galton_reg.dat")
m1 <- apply(dat,1,sum)
dat1 <- cbind(dat,m1)
m2 <- apply(dat1,2,sum)
dat12 <- rbind(dat1,m2)
colnames(dat12)<-c("<61.7",seq(62.2,73.2,by=1),">73.7","Total")
rownames(dat12)<-c(">73.0",seq(72.5,64.5,by=-1),"<64.0","Total")
print(dat12)
d <- scan("regression.dat")
child.factor <- gl(14,1,154)
parent.factor <- gl(11,14,154)
child.height <- 60.2 + as.numeric(child.factor)
parent.height <- 74.5 - as.numeric(parent.factor)
child <- rep(child.height,d)
parent <- rep(parent.height,d)
child.parent <- lm(child ~ parent)
parent.child <- lm(parent ~ child)
plot(child.height,parent.height,type="n")
text(child.height,parent.height,d)
abline(child.parent)
abline(parent.child)
print(summary(child.parent))
print(summary(parent.child))