Rnaseq correlation between two samples
0
0
Entering edit mode
9.5 years ago
jack1 ▴ 10

Dear All,

I have two samples they are extracted from the same source at different times. What I want to do is to understand if inside at unsupervised analysis cluster this samples goes close.

In my hand at the moment both sample are not cluster together but they have a R^2 0f 0.94.

ddHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable=SampleTable,directory="Count/", design= ~1)
dds <- DESeq(ddHTSeq)

 we do blind estimation once, then re-use these for both:
rld <- rlog(ddsblind, blind=FALSE)

## ----euclDist------------------------------------------------------------
sampleDists <- dist(t(assay(rld)))
sampleDists
hc =hclust(sampleDists,labels=metadata$CDT)

What do you think? Could be useful to delete all the counts gene are zero? Or maybe it is wrong the normalization methods I use

Deseq2 RNA-Seq • 2.8k views
ADD COMMENT
1
Entering edit mode

How close/distant any given samples will be after clustering is dependent on the other samples present. Do you only have these two samples in the dataset? If not, it'd be impossible for them to cluster together unless all of their normalized values were identical.

ADD REPLY

Login before adding your answer.

Traffic: 1328 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6