User: Biologist

gravatar for Biologist
Biologist70
Reputation:
70
Status:
Trusted
Location:
Last seen:
4 months, 2 weeks ago
Joined:
3 years, 7 months ago
Email:
v**********@gmail.com

Posts by Biologist

<prev • 158 results • page 1 of 16 • next >
0
votes
1
answer
964
views
1
answers
Comment: C: Use of batch corrected (removeBatchEffect) read counts for further downstream an
... Hi Aaron, Do you think this approach is right for making a clustering heatmap? library(edgeR) dge <- DGEList(counts=count) dge <- calcNormFactors(dge, method = “TMM”) logCPM <- cpm(dge,log=TRUE,prior.count=5) logCPM_bc <- removeBatchEffect(logCPM,batch=batch) ...
written 5 months ago by Biologist70
0
votes
1
answer
253
views
1
answers
Comment: C: edgeR removing batch effect before using expression data for clustering
... Sorry, I got you now. Before your first comment itself I saw the Gordon's comment in that post and thought of asking you whether it is right. Anyways now I got the point, I will use `voom` But If you dont mind could you please help me with the all the syntax I have to use for that. I never did this ...
written 5 months ago by Biologist70
0
votes
1
answer
253
views
1
answers
Comment: C: edgeR removing batch effect before using expression data for clustering
... Yes, James. But in this [Remove batch effect in small RNASeq study (SVA or others?)][1] Gordon recommends to do in the above way and then use it for MDS, PCA or heatmaps. So, I did it in that way. Is this wrong? [1]: https://support.bioconductor.org/p/59195/ ...
written 5 months ago by Biologist70
0
votes
1
answer
253
views
1
answers
Comment: C: edgeR removing batch effect before using expression data for clustering
... Thanks for the reply. So, after calculating logCPM, then I have to use removeBatchEffect function on that and then scale that for clustering analysis. Is it fine? logCPM <- cpm(y, prior.count=2, log=TRUE) logCPM <- removeBatchEffect(logCPM) logCPM <- t(scale(t(logCPM))) I ...
written 5 months ago by Biologist70
0
votes
1
answer
253
views
1
answer
edgeR removing batch effect before using expression data for clustering
... I have gene expression data for 8 samples, 4 controls and 4 FOXCUT OE samples. The raw counts data were in a dataframe `counts`. `coldata` looks like below: Samples Group Experiment Sample1 Control EXP1 Sample2 Control EXP1 Sample3 Co ...
batcheffect clustering edger R removebatcheffect written 5 months ago by Biologist70 • updated 5 months ago by James W. MacDonald51k
0
votes
2
answers
245
views
2
answers
Comment: C: How to get t-statistic values using edgeR?
... Gordon may I know your answer for my previous comment. thanks ...
written 5 months ago by Biologist70
0
votes
2
answers
245
views
2
answers
Comment: C: How to get t-statistic values using edgeR?
... Thanks Gordon. I also gave a try with `camera` now. For `camera` Is the following a right way? pathways.hallmark <- gmtPathways("h.all.v6.2.symbols.gmt") cam <- camera(y, pathways.hallmark, design, contrast=contrast.matrix, inter.gene.cor=0.01) head(cam,5) ...
written 5 months ago by Biologist70
2
votes
2
answers
245
views
2
answers
How to get t-statistic values using edgeR?
... Recently, I have seen a post about `fgsea` for gene set enrichment analysis. In this [DESeq results to pathways in 60 Seconds with the fgsea package][1]. In this post they use `stat` values column from `Deseq2` results to rank the genes for gene set enrichment analysis. With `edgeR` I use the foll ...
edger R gsea fgsea t-statistic written 5 months ago by Biologist70 • updated 5 months ago by Gordon Smyth38k
0
votes
1
answer
556
views
1
answers
Comment: C: EdgeR for differential analysis between two cell lines without replication
... Thank you very much for the reply. ...
written 15 months ago by Biologist70
0
votes
1
answer
556
views
1
answers
Comment: C: EdgeR error for differential analysis between two cell lines
... Hi Gordon, Thank you. I followed the tutorial and did the analysis. I have raw counts data of 72 genes for two cell-lines in dataframe "df". Three columns. First columns has genes and other columns are cell-lines. df <- data.frame(df[,-1], row.names=df[,1]) library(edgeR) y <- DGEList(cou ...
written 15 months ago by Biologist70

Latest awards to Biologist

Popular Question 10 months ago, created a question with more than 1,000 views. For Questions about analysis of Human Gene 2.0 ST Array
Popular Question 10 months ago, created a question with more than 1,000 views. For RNA degradation plot interpretation
Popular Question 10 months ago, created a question with more than 1,000 views. For Difference between DESeq and DESeq2
Epic Question 10 months ago, created a question with more than 10,000 views. For Heatmap with error message: `x' must be a numeric matrix
Popular Question 10 months ago, created a question with more than 1,000 views. For Difference between DESeq and DESeq2
Popular Question 14 months ago, created a question with more than 1,000 views. For Difference between DESeq and DESeq2
Good Question 16 months ago, asked a question that was upvoted at least 5 times. For How to check the read length for all the TCGA data samples in GDC?
Appreciated 2.3 years ago, created a post with more than 5 votes. For How to check the read length for all the TCGA data samples in GDC?
Student 2.3 years ago, asked a question with at least 3 up-votes. For How to check the read length for all the TCGA data samples in GDC?
Supporter 2.3 years ago, voted at least 25 times.
Centurion 2.3 years ago, created 100 posts.
Popular Question 2.3 years ago, created a question with more than 1,000 views. For Questions about analysis of Human Gene 2.0 ST Array
Great Question 2.3 years ago, created a question with more than 5,000 views. For Heatmap with error message: `x' must be a numeric matrix
Popular Question 2.7 years ago, created a question with more than 1,000 views. For Heatmap with error message: `x' must be a numeric matrix
Popular Question 3.5 years ago, created a question with more than 1,000 views. For Heatmap with error message: `x' must be a numeric matrix
Rising Star 3.5 years ago, created 50 posts within first three months of joining.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 249 users visited in the last hour