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