使用 biofam 数据集
library(TraMineR)
data(biofam)
lab <- c("P","L","M","LM","C","LC","LMC","D")
biofam.seq <- seqdef(biofam[,10:25], states=lab)
head(biofam.seq)
Sequence
1167 P-P-P-P-P-P-P-P-P-LM-LMC-LMC-LMC-LMC-LMC-LMC
514 P-L-L-L-L-L-L-L-L-L-L-LM-LMC-LMC-LMC-LMC
1013 P-P-P-P-P-P-P-L-L-L-L-L-LM-LMC-LMC-LMC
275 P-P-P-P-P-L-L-L-L-L-L-L-L-L-L-L
2580 P-P-P-P-P-L-L-L-L-L-L-L-L-LMC-LMC-LMC
773 P-P-P-P-P-P-P-P-P-P-P-P-P-P-P-P
我可以拟合并显示回归树:
seqt <- seqtree(biofam.seq~sex + birthyr, data=biofam)
seqtreedisplay(seqt, type="I", border=NA, withlegend= TRUE, legend.fontsize=2, legendtext = "Biofam Regression Tree")
然后我可以识别叶子成员:
seqt$fitted[,1]
然而,这就是我感到困惑的地方。我怎么知道哪个叶子编号对应于图中的哪个叶子?该图似乎没有显示它,并且运行print(seqt)
似乎也没有给出叶子编号。
我想要实现的是分离出每片叶子中的序列,以便我可以分别在每片叶子上运行描述。我怎样才能做到这一点?