Sempre que vou rodar a análise bifatorial, onde há interação, a função da o seguinte erro:

########################################################################

...

Interacao significativa: desdobrando a interacao
------------------------------------------------------------------------

Desdobrando  F1  dentro de cada nivel de  F2
------------------------------------------------------------------------
------------------------------------------------------------------------
Quadro da analise de variancia
------------------------------------------------------------------------
          GL        SQ       QM    Fc  Pr.Fc
F2         4  194102.2 48525.55 4.768 0.0016
F1:F2      5        NA       NA  <NA>   <NA>
F1:F2 D0   5        NA       NA  <NA>   <NA>
F1:F2 D3   5        NA       NA  <NA>   <NA>
F1:F2 D5   5        NA       NA  <NA>   <NA>
F1:F2 D7   5        NA       NA  <NA>   <NA>
Residuo   87  885431.5 10177.37            
Total    107 3341005.8 31224.35            
------------------------------------------------------------------------

Error in if (des1.tab[(i + 2), 5] <= sigF) { :
  valor ausente onde TRUE/FALSE necessário

Aparentemente é um erro dentro do pacote, contudo, se alguém conseguiu contornar este bug de alguma forma, ou conhecer outro modo de rodar uma análise bifatorial no R, seria grato por receber suas opiniões.

A título de informação, estou chamando a função deste modo: "fat2.dic(data$TratA, data$TratB, data$Resp, quali=c(TRUE, TRUE), mcomp='tukey')"

Desde já, agradeço pela atenção,

Tiago