如何使用序列频率将第二个 y 轴添加到 seqfplot？

Question

我正在与 TraMineR 合作，对教育数据进行序列分析。我可以让 R 使用类似于以下的代码生成数据中最常见的 10 个序列的图：

library(TraMineR)

##Loading the data
data(actcal)

##Creating the labels and defining the sequence object
actcal.lab <- c("> 37 hours", "19-36 hours", "1-18 hours", "no work")
actcal.seq <- seqdef(actcal, 13:24, labels=actcal.lab)

## 10 most frequent sequences in the data
actcal.freq <- seqtab(actcal.seq)
actcal.freq

## Plotting the object
seqfplot(actcal.seq, pbarw=FALSE, yaxis="pct", tlim=10:1, cex.legend=.75, withlegend="right")

但是，我还想在图的右侧显示每个序列的频率（在对象 actcal.freq 中）。例如，上面代码创建的图中的第一个序列代表了 37.9% 的数据（如当前的图所示）。根据seqtab，这是 757 个主题。我希望数字 757 出现在右侧 y 轴上（其他序列以此类推）。

这可能吗？我玩过axis(side=4, ...)，但始终无法让它重现左 y 轴的间距。

Answer 1

好的。这有点乱，但如果您默认包含图例，该功能会重置par 设置，因此您需要将其关闭。然后您可以更轻松地设置轴，然后我们可以返回图例。这应该适用于您上面的测试数据。

#add padding to the right for axis and legend
par("mar"=c(5,4,4,8)+.1)

#plot w/o axis
seqfplot(actcal.seq, pbarw=FALSE, yaxis="pct", tlim=10:1, withlegend=F)

#plot right axis with freqs
axis(4, at = seq(.7, by=1.2, length.out=length(attr(actcal.freq,"freq")$Freq)), 
    labels = rev(attr(actcal.freq,"freq")$Freq), 
    mgp = c(1.5, 0.5, 0), las = 1, tick = FALSE)

#now put the legend on        
legend("right", legend=attr(actcal.seq, "labels"), 
    fill=attr(actcal.seq, "cpal"), 
    inset=-.3, bty="o", xpd=NA, cex=.75)

您可能需要稍微调整一下边距，尤其是legend 的inset= 参数才能正确放置。我希望您的真实数据与此相差不大，因为您确实必须深入研究该函数以查看它如何进行格式化以使事物匹配。