User: nonCodingGene
nonCodingGene • 10
- Reputation:
- 10
- Status:
- New User
- Location:
- Last seen:
- 1 month, 3 weeks ago
- Joined:
- 3 years, 10 months ago
- Email:
- k**********@gmail.com
Posts by nonCodingGene
<prev
• 31 results •
page 1 of 4 •
next >
0
votes
1
answers
68
views
1
answers
... To me I think the two first ways of doing the analysis are the right ones, but I do not understand what the last two comparison mean, so this confuses my a little bit.
~1+Cell.type+Patient
coef = 2
~0+Cell.type+Patient
contrast = c(-1,1,0,0,0,0,0,0)
~0+Cell.type+Patient
coef = 2
~1+Cell.typ ...
written 9 weeks ago by
nonCodingGene • 10
4
votes
1
answer
68
views
1
answer
... I'm trying to perform a DEG analysis of two cell types extracted from 8 different patients (data comes from scRNA).
I'm trying to compare the two cells while adjusting for the patients, there are two possible models:
a) ~0+cell.type+patients
b) ~cell.type.patients
lrt <- glmWeightedF(fit, co ...
written 9 weeks ago by
nonCodingGene • 10
• updated 9 weeks ago by
James W. MacDonald ♦ 47k
0
votes
1
answers
103
views
1
answers
... What can be the source of this and how can I solve it?
I'm using edgeR for scRNA-seq, in order to manage droputs I use zinbwave
Just in case I've checked whether one condition has a gene for which its counts are all 0, and this does not happens.
This is how I calculate weights:
zinb <- zinb ...
written 9 weeks ago by
nonCodingGene • 10
0
votes
1
answers
103
views
1
answers
...
sum(is.na(dge$weights))
outputs 0
...
written 9 weeks ago by
nonCodingGene • 10
0
votes
1
answers
103
views
1
answers
... Sorry, I forgot to show that also tested for max(dge$weights) and min(dge$weights), and values where 1 and 0.
...
written 9 weeks ago by
nonCodingGene • 10
0
votes
1
answer
103
views
1
answer
... I'm using edgeR and an error for which I have not found a solution just appeared.
design <- model.matrix(~Cell.type, data = colData(summ.exp_norm))
dge$weights <- zinb.weights
dge <- estimateDisp(dge, design)
Error in .compressWeights(y, weights) : weights must be finite positive values ...
written 9 weeks ago by
nonCodingGene • 10
• updated 9 weeks ago by
Aaron Lun • 21k
0
votes
1
answers
264
views
1
answers
... Update, I've modify the classes vector such
classes <- factor(c("NEG","NEG","NEG", "MDK","MDK","MDK"), as.character(c("NEG", "MDK")))
Now logFC are quite close to the ones obtained in cuffdiff but none of the top 7 significant genes from cuffdiff are significant in limma.
3157 significant gen ...
written 12 months ago by
nonCodingGene • 10
0
votes
1
answer
264
views
1
answer
... I've the output of a differential expression analysis for RNA done with cuffdiff and I'm trying to check it with limma but I'm not getting the results that I would expect.
w164 <- read.table("COUNTfile.txt", header = T, row.names = 1)
head(w164)
WM164_NEG_1 WM164_NEG_2 WM164_NEG_3 W ...
written 12 months ago by
nonCodingGene • 10
• updated 12 months ago by
Gordon Smyth ♦ 35k
0
votes
1
answers
742
views
1
answers
... The identical() test returns a FALSE, so thats is the problem.
Thanks for the answer.
...
written 2.3 years ago by
nonCodingGene • 10
2
votes
1
answer
742
views
1
answer
... Sorry for this question, but for some reason something is happening that I'm not able to find. I'm having the next problem while creating the expression set.
pd <- data.frame(state = colnames(q@mat))
rownames(pd) <- pd[,1]
pd[,1] <- 1
pheno = new("AnnotatedDataFrame" ...
written 2.3 years ago by
nonCodingGene • 10
• updated 2.3 years ago by
Martin Morgan ♦♦ 22k
Latest awards to nonCodingGene
Popular Question
2.2 years ago,
created a question with more than 1,000 views.
For AnnotatedDataFrame is not a defined class
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.2.0
Traffic: 143 users visited in the last hour