User: Arne Muller

gravatar for Arne Muller
Arne Muller20
Reputation:
20
Status:
New User
Location:
United States
Last seen:
3 years, 6 months ago
Joined:
4 years, 5 months ago
Email:
a**********@gmail.com

Posts by Arne Muller

<prev • 10 results • page 1 of 1 • next >
1
vote
1
answer
495
views
1
answer
slow granges in GAlignmentPairs
... Hello, I'm trying to extract the position and fragment length of paired end reads from bam files using readGAlignmentPairs and granges. There was already a posting (https://support.bioconductor.org/p/60339/) about the performance of readGAlignmentPairs, but another issue is that the granges method ...
genomicalignments written 3.5 years ago by Arne Muller20 • updated 3.5 years ago by Hervé Pagès ♦♦ 14k
3
votes
1
answer
865
views
1
answer
MedIP-Seq: IP over Input with edgeR
... Hello, I'd like to run edgeR on a MeDIP-Seq experiment with 4 treatment groups (each with 5 samples). I'm counting reads in slides over the genome (a few hundred bases with some overlap). I realised that there is a small batch effect as both IP and input separates per treatment group in an MDS plo ...
edge written 4.1 years ago by Arne Muller20 • updated 4.1 years ago by Aaron Lun25k
0
votes
2
answers
2.2k
views
2
answers
Answer: A: Problem installing BSgenome.Hsapiens.UCSC.hg19
... Hi. The BSgenome package has a dependency on the XML package (via rtracklayer). The XML package itself depends on the system library libxml2 which must be installed in your system with version 2.6.3 or higher. On my old Linux redhat system the library were too old. I can give hints how to solve this ...
written 4.2 years ago by Arne Muller20
0
votes
5
answers
1.2k
views
5
answers
Comment: C: readGAlignments + BamViews not working in BioC-Devel
... Hi - I think there's still a problem. Bioc-Release still uses 1.2.9 whereas in the svn there's' already 1.2.18 for release (note, the build/check report says there's a timeout for 1.2.18). ...
written 4.4 years ago by Arne Muller20
0
votes
1
answer
753
views
1
answer
estimateCellCounts: proportions and sex
... Hello, I am using the estimateCellCounts function that implements the Houseman regression to estimate blood cell types in Infinium 450k data derived from blood. As far as I understand this is a quadratic programming method that constraints the estimates to sum up to 1 per sample, but in my dataset ...
minfi written 4.4 years ago by Arne Muller20 • updated 3.9 years ago by Chi Hong0
0
votes
5
answers
1.2k
views
5
answers
Answer: A: readGAlignments + BamViews not working in BioC-Devel
... Hi - thanks for this fix in Rsamtools, but this only fixes the problem in one place, there's another (and possibly more) in GenomicAlignments::summarizeOverlaps: mapq_filter <- function(features, reads, algorithm, ignore.strand, inter.feature) { re ...
written 4.4 years ago by Arne Muller20 • updated 4.4 years ago by Martin Morgan ♦♦ 24k
0
votes
5
answers
1.2k
views
5
answers
Answer: A: readGAlignments + BamViews not working in BioC-Devel
... Hi - the error also occurs with R-3.2.0 and BiocDevel: aln <- readGAlignments(bv) Error in readGAlignments(bv) : Assertion on 'more.args' failed: Vector must be named Enter a frame number, or 0 to exit 1: readGAlignments(bv) 2: readGAlignments(bv) 3: Rsamtools:::.BamViews_delegate("r ...
written 4.4 years ago by Arne Muller20 • updated 4.4 years ago by Martin Morgan ♦♦ 24k
0
votes
5
answers
1.2k
views
5
answers
readGAlignments + BamViews not working in BioC-Devel
... Hi All, I'm wondering why the code below throws an error. I'm using R-devel (3.3.0) + BioC-devel. Is this a bug or did I miss something. From the GenomicAlignments help: ## With a BamViews object: fls <- system.file("extdata", "ex1.bam", package="Rsamtools", mustWork=TRUE) ...
rsamtools genomicalignments written 4.4 years ago by Arne Muller20
1
vote
1
answer
620
views
1
answer
csaw: input coverage RleLists instead of bam files
... Hello, would it be possible+reasonable to enter csaw peak calling with RleLists objects such as those generated by the GenomicAlignments' coverage method? I guess it'd be more performant to store and load the RleList objects compared to bam files. As I understand the coverage vector + the available ...
csaw written 4.4 years ago by Arne Muller20 • updated 4.4 years ago by Aaron Lun25k
1
vote
1
answer
1.2k
views
1
answer
FastqSampler with trimmed paired-end reads
... Hello, I came across a problems with FastqSampler from the ShortRead package. I get non-matching pairs/mates (reads that don;t belong to each other)  when doing: f1 <- FastqSampler("file1.fastq", n=1e6) f2 <- FastqSampler("file2.fastq", n=1e6) set.seed(123L); p1 <- yield(f1) set.seed(12 ...
shortread written 4.4 years ago by Arne Muller20 • updated 4.4 years ago by Martin Morgan ♦♦ 24k

Latest awards to Arne Muller

Popular Question 3.5 years ago, created a question with more than 1,000 views. For readGAlignments + BamViews not working in BioC-Devel
Popular Question 3.5 years ago, created a question with more than 1,000 views. For FastqSampler with trimmed paired-end reads

Help
Access

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