[英]TraMineR graphics not reactive to layout of R
##PACKAGES
library(tidyverse,quietly=TRUE)
library(TraMineR)
library(WeightedCluster, quietly = TRUE)
library(viridis)
library(seqhandbook, quietly = TRUE)
這是我的數據集的表示:
mydata<-structure(list(T1 = structure(c("ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN"
), label = "Type chir. bariat."),
T2 = structure(c("ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN"), label = "Type chir. bariat."),
T3 = structure(c("ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN"), label = "Type chir. bariat."),
T4 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN"), label = "Type chir. bariat."),
T5 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN"), label = "Type chir. bariat."),
T6 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN"), label = "Type chir. bariat."),
T7 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ABL"), label = "Type chir. bariat."),
T8 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ABL"), label = "Type chir. bariat."),
T9 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ABL"), label = "Type chir. bariat."),
T10 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ABL"), label = "Type chir. bariat."),
T11 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ABL",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T12 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ABL",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T13 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "ABL", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T14 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "ABL", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T15 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T16 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T17 = structure(c("ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T18 = structure(c("ANN",
"ABL", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T19 = structure(c("ANN",
"ABL", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat."),
T20 = structure(c("ANN",
"ABL", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN", "ABL", "ANN",
"ANN", "ANN", "ANN", "SLE", "SLE", "ANN", "ANN", "ANN", "BPG",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "SLE", "ANN", "ANN", "ANN", "ANN", "ANN", "ANN",
"ANN", "ANN", "ANN", "SLE"), label = "Type chir. bariat.")),
row.names = c(NA,-50L), class = c("tbl_df", "tbl", "data.frame"))
我想繪制一個圖形和彼此相鄰的圖例,因為顯然, seqheatmap from seqhandbook
中的函數seqheatmap from seqhandbook
不顯示圖例。
我已經分別完成了這兩個工作,並希望使用layout from R base
的函數layout from R base
來做到這一點
我創建序列的代碼
# ALPHABETS AND LABELS
labels <- sort(c("ANN","SLE","BPG","BPD","ABL","DCD"))
Etats_lng<-c("Ablation","Anneau","Derivation bilio-pancréatique","By-pass gastrique",
"Décès","Sleeve")
seq <- seqdef(
mydata[, paste("T",1:20,sep = "")],
alphabet = labels,
states = Etats_lng,
cpal = viridis(6, direction = -1)
)
# distances
couts <- seqsubm(seq, method = "CONSTANT", cval = 2)# Matrice des couts de substitution
seq.dist<- seqdist(seq, method = "OM", indel = 1, sm = couts)# Distance avec optimal-Matching
seq.hclust <- hclust(as.dist(seq.dist), method = "ward.D2") #Clustering des trajectoire
#SETTING OF LAYOUT
layout.matrix <- matrix(c(1,2), nrow = 1, ncol = 2)
layout(mat = layout.matrix,
widths = c(3,1) )# Widths of the two columns
layout.show(2)
在第一幀中,我想放置熱圖,在第二幀中放置圖例
但我不工作
seq_heatmap(seq, seq.hclust)
seqlegend(seq)
函數layout
在您的示例中不起作用,因為layout
不能嵌套並且heatmap
(由seq_heatmap
調用)已經使用layout
來生成熱圖。
我看到的唯一的解決辦法是檢索功能的源代碼heatmap
(從數據),將其重命名說myheatmap
,並對其進行修改以添加顏色圖例的顯示。
檢索heatmap
的代碼
myheatmap <- heatmap
edit(myheatmap)
要使edit
工作,您可能必須指定編輯器。 或者,您可以使用“ View
並將內容復制粘貼到您選擇的編輯器中。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.