Moderator: Steve Lianoglou

gravatar for Steve Lianoglou
Reputation:
11,830
Status:
Trusted
Location:
Genentech
Scholar ID:
Google Scholar Page
Last seen:
1 day, 11 hours ago
Joined:
8 years, 3 months ago
Email:
s*********@gmail.com

I am a Computational Biologist at Genentech working in Cancer Immunology.

In a previous life, I was studying alternative cleavage and polyadenylation in humans.

Posts by Steve Lianoglou

<prev • 1,368 results • page 1 of 137 • next >
0
votes
1
answers
95
views
1
answers
Comment: C: Tac and Bioconductor giving different results
... You are likely right, you shouldn't be seeing such large differences. Again, no one can help you because we can't debug stories -- you need to provide code. You probably want to RMA normalize your data and use limma for analysis. You hopefully have some intuition of the experiment you are analyzing ...
written 16 days ago by Steve Lianoglou12k
0
votes
1
answers
95
views
1
answers
Answer: A: Tac and Bioconductor giving different results
... This is an impossible question to answer. Although I (and likely most people here) have never used "TAC affymetrix software" (I guess you mean this?), I imagine there are a multitude of ways to perform an analysis there just as there are many ways to perform an analysis using various bioconductor t ...
written 17 days ago by Steve Lianoglou12k
0
votes
2
answers
124
views
2
answers
Answer: A: Best DEG tool for datasets with FPKM counts?
... For analysis, you will want to log2 transform the fpkm values and use the limma-trend pipeline, see here: https://support.bioconductor.org/p/56275/ As for the gene_id column, you're only looking in the land of "funny looking" gene ids. Take a random sampling of the gene_id values, ie. if you have ...
written 21 days ago by Steve Lianoglou12k
0
votes
1
answers
148
views
1
answers
Comment: C: edgeR cpm filter
... Everything is working as it should and, further, your filtering strategy should be blind to the groups each sample belongs to. You can either increase your filtering threshold (ie. keep <- rowSums(cpm(y) > 2) >=5) if you want to remove more (any) genes, or you can switch to filtering based ...
written 5 weeks ago by Steve Lianoglou12k
0
votes
1
answers
148
views
1
answers
Comment: C: edgeR cpm filter
... The code itself looks like it should work, so when you say "this doesn't work", what do you mean, exactly? Are you getting an error somewhere, or? Also, if you have 6 replicates per group your keep rowSums threshold should technically be >= 6, but this probably won't change much in the downstrea ...
written 5 weeks ago by Steve Lianoglou12k
0
votes
2
answers
109
views
2
answers
Answer: A: Add new columns to DESeqResults from text file by matching rownames.
... Call as.data.frame on your DESqResults objects before the merge, ie. out <- merge(as.data.frame(result.ASE_Data), geneMap, ...) merge doesn't match by rownames, though, so you will have to specify the correct value for by, or a by.x and by.y combination of the column names in your result.ASE_D ...
written 7 weeks ago by Steve Lianoglou12k
0
votes
1
answers
131
views
1
answers
Comment: C: Differential expression of nomalised and log2 transformed data in limma
... Your design matrix is constructed from a factor (vector) that is 100 elements long. The first 50 are "1" and the second "50" are "2", but you said yourself that you have 25 elements per group, so ... why are you doing rep(..., each=50)? Fix it so that you have as many labels as samples ... and, btw ...
written 8 weeks ago by Steve Lianoglou12k
0
votes
2
answers
135
views
2
answers
Answer: A: DESeq2 Which gene used to normalize?
... Unless you specify a set of control genes to use, every gene which has all-non-zero-observations is used. See the help in ?estimateSizeFactors and ?estimateSizeFactorsForMatrix for more information. ...
written 8 weeks ago by Steve Lianoglou12k
0
votes
2
answers
115
views
2
answers
Answer: A: Expression profile from $E from Voom
... 1) expression profile in Elist has performed log(CPM) and been normalized for library size, right? Yes 2) Has the expression profile in Elist incorporated the precision weights (represented as $weights in v object) yet? No Note that if you want to use the expression data in your voomed "v" objec ...
written 8 weeks ago by Steve Lianoglou12k
0
votes
2
answers
176
views
2
answers
Answer: A: Treat Method - number of DE genes
... I'm not really sure what the issue is that you are trying to chase down is. I would recommend, however, that you don't keep updating your question and changing it slightly after interacting with Gordon, but maybe ask a different question that is specific to the new problem you think you are facing. ...
written 9 weeks ago by Steve Lianoglou12k

Latest awards to Steve Lianoglou

Scholar 6 months ago, created an answer that has been accepted. For C: DEseq2 with limited gene set
Scholar 19 months ago, created an answer that has been accepted. For C: DEseq2 with limited gene set
Scholar 19 months ago, created an answer that has been accepted. For A: unable to find Lumi package in R
Scholar 19 months ago, created an answer that has been accepted. For A: Filtering after a contrast
Teacher 19 months ago, created an answer with at least 3 up-votes. For C: DEseq2 with limited gene set
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: edgeR subsetting DGEList by column/sample
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: Filtering lowly expressed genes in voom-limma analysis
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: Counting ambiguously mapped reads per feature
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: Filtering after a contrast
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: limma experimental design
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: Perform limma on a subset of genes of interest
Scholar 19 months ago, created an answer that has been accepted. For A: unable to find Lumi package in R
Teacher 19 months ago, created an answer with at least 3 up-votes. For A: Perform limma on a subset of genes of interest
Scholar 21 months ago, created an answer that has been accepted. For A: unable to find Lumi package in R
Guru 21 months ago, received more than 100 upvotes.
Scholar 22 months ago, created an answer that has been accepted. For A: "Wobble" patterns for genome searching using Biostrings
Scholar 2.1 years ago, created an answer that has been accepted. For A: "Wobble" patterns for genome searching using Biostrings
Teacher 2.3 years ago, created an answer with at least 3 up-votes. For A: Perform limma on a subset of genes of interest
Teacher 2.3 years ago, created an answer with at least 3 up-votes. For A: Perform limma on a subset of genes of interest
Scholar 2.4 years ago, created an answer that has been accepted. For A: "Wobble" patterns for genome searching using Biostrings
Scholar 2.5 years ago, created an answer that has been accepted. For A: Using RNA-Seq raw count data in Weighted Gene Co-Expression Network Analysis
Scholar 2.7 years ago, created an answer that has been accepted. For A: Using RNA-Seq raw count data in Weighted Gene Co-Expression Network Analysis
Oracle 2.8 years ago, created more than 1,000 posts (questions + answers + comments).
Voter 2.8 years ago, voted more than 100 times.
Teacher 3.0 years ago, created an answer with at least 3 up-votes. For A: Using RNA-Seq raw count data in Weighted Gene Co-Expression Network Analysis

Help
Access

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