User: Paul Harrison

Reputation:
70
Status:
Trusted
Location:
Australia/Melbourne/Monash University Bioinformatics Platform
Website:
http://logarithmic.net...
@paulfharrison
Scholar ID:
Last seen:
1 month, 1 week ago
Joined:
6 years, 4 months ago
Email:
p************@monash.edu

Posts by Paul Harrison

<prev • 13 results • page 1 of 2 • next >
1
183
views
1
... Ok, it's important not to guess wrong. I'll require the user to do something explicitly marking the assays appropriate to use. I think my solution will be to add some metadata fields to a SummarizedExperiment defining the appropriate pair of assays to use. I'll then write my code to produce an erro ...
written 3 months ago by Paul Harrison70
1
183
views
1
... Thanks Davide and Aaron. I had not thought about the distinction between frequencies and precision weights. So if I'm understanding correctly, zinbwave weights can be thought of as (fractional) frequencies, and contribute to the degrees of freedom as such. Also importantly, they are for "counts" ...
written 4 months ago by Paul Harrison70
1
183
views
1
... It looks like SingleCellExperiment has a weights() method, which is hopeful but I'm not clear on which assay the weights are meant to be associated with with. ...
written 5 months ago by Paul Harrison70
1
183
views
1
... I am trying to work out the best way to represent weights in the newer SummarizedExperiment and SingleCellExperiment classes. By weights I mean inverse variance (up to some scaling factor) for example as accepted by the lm() function. Weights are a somewhat simplistic but general way of acco ...
written 5 months ago by Paul Harrison70 • updated 5 months ago by davide risso830
1
2.5k
views
1
... Since RNA-seq is count based, in an RNA-seq analysis at the gene level the noise in individual counts for genes can be assumed to be Poisson (variance equal to the mean). This is an estimate of "technical variation", but does not include "biological varation". This would be the equivalent of the Kal ...
written 21 months ago by Paul Harrison70
1
1.6k
views
1
Comment: C: Limma-voom vs limma-trend
... If you examine the documentation to contrasts.fit, there is a warning that using weights with a non-orthogonal design matrix uses an approximation. Since voom uses weights on inidividual observations, and limma-trend does not, this could be an advantage to using limma-trend. ...
written 23 months ago by Paul Harrison70
1
502
views
1
... limma's barcodeplot function by default draws a pale red rectangle for statistics greater than sqrt(2) and a pale blue rectangle for statistics less than -sqrt(2). What is the reasoning is behind this choice? Is the expectation that statistics follows a standard normal distribution for genes that a ...
written 2.8 years ago by Paul Harrison70 • updated 2.8 years ago by Gordon Smyth37k
1
439
views
1
Comment: C: Behaviour of weights in limma
... On Thu, Jun 5, 2014 at 10:04 AM, Gordon K Smyth wrote: > Dear Paul, > >> Date: Wed, 4 Jun 2014 17:30:59 +1000 >> From: Paul Harrison >> To: Bioconductor mailing list >> Subject: [BioC] Behaviour of weights in limma >> >> Hello, >> >> I have so ...
written 5.0 years ago by Paul Harrison70
1
439
views
1
... Hello, I have some data from a variant of RNA-seq which I am hoping do some moderated t-test differential testing on with limma. In this data, many of the reads have sequenced through into the poly(A) tail, and we believe this gives us information about changes in poly(A) tail length. For each gen ...
written 5.0 years ago by Paul Harrison70 • updated 5.0 years ago by Gordon Smyth37k
2
989
views
2
... On Mon, Feb 4, 2013 at 7:56 PM, Wolfgang Huber wrote: > > Hi Paul > > the model (in the "parametric" case) is > > var(X) = a * E(X) + b * E(X)^2 > > with some a,b>0. And you are right that var(X)=0 <=> E(X)=0 > but since the support of X is positive, the onl ...
written 6.4 years ago by Paul Harrison70

Latest awards to Paul Harrison

No awards yet. Soon to come :-)

Content
Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.