User: Biologist

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

Posts by Biologist

<prev • 158 results • page 1 of 16 • next >
0
votes
1
answer
756
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 7 days ago by Biologist70
0
votes
1
answer
123
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 8 days ago by Biologist70
0
votes
1
answer
123
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 8 days ago by Biologist70
0
votes
1
answer
123
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 8 days ago by Biologist70
0
votes
1
answer
123
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 9 days ago by Biologist70 • updated 9 days ago by James W. MacDonald50k
0
votes
2
answers
112
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 15 days ago by Biologist70
0
votes
2
answers
112
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 16 days ago by Biologist70
2
votes
2
answers
112
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 16 days ago by Biologist70 • updated 16 days ago by Gordon Smyth37k
0
votes
1
answer
428
views
1
answers
Comment: C: EdgeR for differential analysis between two cell lines without replication
... Thank you very much for the reply. ...
written 10 months ago by Biologist70
0
votes
1
answer
428
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 10 months ago by Biologist70

Latest awards to Biologist

Popular Question 5 months ago, created a question with more than 1,000 views. For Difference between DESeq and DESeq2
Popular Question 9 months ago, created a question with more than 1,000 views. For Difference between DESeq and DESeq2
Good Question 11 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?
Centurion 23 months ago, created 100 posts.
Great Question 23 months ago, created a question with more than 5,000 views. For Heatmap with error message: `x' must be a numeric matrix
Appreciated 23 months ago, created a post with more than 5 votes. For How to check the read length for all the TCGA data samples in GDC?
Popular Question 23 months ago, created a question with more than 1,000 views. For Questions about analysis of Human Gene 2.0 ST Array
Student 23 months 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 23 months ago, voted at least 25 times.
Popular Question 2.3 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.1 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.1 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: 171 users visited in the last hour