0

R
を使用してヒートマップを生成しています。hclust関数を使用したところ、エラーが発生しました: Colv dendrogram doesn't match size of x.

以下は私のコードです。正方行列で完全に機能します。今、私の入力行列は400x2000です。このエラーが発生する理由を教えてもらえますか?

data <- read.table("mydata.txt",sep="\t",header=TRUE,row.names=1)
data_matrix <- data.matrix(data)
library(gplots)
library("RColorBrewer")
colors <- colorRampPalette(rev(brewer.pal(9, "Blues")) )(255)
hc_dist= dist(data_matrix)
hc_clust= hclust(hc_dist)
hr_dist= dist(t(data_matrix))
hr_clust= hclust(hr_dist)
heatmap.2(data_matrix, col=brewer.pal(11,"RdBu"), labRow=NA, density.info="none", scale="row",trace="none",Colv=as.dendrogram(hc_clust),     Rowv=as.dendrogram(hr_clust))

 Error in heatmap.2(data_matrix, col = brewer.pal(11, "RdBu"), labRow = NA,:Colv dendrogram doesn't match size of x
4

1 に答える 1

0

私のコメントを作成するには：

次の行列を想像してください。

matex<-matrix(1:50,ncol=5)
colnames(matex)<-paste0("col",1:5)
rownames(matex)<-paste0("row",1:10)
head(matex)
#      col1 col2 col3 col4 col5
#row1     1   11   21   31   41
#row2     2   12   22   32   42
#row3     3   13   23   33   43
#row4     4   14   24   34   44
#row5     5   15   25   35   45
#row6     6   16   26   36   46

今、あなたがしたように距離を計算します：

hc_dist= dist(matex) # what you think are between columns distances
hr_dist= dist(t(matex)) # what you think are between rows distances

そしてオブジェクトhc_distをチェックしてくださいhr_dist：

hc_dist
           # row1      row2      row3      row4      row5      row6      row7      row8      row9
# row2   2.236068                                                                                
# row3   4.472136  2.236068                                                                      
# row4   6.708204  4.472136  2.236068                                                            
# row5   8.944272  6.708204  4.472136  2.236068                                                  
# row6  11.180340  8.944272  6.708204  4.472136  2.236068                                        
# row7  13.416408 11.180340  8.944272  6.708204  4.472136  2.236068                              
# row8  15.652476 13.416408 11.180340  8.944272  6.708204  4.472136  2.236068                    
# row9  17.888544 15.652476 13.416408 11.180340  8.944272  6.708204  4.472136  2.236068          
# row10 20.124612 17.888544 15.652476 13.416408 11.180340  8.944272  6.708204  4.472136  2.236068

hr_dist
          # col1      col2      col3      col4
# col2  31.62278                              
# col3  63.24555  31.62278                    
# col4  94.86833  63.24555  31.62278          
# col5 126.49111  94.86833  63.24555  31.62278

実際にhc_distは、行間の距離とhr_dist列間の距離が含まれているため、エラーが発生します（行と列の数が同じであるため、正方行列では得られません）

于 2015-02-06T13:51:37.963 に答える