Pheatmap: ошибка в annotation_colors [[colnames (annotation) [i]]]: нижний индекс вне границ - PullRequest
1 голос
/ 30 марта 2019
library(pheatmap)

another<-read.table("~/Desktop/tcga3dcategorized.csv",sep=',',header=FALSE)
test<-another[1:3]
labels<-another[4]
colnames(test)=c("Dim 1","Dim 2","Dim 3")

annotation_row=data.frame(CancerType=factor(labels))
ann_colors=list(CancerType=c(KIRC="red", UCEC="#B2509E",LUSC="#D49DC7",
    LGG="#C1A72F",PCPG="#E8C51D",HNSC="#F9ED32",CESC="#104A7F",ESCA="#9EDDF9",
    BRCA="#007EB5",THCA="#CACCDB",PRAD="#6E7BA2",MESO="#DAF1FC",PAAD="#00AEEF",
    LUAD="#F6B667",LAML="#D97D25",UVM="#FBE3C7",GBM="#F89420",READ="#97D1A9",
    SKCM="#009444",KIRP="#754C29",COAD="#CEAC8F",LIHC="#3953A4",SARC="#BBD642",
    OV="#00A99D",BLCA="#D3C3E0",STAD="#A084BD",TGCT="#542C88",ACC="#FAD2D9",
    THYM="#ED1C24",KICH="#F8AFB3",DLBC="#EA7075",UCS="#7E1918",CHOL="#BE1E2D"))

pheatmap(test,annotation_row=annotation_row,annotation_colors=ann_colors,main="title")
```

Ошибка:

Error in annotation_colors[[colnames(annotation)[i]]] : 
  subscript out of bounds

Ответы [ 2 ]

1 голос
/ 08 июля 2019

Вы должны заставить свой фрейм данных иметь имена строк.

# after read.table
rownames(another) <- paste0("row_", seq(nrow(another)))

Посмотрите на Что-то странное в pheatmap (ошибка?)

0 голосов
/ 30 марта 2019

Проверьте ваши структуры данных - они не верны. Ваша структура another, созданная read.table и test, которая является ее подмножеством, dataframes. И вы пытаетесь передать их в pheatmap функцию, которая использует matrix - которая отличается от dataframe. Вы можете преобразовать ваш объект в матрицу следующим образом:

test_matrix <- data.matrix(test, rownames.force=NA) # change the rownames.force -argument if you want rownames

Тогда вы можете скормить его pheatmap. Если это не работает, вы можете проверить размеры объектов, подаваемых на pheatmap.

...