bamp (version 2.1.3)

apcSimulate: Simulate from age-period-cohort model

Description

This functions simulates a data set of cases on the Lexis diagram from given age, period and cohort effects. Population numbers have to be given; can be one number for all age group/period combinations.

Usage

apcSimulate(intercept, age, period, cohort, periods_per_agegroup, population)

Value

List with number of cases (matrix) and population numbers (matrix).

Arguments

intercept

Intercept

age

Vector of effect for age groups

period

Vector of effects for periods

cohort

Vector of effect for cohorts

periods_per_agegroup

Periods per age group

population

Population number. Either a matrix or a scalar.

See Also

vignette("simulation", package = "bamp")

Examples

age=sqrt(seq(5,0,length=10)); age<-1-age-mean(age)
period=15:1; period[8:15]<-8:15; period<-period/6; period<-period-mean(period)
periods_per_agegroup=5; number_of_cohorts <- periods_per_agegroup*(10-1)+15
cohort<-rep(0,60); cohort[1:10]<-10:1; cohort[41:60]<- -(1:20)/2; cohort<-cohort/10;
cohort<-cohort-mean(cohort)
simdata<-apcSimulate(-5, age, period, cohort, periods_per_agegroup, 1e6)
par(mfrow=c(3,1))
plot(age, type="l")
plot(period, type="l")
plot(cohort, type="l")
if (FALSE) {
simmod <- bamp(cases = simdata$cases, population = simdata$population, age = "rw1", 
period = "rw1", cohort = "rw1", periods_per_agegroup =periods_per_agegroup)
plot(simmod)
}