Hi,

I am using Amelia to simulate missing vote values for French elections (before performing multiparty electoral data analysis using Clarify).

I need to make sure that after the simulation the sum of the votes for the different parties (vFN + vPC + vPS + vUMP + vVerts) is below one (these are the main parties and some much smaller parties are not included).

I thought the prior I generated (see below) would make the trick but it does not seem to work.

And the thing is that in the imputed data I generate now, the sum of the votes is very often above 1 which makes no sense (and generate issue then on Clarify with the logistic transformation).

Any idea of how I could handle that?

Many thanks in advance,

Best,

Julia

Here is my code:

database <- read.dta13("rall.dta")

prior <- matrix(NA, nrow=nrow(database),ncol=5)
for (i in 1:nrow(database)){
  v3 <- database$vFN[i]
  v5 <- database$vPC[i]
  v7 <- database$vPS[i]
  v9 <- database$vUMP[i]
  v11 <- database$vVerts[i]
 
  prior[i,] <-  c(i, 3, 0, 1 - v5+v7+v9+v11, 0.999999)
}

prior <- prior[!is.na(prior[,4]),]

a.out <- amelia(database, m = 5,ts="year", cs = "district",priors=prior,lgstc=c("vFN","vPC","vPS","vUMP","vVerts"),bound=rbind(c(4,0,Inf),c(6,0,Inf),c(8,0,Inf),c(10,0,Inf),c(12,0,Inf)))
write.amelia(obj=a.out, file.stem = "R19932012/outdata", format = "dta")

--
Amelia mailing list served by HUIT
[Un]Subscribe/View Archive: http://lists.gking.harvard.edu/?info=amelia
More info about Amelia: http://gking.harvard.edu/amelia
Amelia mailing list
Amelia@lists.gking.harvard.edu

To unsubscribe from this list or get other information:

https://lists.gking.harvard.edu/mailman/listinfo/amelia