install.packages("Amelia")
update.packages()
require(Amelia)
summary(data)
a.out <- amelia(data, m=30, ts = "visit", cs = "id", boot.type = "none", noms=c("mar", "abst", "famalc1", "spousalc", "gender", "race", "income", "A", "D"),
logs=c("ASIalc", "ASIdrg", "ASIemp", "ASIleg", "ASImed", "ASIpsy", "ASIsoc", "ASIpdgt", "pcntbpov"),
bound = rbind(c(12, 0, Inf), c(13, 0, Inf), c(14, 0, Inf), c(15, 0, Inf), c(16, 0, Inf), c(17, 0, Inf), c(18, 0, Inf), c(19, 0, Inf), c(8, 0, Inf)),
polytime = 1, intercs = TRUE, p2s = 2, empri = .1*nrow(data))
Deysia Levin, MPH
Epidemiology Doctoral Candidate
Department of Epidemiology
UC Berkeley School of Public Health
101 Haviland Hall
Berkeley, CA. 94720
Phone: (510) 926-2496
Email: Deysia@gmail.com