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 analyze
but now it
didn't work at
all.
fs_write_Y
worked fine
for Beta,
also
this
time I used
the mri1
=
mri cmd
-
so
the reshape
error
shouldn't be
associated
with this.
Please help !
best
regards,
Alex.
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
. 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.