Coder Social home page Coder Social logo

gerstung-lab / deepsnv Goto Github PK

View Code? Open in Web Editor NEW
14.0 14.0 6.0 8.95 MB

The official repository for the deepSNV bioconductor packages. Branch master tracks the official Bioconductor development branch.

Home Page: https://master.bioconductor.org/packages/deepSNV/

R 90.93% C++ 5.98% C 1.74% Dockerfile 1.34%

deepsnv's People

Contributors

dtenenba avatar hpages avatar im3sanger avatar jmarshall avatar mg14 avatar nhayden avatar nturaga avatar sb43 avatar vobencha avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar

deepsnv's Issues

Create VCF-class object

Hi,

I'm trying to create a VCF-class object using:

summary(
object,
sig.level = 0.05,
adjust.method = "bonferroni",
fold.change = 1,
value = c("VCF")
)

but keep getting the following error: "Error in metadata(v)$header@header$META["date", 1] <- paste(Sys.time()) :
incorrect number of subscripts on matrix"

Creating the data.frame works... Please could you let me know what this error means and how I can solve it?

Thanks!

problem with loadAllData in R 4.0

I've recently upgraded to R 4.0

R version 4.0.3 (2020-10-10)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Mojave 10.14.6

when I run the loadAllData function in deepSNV

I get this error message:

Error in counts[i, , ] <- c[[i]] : 
  number of items to replace is not a multiple of replacement length
In addition: Warning message:
In mclapply(files, function(f) { :
  scheduled cores 1, 2 did not deliver results, all values of the jobs will be affected

The function works fine in R3.6

loadAllData seg faults

The test example works fine, but when I run it on some small test bams of my own I get a segfault:

> library(deepSNV)
> regions <- GRanges("chr6", IRanges(start = 45422511, end=45422911))
> files <- c("/home/nicholas_calicolabs_com/test_ref.bam", "/home/nicholas_calicolabs_com/test_test.bam")
> loadAllData(files, regions, q=30)

 *** caught segfault ***
address (nil), cause 'memory not mapped'

Traceback:
 1: bam2R(f, regions$chr[j], regions$start[j], regions$stop[j], ...)
 2: FUN(X[[i]], ...)
 3: lapply(X = X, FUN = FUN, ...)
 4: mclapply(files, function(f) {    test.matrix <- matrix(0, ncol = length(nucleotides), nrow = rows)    for (j in 1:nrow(regions)) {        test.matrix[beg[j]:end[j], ] = bam2R(f, regions$chr[j],             regions$start[j], regions$stop[j], ...)[, nucleotides]    }    mode(test.matrix) <- "integer"    test.matrix}, mc.cores = mc.cores)
 5: loadAllData(files, regions, q = 30)

The same thing happens when I run it on a single file in bam2R. These are sorted indexed bam files generated from subsetting a larger bam file using samtools.

Also this files url no longer exists "http://www.bsse.ethz.ch/cbg/software/deepSNV/data/test.bam"

Error building vignettes

Bioconductor currently fails to build the current release due to an issue with reading the VCF headers when building the vignettes.

##############################################################################
##############################################################################
###
### Running command:
###
###   /home/biocbuild/bbs-3.16-bioc/R/bin/R CMD build --keep-empty-dirs --no-resave-data deepSNV
###
##############################################################################
##############################################################################


* checking for file ‘deepSNV/DESCRIPTION’ ... OK
* preparing ‘deepSNV’:
* checking DESCRIPTION meta-information ... OK
* cleaning src
* installing the package to build vignettes
* creating vignettes ... ERROR
--- re-building ‘shearwaterML.Rmd’ using rmarkdown
[E::bcf_hdr_parse_sample_line] Could not parse the "#CHROM.." line, either FORMAT is missing or spaces are present instead of tabs:
	#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT

Quitting from lines 35-39 (shearwaterML.Rmd) 
Error: processing vignette 'shearwaterML.Rmd' failed with diagnostics:
error in evaluating the argument 'x' in selecting a method for function 'as.list': error in evaluating the argument 'x' in selecting a method for function 'as.list': no 'header' line "#CHROM POS ID..."?
--- failed re-building ‘shearwaterML.Rmd’

--- re-building ‘deepSNV.Rnw’ using knitr
--- finished re-building ‘deepSNV.Rnw’

--- re-building ‘shearwater.Rnw’ using knitr
[E::bcf_hdr_parse_sample_line] Could not parse the "#CHROM.." line, either FORMAT is missing or spaces are present instead of tabs:
	#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT

Quitting from lines 348-353 (shearwater.Rnw) 
Error: processing vignette 'shearwater.Rnw' failed with diagnostics:
error in evaluating the argument 'x' in selecting a method for function 'as.list': error in evaluating the argument 'x' in selecting a method for function 'as.list': no 'header' line "#CHROM POS ID..."?
--- failed re-building ‘shearwater.Rnw’

SUMMARY: processing the following files failed:
  ‘shearwaterML.Rmd’ ‘shearwater.Rnw’

Error: Vignette re-building failed.
Execution halted

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.