ageDistribution.R

Nathan — Apr 7, 2014, 12:21 PM

setwd("C:/Users/Nathan/Desktop/Project/")
# We know that there are 1000 policyholders
# 20% are between ages 30-40
# 30% are between ages 40-50
# 50$ are between ages 50-60
# Assume that as age increases, number of policyholders at that age increases
# So amount_Policyholders is a function of age, and weighted by 1000
# amount_Policyholders must integrate to 1 to 30 to 60

#group1 = floor(sort(runif(2000,30,40))) + 1
#group2 = floor(sort(runif(3000,40,50))) + 1
#group3 = floor(sort(runif(5000,50,60))) + 1
group1 = floor(sort(runif(2000,30,40)))
group2 = floor(sort(runif(3000,40,51)))
group3 = floor(sort(runif(5000,50,60))) + 1

group = c(group1,group2,group3)
distribution = as.data.frame(table(group))

# exoort a column of values
count = distribution[,2]
write.csv(count, "ageAmounts.csv",row.names=FALSE)