FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_babies.wasp
Title produced by softwareExercise 1.13
Date of computationSat, 11 Oct 2008 16:00:15 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Oct/12/t12237624768ozdgdk92i81kqw.htm/, Retrieved Sun, 19 May 2024 13:06:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=15372, Retrieved Sun, 19 May 2024 13:06:00 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact223

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Exercise 1.13] [Exercise 1.13 (Wo...] [2008-10-01 13:28:34] [b98453cac15ba1066b407e146608df68]
F   P     [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:00:15] [e08fee3874f3333d6b7a377a061b860d] [Current]
-   P       [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:02:21] [819b576fab25b35cfda70f80599828ec]
F   P         [Exercise 1.13] [Taak 1 Compuation...] [2008-10-11 22:03:56] [819b576fab25b35cfda70f80599828ec]
-   P           [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:05:14] [819b576fab25b35cfda70f80599828ec]
-                 [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:06:41] [819b576fab25b35cfda70f80599828ec]
-                   [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:08:00] [819b576fab25b35cfda70f80599828ec]
F                     [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:09:26] [819b576fab25b35cfda70f80599828ec]
-                       [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:10:48] [819b576fab25b35cfda70f80599828ec]
F                         [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:13:00] [819b576fab25b35cfda70f80599828ec]
-                           [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:18:42] [819b576fab25b35cfda70f80599828ec]
-                             [Exercise 1.13] [Taak 1 Computatio...] [2008-10-11 22:20:24] [819b576fab25b35cfda70f80599828ec]
-   P                           [Exercise 1.13] [Taak 1 Computatio...] [2008-10-13 11:01:01] [819b576fab25b35cfda70f80599828ec]
-   P                           [Exercise 1.13] [Taak 1 Computatio...] [2008-10-13 11:03:34] [819b576fab25b35cfda70f80599828ec]
-   P                           [Exercise 1.13] [Verbetering: 80%] [2008-10-18 15:10:04] [b85eb1eb4b13b870c6e7ebbba3e34fcc]
-   P                 [Exercise 1.13] [Taak 1 Computatio...] [2008-10-13 10:42:41] [a18c43c8b63fa6800a53bb187b9ddd45]
-   P                 [Exercise 1.13] [Taak 1 Computatio...] [2008-10-13 10:46:18] [819b576fab25b35cfda70f80599828ec]
Feedback Forum
2008-10-17 17:43:41 [Gregory Van Overmeiren] [reply] 
Je had in vraag 1b in je werk,enkel de parameter aantal dagen moeten veranderen. Hier heb je dat gedaan maar je had beter het maximum aantal dagen gepakt (3650 dagen). Zo wordt de berekening gedaan over een grotere tijdspanne en wordt je resultaat daardoor nauwkeuriger (=> wet van de grote getallen )
2008-10-18 15:05:56 [Ellen Smolders] [reply] 
De student heeft proefsgewijs telkens andere parameters veranderd, om te zien wat het resultaat was. Uiteindelijk beslist de student dat de parameter 'percentage mannelijke geboortes' de parameter is die de waarschijnlijkheidsgraad het hardst doet dalen. Dit is niet correct. Uit berekeningen blijkt dat de parameter 'aantal dagen' moet worden verhoogt, zodat er meer trekkingen plaatsvinden en men een stabieler resultaat krijgt.
Volgende link: http://www.freestatistics.org/blog/index.php?v=date/2008/Oct/18/t1224342342fwx6y5sswcwbk42.htm
bewijst dat met het verhogen van het aantal dagen, de waarschijnlijkheidsgraad het laagst is (door de wet van grote getallen).
2008-10-18 19:09:32 [Astrid Sniekers] [reply] 
Uitleg oplossing vraag 1b:
Opnieuw is de student niet duidelijk in zijn antwoord. Hij begrijpt ook de vraag niet. Hij verandert constant de opgave van de oefening. Hij heeft wel de tijdspanne vergroot, maar zegt niet expliciet dat hierdoor de resultaten nauwkeuriger worden.
2008-10-20 19:22:27 [Steven Symons] [reply] 
de student doet in mijn opinie teveel werk waaruit weinig voorkomt. Hij moet gewoon de accuraatheid verbeteren door de parameter 'aantal dagen' te veranderen. hij doet dit in dit voorbeeld met 730 dagen maar het maximum ligt op 3650, dan zou hij deze optie moeten kiezen. Want de regel is hoe groter het aantal dagen, hoe neuwkeuriger de oplossing zal zijn. Als je dit dan nog eens extra neuwkeurig wil uitzoeken kan je de berekening met 3650 dagen enkele malen herhalen en dan krijg je een maximum en minimum.

Post a new message

Dataset

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 4 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=15372&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=15372&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=15372&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001






Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.)
Number of simulated days730
Expected number of births in Large Hospital45
Expected number of births in Small Hospital15
Percentage of Male births per day(for which the probability is computed)0.6
#Females births in Large Hospital16459
#Males births in Large Hospital16391
#Female births in Small Hospital5393
#Male births in Small Hospital5557
Probability of more than 60 % of male births in Large Hospital0.0698630136986301
Probability of more than 60 % of male births in Small Hospital0.164383561643836
#Days per Year when more than 60 % of male births occur in Large Hospital25.5
#Days per Year when more than 60 % of male births occur in Small Hospital60

\begin{tabular}{lllllllll}
\hline
Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.) \tabularnewline
Number of simulated days & 730 \tabularnewline
Expected number of births in Large Hospital & 45 \tabularnewline
Expected number of births in Small Hospital & 15 \tabularnewline
Percentage of Male births per day(for which the probability is computed) & 0.6 \tabularnewline
#Females births in Large Hospital & 16459 \tabularnewline
#Males births in Large Hospital & 16391 \tabularnewline
#Female births in Small Hospital & 5393 \tabularnewline
#Male births in Small Hospital & 5557 \tabularnewline
Probability of more than 60 % of male births in Large Hospital & 0.0698630136986301 \tabularnewline
Probability of more than 60 % of male births in Small Hospital & 0.164383561643836 \tabularnewline
#Days per Year when more than 60 % of male births occur in Large Hospital & 25.5 \tabularnewline
#Days per Year when more than 60 % of male births occur in Small Hospital & 60 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=15372&T=1

[TABLE]
[ROW][C]Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.)[/C][/ROW]
[ROW][C]Number of simulated days[/C][C]730[/C][/ROW]
[ROW][C]Expected number of births in Large Hospital[/C][C]45[/C][/ROW]
[ROW][C]Expected number of births in Small Hospital[/C][C]15[/C][/ROW]
[ROW][C]Percentage of Male births per day(for which the probability is computed)[/C][C]0.6[/C][/ROW]
[ROW][C]#Females births in Large Hospital[/C][C]16459[/C][/ROW]
[ROW][C]#Males births in Large Hospital[/C][C]16391[/C][/ROW]
[ROW][C]#Female births in Small Hospital[/C][C]5393[/C][/ROW]
[ROW][C]#Male births in Small Hospital[/C][C]5557[/C][/ROW]
[ROW][C]Probability of more than 60 % of male births in Large Hospital[/C][C]0.0698630136986301[/C][/ROW]
[C]Probability of more than 60 % of male births in Small Hospital[/C][C]0.164383561643836[/C][/ROW]
[ROW][C]#Days per Year when more than 60 % of male births occur in Large Hospital[/C][C]25.5[/C][/ROW]
[C]#Days per Year when more than 60 % of male births occur in Small Hospital[/C][C]60[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=15372&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=15372&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.)
Number of simulated days730
Expected number of births in Large Hospital45
Expected number of births in Small Hospital15
Percentage of Male births per day(for which the probability is computed)0.6
#Females births in Large Hospital16459
#Males births in Large Hospital16391
#Female births in Small Hospital5393
#Male births in Small Hospital5557
Probability of more than 60 % of male births in Large Hospital0.0698630136986301
Probability of more than 60 % of male births in Small Hospital0.164383561643836
#Days per Year when more than 60 % of male births occur in Large Hospital25.5
#Days per Year when more than 60 % of male births occur in Small Hospital60


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 730 ; par2 = 45 ; par3 = 15 ; par4 = 0.6 ;
Parameters (R input):
par1 = 730 ; par2 = 45 ; par3 = 15 ; par4 = 0.6 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
numsuccessbig <- 0
numsuccesssmall <- 0
bighospital <- array(NA,dim=c(par1,par2))
smallhospital <- array(NA,dim=c(par1,par3))
bigprob <- array(NA,dim=par1)
smallprob <- array(NA,dim=par1)
for (i in 1:par1) {
bighospital[i,] <- sample(c('F','M'),par2,replace=TRUE)
if (as.matrix(table(bighospital[i,]))[2] > par4*par2) numsuccessbig = numsuccessbig + 1
bigprob[i] <- numsuccessbig/i
smallhospital[i,] <- sample(c('F','M'),par3,replace=TRUE)
if (as.matrix(table(smallhospital[i,]))[2] > par4*par3) numsuccesssmall = numsuccesssmall + 1
smallprob[i] <- numsuccesssmall/i
}
tbig <- as.matrix(table(bighospital))
tsmall <- as.matrix(table(smallhospital))
tbig
tsmall
numsuccessbig/par1
bigprob[par1]
numsuccesssmall/par1
smallprob[par1]
numsuccessbig/par1*365
bigprob[par1]*365
numsuccesssmall/par1*365
smallprob[par1]*365
bitmap(file='test1.png')
plot(bigprob,col=2,main='Probability in Large Hospital',xlab='#simulated days',ylab='probability')
dev.off()
bitmap(file='test2.png')
plot(smallprob,col=2,main='Probability in Small Hospital',xlab='#simulated days',ylab='probability')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Number of simulated days',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Expected number of births in Large Hospital',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Expected number of births in Small Hospital',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Percentage of Male births per day
(for which the probability is computed)',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#Females births in Large Hospital',header=TRUE)
a<-table.element(a,tbig[1])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#Males births in Large Hospital',header=TRUE)
a<-table.element(a,tbig[2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#Female births in Small Hospital',header=TRUE)
a<-table.element(a,tsmall[1])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#Male births in Small Hospital',header=TRUE)
a<-table.element(a,tsmall[2])
a<-table.row.end(a)
a<-table.row.start(a)
dum1 <- paste('Probability of more than', par4*100, sep=' ')
dum <- paste(dum1, '% of male births in Large Hospital', sep=' ')
a<-table.element(a, dum, header=TRUE)
a<-table.element(a, bigprob[par1])
a<-table.row.end(a)
dum <- paste(dum1, '% of male births in Small Hospital', sep=' ')
a<-table.element(a, dum, header=TRUE)
a<-table.element(a, smallprob[par1])
a<-table.row.end(a)
a<-table.row.start(a)
dum1 <- paste('#Days per Year when more than', par4*100, sep=' ')
dum <- paste(dum1, '% of male births occur in Large Hospital', sep=' ')
a<-table.element(a, dum, header=TRUE)
a<-table.element(a, bigprob[par1]*365)
a<-table.row.end(a)
dum <- paste(dum1, '% of male births occur in Small Hospital', sep=' ')
a<-table.element(a, dum, header=TRUE)
a<-table.element(a, smallprob[par1]*365)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')