如何提取adipart/oecosimu对象中的2.5%与97.5%置信区间？

方法说明 #

adipart()返回的是oecosimu类对象，其2.5%和97.5%置信区间并非直接存储为对象的顶级元素，而是基于内置的模拟结果计算得到的。你可以通过以下步骤提取目标向量：

步骤1：提取模拟数据集 #

adipart的模拟结果存储在对象的$simulations属性中，这是一个矩阵——每行对应一次模拟，每列对应一个多样性分区组分：

sim_data <- ADP.mite$simulations

步骤2：计算置信区间向量 #

使用apply()函数对每个组分的模拟结果计算分位数，得到2.5%和97.5%置信区间的上下限向量：

# 提取2.5%置信区间下限
ci_lower <- apply(sim_data, 2, quantile, probs = 0.025, na.rm = TRUE)

# 提取97.5%置信区间上限
ci_upper <- apply(sim_data, 2, quantile, probs = 0.975, na.rm = TRUE)

验证一致性 #

你可以将计算结果与ADP.mite的打印输出对比，确认置信区间完全匹配：

# 查看原始adipart输出
ADP.mite

# 查看提取的置信区间向量
ci_lower
ci_upper

整合到其他数据框 #

直接将这两个向量作为新列添加到目标数据框即可：

# 假设已有数据框target_df，包含对应多样性组分的列
target_df$ci_2.5 <- ci_lower
target_df$ci_97.5 <- ci_upper

内容的提问来源于stack exchange，提问作者AugieRockDoc