地上と地下のバイオマスを 3 つの要因に対してバープロットする必要があります。各バイオマスについて、格子からの棒グラフを使用してプロットできますが、そのうちの 2 つ、つまり Y 軸の正の部分に地上のバイオマス、負の部分に地下のバイオマスをプロットすることについてはわかりません。これが私のデータのテスト コードです。ここで、A、B、C は完全な要因因子です。
library(lattice)
above <- runif(108)
below <- runif(108)
A <- rep(c("ab", "bc", "cd"), each=36)
B <- rep(1:3, 36)
C <- gl(3, 4, 108, c(30, 60, 90))
barchart(above~A|B+C)
barchart(below~A|B+C)
#This is what I used to do for two factors using barplot
par(mfrow=c(2, 1), mai=c(0, 1, 0.5, 0.5))
agg.abv <- aggregate(above, by=list(A, B), mean)
abv <- matrix(agg.abv[, 3], ncol=nlevels(A), dimnames=list(levels(A), levels(B)))
agg.bel <- aggregate(below, by=list(A, B), mean)
bel <- matrix(agg.bel[, 3], ncol=nlevels(A), dimnames=list(levels(A), levels(B)))
barplot(abv, beside=T, ylim=c(0, 1))
barplot(abv, beside=T, ylim=c(1, 0))
しかし、3つの要因はどうですか?
私はトレリスプロットシステムにまったく慣れていないので、どんなアイデアでも大歓迎です
ではごきげんよう
マルコ