28089

example <- read.table("dane_kosiarki.txt",header=TRUE) example

s.wl <- example$Y==l s.n.wl <- example$Y==2 wl <- example[s.wl,] n.wl <- example(s.n.wl,]

plot(wl$Xl,wl$X2,xlim=c(20,120),ylim=c(10,25), xlab="Dochi»d w tysiącach dolarów",

ylab="Wielko?6 dziadki w tysifcach metr&w kwadratowych") points(n.wl$Xl,n.wl$X2,pch=19)

wl.x <- data.framę(wl$Xl,wl$X2) n.wl.x <• data.frame(n.wl$Xl,n.wl$X2)

ml <- as.vector(mean(wl.x)) m2 <- as.vector(mean(n.wl.x)) sl <- cov(wl.x) s2 <- cov(n.wl.x)

n<-nrow(example)

nl<-nrow(wl)

n2<-nrow(n.wl)

W<-(l/(n-2))*((nl*1)*sl+(n2-l)*s2)

Wi<- solve(W)

a <- t(ml-m2)S*%Wi m <- t(ml-m2)V%WiVMml+m2)/2 plot(wl$Xl,wl$X2,xlim=c(20,120),ylim=c(10,25), xlab="Dochfid w tysiącach dolarów",

ylab="Wielko?6 dziadki w tysiącach metr&w kwadratowych")

points(n.wl$X1,n.wl$X2,pch=19)

text(110,22,"Rl",col="blue")

text(30,15,"R2",col="blue")

abline(m/a[2],-a[1]/a[2],lwd=2,col="blue")