I am using the test data with the commands outlined in the manual and am getting an error when using the epistatic.correlation function. Any suggestions?
library(WISH)
library(data.table)
data.table 1.12.2 using 16 threads (see ?getDTthreads). Latest news: r-datatable.com
ped <- fread("./test.ped", data.table=F)
tped <- fread("./test.tped", data.table=F)
genotype <- generate.genotype(ped, tped)
1103 passed QC
phenotype <- as.data.frame(fread("./test_pheno.txt"))
epistatic.correlation(phenotype, genotype, threads=16, test=T)
Running Test
Estimating run time based on ~100.000 models
Error in { : task 1 failed - "replacement has 406 rows, data has 203"
I am using the test data with the commands outlined in the manual and am getting an error when using the epistatic.correlation function. Any suggestions?