Hi Bill,
Thanks, I did try out what you suggested but it doesn't seem to work.
I get the same error again.
There's obviously something here that I don't get. Need to figure it out.
Aditi
--On 18 September 2009 08:34 -0700 William Dunlap <wdun...@tibco.com> wrote:
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of A Singh
Sent: Friday, September 18, 2009 4:42 AM
To: r-help@r-project.org
Subject: [R] Error: length(f1) == length(f2) is not TRUE
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)
Note that you did not change columns 4:13 of vcdf to factors:
you changed copies of them to factors. Do
vcdf[4:13] <- lapply(vcdf[4:13], factor)
and lmer should be happier.
Bill Dunlap
TIBCO Software Inc - Spotfire Division
wdunlap tibco.com
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.ht
ml> (excel)
http://www.4shared.com/file/131980512/dc7308b/Testvcomp10.html
(txt)
Cheers,
Aditi
----------------------
A Singh
aditi.si...@bristol.ac.uk
School of Biological Sciences
University of Bristol
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
----------------------
A Singh
aditi.si...@bristol.ac.uk
School of Biological Sciences
University of Bristol
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
