Dear all,

typing

fitlm(rand(100,12),rand(100,1),'interactions')

yields an error message

error: anovan: the model terms matrix must list main effects above/before
 interactions
error: called from
    anovan at line 547 column 7
    fitlm at line 272 column 23

It seems that the index calculation for filling the TERM matrix is broken in this case. Here's a crude fix.

anovan.m+501:

        case 2
          ## Create term definitions for a model with two factor interactions
          if (N > 1)
            Nx = nchoosek (N, 2);
          else
            Nx = 0;
          endif
          TERMS = zeros (N + Nx, N);
          TERMS(1:N,:) = eye (N);
          cnt = N+1;
          for j = 1:N
            for i = j:N-1
              TERMS(cnt,j) = 1;
              TERMS(cnt,i+1) = 1;
              cnt = cnt+1;
            endfor
          endfor