Closed AbubakariAbdulwasid closed 7 months ago
I think this error might be due to the way your vcf file was compressed/indexed. What happens if you decompress it and try again?
Please I use Fastree to generate the SNP.VCF and converted it with gzip command
when I use the decompressed file, this is what happened Read 98 phenotypes Detected binary phenotype Writing lineage effects to lineage_effects.txt Traceback (most recent call last): File "/home/public/PublicSoftware/pyseer", line 10, in
ValueError: zero-size array to reduction operation maximum which has no identity
Can you provide the commands you are using? I think for the purpose of testing what is going on with your VCF file you can probably omit all the --lineage
arguments
Perhaps also try compressing the vcf with bgzip
(from the samtools suite) rather than gzip
Please it has worked with "pyseer --phenotypes input.pheno --vcf snps.vcf --print-samples --no-distances> output_SNPs.txt"
Please I am running GWAS with pyseer using the SNPs variant with pyseer --phenotypes snps.pheno --vcf snps.vcf.gz --lineage --print-samples --no-distances> SNPs.txt and had this output Read 98 phenotypes Detected binary phenotype Writing lineage effects to lineage_effects.txt Traceback (most recent call last): File "/home/public/PublicSoftware/pyseer", line 10, in sys.exit(main()) File "/home/public/anaconda3/lib/python3.8/site-packages/pyseer/main.py", line 465, in main infile, sample_order = open_variant_file(var_type, var_file, options.burden, burden_regions, options.uncompressed) File "/home/public/anaconda3/lib/python3.8/site-packages/pyseer/input.py", line 278, in open_variant_file infile = VariantFile(var_file) File "pysam/libcbcf.pyx", line 4036, in pysam.libcbcf.VariantFile.init File "pysam/libcbcf.pyx", line 4299, in pysam.libcbcf.VariantFile.open File "pysam/libchtslib.pyx", line 516, in pysam.libchtslib.HTSFile.tell
NotImplementedError: seek not implemented in files compressed by method 1
Please is there a way I can resolved it?