Ashley, its better to reparameterize the variable than to use this
option, if its possible. you can use the same tricks we used for
reparameterizing parameters. so y>0 can be fed into amelia as ln(y),
or 0<y<1 can be fed in as ln(y/(1-y)), etc. and then
un-reparameterize when it comes out.
Gary
---
On Sun, Apr 25, 2010 at 1:20 PM, Ashley Anderson
<aaanders at fas.harvard.edu> wrote:
I figured out the first question but now I'm
getting an error with my
bounding matrix. This is what I have:
row1 <- c(5,0,100000)
row2 <- c(6,0,100000)
row3 <- c(7,0,100000)
row4 <- c(8,0,100000)
row5 <- c(9,0,100000)
grr <- matrix(rbind(row1, row2, row3, row4, row5), ncol=3, nrow=5)
bds <- grr
a.out <- amelia(data, m=5, ts="year", cs="cty", idvars =
c("id", "id2"),
bounds=bds, max.resample = 100000)
I was pretty sure that this was how you set up a bounds matrix (it tells
amelia that columns 5,6,7,8, and 9 should be bounded from 0 to 100000), but
I keep getting this error:
Error in matrix(NA, nrow = AMn.ss, ncol = AMp) :
??non-numeric matrix extent
Does anyone know what I could be doing wrong?
-Ashley Anderson
_______________________________________________
gov2001-l mailing list
gov2001-l at
lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l