Re: [R] PAN: Need Help for Multiple Imputation Package

From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>
Date: Fri 27 May 2005 - 00:47:17 EST

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@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


R-help@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 Received on Fri May 27 01:03:45 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:07 EST