User: zhengx

gravatar for zhengx
zhengx30
Reputation:
30
Status:
New User
Location:
United States
Last seen:
2 months, 1 week ago
Joined:
4 years, 4 months ago
Email:
z*****@u.washington.edu

Posts by zhengx

<prev • 18 results • page 1 of 2 • next >
0
votes
2
answers
82
views
2
answers
Answer: A: seqVCF2GDS Error Converting VCF to GDS file
... You can directly modify the header in R: hdr <- seqVCF_Header("data/72S1.vcf.gz") hdr$info$Number[hdr$info$ID == "GMAF"] <- "." gdsfile <- seqVCF2GDS(vcffile, gdsfile, header=hdr) ...
written 10 weeks ago by zhengx30
1
vote
1
answer
85
views
1
answers
Answer: A: What does SNPRelate::snpgdsDiss calculate?
... See: Weir BS, Goudet J. A Unified Characterization of Population Structure and Relatedness. Genetics. 2017 Aug;206(4):2085-2103. doi: 10.1534/genetics.116.198424. `snpgdsDiss()` returns 1 - beta_ij. ...
written 4 months ago by zhengx30
0
votes
1
answer
105
views
1
answers
Answer: A: Different order of samples in output of snpgdsSlidingWindow to calculate Fst and
... SNPRelate re-orders "population" internally according to the order of sample IDs. res$sample.id is the sample order in the GDS file. If you are not sure whether the order of population is correct, you could order your input sample IDs as the order in the GDS file and provide population information ...
written 5 months ago by zhengx30
0
votes
1
answer
193
views
1
answers
Answer: A: SeqArray problem in seqVCF2GDS and convert back to vcf
... SeqArray_1.22.6 fixes the bug in the merged GDS file: merged <- seqOpen("merged.gds") seqAsVCF(merged) "Error in .Call2("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") : solving row 1: range cannot be determined from the supplied arguments (too many NAs)" ...
written 7 months ago by zhengx30
0
votes
1
answer
364
views
1
answers
Answer: A: SNPrelate, snpgdsFst is excluding all snps on non-autosomes, how can I get it to
... Try: snpgdsFst(genofile, sample.id = samp.sel, population = as.factor(pop.sel), method = "W&C84", autosome.only=FALSE)   ...
written 12 months ago by zhengx30
0
votes
1
answer
206
views
1
answers
Answer: A: Unable to build GDS with more than 255 scaffolds?
... Please use cvt.chr="char" in snpgdsBED2GDS(). See: ?SNPRelate::snpgdsBED2GDS​ ...
written 13 months ago by zhengx30
0
votes
2
answers
244
views
2
answers
Answer: A: Seq Array converting VCF to GDS using seqVCF2GDS, unable to find file problems w
... Please make sure your file exists. You might try file.size(vcf.fn) to see whether it returns the correct size of that VCF file.   ...
written 14 months ago by zhengx30
0
votes
0
answers
507
views
0
answers
Comment: C: Avoid 'dStr8' format in the GDS files
... Please downgrade the gdsfmt package. snpgdsCreateGeno() has no option to control the type of characters. ...
written 17 months ago by zhengx30
0
votes
1
answer
312
views
1
answers
Answer: A: Write gds in parallel
... The answer is simple: I have not implemented GDS writing with multiple processes at the same time so far. We usually create multiple independent GDS files, and each process writes data to each separate GDS file. After all child processes are done, multiple GDS files are merged in the main process. ...
written 17 months ago by zhengx30
0
votes
2
answers
451
views
2
answers
Answer: A: Error converting vcf to gds
... You might try SeqArray::seqVCF2GDS(), which could provide you more error information.   ...
written 19 months ago by zhengx30

Latest awards to zhengx

No awards yet. Soon to come :-)

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 354 users visited in the last hour