Best way to normalize transcript reads data for PCA & correlation and further biostat analysis
1
1
Entering edit mode
12 weeks ago
Faith ▴ 40

Hello all,

I would like your advice on how to normalize, a Log2 normalized values of transcripts from DEGs flowerets in eight consecutive stages of development across, each stage has 3 biological replicates with their reads data.

I know they're logged, but my supervisor suggested we scale(scale=T, center=T) is this the best way to scale them? I tried it but the heatmap with scaling and without is nearly the same.

I'm planning to do PCA &Correlation on the data and further biostat analysis. Any advice is appreciated.

Snippet of the data

transcript_reads • 260 views
ADD COMMENT
1
Entering edit mode
12 weeks ago
Jeremy ▴ 870

Scaling the data will give you Z-scores. This is a standard way of preparing data for PCA. I would also suggest trying a center-log-ratio transformation, which sometimes does a better job of separating the data for PCA. You can use the clr() function in the compositions package in R to do that.

ADD COMMENT

Login before adding your answer.

Traffic: 2709 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