Question: get different heatmap from website and R
0
gravatar for ebrahimi.re67
12 months ago by
ebrahimi.re670 wrote:

Hi all,

enter image description here

I have differntial gene expression and I want to make a heat map. I put their expression on the clustvis website, and I got the different heat map from R. I have 142 genes and 4 samples. I used the R

this code:

####
file <- read.csv ("")
file <- file[order(),]
row.names() <- file$B
file <- [,2:5]
file_matrix <- data.matrix()
file_heatmap <- heatmap(_matrix, Rowv=NA, Colv=NA, col = cm.colors(256), sccale="column", margins=c(5,10))
###

the attached image I got (https://ibb.co/kyZ0af)

and the expression of gene 1 is

genes           S1                  S2              S3            S4
gene 1  174136.1093 160236.9196 145866.516  181905.6098
gene 18 52.56674663 79.04848213 55.4397647  1883.182956

the expression for gene18 in heatmap for S4 is white, but its expression in this sample is higher than others. the code and my heatmap are correct?

heatmap clustering excel R • 503 views
ADD COMMENTlink modified 10 months ago by Biostar ♦♦ 20 • written 12 months ago by ebrahimi.re670

Your 4th column seems to be empty here, but the heatmap doesn't. Your code seems somewhat correct, but as its partially given can't say much.

Please describe correctly what you want to know?

ADD REPLYlink written 12 months ago by sangram_keshari220

I want to make a heatmap and cluster these four stages, and S4 is not empty. in heatmap gene 18 in the S4 is white, but as you could see its expression is higher than others. the similar pattern could be seen for gene 1, I want to know my heatmap is correct or not. As I am a beginner in R, I tried a lot to make a matrix for using pheatmap, but I couldn't so I used this easy code.

genes S1 S2 S3 S4 gene1 174136.1093 160236.9196 145866.516 181905.6098 gene18 52.56674663 79.04848213 55.4397647 1883.182956

ADD REPLYlink modified 12 months ago • written 12 months ago by ebrahimi.re670

Okay got your point. So basically you saying your 18th gene have a high degree of expression but in heatmap its showing white.

Let me tell you the heatmap you generated using the webtool is clustered one, So it may possible that the order of genes may vary. Thats why you seeing on 18th row a white scale which may belong to a different gene.

ADD REPLYlink written 12 months ago by sangram_keshari220

actually, this heatmap is generated by R, the heatmap based on clustvis website show different cluster and for gene 1 in S1 shows the lowest expression.

The image from the website: (https://ibb.co/fVRgY0) So I don't know which one is correct!! image from R or website

ADD REPLYlink written 12 months ago by ebrahimi.re670

actually, this heatmap is generated by R, the heatmap based on clustvis website show different cluster and for gene 1 in S1 shows the lowest expression.

The image from the website: (https://ibb.co/fVRgY0) So I don't know which one is correct!! image from R or website

ADD REPLYlink written 12 months ago by ebrahimi.re670
1
gravatar for cpad0112
12 months ago by
cpad011212k
India
cpad011212k wrote:

OP code does column wide scaling (incorrectly typed sccale in OP code). Try drawing heatmap with column wise scaled values (do not scale scaled values).

ADD COMMENTlink modified 12 months ago • written 12 months ago by cpad011212k

Is it possible to explain more detailed about "do not scale scaled values"?

ADD REPLYlink written 12 months ago by ebrahimi.re670

It makes little sense to scale values that are already scaled. Can you confirm that changing sccale to scale improves the situation?

ADD REPLYlink written 10 months ago by Kevin Blighe49k
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: 723 users visited in the last hour