I was finally able to run linear mixed models with my data - thank you very much Martin for your precious help!
No I would need some help/advice for visualization and description of my results
The command fs_write_Y(sided_pval,mri,'spval.mgh') is not working on my matlab. This is the message error: Error using reshape To RESHAPE the number of elements must not change. Error in fs_write_Y (line 21) save_mgh(reshape(Y',mri.volsz),fname,mri.M,mri.mr_parms);
so, I have used the command: fs_write_fstats(F_rhstats_lin,mri,'sig.mgh','sig') to export the findings, as described in the wiki. However, how can I describe my data, in terms of measure of the clusters, MNI coordinates and statistical values? and is it possible to smooth the image?
In addition, this map is 'uncorrected' for multiple testing. How can I perform this correction?
I am testing null hypothesis of no group differences in the rate of change over time between two groups and I found significant differences both using quadratic effects (time per time per group) and linear effects (time per group) in the same model (with one random effect, that appeared to perform better than the model with 2 RF). The map of the first analysis appeared in blue in Freeview, whereas the map for linear effects appeared in red. Does sig.mgh include the sign of the difference? and is it possible that the two analyses have opposite sign? am I doing something wrong?
Thank you!
Angela
Martin Reuter mreuter@nmr.mgh.harvard.edu ha scritto:
Actually I am just adding a fix for this submitted by a user (thanks Knut) , so it should work with modern versions soon.
Best, Martin
On 04/06/2016 11:02 AM, Martin Reuter wrote:
Hi Angela,
newer versions of matlab don't support the old matlabpool comands any longer. You need to use an older matlab version (less than 8.2.0.29).
Best, Martin
On 04/05/2016 06:06 PM, Angela Favaro wrote:
Hi Martin, thank you for your help. I still have problems in running mixed models analyses. Do I need some additional toolbox in matlab? (I did not find this information in the web) I tried to run the module in two versions of matlab, but the same problem occurr.
I used these commands: [M,Y,ni] = sortData(M,1,Y,sID) X = [ones(length(M),1) M M(:,1).*M(:,2)]
[lhTh0,lhRe] = lme_mass_fit_EMinit(X,[1 2],Y,ni,lhcortex,3)
this is the error message: Undefined function 'matlabpool' for input arguments of type 'char'.
Error in lme_mass_fit_init (line 74) if (prs==1) || (matlabpool('size') ~= prs)
Error in lme_mass_fit_EMinit (line 69) [Theta1,Re1] = lme_mass_fit_init(X,Zcols,Y,ni,[],prs);
What can it be the problem? the matrices seem to be OK (and M has the same number of rows as Y (and as the sID)).
thank you for any help!
Angela
Martin Reuter mreuter@nmr.mgh.harvard.edu ha scritto:
Hi Angela,
the smooth function seems to be part of the curve fitting toolbox http://www.mathworks.com/help/curvefit/smooth.html so either that is missing, or your version is too old. You can skip the lowess Plot (it is just for initially checking the data).
Make sure your M has the same number of rows as your Y (and as the sID). Also make sure you are computing the ni (and sorted M and Y) like this (if your first column in M is the time, else specify the column with the time as the second parameter):
[M,Y,ni] = sortData(M,1,Y,sID); # (sorts the data)
Best, Martin
On 03/20/2016 04:38 PM, Angela Favaro wrote:
Dear Freesurfer experts, I need some help with the longitudinal streaming and the use of mixed models in matlab. I followed all the instructions and tried to use the following commands in matlab:
lme_timePlot(M(:,1),Y,ni) and this worked well
However, when I tried this: lme_lowessPlot(M(:,1),Y,0.70,M(:,2))
matlab gave the following error message:
Undefined function 'smooth' for input arguments of type 'double'. Error in lme_lowessPlot (line 49) ybw = smooth(st,dat(ix),bw,'rlowess')
I checked in my version of matlab (R2012b) and the command 'smooth' does not exist
Further, when I tried this:
lhstats = lme_fit_FS(X,[1 2],Y,ni)
Error using - Matrix dimensions must agree.
Error in lme_fit_init (line 42) ri = yi-Xi*Bhat;
Error in lme_fit_FS (line 71) [D,phisq] = lme_fit_init(X,Zcols,y,ni);
Looking into the matlab code, I observed that in line 42 there is no command..... and I am sure that matrix dimensions are correct.
So, I have problems in solving these issues... Any help is greatly appreciated! Is it a problem of Matlab version? Thank you!
Angela
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
-- Martin Reuter, PhD Assistant Professor of Radiology, Harvard Medical School Assistant Professor of Neurology, Harvard Medical School A.A.Martinos Center for Biomedical Imaging Massachusetts General Hospital Research Affiliate, CSAIL, MIT Phone: +1-617-724-5652 Web : http://reuter.mit.edu
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
-- Martin Reuter, PhD Assistant Professor of Radiology, Harvard Medical School Assistant Professor of Neurology, Harvard Medical School A.A.Martinos Center for Biomedical Imaging Massachusetts General Hospital Research Affiliate, CSAIL, MIT Phone: +1-617-724-5652 Web : http://reuter.mit.edu
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer