Statistik-Laborstepfunts1par(mfrow=c(1,2))
plot(c(0,Anzahl),rep(0,2),type="n",main = "Histogramm von Saldo")
plot(c(0,Anzahl),rep(0,2),type="n",main = "Histogramm von Serienlänge")
Saldo<-rep(0,AnzahlLaeufe)
Laenge<-rep(0,AnzahlLaeufe)
pair<-seq(2,36,by=2)
impair<-seq(1,35,by=2)
rouge<-c(1,3,5,7,9,12,14,16,18,19,21,23,25,27,30,32,34,35)
noir<-c(2,4,6,8,10,11,13,15,17,20,22,24,26,28,29,31,33,35)
manque<-c(1:18)
passe<-c(19:36)
# Progressionsmatrix
gmat<-matrix(c(2,3,8,20,2,3,8,20,0,4,12,30,0,4,12,30,0,6,16,40,0,6,16,40),4,6)
# Initialisierung des Setzens
a<-Coup
# Zuordnen der gesetzten Chance
{
if(a=="pair") { satz<-pair }
else if (a=="impair") {satz<-impair }
else if (a=="rouge") {satz<-rouge}
else if (a=="noir") {satz<-noir }
else if (a=="manque") {satz<-manque}
else if (a=="passe") {satz<-passe}
else stop(message="Nur einfache Chancen sind zugelassen")
}
for (wieder in 1:AnzahlLaeufe)
{
satzhoehe <- Einsatz*2 # Begonnen wird mit 2 Satzeinheiten
it<-0
saldo<-0
hstufe<-1
ustufe<-1
while(satzhoehe>0 & satzhoehe<=Limit & ((Kapital+saldo)>=satzhoehe) & it<Anzahl)
{
it<-it+1
coup<-floor(runif(1,min=0,max=37)) # Werfen der Kugel
gewinn<-0
if ( coup%in%satz ) # Progression bei Gewinn
{
gewinn <- 1
saldo <- saldo+satzhoehe
{
if ((hstufe == 1) & (ustufe == 1)){ ustufe <- 2}
if ((hstufe == 1) & (ustufe == 2)){ hstufe <- 2 ; ustufe <- 1}
else if ((hstufe > 1) & (hstufe <= 4) & (ustufe < 6)){ ustufe <- ustufe+1 }
else if ((hstufe > 1) & (hstufe < 4) & (ustufe == 6)) { hstufe <- hstufe+1; ustufe<-1 }
else { hstufe<-1 ; ustufe <- 1 }
}
satzhoehe <- Einsatz*gmat[hstufe,ustufe]
}
if ( coup==0 ) # Bei Zero: 1/2 Einsatz als Verlust
{
saldo <- saldo-0.5*satzhoehe
}
if (( coup!=0 ) & ( gewinn==0 ))
{
saldo <- saldo-satzhoehe
}
if ( gewinn==0 ) # Progression bei Verlust (egal ob halb oder ganz)
{
if (ustufe %in% c(1,3,5))
{
ustufe<-1;
hstufe <- max(1,hstufe-1)
satzhoehe <- gmat[hstufe,ustufe]
}
else
{
ustufe<-ustufe-1
satzhoehe <- Einsatz*gmat[hstufe,ustufe]
}
}
}
Saldo[wieder]<-saldo
Laenge[wieder]<-it
}
hist(Saldo,breaks=floor(10*log10(AnzahlLaeufe)),freq=T,main = "Histogramm von Saldo",border="blue",density = 20, angle = 45, col = "blue")
hist(Laenge,breaks=floor(10*log10(AnzahlLaeufe)),freq=T,main = "Histogramm von Serienlänge",border="blue",density = 20, angle = 45, col = "blue")
1{\rtf1\ansi\ansicpg1252\deff0\deflang1031{\fonttbl{\f0\fswiss\fprq2\fcharset0 Verdana;}}
{\colortbl ;\red0\green0\blue0;}
\viewkind4\uc1\pard\cf1\f0\fs24 Mittlerer Gewinn : @(print(mean(Saldo)))
\par Standardabweichung Gewinn: @(print(sd(Saldo)))
\par Mittlerer Serienl\'e4nge : @(print(mean(Laenge)))
\par Standardabweichung Serienl\'e4nge: @(print(sd(Laenge)))
\par }
1Datensatz(D11)1"pair"10000101pair1000010100100200