5

我正在尝试使用“for”的“by”来使用一个或两个组变量创建许多子图。两个组变量都是因子变量（性别是一个虚拟变量，父亲的社会地位有多个层次）。如何在图例或图表标题中添加组的级别（也称为名称）？

这是我正在使用的代码。

library(TraMineR)  
library(Hmisc)
data(biofam)
biofam.lab <- c("Parent", "Left", "Married", "Left+Marr",
            "Child", "Left+Child", "Left+Marr+Child", "Divorced")
biofam.seq <- seqdef(biofam, 10:25, labels=biofam.lab)

class(biofam$sex)
levels(biofam$sex)
describe(biofam$sex)

class(biofam$cspfaj)
levels(biofam$cspfaj)
describe(biofam$cspfaj)

### Simple plots
seqdplot(biofam.seq)
seqdplot(biofam.seq, group=biofam$sex, title="Marital status by gender")

### Plot with automatic title using "by"
by(biofam.seq, biofam$sex, function(X) seqdplot(X, title="X$sex[1]"))
by(biofam.seq, biofam$sex, function(X) seqdplot(X, title=X$sex[1]))

### Plot with automatic title and multiple-grouping using "for"
for(n in c(1, 2, 3)) {
  seqdplot(subset(biofam.seq, subset=biofam$cspfaj==(n)), title="(n)")
}
for(n in c(1, 2, 3)) {
  seqdplot(subset(biofam.seq, subset=biofam$cspfaj==(n)), group=biofam$sex, title="(n)")
}
4

2 回答 2

2

您可以通过将 data.frame 子集传递为function(X)：

by(data.frame, data.frame$group, function(X) seqdplot(X, title=X$group[1]))
于 2012-12-07T02:14:37.373 回答
2

seqdplot 函数可以使用 group 参数自动为您完成：

seqdplot(seqobject, group=data.frame$group)

其中 seqobject 是使用 seqdef 函数创建的对象。

于 2012-12-07T12:30:56.843 回答