使用 vegan 包更改颜色并向 nestedtemp（） 图添加标签

问：

我使用 R 包 vegan 来检查数据集中的嵌套性：

bc_rp <- read.csv("BC_RPc.csv", na.strings = c("","NA"), header=TRUE, check.names = FALSE)

bc_rp

bc <- data.frame(
  cell_layer = c("H", "B", "A", "G", "O"),
  Np = c(0, 1, 1, 0, 0),
  Pn = c(1, 1, 0, 0, 1),
  Npn = c(0, 1, 0, 0, 0),
  Npnp = c(0, 1, 0, 0, 0),
  Pnp = c(1, 0, 0, 0, 0),
  Pnpn = c(0, 1, 0, 0, 0)
)


# Renamed the columns to be easier to manipulate and have no special characters
# colnames(hc) <-c('species', 'a','b','c','d','e','f','g')

# Transpose the rows and columns in data. 
bc_t <-t(bc)

out_bc_t <- nestedtemp(bc_t)
out_bc_t
#nestedness temperature: 10.24168 
#with matrix fill 0.4 

plot(out_bc_t, kind="incid")

此代码按预期工作，以生成嵌套温度和绘图。

嵌套温度：10.24168 矩阵填充 0.4

但是，我希望我可以将绘图的颜色更改为红色以外的颜色，并为每种单元格类型和类别（np、pn 等）添加标签我该怎么做？