User: maltethodberg

gravatar for maltethodberg
Reputation:
70
Status:
Trusted
Location:
UCPH
Last seen:
20 hours ago
Joined:
2 years, 5 months ago
Email:
m************@gmail.com

Posts by maltethodberg

<prev • 52 results • page 1 of 6 • next >
4
votes
1
answer
32
views
1
answer
bpaggregate with formula returns strange error message
... Is bpaggregate from the BiocParallel package intended to be a drop-in replacement for aggregate? I am getting a somewhat cryptic error message when trying to use it: library(BiocParallel) aggregate(Sepal.Length~Species, data=iris, FUN=sum) # Works bpaggregate(Sepal.Length~Species, data=iris, FUN=s ...
biocparallel bpaggregate aggregate written 11 days ago by maltethodberg70 • updated 9 days ago by Martin Morgan ♦♦ 22k
5
votes
1
answer
100
views
1
answer
edgeR normalisation factors different between experimental groups
... I'm analysing an RNA-Seq datasets with edgeR and using calcNormFactors to estimate normalisation factors. In one of my experimental groups all samples have normalisation factors above 1,  while in another group all samples have normalisation factors below 1. Library sizes do not vary systematically ...
edger written 3 months ago by maltethodberg70 • updated 3 months ago by Aaron Lun20k
0
votes
1
answers
204
views
1
answers
Comment: C: Future plans for the new GPos class
... Just had look at this again: `coverage` now works on GPos-objects, but the "weight" argument behaves differently. It fails when trying to use a column as weights: pos_runs <- GRanges(c("chr1", "chr1", "chr2"), IRanges(c(1, 5, 9), c(10, 8, 15))) score(pos ...
written 3 months ago by maltethodberg70
3
votes
1
answer
204
views
1
answer
Future plans for the new GPos class
... In continuation of my previous question (https://support.bioconductor.org/p/104597/), I was wondering about the the plans for the new GPos class. I can see it has already been added to vignette, but also that it is currently being actively worked on (along with the underlying IPos object), with one ...
genomicranges gpos written 6 months ago by maltethodberg70 • updated 5 months ago by Hervé Pagès ♦♦ 13k
0
votes
0
answers
165
views
0
answers
coverage no longer works on GPos objects?
... After the recent update, it seems coverage no longer works on GPos objects?  library(GenomicRanges) ## Example 2: pos_runs <- GRanges("chrI", IRanges(c(1, 6, 12, 17), c(5, 10, 16, 20)), strand=c("+", "-", "-", "+")) gpos2 <- GPos(pos_runs) coverage(a ...
genomicranges gpos written 6 months ago by maltethodberg70
0
votes
1
answers
221
views
1
answers
Comment: C: breakTies for a Hits-object based on value in mcols?
... Thanks! So I guess that means there is no solution using breakTies - although as you write in the edited post, that seems like an obvious functionality for that function to have (The current documentation for breakTies is also a bit unclear). As was unaware of the global=TRUE argument to which.max/ ...
written 6 months ago by maltethodberg70
0
votes
1
answers
221
views
1
answers
Comment: C: breakTies for a Hits-object based on value in mcols?
... I get the following error when running this code: ​ > sort(hits, by = ~val, decreasing=TRUE) Error in as.vector(x) : no method for coercing this S4 class to a vector  In any case, doesn't breakTies always choose based on the index value, rather than row order? hits2 <- as(hits, "Hits") bre ...
written 6 months ago by maltethodberg70
4
votes
1
answer
221
views
1
answer
breakTies for a Hits-object based on value in mcols?
... Say I have Hits-object along these lines: from <- c(5, 2, 3, 3, 3, 2) to <- c(11, 15, 5, 4, 6, 11) hits <- Hits(from, to, 7, 15, sort.by.query=TRUE) For every hit, I can assign some value: mcols(hits)$val <- c(10, 11, 15, 12, 10, 10) I then want to break all ties, similar to brea ...
s4vectors hits written 7 months ago by maltethodberg70 • updated 7 months ago by Michael Lawrence10k
0
votes
0
answers
213
views
0
answers
interGeneCorrelation and voom
... I'm looking at using interGeneCorrelation to calculate correlation between sets of genes for RNA-Seq data. camera can take weights obtained via voom into account, but I was surprised to see that camera does not actually call interGeneCorrelation to do this. As the interGeneCorrelation function itsel ...
limma voom intergenecorrelation written 8 months ago by maltethodberg70
0
votes
2
answers
333
views
2
answers
Comment: C: limma-trend with random effects and/or array weights
... In line with the original question, which sequence of function calls would you use in that case for combining duplicateCorrelation with voom/voomWithArrayWeights? ...
written 9 months ago by maltethodberg70

Latest awards to maltethodberg

Popular Question 13 months ago, created a question with more than 1,000 views. For Correct sequence of functions for edgeR, edgeR robust and edgeR QL.
Popular Question 13 months ago, created a question with more than 1,000 views. For Correct sequence of functions for edgeR, edgeR robust and edgeR QL.
Supporter 2.2 years ago, voted at least 25 times.

Help
Access

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