Hi Emma,
yes, LME is Matlab, but the wiki page describe the use step-by-step.
LME makes sense in your case because some subjects have less time points than others.
The alternative (if only a very small minority of subjects differs from the rest) would be the 2-stage approach, where you first reduce the longitudinal data to a single estimate of change (the slope of a linear fit within-subject) and then compare these slopes across groups. However, - this is not a good idea in general as it does not consider the fact that some subjects have less time points than others, which is especially problematic, if there is a bias across groups (if you did a group analysis) - also, while QDEC could be used to compare this slope/atrophy measure across groups, it currently does not allow 'one sample group mean' to check if this slope is different from zero. You'd have to use mri_glmfit for that.
You can upgrade to 5.3 to use a more recent freeview and to run your post-processing (LME etc), but don't mix versions for the image processing part (recon-all).
And finally about the design in general. Not having a control group can be a big problem. For example just looking at atrophy rates in your treatment group will not tell you anything. You'll probably find some atrophy here and there, but is it different from no-treatment? Is it different from normal ageing? You would not be able to tell. But maybe your design is different, e.g. you could look for a correlation with drug dose etc.
I hope that helps. Best, Martin
On 09/23/2014 05:55 PM, Emma Thompson wrote:
Hi FS, I want to conduct a within-subjects longitudinal analysis, I thought I would be able to run a LME model in QDEC but this doesn't seem to be the case at all, is this correct? It seems I have to do everything using Matlab (big sigh).
https://surfer.nmr.mgh.harvard.edu/fswiki/LinearMixedEffectsModels
I have already preprocessed the data: 1) cross for all time points, 2) base, and then 3) long according to the awesome tutorial provided by FS.
I have only one group, a patient population that was scanned prior to (bseline) and at various time points (3 time points post) following treatment. Would you agree that the best analytical approach for me would be the LME model, especially since I have a few subjects with a couple of missing post-treatment time points?
Is there any way you would recommend another approach using QDEC?
Lastly, I have done everything in FS version 5.1, I'm considering upgrading to 5.3, since my freeview software seems to be out of date. It also seems that the Matlab tools I need to run LME are only available in FS version 5.2. At this point a little concerned since I'm hoping I didn't waste a bunch of time just realizing all this now, would upgrading to 5.3 negatively impact all the work I've done thus far using version 5.1?
Thanks for the help!!
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer