User: evocanres

gravatar for evocanres
evocanres • 0
Reputation:
0
Status:
New User
Last seen:
4 months ago
Joined:
11 months, 3 weeks ago
Email:
e********@gmail.com

Profile information, website and location are not shown for new users.

This helps us discourage the inappropriate use of our site.

Posts by evocanres

<prev • 11 results • page 1 of 2 • next >
0
votes
1
answer
229
views
1
answers
Comment: C: Genes that have NOT changed
... resLA <- results(dds, lfcThreshold=.5, altHypothesis="lessAbs")  sum(!is.na(resLA$padj) < 0.1) can we say, this is the number of genes which are NOT different within logfoldchange of 0.5 ? ...
written 9 months ago by evocanres • 0
0
votes
1
answer
229
views
1
answers
Comment: C: Genes that have NOT changed
... Thank you Micheal, I did read it would be incorporated in the package, but didn't read the vignette carefully I guess.  Am I correct in reading that "lessAbs - |β| < x - p values are the maximum of the upper and lower tests " is the replacement for the that code for TOST (test for equivalency)?  ...
written 9 months ago by evocanres • 0
0
votes
1
answer
229
views
1
answers
Comment: A: Genes that have NOT changed
... I am trying to put the code that was run :) running into formatting issues ...
written 9 months ago by evocanres • 0
0
votes
1
answer
229
views
1
answers
Comment: C: Genes that have NOT changed
... 2. If I can replace that with betaSE <- rowData(dse)$SE_condition_untreated_vs_treated` because all.equal(mcols(dse), rowData(dse)) gives >TRUE > pvalDE <- 2 * pnorm( abs( beta ), sd = betaSE, lower.tail=FALSE ) > all( abs( pvalDE - results(dse)$pvalue ) < 1e-15, na.rm=TRUE ) ...
written 9 months ago by evocanres • 0
0
votes
1
answer
229
views
1
answers
Comment: C: Genes that have NOT changed
... I tried to run the exact code ```{r} > # Let's use the example data from the pasilla package > library( DESeq2 ) > library( pasilla ) > data( "pasillaGenes" ) > > # Create a DESeq2 data object from the pasilla data > dse <- DESeqSummarizedExperimentFromMatrix( + counts(pasill ...
written 9 months ago by evocanres • 0
0
votes
1
answer
229
views
1
answer
Genes that have NOT changed
... Hi I am trying to follow through this post as an exercise to find genes that are not differentially expressed. Few things in this post does not work [anymore]. Can somehow help ? 1 Specifically, 1. DESeqSummarizedExperimentFromMatrix. #does not work #Instead I use the `DESeqDataSetFromMatrix` t ...
deseq2 written 9 months ago by evocanres • 0
0
votes
1
answer
220
views
1
answers
Comment: C: DESeq results gives all genes
... Thank you. That helped actually. ...
written 10 months ago by evocanres • 0
0
votes
1
answer
168
views
1
answer
DGE Analysis between 3 conditions
... Hi, I have datasets for 3 different conditions["WT", "KO", "KO1"] (each with 3 replicates). [and two time points["0" and "12"]]. To use in DESeq, I have constructed a DSElist, with design ~ group where group has 6 levels. dds <- DESeqDataSetFromMatrix(countData = readCount, ...
rnaseq deseq2 written 10 months ago by evocanres • 0 • updated 10 months ago by Michael Love ♦ 25k
2
votes
1
answer
220
views
1
answer
DESeq results gives all genes
... Hi, I would appreciate any suggestions. I have three genotypes and two treatments. kable(readCountColData) sampleList |genotype |timePoint |group | |:----------|:--------|:---------|:------| |WT_0_1 |WT |0 |WT_0 | |WT_0_2 |WT |0 |WT_0 | |WT_0_3 |WT ...
deseq2 design matrix results written 10 months ago by evocanres • 0 • updated 10 months ago by Michael Love ♦ 25k
0
votes
1
answer
3.9k
views
1
answers
Comment: C: DESeq2 likelihood ratio test (LRT) design - 2 genotypes, 4 time points
... I figured out. In my samplelist i.e readCountcolData, I should have the group as a column.  ...
written 10 months ago by evocanres • 0

Latest awards to evocanres

No awards yet. Soon to come :-)
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 322 users visited in the last hour