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]]

