根据R中的功能选择boxplot的颜色

Question

我正按照中间生物量的顺序绘制礁石名称的鱼类生物量的箱线图。所有的珊瑚礁名称（网站）都在MPA之内或之外，例如MPA="1"或MPA=="0"。目前所有地块都显示绿色。

如何将MPA=="0"网站显示为蓝色，将MPA=="1"显示为绿色。同时保持鱼类生物量的顺序。

MPA <- factor(Fish$MPA)
bymedian <- with(Fish, reorder(ReefName, log10(Biomassm+1)), median)
boxplot(log10(Biomassm+1) ~ bymedian, data = Fish,
        xlab = "ReefName", ylab = "Biomassm",
        main = "Biomassm in Caribbean", varwidth = TRUE, 
        col=(c("darkgreen")), las=3, cex.axis=0.3)

谢谢

Answer 1

使用ggplot2包可能更好一点。您的代码将如下所示：

ggplot(data=Fish, aes(x=reorder(ReefName, log10(Biomassm+1)), median), y=Biomassm, fill=MPA)) +
  geom_boxplot() +
  scale_y_log10("Biomassm") +
  xlab("ReefName") +
  scale_fill_manual(values=c("blue", "green")) +
  ggtitle("Biomassm in Caribbean")

Answer 2

这里有一组根据MPA值着色的箱形图：

# generate some data
set.seed(1)
X = matrix(rnorm(100), ncol=10)

# order by median
X = X[,order(apply(X, 2, median))]

# some fake MPA values  
MPA = round(runif(n=10, min=0, max=1))

# generate boxplots and check if MPA==1
boxplot(X, col=ifelse(test=MPA==1, yes='green', no='blue'))

# add legend
legend(x='bottomleft', fill=c('green','blue'), legend=c('MPA=1', 'MPA=0'), inset=c(0.01))

ifelse的输出是根据MPA值的颜色矢量，它们用于为方框着色：

[1] "blue"  "blue"  "green" "blue"  "blue"  "green" "green" "blue"  "blue"  "green"