2

我正在尝试使用 R-package 配置序列索引图的 y 轴TraMineR，我想出了如何为 y 轴设置全局限制，如果您想比较两个或多个之间的某些序列的数量，这将很有帮助分组，因为它均衡了规模。但我没有设法设置 y 轴刻度（如在 xtstep 中）。也许您可以使用以下示例代码快速帮助我：

library(TraMineR) 
data(mvad)
mvad.alphabet <- c("employment", "FE", "HE", "joblessness", "school",
                      "training")
mvad.labels <- c("Employment", "Further Education", "Higher Education",
                    "Joblessness", "School", "Training")
mvad.scodes <- c("EM", "FE", "HE", "JL", "SC", "TR")

## Define sequence objects
mvad.seq <- seqdef(mvad[, 17:86], alphabet = mvad.alphabet,
                     states = mvad.scodes, labels = mvad.labels, weights = mvad$weight, xtstep = 6)

## Plots
seqIplot(mvad.seq, group=mvad$gcse5eq, withlegend=TRUE, border=NA, xtstep=3, sortv="from.start") ## Default plot
seqIplot(mvad.seq, group=mvad$gcse5eq, withlegend=TRUE, border=NA, xtstep=3, sortv="from.start", ylim=c(0, 400)) ## Plot with custom ylim to compare the number of sequences between groups

默认的序列索引图如下所示，很难比较两组： 默认序列索引图

4

1 回答 1

2

当 a ylim在 中给出时seqIplot，它用于所有组。为了使绘图高度与每组中的加权序列数成比例，ylim应将上限值设置为最频繁组的值。

group <- mvad$gcse5eq
(nseq <- xtabs(mvad$weight ~ group))
(nmax <- max(nseq))
seqIplot(mvad.seq, group=group, withlegend=TRUE,
         border=NA, xtstep=3, sortv="from.start",
         ylim=c(0, nmax) )

y 轴上的刻度标签是序列索引。你可以通过给yaxis = FALSE来压制他们seqIplot。要显示您自己的标签，您可以发出类似的内容（有关详细信息，请参阅该axis功能的帮助）

axis(2, at = c(1, nseq[1]))

但在这种情况下，您应该使用 或seqIplot为每个组单独生成图，withlegend=FALSE并自己将图组织在单个图形layout中par(mfrow=...)。

于 2013-08-14T06:01:28.053 回答