search for: p1l55

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&...

...trange that it shouldn't do the action. Has anybody come across this before? Any input on this strange issue much appreciated.. Hope I haven't missed something obvious. Thanks a lot, Aditi (P.s.- I've tried converting columns one by one to factors this time, and that works. >P1L55<-factor(P1L55) >levels(P1L55) >[1] "0" "1" Code: prm<-read.table("P:\\..... .csv", header=T, ...sep=",", ...) prmdf<-data.frame(prm) prmdf[2:13]<-lapply(prmdf[2:13], factor) ## action performed, no error message ##I tried to pick rand...

...E: [R] Error: length(f1) == length(f2) is not TRUE Yup, they are all factors- and its still doesn't work. Getting to the stage where I can use 'summary()' is the problem- the error stalls the process before a summary can be output. I tried the same code, but using 'family' and P1L55 as independent random factors, and that works. (I actually tried this with a new marker P1L74) > try(fit<-lmer(peg.no~1 + (1|family) + (1|P1L74), na.action=na.exclude)) > summary(fit) Linear mixed model fit by REML Formula: peg.no ~ 1 + (1 | family) + (1 | P1L74) AIC BIC logLik devia...

...w) that have two levels- 1 (presence) or 0 (absence). I need to estimate Variance Components for 'peg.no' (family wise, which is why this is the blocking factor), explained by each PxLyy, and by 'family' as well. I have managed to run the lmer model for the first marker column (P1L55), and have been able to block it according to 'family' as follows: vc<-read.table("P:\\R\\Testvcomp10.txt",header=T) attach(vc) family<-factor(family) colms<-(vc)[4:13] ## these are the marker columns lapply(colms,factor) library(lme4) for (f in levels(family)) { print...

...P:\\R\\Text Files\\Testvcomp10.txt",header=T,sep="",dec=".",na.strings=NA,strip.white=T) > attach(vc) > > family<-factor(family) > colms<-(vc)[,4:13] > > names(vc) [1] "male.parent" "family" "offspring.id" "P1L55" "P1L73" [6] "P1L74" "P1L77" "P1L91" "P1L96" "P1L98" [11] "P1L100" "P1L114" "P1L118" "peg.no" "ec.length" [16] "syll.len...