1つまたは2つのグループ変数を使用して多数のサブグラフを作成するために "by"を使用しようとしています。両方のグループ変数は因子変数である(性別はダミーであり、父親の社会的地位は複数のレベルを持つ)。凡例またはグラフのタイトルにグループのレベル(別名)を追加するにはどうすればよいですか?"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)")
}
切り取りと貼り付けが可能なコードを投稿してください。そのページにdata.frameというデータオブジェクトはありません。 –
おそらく、私のサンプルコードは私の質問をより明確にしています。 –
はい、より明確ですが、forループ内の 'seqdplot'呼び出しでエラーが発生します。 –