[MINC-users] R/RMINC Linear Mixed-Effect Models with mincApply

[nicolas guizard]

Hello Minc-Users,

I am trying to use R with Rminc to compute a mixed-effect model on
longitudinal Jacobians of multiple subjects with multiple visits.

I didn't find a mincfunction in the Rminc library to do so, therefore, I
used another library (http://cran.r-project.org/web/packages/nlme/index.html
).
But I am not sure how to use the mincApply function to perform the analyses
with the different variables.

Here is what I do but doesn't work:

> f <- function(x) {
> lme(x ~  1 + gf$Visit + gf$Treatment, random = ~ gf$ID)
> }

> lmemodel  <- mincApply(gf$File, quote(f(x)), mask)

> Error in eval(expr, envir, enclos) : object 'ID' not found


Where gf is my data table.

Regards,
Nicolas