Traminer seqtree和seqtreedisplay的错误消息

Question

我在Mac上使用最新更新的R和Traminer来分析序列数据，并且在让seqtree和seqtreedisplay运行回归树时遇到了问题。使用TraMineR附带的biofam数据，R使用简单树和回归树给出了三个错误消息。

library(TraMineR)
library(TraMineRextras)

data(biofam)
bio.seq <- seqdef(biofam, 10:25, weights=bio$wp00tbgs, xtstep = 2)

## Create Ward clusters using an OM distance matrix

properties <- matrix(c(# left, married, child, divorced
  0, 0, 0, 0, # parent
  1, 0, 0, 0, # left
  0, 1, .5, 0,# married
  1, 1, 0, 0, # left+married
  0, 0, 1, 0, # child
  1, 0, 1, 0, # left+child
  1, 1, 1, 0, # left+married+child
  .5, 1, .5, 1 #divorced
), 8, 4, byrow=TRUE)

sm <- as.matrix(dist(properties))
indel <- .5 * max(sm)
bio.dist <- seqdist(bio.seq, method="OM", indel=indel, sm=sm, full.matrix=FALSE)

weight <- attr(bio.seq, "weights")
ward <- hclust(bio.dist, method="ward", members=weight)

## Tree display of the cluster solution and growing a regression tree

tree2 <- as.seqtree(ward, seqdata=bio.seq, diss=bio.dist, ncluster=2)
seqtreedisplay(tree2, type="I", border=NA, sortv="from.start", showdepth=TRUE) 

treereg <- seqtree(bio.seq ~ sex + plingu02, data=biofam, 
               diss=bio.dist)
seqtreedisplay(treereg, type="I", border=NA, sortv="from.start")

错误消息为

> seqtreedisplay(tree2, type="I", border=NA, sortv="from.start", showdepth=TRUE)
sh: dot: command not found
The file /private/var/folders/YY/YYp18zjXELKZUAbIpHvtIk+++TI/-Tmp-
/RtmpVNNMlt/tmpseqtree3745b091371.png does not exist.

> treereg <- seqtree(bio.seq ~ sex + cohort2 + religion + plingu02, data=bio,
 diss=bio.dist)
Error in DTNdisstree(dissmatrix = dissmatrix, predictor = predictor, terms = tterms,  : 
[!] To permute replicate, you should specify integer weights

> seqtreedisplay(treereg, type="I", border=NA, sortv="from.start")
Error in inherits(seqdata, "stslist") : object 'treereg' not found

Answer 1

您有以下错误：

• as.seqtree是由library(WeightedCluster)).
• You软件包提供的(因此你需要用GraphViz)
• Since加载它，需要安装http://www.graphviz.org/才能正常工作:请参阅weight.permutation = "diss" (安装GraphViz)
• Since后，你需要重新启动R你的序列对象是加权的，你需要指定权重应该如何排列(通常，我们使用weight.permutation = "diss"来排列距离矩阵)。

希望这能有所帮助。