0

応答変数 ('Y' と呼びます) と 2 つの因子 (水準 A1 と A2 の因子 A、および水準 B1 と B2 の因子 B) があるとします。

ラティス パッケージの「ヒストグラム」関数を使用して、(A2 および B2) の応答に対して (A1 および B1) の応答をプロットできますか?

知っている

ヒストグラム(~y|因子A*因子B)

4つの組み合わせすべてをプロットします。しかし、その 2 つだけが必要な場合はどうすればよいでしょうか。

このサイトのすべての助けに感謝したいと思います!

4

2 に答える 2

2

1 つの方法は、4 つの組み合わせで新しい変数を作成し、subsetコマンドを使用することです。

FactorAB <- factor(paste(FactorA, FactorB, sep=""))
histogram(~y|FactorAB, subset=FactorAB %in% c("A1B1", "A2B2"))
于 2011-02-15T01:59:04.343 に答える
0

これを行うにはもっと簡単な方法があるはずですが、ここにあなたが望むものにかなり近づくはずのハックがあります.

# sample data
dat <- data.frame(Y=rpois(100,20),A=factor(c(rep("A1",50),rep("A2",50))),B=factor(c(rep("B1",50),rep("B2",50))))
dat$B <- sample(dat$B)

# create blank (colourless) histogram
p <- histogram(~Y|A*B,dat,col=0,border=0)
# subset and print blank panels
p[1,]
# draw data from desired panels onto blank "template"
trellis.focus("panel",1,1)
do.call("panel.histogram",trellis.panelArgs(p,1)[1:5])
trellis.focus("panel",1,2)
do.call("panel.histogram",trellis.panelArgs(p,4)[1:5])
trellis.unfocus()

残っているのは、上のプロットの下のストリップを変更することだけです。呼び出すtrellis.focus("strip",1,2)とフォーカスが置かれ、適切な他の呼び出しで変更できるはずですが、それが何であるかがわかりません。ただし、オーバープロットする可能性があります。

于 2011-02-15T14:35:46.990 に答える