Description Usage Arguments Details Value Author(s) See Also Examples
To merge GDS files of SNP genotypes into a single GDS file
1 2 3 | snpgdsCombineGeno(gds.fn, out.fn, method=c("position", "exact"),
compress.annotation="ZIP_RA.MAX", compress.geno="ZIP_RA",
same.strand=FALSE, snpfirstdim=FALSE, verbose=TRUE)
|
gds.fn |
a character vector of GDS file names to be merged |
out.fn |
the name of output GDS file |
method |
|
compress.annotation |
the compression method for the variables except
|
compress.geno |
the compression method for the variable
|
same.strand |
if TRUE, assuming the alleles on the same strand |
snpfirstdim |
if TRUE, genotypes are stored in the individual-major mode, (i.e, list all SNPs for the first individual, and then list all SNPs for the second individual, etc) |
verbose |
if TRUE, show information |
This function calls snpgdsSNPListIntersect
internally to
determine the common SNPs. Allele definitions are taken from the first GDS file.
None.
Xiuwen Zheng
snpgdsCreateGeno
, snpgdsCreateGenoSet
,
snpgdsSNPList
, snpgdsSNPListIntersect
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 | # get the file name of a gds file
fn <- snpgdsExampleFileName()
f <- snpgdsOpen(fn)
samp_id <- read.gdsn(index.gdsn(f, "sample.id"))
snp_id <- read.gdsn(index.gdsn(f, "snp.id"))
geno <- read.gdsn(index.gdsn(f, "genotype"), start=c(1,1), count=c(-1, 3000))
snpgdsClose(f)
# split the GDS file with different samples
snpgdsCreateGenoSet(fn, "t1.gds", sample.id=samp_id[1:10],
snp.id=snp_id[1:3000])
snpgdsCreateGenoSet(fn, "t2.gds", sample.id=samp_id[11:30],
snp.id=snp_id[1:3000])
# combine with different samples
snpgdsCombineGeno(c("t1.gds", "t2.gds"), "test.gds", same.strand=TRUE)
f <- snpgdsOpen("test.gds")
g <- read.gdsn(index.gdsn(f, "genotype"))
snpgdsClose(f)
identical(geno[1:30, ], g) # TRUE
# split the GDS file with different SNPs
snpgdsCreateGenoSet(fn, "t1.gds", snp.id=snp_id[1:100])
snpgdsCreateGenoSet(fn, "t2.gds", snp.id=snp_id[101:300])
# combine with different SNPs
snpgdsCombineGeno(c("t1.gds", "t2.gds"), "test.gds")
f <- snpgdsOpen("test.gds")
g <- read.gdsn(index.gdsn(f, "genotype"))
snpgdsClose(f)
identical(geno[, 1:300], g) # TRUE
# delete the temporary files
unlink(c("t1.gds", "t2.gds", "t3.gds", "t4.gds", "test.gds"), force=TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.