我正在按珊瑚礁名称绘制鱼类生物量的箱形图,按生物量中位数的顺序排列。所有reefnames (站点)都在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)谢谢
发布于 2014-05-22 23:34:36
使用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")发布于 2014-05-23 03:16:39
下面是一组根据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"

https://stackoverflow.com/questions/23811077
复制相似问题