Unable to build GDS with more than 255 scaffolds?
1
0
Entering edit mode
naglemi ▴ 10
@naglemi-16588
Last seen 5.4 years ago

I'm preparing to run assocTestMM (GENESIS) but first need to have a GDS file prepared by snpgdsBED2GDS. I'm working with a nonmodel organism that has 19 autosomes and over 1000 scaffolds which are not assigned to any autosome. Interestingly, my GDS files seem to be capped at 255 scaffolds and everything else is being lumped into scaffold "0". How can I resolve this issue so that I can run GWAS and detect associations on all scaffolds?

snpgdsBED2GDS(bedfile, famfile, bimfile, "gatk_882_WG_genotypes_biallelic_snps_VQSR_0.05maf_bis_pruned_totped_recoded_0.1geno.gds")

Some values of 'snp.chromosome' are invalid (should be finite and >= 1)!
Hint: specifying 'autosome.only=FALSE' in the analysis could avoid detecting chromosome coding.

 

> max(unique(read.gdsn(index.gdsn(genofile, "snp.chromosome"))))
[1] 255
> min(unique(read.gdsn(index.gdsn(genofile, "snp.chromosome"))))
[1] 0

 

 

 

 

GENESIS GWAS GDS • 1.1k views
ADD COMMENT
0
Entering edit mode
zhengx ▴ 30
@zhengx-7950
Last seen 4.9 years ago
United States

Please use cvt.chr="char" in snpgdsBED2GDS().

See:

?SNPRelate::snpgdsBED2GDS​

ADD COMMENT

Login before adding your answer.

Traffic: 1306 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6