Question: Rnaseq correlation between two samples
0
gravatar for jack1
4.4 years ago by
jack110
Italy
jack110 wrote:

Dear All,

 

I have two sample they  are extracted from the same source at  different time. 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 togheter 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

rna-seq deseq2 • 1.6k views
ADD COMMENTlink modified 4.4 years ago by Devon Ryan89k • written 4.4 years ago by jack110
1

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 REPLYlink modified 4.4 years ago • written 4.4 years ago by Devon Ryan89k
Please log in to add an answer.

Help
Access

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