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

[Marc BOUFFARD]

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
>