当前位置: 代码迷 >> 综合 >> Review: R.MRPP分析
  详细解决方案

Review: R.MRPP分析

热度:97   发布时间:2023-12-16 23:24:15.0

MRPP(Multi ResponsePermutation Procedure)分析,类似于ANOSIM 分析,主要是用于分析高维度数据组间相似性的统计方法。

目的:用于评估两组实验数据的整体相似性,以及相似性是否显著。

执行:R中vegan包提供了MRPP分析的函数mrpp

结果:得到A值、P 值

A > 0,表示组间差异大于组内差异

A < 0,表示组内差异大于组间差异

P值评价组内差异与组间差异的显著性

### 输入文件的格式处理
>head(data)sp1    sp2    sp3    sp4    ... group
sample1    相对丰度...                        A
sample2    ...                               B
sample3    ...                               Alibrary(vegan)data.dist = vegdist(subset(data, select = -group))  ### bray Curtis 计算相似性矩阵
data.mrpp = mrpp(data.dist, grouping = data$group, permutations = 999) ### 得出A值、P值###
data.df = subset(data, select = -group)
data.mrpp = mrpp(data.df, grouping = data$group, permutations = 999,  distance = "bray") ### data 是data.frame,指定distance 参数用于计算相似性矩阵的方法