sczhang at mail.ustc.edu.cn wrote:> Subject: [R] Ask for your help
> I got following error information when I run the haplo.glm program.Could
> you tell me which wrong was happened in my program?Many thanks. Shanchun
> 
> fit.gaus <- haplo.glm(y ~ SEX+geno, family = gaussian,
> + data=my.data, locus.label=label, control > +
haplo.glm.control(haplo.freq.min = 0.02))
> 
> Error in haplo.model.frame(m, locus.label = locus.label, allele.lev >
allele.lev, :
> Missing allele.lev = list of vectors for labels of alleles
> Check par list for haplo.glm
> 
allele.lev needs to be correctly specified, usually as:
allele.lev=attributes(geno)$unique.alleles
As documented, this is R specific.
| David Duffy (MBBS PhD)   email: David.Duffy at qimr.edu.au              ,-_|\
| ph: INT+61+7+3362-0217 fax: -0101                                   /     *
| Genetic Epidemiology Lab, Queensland Institute of Medical Research  \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A       v