## User: Paul Harrison

Paul Harrison •

**70**- Reputation:
**70**- Status:
- Trusted
- Location:
- Australia/Melbourne/Monash University Bioinformatics Platform
- Website:
- http://logarithmic.net...
- Twitter:
- @paulfharrison
- Scholar ID:
- Google Scholar Page
- Last seen:
- 2 days ago
- Joined:
- 6 years ago
- Email:
- p************@monash.edu

#### Posts by Paul Harrison

<prev
• 12 results •
page 1 of 2 •
next >

0

votes

1

answer

100

views

1

answers

... 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 weeks ago by
Paul Harrison •

**70**1

vote

1

answer

100

views

1

answers

... 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 4 weeks ago by
Paul Harrison •

**70**4

votes

1

answer

100

views

1

answer

...
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 4 weeks ago by
Paul Harrison •

**70**• updated 4 weeks ago by davide risso •**810**0

votes

1

answer

2.3k

views

1

answers

... 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 17 months ago by
Paul Harrison •

**70**0

votes

1

answer

1.3k

views

1

answers

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 19 months ago by
Paul Harrison •

**70**3

votes

1

answer

459

views

1

answer

... 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.4 years ago by
Paul Harrison •

**70**• updated 2.4 years ago by Gordon Smyth ♦**36k**0

votes

1

answer

404

views

1

answers

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 4.7 years ago by
Paul Harrison •

**70**0

votes

1

answer

404

views

1

answer

... 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 4.7 years ago by
Paul Harrison •

**70**• updated 4.7 years ago by Gordon Smyth ♦**36k**0

votes

2

answers

931

views

2

answers

... 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.0 years ago by
Paul Harrison •

**70**0

votes

2

answers

931

views

2

answers

... On Tue, Feb 5, 2013 at 11:23 AM, Gordon K Smyth
wrote:
> Dear Paul,
>
> The transformation that you propose is the same transformation that
is done
> by predFC(y) in the edgeR package, or by cpm(y,log=TRUE) in the
> developmental version of the edgeR package. The argument
prior.coun ...

written 6.0 years ago by
Paul Harrison •

**70**#### Latest awards to Paul Harrison

No awards yet. Soon to come :-)

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

Powered by Biostar
version 16.09

Traffic: 264 users visited in the last hour