Message-ID: <4295E175.1010205@statistik.uni-dortmund.de>
Date: 2005-05-26T14:47:17Z
From: Uwe Ligges
Subject: PAN: Need Help for Multiple Imputation Package
In-Reply-To: <20050526002314.1372.qmail@web53709.mail.yahoo.com>
Please read the posting guide.
It tells you to provide (a) small reproducible examples and (b)
'If the question relates to a package that is downloaded from CRAN try
contacting the package maintainers first. You can also use
find("functionname") and packageDescription("packagename") to find this
information.'
Uwe Ligges
AC wrote:
> Hello all. I am trying to run PAN, multilevel
> multiple imputation program, in R to impute missing
> data in a longitudinal dataset. I could successfully
> run the multiple imputation when I only imputed one
> variable. However, when I tried to impute a
> time-varying covariate as well as a response variable,
> I received an error message, ?Error: subscript out of
> bounds.? Can anyone tell if my commands contain any
> mistakes?
>
> First I imported SAS dataset ?sim? which includes a
> response variable ?MIY1?, a time-varying covariate
> ?TCOV1?, TIME, GROUP (0 or 1), and ID. 200
> participants were included and measurement occurred
> six times. Approximately 25% of participants dropped
> out at end.
>
>
>>sim <- read.xport('c:\\xptds.dat')
>>
>>int <- rep(1,1200)
>>y <- cbind(sim$MIY1,sim$TCOV1)
>>subj <- sim$ID
>>pred <- cbind(int, sim$TIME, sim$GROUP)
>>
>>xcol <- 1:3
>>zcol <- 1
>
>
>>prior <- list(a=2,Binv=4,c=2,Dinv=4)
>
>
>>result <-
>
> pan(y,subj,pred,xcol,zcol,prior,seed=13579,iter=1000)
> Error: subscript out of bounds
>
>
> By the way, I also received the same error message
> when I tried to include intercept and time in Zcol, a
> matrix for random effect specification. I used
> command ? zcol <- 1:2?. Does anybody know if this
> error is due to sample size/proportion of missing data
> or due to command mistake?
>
> I truly appreciate any feedbacks.
> Best regards,
> Eishi
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html