どうすれば次のようになりますか?
単一のプロットで、複数の箱ひげ図を作成したいと思います。各 X 変数は、以下に示すカテゴリ変数の組み合わせです。
data(mtcars)
# y variable is mpg
mtcars$vs = as.factor(mtcars$vs)
mtcars$cyl = as.factor(mtcars$cyl)
どうすれば次のようになりますか?
単一のプロットで、複数の箱ひげ図を作成したいと思います。各 X 変数は、以下に示すカテゴリ変数の組み合わせです。
data(mtcars)
# y variable is mpg
mtcars$vs = as.factor(mtcars$vs)
mtcars$cyl = as.factor(mtcars$cyl)
2 つの要因 (エンジン形状 (V/S) とシリンダー) のすべての組み合わせに対する単一のプロットを意味する場合、次のようになります。
with(mtcars,boxplot(mpg~interaction(as.factor(ifelse(vs,"S","V")),as.factor(cyl))))
abline(v=c(2.5,4.5),col=8)
(私が「S」と「V」を持っていると仮定すると-「ストレート」と「ヴィー」エンジン構成の場合-正しい方法; R 実装がホッキングの 1976 年の論文と一致している場合-ここの最初のページで説明されているように-私はそれだと思います正しいはずです)
次のような意味です。
data(mtcars)
y <- mtcars$mpg
vs <- as.factor(mtcars$vs)
cyl <- as.factor(mtcars$cyl)
par(mfrow=c(1,2))
plot(formula = y ~ cyl + vs)
lattice
同様に、次のようにパッケージを使用できます。
require(lattice)
bwplot( ~ y | vs + cyl)
注: このスタイルのプロットの唯一の問題は、結果が実際には 1 つに 6 つのプロットになることです ("vs" には 2 つのレベル * "cyl" には 3 つのレベルがあるため)。したがって、可能な組み合わせごとに十分なデータが得られない限り、結果のプロットは見栄えがよくないかもしれません...