User: wunderl
wunderl • 20
- Reputation:
- 20
- Status:
- New User
- Location:
- Last seen:
- 2 months ago
- Joined:
- 5 months, 1 week ago
- Email:
- w******@cshl.edu
Posts by wunderl
<prev
• 10 results •
page 1 of 1 •
next >
0
votes
1
answer
137
views
1
answers
... Thanks, that makes sense.
I can confirm from my data that the 'outlier' sample from the raw count plot had both the largest library size and highest variance.
What fully convinced me was doing a simple normalization by library size. This ended up producing the same pattern of sample distances tha ...
written 9 weeks ago by
wunderl • 20
0
votes
2
answers
182
views
2
answers
... You are likely getting this error because of `log=TRUE` in your call to `cpm()`. The `log` of anything less than 1 will be negative, so if after normalization you end up with any fractional counts they will produce a negative number when you take the `log`. This issue also occurs with DESeq2's `rl ...
written 10 weeks ago by
wunderl • 20
0
votes
2
answers
182
views
2
answers
... If you want to truly compare pipeline performance, you should let each method do the normalization it's own way. Each method has been designed with its own philosophy and underlying assumptions, so mixing parts of one method with another is likely going to give you sub-optimal performance. (or, if ...
written 10 weeks ago by
wunderl • 20
3
votes
1
answer
137
views
1
answer
... ::Edit:: The embedded images are not formatting well, so I have changed them to external links.
I have been following the [DESeq2 vignette][1] for bulk RNA-Seq data and got the to section where it is recommended you make a [heatmap of the euclidian distances between samples][2].
In addition to the ...
written 12 weeks ago by
wunderl • 20
• updated
12 weeks ago by
Michael Love ♦ 25k
0
votes
1
answer
136
views
1
answers
... Thank you for the elaboration. I didn't realize `removeBatchEffect` operated at the log-transformed level, I can see now why it is distinct from the suite of tools `DESeq2` provides. The point about GC bias was very interesting as well. I am very new to this sort of space so I really appreciate th ...
written 4 months ago by
wunderl • 20
2
votes
2
answers
161
views
2
answers
... I am new to `DESeq2` myself so take this answer with a grain of salt.
I was initially confused by the difference between a `Wald test` and a `t-test` as well. To explain why a Wald test is used we first need to explain why we are fitting a generalized linear model (`GLM`).
A **lot** is going on ...
written 4 months ago by
wunderl • 20
0
votes
1
answer
136
views
1
answers
... Thanks Michael, I really appreciate the help!!
I will look into Salmon. So far my pipeline has been `STAR` --> `featureCounts` --> `DESeq2`. Would you say that not correcting for GC bias is a significant oversight?
Also, is there a reason `DESeq2` does not have anything analogous to `limma ...
written 4 months ago by
wunderl • 20
1
vote
1
answer
122
views
1
answer
... I have the following design matrix:
```
batch treatment
A 1
A 2
A 3
B 1
B 2
B 3
C 1
C 2
C 3
> batch
[1] A A A B B B C C C
Levels: A B C
> treatment
[1] 1 2 3 1 2 3 1 2 3
Levels: 1 2 3
```
I wo ...
written 4 months ago by
wunderl • 20
• updated
4 months ago by
Michael Love ♦ 25k
4
votes
1
answer
136
views
1
answer
... I am analyzing bulk RNA data from neurons derived from iPSC cells. The experimental design matrix is as follows:
| Lineage | | Method
------------- | ------------
A | |1
A | | 2
A | | 3
B | | 1
B | | 2
B | | 3
C | | 1
C | | 2
C | | 3
Here `lineage` indicates what iPSC line the neuronal ...
written 4 months ago by
wunderl • 20
• updated
4 months ago by
Michael Love ♦ 25k
1
vote
1
answer
117
views
1
answer
... When using Rsubread v 1.28.1 it sucessfully parsed my input file names by removing the file path and just keeping the names. As of the current version, however, the names contain the entire file path, this makes for very long and ugly column names.
Example
----
Setup
---
files: /path/to/file/here/ ...
written 5 months ago by
wunderl • 20
• updated
5 months ago by
Gordon Smyth ♦ 39k
Latest awards to wunderl
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: 390 users visited in the last hour