[MINC-users] Multistat error - ??? Maximum variable size allowed by the program is exceeded.

[Penhune Lab]

Thank you for your input. I also sent this to Dr. Worsley and he replied
with basically the same answer that Marc has given here - which works just
fine :). This has been corrected in multistat by Dr. Worsley.

Cheers,
Chris

On Thu, Jun 5, 2008 at 9:00 AM, Marc BOUFFARD <bouffard at bic.mni.mcgill.ca>
wrote:

> Hi Chris,
>
> I have encountered this error before:
>
> try
>
> fwhm_varatio = inf
>
> instead of -100
>
> which gives a fixed effects analysis which seems to be what the error is
> telling you in this case.  I think (to be confirmed) that this is a bug
> whithin multistat.  It is related to this section of the code in
> multistat:
> *******************************************************************
>
>      if fwhm_varatio<0
>         % find fwhm to achieve target df:
>         df_target=-fwhm_varatio;
>         fwhm_varatio_limit=50;
>         ff=0:1:100;
>         dfs=[];
>         for f=ff;
>            df=regularized_df(f,D,Steps,numslices,df,fwhm_data);
>            dfs=[dfs df.t];
>         end
>         fwhm_varatio=interp1(dfs,ff,df_target);
>         if isnan(fwhm_varatio) | fwhm_varatio > fwhm_varatio_limit
>            fwhm_varatio=Inf;
>            varatio_vol=zeros(numpix,numslices);
>         end
>         fwhm_varatio
>      end
>
>      [df, ker_x, ker_y,
> K]=regularized_df(fwhm_varatio,D,Steps,numslices,df,fwhm_data);
>      df.rfx=round(df.rfx);
>      df.t=round(df.t);
>
>      if fwhm_varatio>0 & fwhm_varatio<Inf
>         % Smoothing varatio in slice is done using conv2 with a kernel
> ker_xy.
>         for slice=1:numslices
>            varatio_slice=reshape(varatio_vol(:,slice),numxs,numys);
>
>
> varatio_vol(:,slice)=reshape(conv2(ker_x,ker_y,varatio_slice,'same'),numpix,1);
>         end
>         % Smoothing betwen slices is done by straight matrix
> multiplication
>         % by a toeplitz matrix K normalized so that the column sums are
> 1.
>         varatio_vol=varatio_vol*K;
>      end
>
>   else
>      df.rfx=Inf;
>      df.t=df.fixed;
>   end
>
> *************************************************************
> In the case where fwhm_variatio > fwhm_varatio_limit it tries to calculate
> regularized_df which gives an error.  Instead it should go directly to:
>
>  else
>      df.rfx=Inf;
>      df.t=df.fixed;
>   end
>
>
>
>
> On Wed, 4 Jun 2008, Michael Ferreira wrote:
>
> > Hi Chris,
> >
> > Ditto to what Andrew said.
> >
> > I don't see anything wrong w/ your code.  And I've never come across
> > that error before.  So this is certainly one for Prof. Worsley.
> >
> > mike
> >
> >
> > Andrew Janke wrote:
> > > Hi,
> > >
> > > This really is one for Keith Worsley given that fmristat is his baby,
> > > I am unsure if he still reads this list so perhaps send him an email
> > > directly.
> > >
> > >
> > > a
> > >
> > > On Wed, Jun 4, 2008 at 5:54 AM, Penhune Lab <penhunelab at gmail.com>
> wrote:
> > >
> > >>Greetings!
> > >>
> > >>I am getting the following error from a multistat script that I have
> > >>successfully used for ... about 15 other analyses. Any ideas?
> > >>
> > >>Cheers,
> > >>Chris
> > >>
> > >>??? Maximum variable size allowed by the program is exceeded.
> > >>
> > >>Error in ==> multistat>regularized_df at 809
> > >>  ker_x=exp(-(-ceil(fwhm_x):ceil(fwhm_x)).^2*4*log(2)/fwhm_x^2);
> > >>
> > >>Error in ==> multistat at 418
> > >>     [df, ker_x, ker_y,
> > >>K]=regularized_df(fwhm_varatio,D,Steps,numslices,df,fwhm_data);
> > >>
> > >>All I am doing is averaging three files with the following code: (each
> of
> > >>the three files is the result of a covariation analysis of 15 Ss with
> their
> > >>Bx scores)
> > >>
> > >>input_files_effect =
> >
> >>['Lonlycov_PCOR_d1r1_ef.mnc';'Lonlycov_PCOR_d1r2_ef.mnc';'Lonlycov_PCOR_d1r3_ef.mnc'];
> > >>input_files_sdeffect =
> >
> >>['Lonlycov_PCOR_d1r1_sd.mnc';'Lonlycov_PCOR_d1r2_sd.mnc';'Lonlycov_PCOR_d1r3_sd.mnc'];
> > >>
> > >>x = [ones(size(input_files_effect,1),1)]
> > >>contrast = [1];
> > >>input_files_df = [];
> > >>input_files_fwhm = [];
> > >>which_stats = '_t _ef _sd _fwhm';
> > >>fwhm_varatio = -100;
> > >>df = multistat(input_files_effect, input_files_sdeffect,
> input_files_df,
> > >>input_files_fwhm, x, contrast, output_file_base, which_stats,
> fwhm_varatio);
> > >>
> > >>--
> > >>Laboratory for Motor Learning and Neural Plasticity
> > >>Directed by Dr. Virginia Penhune
> > >>Department of Psychology, Concordia University
> > >>SP-A 244, 7141 Sherbrooke St. W
> > >>Montreal, QC H4B 1R6 Canada
> > >>(514)848-2424 ext. 7567
> > >>http://psychology.concordia.ca/fac/penhune/winindex.html
> > >>_______________________________________________
> > >>MINC-users at bic.mni.mcgill.ca
> > >>http://www2.bic.mni.mcgill.ca/mailman/listinfo/minc-users
> > >>
> > >
> > >
> > >
> > >
> > _______________________________________________
> > MINC-users at bic.mni.mcgill.ca
> > http://www2.bic.mni.mcgill.ca/mailman/listinfo/minc-users
> >
> _______________________________________________
> MINC-users at bic.mni.mcgill.ca
> http://www2.bic.mni.mcgill.ca/mailman/listinfo/minc-users
>



-- 
Laboratory for Motor Learning and Neural Plasticity
Directed by Dr. Virginia Penhune
Department of Psychology, Concordia University
SP-A 244, 7141 Sherbrooke St. W
Montreal, QC H4B 1R6 Canada
(514)848-2424 ext. 7567
http://psychology.concordia.ca/fac/penhune/winindex.html