Dear R Users,
I am having trouble with lmer. I am looking at recombinant versus non
recombinant individuals. In the response variable recombinant
individuals are coded as 1's and non-recombinant as 0's. I built a model
with 2 fixed factors and 1 random effect. Sex (males/females) is the
first fixed effect and sexual genotype (XY, YY, WX and WY) the second
one. Sexual Genotype is nested in sex. XY and YY individuals are males
and WX and WY females. I crossed 8 XY males with 8 WY females and 8 YY
males with 8 WX females. Each cross corresponds to a family (i.e. family
1 to 8 for XY x WY and family 9 to 16 for YY WX). For each family I have
20 offspring. Family is nested in sexual genotype as a random factor.
My data looks like:
reps<-factor(sample(c(0,1),640,replace=T,prob=c(0.95,0.05)))
sex<-factor(rep(c("M","F"),each=320))
geno<-factor(rep(c("XY","YY","WY","WX"),each=160))
fam<-factor(rep(c(1:8,9:16,1:8,9:16),each=20))
dat<-data.frame(reps,sex,geno,fam)
and I built the following model:
lmer(reps~sex+geno+(1|fam),data=dat,family=binomial)
and tried also
lmer(reps~sex/geno+(1|fam),data=dat,family=binomial)
but I keep getting this error:
asMethod(object) : matrix is not symmetric [1,2]
Does someone have an idea where the error might come from?
Thank you very much in advance!
Olivier
Sex-chromosome turnovers induced by deleterious mutation load
[[alternative HTML version deleted]]
______________________________________________
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.
