## User: Yuya Liang

Reputation:
0
Status:
New User
Location:
United States
Last seen:
3 months, 4 weeks ago
Joined:
1 year, 1 month ago
Email:
r********@tamu.edu

#### Posts by Yuya Liang

<prev • 26 results • page 1 of 3 • next >
1
223
views
1
... I will read them again more carefully. Thank you for explaining everything! ...
written 5 months ago by Yuya Liang0
1
223
views
1
... Thank you for bringing this out again! I didn't express this clearly. I guess "Not recommend to use vst or rlog for statistical analysis of LFCs per gene" meaning not use the rlog counts for comparing control and treatment.  If rlog is not helpful for solving βir, why formula (2) we put rlog inste ...
written 5 months ago by Yuya Liang0
1
223
views
1
... Thank you so much!! So I can just simply say that this difference between raw count and rlog count due to "counts in all samples were moved toward the center(middle)" and this will help to calculate distances, correct?  (also help to calculate coefficient as well?) ...
written 5 months ago by Yuya Liang0
1
223
views
1
... Can you give me some information for why rlog makes raw counts from zero to non-zero? I expected to see huge different in this kind of extremely opposite (i.e. zero vs. some counts) but not. Even the rlog counts were not that different, they still got large LFC and significant pdaj. I'm just curiou ...
written 5 months ago by Yuya Liang0
1
223
views
1
... Got it!  I found a previous reply https://support.bioconductor.org/p/64014/#64016. Is this the reason why sample with 0 raw counts will have the value greater than 0?  ...
written 5 months ago by Yuya Liang0
1
223
views
1
... Thank you for replying! It helps a lot. I have two more questions:  This is the formula I cropped from DESeq2 paper:                          The log_qij for solving βir is the value from rlog(), is that correct? (I'm not sure is log2 or log10.)       2.   I also found a formula in the above ...
written 5 months ago by Yuya Liang0
1
223
views
1
... I was wondering (1) how does rlog( )calculate values from raw counts, (2) the relationship between rlog counts and the effect estimator (βir), and (3) what information rlog counts tell.   In vignette:                     log2(qij) = Xj βi Q1: Is log2(qij) here represent the values from rlog  ...
written 5 months ago by Yuya Liang0 • updated 5 months ago by Michael Love20k
1
145
views
1
... I saw it... Thank you!! ...
written 5 months ago by Yuya Liang0
1
145
views
1
... It seems not too much different. Is that normal? This is the code I used: ​ DESeq.A <- DESeq(DESeq.ds, test = "Wald") a <- results(DESeq.A) resLFC <- lfcShrink(DESeq.A, coef = 2, res = a, type = "normal") resLFC.ape <- lfcShrink(DESeq.A, coef = 2, res = a, type = "apeglm") par(mfrow ...
written 5 months ago by Yuya Liang0
1
145
views
1
... Thank you! It works. ...
written 5 months ago by Yuya Liang0

#### Latest awards to Yuya Liang

No awards yet. Soon to come :-)

Content
Help
Access

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