rwdavies / STITCH

STITCH - Sequencing To Imputation Through Constructing Haplotypes
http://www.nature.com/ng/journal/v48/n8/abs/ng.3594.html
GNU General Public License v3.0
76 stars 17 forks source link

Error with reference panel usage in input to eigen #43

Open rwdavies opened 3 years ago

rwdavies commented 3 years ago

From an email to the mailing list, seems the error is in the reference panel subsetting

Hi : when I use STITCH with reference panel, an error occurred: Error in eigen(h2, symmetric = TRUE) : error code 1 from Lapack routine 'dsyevr' and my command line: Rscript STITCH.R --outputdir=./ --method="diploid" --chr=22 --posfile=s5_22.allele --K=5 --nGen=2600 --nCores=6 --bamlist=bam06x.list --reference_populations='c("CHB","CHS","CHD")' --reference_haplotype_file=1000GP_Phase3_chr22.hap.gz --reference_legend_file=1000GP_Phase3_chr22.legend.gz --reference_sample_file=1000GP_Phase3.sample --S=10 I don't know how to fix it Thankes ! Sincerely , Lee