Dear colleagues, I face a problem doing haplotype analyses with haplostats: when I use the haplo.em function, the programme gives an error message because of 'recursive default argument reference.' I am not able to figure out what this means. Could you perhaps help me?
The full output is the following: > library(haplo.stats) > datafile.dat<-read.table("datafile.dat",header=TRUE) > attach(datafile.dat) > names(datafile.dat) [1] "SEXE" "AGE" "OUTCOME1" "OUTCOME2" "OUTCOME3" "OUTCOME4" [7] "OUTCOME5" "SNP1X" "SNP1Y" "SNP2X" "SNP2Y" "SNP3X" [13] "SNP3Y" "SNP4X" "SNP4Y" "SNP5X" "SNP5Y" "SNP6X" [19] "SNP6Y" "SNP7X" "SNP7Y" > block1<-datafile.dat[,c(8:21)] > loci<-c("1","2","3","4","5","6","7") > em<-haplo.em(geno=block1,locus.label=loci,miss.val=c(0,NA)) Error in .Call("R_lazyLoadDBfetch", key, file, compressed, hook, PACKAGE = "base") : recursive default argument reference Thank you very much in advance! Michiel Bos -- -------------------------------------------- Michiel J. Bos, MD MSc PhD-student Neuro-epidemiology Erasmus MC Dept. of Epidemiology and Biostatistics Room Ee2130 PO Box 1738 3000 DR Rotterdam tel: +31 (0)10-4087478 fax: +31 (0)10-4089382 E-mail: [EMAIL PROTECTED] ______________________________________________ R-help@stat.math.ethz.ch 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.