Home » date » 2009 » Oct » 06 »

*The author of this computation has been verified* R Software Module: rwasp_babies.wasp (opens new window with default values) Title produced by software: Exercise 1.13 Date of computation: Tue, 06 Oct 2009 14:01:06 -0600   Cite this page as follows: Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Oct/06/t1254859306bltokhx5xgfbgbo.htm/, Retrieved Sat, 25 May 2013 00:01:18 +0000   Original text written by user:   IsPrivate? No (this computation is public)   User-defined keywords: Vraag 2 test 1   System-generated keywords (parent): (pk = 0) Estimated Impact 36   Output produced by software: Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 3 seconds R Server 'George Udny Yule' @ 72.249.76.132 Exercise 1.13 p. 14 (Introduction to Probability, 2nd ed.) Number of simulated days 1825 Expected number of births in Large Hospital 45 Expected number of births in Small Hospital 15 Percentage of Male births per day (for which the probability is computed) 0.6 #Females births in Large Hospital 40990 #Males births in Large Hospital 41135 #Female births in Small Hospital 13599 #Male births in Small Hospital 13776 Probability of more than 60 % of male births in Large Hospital 0.0706849315068493 Probability of more than 60 % of male births in Small Hospital 0.158904109589041 #Days per Year when more than 60 % of male births occur in Large Hospital 25.8 #Days per Year when more than 60 % of male births occur in Small Hospital 58   Charts produced by software: http://www.freestatistics.org/blog/date/2009/Oct/06/t1254859306bltokhx5xgfbgbo/1suya1254859263.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Oct/06/t1254859306bltokhx5xgfbgbo/1suya1254859263.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Oct/06/t1254859306bltokhx5xgfbgbo/2q93d1254859263.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Oct/06/t1254859306bltokhx5xgfbgbo/2q93d1254859263.ps (opens in new window) Click here to open pdf file.   Parameters (Session): par1 = 1825 ; par2 = 45 ; par3 = 15 ; par4 = 0.6 ;   Parameters (R input): par1 = 1825 ; 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<br />(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')