# regrese rm(list=ls()) # vyčištění pracovního prostoru data(Melanom) # podá data summary(Melanom) # přehled proměnných plot(mort~lat,data=Melanom) plot(mort~long,data=Melanom) plot(mort~Ocean,data=Melanom) # plot(mort~lat,data=Melanom) summary(lm(mort~lat,data=Melanom)) abline(lm(mort~lat,data=Melanom)) # anova(lm(mort~lat,data=Melanom)) anova(lm(mort~1,data=Melanom)) # plot(mort~long,data=Melanom) abline(lm(mort~long,data=Melanom)) summary(lm(mort~long,data=Melanom)) # summary(lm(mort~lat+long,data=Melanom)) summary(lm(mort~long+lat,data=Melanom)) # ale pozor: summary(lm(mort~lat+long+I(long^2),data=Melanom)) # plot(mort~lat,data=Melanom,col=ifelse(Ocean==0,"red","blue")) abline(lm(mort~lat,data=Melanom,subset=Ocean==0),col="red") summary(lm(mort~lat,data=Melanom,subset=Ocean==0)) abline(lm(mort~lat,data=Melanom,subset=Ocean==1),col="blue") summary(lm(mort~lat,data=Melanom,subset=Ocean==1)) # summary(lm(mort~Ocean+lat,data=Melanom)) # # data(Policie) summary(Policie) plot(fat~height,data=Policie) abline(lm(fat~height,data=Policie)) summary(lm(fat~height,data=Policie)) summary(lm(fat~height+weight,data=Policie))