# First example: generate some raw data
data <- data.frame(
i=i <- rep(1, len = 100),
t=t<-sample(c("2008", "2012"), 100, replace=TRUE),
p=p <- sample(LETTERS[1:3], 100, replace=TRUE),
j=j <- sample(letters[1:13], 100, replace=TRUE) )
## IPC for each party nested within districts over time.
data$tcan = with(data, ave(i, t, p, j, FUN = sum))
get.ipc(data$tcan, mag = 2)
# Second example: generate some aggregate data
data2 <- data.frame(
a_i=i <- rep(1:5, c(1, 1, 4, 4, 3)),
t=t<-sample(c("2008", "2012"), 13, replace=TRUE),
p=p <- sample(LETTERS[1:3], 13, replace=TRUE),
j=j <- sample(letters[1:13], 13, replace=FALSE) )
## IPC for each party nested with districts and over time.
get.ipc(data2$a_i, mag = 1)
Run the code above in your browser using DataLab