# Jedna sestistenna kostka
kostka1<-sample(1:6,200,replace=TRUE)


# Histrogram
hist(kostka1,breaks=c(0.5,1.5,2.5,3.5,4.5,5.5,6.5))
hist(kostka1,breaks=((0:6)+0.5))


# Soucet na sedmi kostkach
kostky1az7<-sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)+
            sample(1:6,200,replace=TRUE)

  
# Histogram
hist(kostky1az7,breaks=12)


# Test normality
ks.test(kostky1az7,"pnorm",mean(kostky1az7),sd(kostky1az7))


# Hod dvanactistennou kostkou
kostka2<-sample(1:12,100,replace=TRUE)


# Vypocet prumeru
prumer<-c()
for (i in 1:length(kostka2)) {
prumer[i]<-mean(kostka2[1:i])  
}


# Graf konvergence
plot(1:length(kostka2),prumer,cex=0.2,xlab="Pocet hodu",ylab="Prumer")
lines(c(1,length(kostka2)),c(6.5,6.5),lwd=3,col="blue")
lines(1:length(kostka2),prumer,lwd=3,col="darkred")