par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
par5 <- as.numeric(par5)
if (par6 == 'TRUE') par6 <- TRUE
if (par6 == 'FALSE') par6 <- FALSE
reps = par1
policies = par2
prem = par3
cost = par4
prob = par5
Poisson = par6
gross.profit <- 1:policies
gross.profit.rate <- 1:reps
if (Poisson==T) {
for (j in 1:reps) {
totalcost=cost*rpois(policies,prob)
gross.profit=prem-totalcost
gross.profit.rate[j]=sum(gross.profit)/policies
}
}
if (Poisson==F){
for (j in 1:reps) {
totalcost=cost*rbinom(policies,1,prob)
gross.profit = prem-totalcost
gross.profit.rate[j] = sum(gross.profit)/policies
}
}
risk <- 100 * (length (gross.profit.rate[gross.profit.rate < 0]))/reps
print(risk)
avgprofit <- format(mean(gross.profit.rate), digits = 5)
sdprofit <- format(sd(gross.profit.rate), digits = 5)
bitmap(file='test1.png')
stripchart(gross.profit.rate, method = 'stack', main = main, xlab = xlab, ylab = ylab, col = 'blue',at=0.1)
stripchart(gross.profit.rate[gross.profit.rate<0], method = 'stack', col = 'red', add = TRUE, at=0.1)
dev.off()
x<-rnorm(100)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Summary Statistics',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Number of policies',header=TRUE)
a<-table.element(a,policies)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Number of simulations',header=TRUE)
a<-table.element(a,reps)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Mean of gross profit per policy in USD',header=TRUE)
a<-table.element(a,avgprofit)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Risk of loss in percent',header=TRUE)
a<-table.element(a, round(risk,0) )
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'SD of gross profit per policy',header=TRUE)
a<-table.element(a,sdprofit)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
|
