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