User: riccardo

gravatar for riccardo
riccardo70
Reputation:
70
Status:
New User
Location:
Last seen:
7 months, 2 weeks ago
Joined:
3 years, 6 months ago
Email:
r**************@gmail.com

Posts by riccardo

<prev • 93 results • page 1 of 10 • next >
0
votes
1
answer
163
views
1
answer
Bad replicate DESeq2
... Hello, Here: https://imgur.com/z7Fkyjv you can find the PCA of the experiment that I am analyzing. Cond1 e cond2 were designed to have an over expression of two different genes. The sample cond1 near to controls has an expression of that gene much less than the others biological replicates. In this ...
deseq2 written 12 months ago by riccardo70 • updated 12 months ago by Michael Love24k
0
votes
1
answer
663
views
1
answers
Comment: C: DESeq2 condition table
... I just realized that if I put "control" in the condition I get the message and checking the code I saw this: controlSynonyms <- c("control", "Control", "CONTROL") for (cSyn in controlSynonyms) { if (cSyn %in% lastDVLvls) { if (cSyn != lastDVLvls[1]) { message(past ...
written 12 months ago by riccardo70
0
votes
1
answer
663
views
1
answers
Comment: C: DESeq2 condition table
... Hello, I did again the same thing that I posted before but this time using the same approach but with different values in the condition table I got this message: it appears that the last variable in the design formula, 'condition', has a factor level, 'control', which is not the reference lev ...
written 12 months ago by riccardo70 • updated 12 months ago by Michael Love24k
0
votes
1
answer
510
views
1
answers
Comment: C: find fold change in DESeq2 results
... Thank you to all for the information. ...
written 2.0 years ago by riccardo70
1
vote
1
answer
510
views
1
answer
find fold change in DESeq2 results
... Hi, maybe I am doing something wrong but I cannot find a log2FoldChange inside the results of DESeq2: res = results(dds, contrast=c("cond", test, control)) dRes = as.data.frame(res) dRes$log2FoldChange[1] [1] 12.62795 dRes[dRes$log2FoldChange == 12.62795,] [1] baseMean log2FoldChange lfc ...
deseq2 written 2.0 years ago by riccardo70 • updated 2.0 years ago by Michael Love24k
0
votes
1
answer
826
views
1
answers
Comment: A: deseq2 course experiment
... Hi Michael, thank you for your help. My WT2 and WT5 are my controls at different times but what if later I will be only interested in the comparison of, for example, A2 vs B2 and A5 vs B5? Should I perform again the analysis excluding WT2 and WT5 from the table of counts used by DESeq2? Or is it be ...
written 2.1 years ago by riccardo70
0
votes
1
answer
826
views
1
answers
Comment: C: deseq2 course experiment
... Hi Michael, I think that I found the problem, removeBatchEffect wants the log counts so using the log2 of the counts normalized by DESeq2 it seems work, do you agree?   ...
written 2.1 years ago by riccardo70
0
votes
1
answer
826
views
1
answers
Comment: C: deseq2 course experiment
... Hi I think that there is a problem because some values are negative, what could be the problem?  I have also another question regarding this code: dat <- counts(dds, normalized = TRUE) idx <- rowMeans(dat) > 1 dat <- dat[idx, ] mod <- model.matrix(~cond, colData(dds)) mod0 < ...
written 2.1 years ago by riccardo70
0
votes
1
answer
826
views
1
answers
Comment: C: deseq2 course experiment
... Thank you! If I wanted the read counts normalized without batch effects could I use this command? removeBatchEffect(counts(dds, normalized = TRUE), covariates = as.matrix(ddssva@colData$SV1, ddssva@colData$SV2, ddssva@colData$SV3)) Thank you ...
written 2.1 years ago by riccardo70
0
votes
1
answer
826
views
1
answers
Comment: C: deseq2 course experiment
... Sorry I think that this is more correct: rldBlindTrueRemBatch = removeBatchEffect(rldBlindTrueCounts, covariates = as.matrix(ddssva@colData$SV1, ddssva@colData$SV2, ddssva@colData$SV3)) Since my ddsva@colData is:                cond          SV1         SV2         SV3            <factor> ...
written 2.1 years ago by riccardo70

Latest awards to riccardo

Popular Question 12 months ago, created a question with more than 1,000 views. For edgeR TMM values
Popular Question 15 months ago, created a question with more than 1,000 views. For edgeR TMM values
Great Question 15 months ago, created a question with more than 5,000 views. For edgeR TMM values
Popular Question 2.0 years ago, created a question with more than 1,000 views. For edgeR TMM values
Popular Question 2.0 years ago, created a question with more than 1,000 views. For DESeq2 old versions
Great Question 2.0 years ago, created a question with more than 5,000 views. For DESeq2 FPKM normalization
Great Question 2.0 years ago, created a question with more than 5,000 views. For edgeR TMM values
Popular Question 2.2 years ago, created a question with more than 1,000 views. For DESeq2 and shrinkage of log2 fold changes
Popular Question 2.4 years ago, created a question with more than 1,000 views. For DESeq2 and ComBat
Popular Question 3.2 years ago, created a question with more than 1,000 views. For DESeq2 FPKM normalization
Popular Question 3.2 years ago, created a question with more than 1,000 views. For DESeq2 and ComBat

Help
Access

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