[R] Error: length(f1) == length(f2) is not TRUE

[A Singh]

I just wanted to know whether the chances of getting help with this would 
be greater if I posted this on the mixed models list, even though I would 
be cross-posting.. ?

--On 18 September 2009 12:41 +0100 A Singh <bzwas at bristol.ac.uk> wrote:

>
> Dear R users,
>
> I am trying to fit an lmer model with only random effects which is giving
> me the following error:
>
> Error : length(f1) == length(f2) is not TRUE
> In addition: Warning messages:
> 1: In P1L55:family :
>   numerical expression has 390 elements: only the first used
> 2: In P1L55:family :
>   numerical expression has 390 elements: only the first used
>
>
> I am trying to extract variance components for a phenotype 'peg.no',
> using the variable 'family' and a marker column 'P1L55' (which is
> categorical and has 2 levels- 0 and 1), as random effects. There are no
> fixed effects.
>
> The code I used is as follows:
>
> vc<-read.table(...)
>
> vcdf<-data.frame(vc)
> colms<-(vc)[4:13] ##these are the markers
> lapply(colms,factor)
>
> try(fit<-lmer(peg.no~1 + (1|family/P1L55), na.action=na.include,
> data=vcdf))
>
>
> I thought that putting the data into a dataframe would help, along with
> the na.exclude command, because there is a lot of missing data in patches
> which I have replaced with NA's, but I  don't know how to fix this error
> at all.
>
> Its a bit urgent, and any help is hugely appreciated.
>
> The data files are at:
>
>
> <http://www.4shared.com/file/131980362/460bdafe/Testvcomp10.html> (excel)
> http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html
> (txt)
>
>
> Cheers,
> Aditi
>
> ----------------------
> A Singh
> Aditi.Singh at bristol.ac.uk
> School of Biological Sciences
> University of Bristol
>
>
>



----------------------
A Singh
Aditi.Singh at bristol.ac.uk
School of Biological Sciences
University of Bristol