【问题标题】:Can I extract certain dissimilarity values from a bray curtis matrix in R?我可以从 R 中的布雷柯蒂斯矩阵中提取某些相异值吗?
【发布时间】:2019-05-09 15:54:00
【问题描述】:

我正在尝试评估入侵物种清除后河流群落的变化持续多长时间。在移除之前和移除之后的几个时间点对社区进行了一系列评估。我想计算社区在移除之前和移除后的每个时间点(即前后 1、前后 2、前后 3 等)之间的差异,然后跨时间分析这些值以确定多长时间溪流社区需要恢复到搬迁前的条件。预计流参数的差异会影响该过程,因此我需要按流来执行此操作。我在 x 轴上描绘了一个时间图,在 y 轴上描绘了不同程度,并且每个流都有不同的函数。有没有办法从布雷柯蒂斯矩阵中提取这些差异值?

这是我的第一个问题,对于任何遗漏的信息或不清楚的地方,我深表歉意。

【问题讨论】:

  • 你能告诉你你做了什么吗?在谷歌上快速搜索也会返回这个 - cc.oulu.fi/~jarioksa/softhelp/vegan/html/vegdist.html 这是你想要的吗?
  • 我已经在数据上运行了 vegdist 以获取 Bray Curtis 值,下面的答案解决了我如何获取特定值的问题。我只是不确定如何确定矩阵中的哪些相异值以知道要提取哪些相异值。

标签: r vegan


【解决方案1】:

如果您将输出从 vegdist() 转换为矩阵,那么您可以轻松地提取样本对之间的差异,使用行名或矩阵的索引来提取您想要的内容。

要转换为矩阵使用

distmat <- as.matrix(bc_obj)

其中bc_obj 是从vegdist() 返回的包含布雷柯蒂斯距离的对象。

【讨论】:

  • 谢谢,这非常有帮助。但是,由于标签不包含在矩阵中,因此我无法确定要提取哪些值。您是否有关于如何构造矩阵的参考资料(即比较日期和地点的哪些组合以判断 1 和 3 之间的差异)或知道如何标记矩阵本身?
  • @AveryScherer 如果您将数据框中的行名传递给vegdist(),当您转换为矩阵时,它们应该被保留并呈现。
猜你喜欢
  • 2021-12-16
  • 1970-01-01
  • 1970-01-01
  • 2021-05-04
  • 1970-01-01
  • 1970-01-01
  • 2021-04-26
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多