#Teht 1.

birthweight<-read.table("http://users.jyu.fi/~junyblom/birthweight.dat", header=TRUE)
attach(birthweight)

malli<-lm(weight ~ gestation + sex)
kert<-coef(malli)
# Jos raskauden kesto on viikon pidempi,
# niin ennustettu keskimääräinen paino on
# 120 g korkeampi
# Tytöt ovat keskimäärin 163 g kevyempiä
# Vakio ei tulkinnallinen suoraan

plot(gestation[sex==0],weight[sex==0],pch=".",xlab="kesto",ylab="paino",col=1,cex=6)
points(gestation[sex==1],weight[sex==1],pch=".",col=2,cex=6)
# Sex=0
abline(kert[1],kert[2])
# Sex=1
abline(kert[1]+kert[3],kert[2],col=2)

Teht 2.

peking=read.table("http://users.jyu.fi/~junyblom/peking.dat",header=TRUE)
attach(peking)
matka.log2=log2(matka)
aika.log2=log2(aika)
#peking.uusi=cbind(peking,matka.log2,aika.log2)
plot(matka.log2[sp==0],aika.log2[sp==0],col="blue",xlab="log2(matka)",ylab="log2(aika)")
points(matka.log2[sp==1],aika.log2[sp==1],col="red")

pojat=lm(aika.log2~matka.log2,subset=(sp==0),data=peking.uusi)
tytot=lm(aika.log2~matka.log2,subset=(sp==1),data=peking.uusi)
pojat
tytot  
#Koko data:log2(matka) kasvaa yhdellä=matka kaksinkertaistuu ---> log2(aika)
#kasvaa 1.12:lla = aika kasvaa
# 2^1.12=2.17 -kertaiseksi

#sp:t vastaavasti
abline(fit)