#Így készítettem a kuiper.xls beolvasását. Kimentettem csv-be, majd a vesszőt pontra
# cseréltem szövegszerkesztővel, aztán futtattam a lenti scriptet

teszt<-read.table("kuipermod.csv",header=TRUE,sep=";",colClasses=c("numeric","numeric", rep("character",4),"numeric","numeric","numeric",rep("character",15)))
head(teszt)
str(teszt)
kruiser<-teszt[1:12]
head(kruiser)
save(kruiser,file="kruiser.RData")

# Az órai anyag innen kezdődik

load("kuiper.RData")
head(kruiser)

str(kruiser)

modellunkkezdet<-with(kruiser,lm(Price ~ Mileage + Cylinder + Sound + Leather))

modellunkkezdet

summary(modellunkkezdet)

fitted(modellunkkezdet)


str(kruiser)

kruiser$Cylinder<-as.factor(kruiser$Cylinder)

str(kruiser)

kruiser$Cylinder

modellunkkezdet2<-with(kruiser,lm(Price ~ Mileage + Cylinder + Sound + Leather))

summary(modellunkkezdet2)

unique(kruiser$Cylinder)

kruiserepit<-kruiser[1:402,]

kruiserteszt<-kruiser[403:804,]

modellunk<-with(kruiserepit,lm(Price ~ Mileage + Cylinder + Sound + Leather))

kruiserteszt$Priceb<-predict(modellunk,kruiserteszt)

View(kruiserteszt[-c(2:12)])

dim(kruiser)

plot(kruiserteszt[[1]],kruiserteszt[[13]])

sapply(kruiserteszt,mean,na.rm=TRUE)

by(kruiserteszt[,1:2], kruiserteszt$Make, summary)

library(rgl)

plot3d(kruiserteszt[[1]],kruiserteszt[[2]],kruiserteszt[[13]], col="red", size=3)