User: maltethodberg

gravatar for maltethodberg
maltethodberg100
Reputation:
100
Status:
Trusted
Location:
UCPH
Last seen:
2 days, 1 hour ago
Joined:
2 years, 7 months ago
Email:
m************@gmail.com

Posts by maltethodberg

<prev • 53 results • page 1 of 6 • next >
0
votes
1
answer
60
views
1
answer
mroast/fry with gene weights for most or all genes
... I’m experimenting with the ability of fry/mroast functions to include gene weights. I have two use cases:   1) Comparing the result of a previous experiment in mouse, with a current experiment in human: Basically I want to see if a previous observed response in mouse is observed in a similar human ...
limma mroast fry written 9 days ago by maltethodberg100 • updated 3 days ago by Gordon Smyth35k
4
votes
1
answer
55
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 10 weeks ago by maltethodberg100 • updated 10 weeks ago by Martin Morgan ♦♦ 22k
5
votes
1
answer
127
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 5 months ago by maltethodberg100 • updated 5 months ago by Aaron Lun21k
0
votes
1
answers
241
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 5 months ago by maltethodberg100
3
votes
1
answer
241
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 8 months ago by maltethodberg100 • updated 8 months ago by Hervé Pagès ♦♦ 13k
0
votes
0
answers
195
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 8 months ago by maltethodberg100
0
votes
1
answers
251
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 9 months ago by maltethodberg100
0
votes
1
answers
251
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 9 months ago by maltethodberg100
4
votes
1
answer
251
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 9 months ago by maltethodberg100 • updated 9 months ago by Michael Lawrence10k
0
votes
0
answers
231
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 10 months ago by maltethodberg100

Latest awards to maltethodberg

Student 7 weeks ago, asked a question with at least 3 up-votes. For Reference paper or resource for limma::diffSplice and edgeR::diffSpliceDGE methods?
Popular Question 15 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 15 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.4 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: 344 users visited in the last hour