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'#
|