tests/testthat/test_02_scoreInvHap.R
In isglobal-brge/scoreInvHap: Get inversion status in predefined regions
context("Test scoreInvHap")
library(VariantAnnotation)
vcf_file <- system.file("extdata", "example.vcf", package = "scoreInvHap")
vcf <- readVcf(vcf_file, "hg19")
res1 <- scoreInvHap(SNPlist = vcf, inv = "inv7_005")
res1
res2 <- scoreInvHap(SNPlist = vcf, inv = "inv7_005", probs=TRUE)
res2
# microbenchmark::microbenchmark(res1 = scoreInvHap(SNPlist = vcf,
# inv = "inv7_005"),
# res2 = scoreInvHap(SNPlist = vcf,
# imputed=TRUE,
# inv = "inv7_005"))
isglobal-brge/scoreInvHap documentation built on Feb. 7, 2021, 8:43 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
context("Test scoreInvHap")
library(VariantAnnotation)
vcf_file <- system.file("extdata", "example.vcf", package = "scoreInvHap")
vcf <- readVcf(vcf_file, "hg19")
res1 <- scoreInvHap(SNPlist = vcf, inv = "inv7_005")
res1
res2 <- scoreInvHap(SNPlist = vcf, inv = "inv7_005", probs=TRUE)
res2
# microbenchmark::microbenchmark(res1 = scoreInvHap(SNPlist = vcf,
# inv = "inv7_005"),
# res2 = scoreInvHap(SNPlist = vcf,
# imputed=TRUE,
# inv = "inv7_005"))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.