Hi,
When I used the data data.zip to run the code below, It reported the error Error in matexpo(Q * EL[i]) : NA/NaN/Inf in foreign function call (arg 1).
tree <- ape::read.tree("J:/MTB/gwas/tree0909.txt")
genotype1 <- read.table("J:MTB/gwas/test2",header = T,check.names = F)
rownames(genotype1) <- genotype1$sample_id
genotype1$sample_id <- NULL
rownames(genotype1)
phenotype1 <- read.table("J:MTB/gwas/test")
rownames(phenotype1) <- phenotype1$V1
phenotype1$V1 <- NULL
hogwash(pheno = as.matrix(phenotype1), geno = as.matrix(genotype1), tree = tree)
What happend? Can you give me some advice?
Thanks,
sun,