#part1
data(galtonpeas)
attributes(galtonpeas)
plot(galtonpeas$Parent, galtonpeas$Progeny)

#part2
fit1 = lm(Progeny~Parent, galtonpeas,weights=1/(SD*SD))
summary(fit1)
plot(galtonpeas$Parent, galtonpeas$Progeny)
abline(fit1)
par(mfrow=c(2,2))
plot(fit1)



Progenynew = galtonpeas$Progeny-galtonpeas$Parent
galtonpeasnew = data.frame(cbind(Progenynew , galtonpeas$Parent,  galtonpeas$SD))
galtonpeas$SD
fit2 = lm(Progenynew ~+1, galtonpeasnew ,weights=1/(galtonpeas$SD*galtonpeas$SD))
summary(fit2)
anova(fit2)