


kockaseged<-runif(500,0.5,6.5)
kockaseged
kocka<-round(kockaseged)
kocka
hist(kocka)
table(kocka)
eredmeny<-chisq.test(table(kocka))
eredmeny

load("kolegyes.RData")

ab<-mean(subset(kolegyes,group==1)$cholest14,na.rm=TRUE)

szb<-sd(subset(kolegyes,group==1)$cholest14,na.rm=TRUE)

ks.test(subset(kolegyes,group==1)$cholest14,"pnorm",ab,szb)

ae<-mean(subset(kolegyes,group==2)$cholest14,na.rm=TRUE)

sze<-sd(subset(kolegyes,group==2)$cholest14,na.rm=TRUE)

ks.test(subset(kolegyes,group==2)$cholest14,"pnorm",ae,sze)








#Í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 itt folytatódik

load("kuiper.RData")
head(kruiser)


chisq.test(kruiser$Sound,kruiser$Leather)

str(kruiser)

?cor.test

with(kruiser, cor.test(Price,Mileage))

kruiser$dPrice<-kruiser$Price
kruiser$dMileage<-kruiser$Mileage

kruiser$dPrice[kruiser$Price<15000]<-1
kruiser$dPrice[kruiser$Price>=15000 & kruiser$dPrice<=30000]<-2
kruiser$dPrice[kruiser$Price>30000]<-3

kruiser$dMileage[kruiser$Mileage<10000]<-1
kruiser$dMileage[kruiser$Mileage>=10000 & kruiser$Mileage<=20000]<-2
kruiser$dMileage[kruiser$Mileage>20000]<-3

teszt2<-table(kruiser$dPrice, kruiser$dMileage)

chisq.test(teszt2)



#McNemar teszt

teszt<-read.table("mcnemar.csv",header=FALSE,sep=";")

teszt[,1]<-as.numeric(teszt[,1])
teszt[,2]<-as.numeric(teszt[,2])


str(teszt)
table(teszt)

mean(teszt[,1])
mean(teszt[,2])


mcnemar.test(table(teszt[,1],teszt[,2]))







