# Code von Franziska Wandtner, 2011, wandtner@math.uni-frankfurt.de

##########################################################################
# Wir ziehen mit Zurücklegen aus 10 roten und 20 blauen Kugeln 12 Mal.
# Wie oft wird eine blaue Kugeln gezogen?
# Das Ganze wird w Mal wiederholt.

##########################################################################
rm(list=ls())

############################
# Anzahl der Wiederholungen
w<-1000

###########################################################################
# Graphikparameter werden gesetzt     
# Für Windows
windows(width=9, height=6)
par(cex.main=1.3,cex.axis=1.0,cex.lab=1.0)
par(font.main=1,mar=c(2,2,5,2)) 

###########################################################################
plot(c(0,31),c(0,10),type="n",main="Ziehe 12 Kugeln mit Zurücklegen",ylab="",xlab="",axes=F)
points(c(1:10),rep(8,10),cex=3,col="red",pch=16)
points(c(11:30),rep(8,20),cex=3,col="blue",pch=16)

readline("Greife 12 Kugeln mit Zurücklegen raus")

s<-sample(30,12,replace=T)

for(i in 1:30)
{
	points(c(1:10),rep(8,10),cex=3,col="red",pch=16)
	points(c(11:30),rep(8,20),cex=3,col="blue",pch=16)
}


sn<-as.data.frame(table(s))
st<-as.matrix(sn)



b<-0
freq<-0
help<-0
ll<-length(st)/2

for(i in 1:12)
{
	if(s[i]<=10)	{
				points(s[i],5,cex=3,col="red",pch=16)
				points(s[i],8,cex=3,col="white",pch=16)
			   	}
			   else	{
			   		points(s[i],5,cex=3,col="blue",pch=16)
			   		points(s[i],8,cex=3,col="white",pch=16)
			   		b<-b+1
			   	   	}
			   	arrows(s[i]-0.2,7.5,s[i]-0.2,5.4,length=0.1)
			   	arrows(s[i]+0.2,5.4,s[i]+0.2,7.5,length=0.1)
			   	
			   	for(j in 1:ll)
			   	{	help<-as.numeric(st[j,1])
			   		if(s[i]==help){freq<-as.numeric(st[j,2])}
			   	}
			   	
			   	text(s[i],4,paste(freq))
			   	
}




text(15,2,paste("Anzahl blaue Kugeln:",b),cex=1.2,font=2)

###########################################################################
# Wiederhole w Mal

readline("Wiederhole w Mal")
par(font.main=1,mar=c(5,5,5,5)) 

smatrix<-matrix(rep(0,w*12),nrow=w)
b<-rep(0,w)

for(i in 1:w)
{
	smatrix[i,]<-sample(30,12)
	for(j in 1:12)
	{
		if(smatrix[i,j]>10){b[i]<-b[i]+1}}
			   	   	
}

breaks<-seq(0,12,by=0.1)

hist(b,breaks=breaks,axes=F,xlab="Ziehen mit Zurücklegen: Anzahl blaue Kugeln",ylab="Häufigkeit",freq=T,main=paste("Anzahl blaue Kugeln (",w, "Wiederholungen)"),col="black")
axis(2)
axis(1,at=c(-0.05,0.95,1.95,2.95,3.95,4.95,5.95,6.95,7.95,8.95,9.95,10.95,11.95),labels=c("0","1","2","3","4","5","6","7","8","9","10","11","12"))

###########################################################################
#Ende

readline("Ende")
graphics.off()