Hi Alex
The cmd
[detvtx,sided_pval,pth] = lme_mass_FDR2(fstats.pval,fstats.sgn,lhcortex,0.05,0);
only takes a few seconds to finish. Please check if Y, Beta2, fstats.pval,fstats.sgn, sided_pval
all have the same number of columns.
Best -Jorge
De: Alex Hanganu al.hanganu@yahoo.ca Para: Jorge L. Bernal-Rusiel jbernal@nmr.mgh.harvard.edu; FS Mailing List Freesurfer@nmr.mgh.harvard.edu Enviado: Lunes 11 de marzo de 2013 10:53 Asunto: [Freesurfer] LME - lme_mass_FDR2 correction doesn't work
Hi Jorge,
I got an error in the correction, and I can't find a way out. this is the cmd I use:
[Y,mri] = fs_read_Y('.../thickness-sm10.mgh'); Qdec = fReadQdec('.../LME.dat'); Qdec = rmQdecCol(Qdec,1); sID = Qdec(2:end,1); Qdec = rmQdecCol(Qdec,1); M = Qdec2num(Qdec); ni = [2 2... 2]; X = [ones(length(M),1) M M(:,1).*M(:,2)]; [stats,st] =
lme_mass_fit_vw1(X,[1],Y,ni,lhcortex);
CM.C = [0 0 0 1]; fstats = lme_mass_F1(stats,CM); fs_write_fstats(fstats,mri,'.../lh.sig.mgh','sig'); nv = length(stats); Beta2 = zeros(1,nv); Beta4= zeros(1,nv); for .... . . end; Beta2_Beta4_sum = Beta2+Beta4; mri1 = mri; mri1.volsz(4) = 1; fs_write_Y(Beta2,mri1,'..../Beta2.mgh'); [detvtx,sided_pval,pth] = lme_mass_FDR2(fstats.pval,fstats.sgn,lhcortex,0.05,0); fs_write_Y(sided_pval,mri1,'.../spval.mgh');
??? Error using ===> reshape To RESHAPE the number of elements must not change.
Error in ==> fs_write_Y at 21 save_mgh(reshape(Y',mri.volsz),fname,mri.M,mri.mr_parms);
the cmd: [detvtx ....] didn't work, but there was no error message. I remember in the past this cmd took around 1 day to analyzebut now it didn't work at all.
fs_write_Y worked fine for Beta, alsothis time I used the mri1 = mri cmd- so the reshape error shouldn't be associated with this.
Please help !
best regards,
Alex.
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.