Sorry to pester the list with a simple question, but from googling it, I've
seen questions on the same error message several times but no answer.

I've run across this error using lme().  Seems like it's got to be something
small and format-y.  In trying to set up some reproducible data, I found a
slight difference that changes whether or not you get the error message.

So, originally the dataframe that I was working with has these attributes:
> str(tmp)
'data.frame':   22 obs. of  3 variables:
 $ V1: int  2 2 2 2 2 3 3 3 3 3 ...
 $ V2: num [1:22(1d)] 1.01 1.04 1.07 1.02 1.03 ...
 $ V3: num  4.21 3.68 7.12 3.32 1.82 ...

And if I run even a really simple model, I get the non-conformable arrays
error:

> lme(V2~V3,random=~1|V1, data=tmp)
Error in y[revOrder] - Fitted : non-conformable arrays

But, if you copy and paste the data back in, it drops the dimnames
attribute. Here's the output of str():
> str(tmp)
'data.frame':   22 obs. of  3 variables:
 $ species: int  2 2 2 2 2 3 3 3 3 3 ...
 $ hold   : num  1.01 1.04 1.07 1.02 1.03 ...
 $ react  : num  4.21 3.68 7.12 3.32 1.82 ...

And now it runs, no problem:
> lme(V2~V3,random=~1|V1, data=tmp)
Linear mixed-effects model fit by REML
  Data: tmp
  Log-restricted-likelihood: 17.31656
  Fixed: V2 ~ V3
(Intercept)          V3
 0.98104294  0.02465302

Random effects:
 Formula: ~1 | V1
        (Intercept)   Residual
StdDev:  0.03789563 0.08126392

Number of Observations: 22
Number of Groups: 5

Good to know that I'm not just crazy and something is different, but I still
can't find a way to get rid of that attribute.
 attributes(tmp)$dimnames<-NULL doesn't seem to do anything if you call it
again with str(), you get the same output.  Besides that, I don't understand
why the dimnames should cause problems with lme().

Any suggestions? Thanks in advance for the help.

	[[alternative HTML version deleted]]


