香农指数和辛普森指数的意义(辛普森指数和香农指数)

Alpha多样性指数求标准差及平均值

现在使用高通量测序技术来对微生物群落来进行测序分析已经发展很多年了,在论文中经常有使用到Alpha多样性来对微生物群落进行组内差异分析,所以今天先来给大家介绍一下如何使用单因素方差来对Alpha多样性指数求显著性差异。 今天内容比较简单,但是对我来说比较有用。

1. Alpha多样性指数概述

Alpha多样性(α多样性):指某个群落或生境内部的物种多样性,主要关注群落内的物种多样性。

Alpha多样性主要与两个因素有关:

样本中物种种类数目,即丰富度(Richness)。群落丰富度(Community richness)的指数主要包括Chao1指数和ACE指数​ 样本中各个种的相对密度,即群落中个体分配上的均匀度(Evenness)。

群落多样性(Community diversity)的指数,包括Shannon指数和Simpson指数。ACE指数:利用稀有物种估算物种多样性的指数,其值越高代表群落物种种类越丰富。

ACE是用来估计群落中含有OTU 数目的指数。

Chao1 指数:在生态学中作为度量物种丰富度的指标,其值越高代表群落物种越丰富。

Simpson指数:用来估算样本中微生物多样性指数之一,指数值越大,说明群落多样性越低。但是,现在常用演变而来的Gini-Simpson指数代表原来的辛普森指数,即用1减去经典辛普森指数的数值后得到,此时Gini-Simpson指数随着丰富度的增加而增加。

香农指数(Shannon index):用来估算样本中微生物多样性指数之一。值越大,说明群落多样性越高

今天就先简单说一下这4个指数,以后应该会细说,网上也很多资源。

2. Alpha指数的数据求平均值和方差

2.1数据样本解释说明

求平均值和方差很简单且有很多的方法,感觉接触最多的应该就是EXCEL来处理,那么用R语言应该怎么做呢?

这里用到的是一组Alpha多样性指数,这组数据有6个样本(sample1 到sample6),每个样本具有3个重复,一共有18个样本进行高通量测序,高通量测序样本名称为(A1-A18)。

Alpha多样性指数:

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

group分组信息表:

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

2.2 加载数据及样本预处理

setwd("C:/Users/shanpengloveforever/Desktop/图/微信/alpha") #设置工作目录
data<-read.table("Alpha.txt",header=T,sep="t",row.names=1) 加载Alpha数据表
data <- data[,2:5] #  取出四种指数
#加载group分组信息表
group  <- read.table('group.txt', sep = 't', header = TRUE, stringsAsFactors = FALSE, check.names = FALSE)
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

#添加样本名
data$sample<- factor(rownames(data), levels = rev(rownames(data)))
#合并两个表格
data <- merge(data, group, by = 'sample')
#将site转化为因子变量
data$site <- factor(data$site)
str(data) #查看数据类型
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

2.3 计算ACE指数的标准差及平均值

因为本人目前能力有限,所以大部分代码还是使用重复手段来解决,这里以ACE指数为例子来讲解使用R语音如何计算ACE指数的标准差及平均值。如果有小伙伴会简化流程的,欢迎交流。

#计算ACE平均值
ACE_mean <- aggregate(data$ACE, by = list(data$site), FUN = mean)
names(ACE_mean) <- c('sample','ACE_mean' ) #修改列名
#保留两位小数
ACE_mean$ACE_mean <- sprintf("%0.2f", ACE_mean$ACE_mean)
#计算ACE标准差
ACE_sd <- aggregate(data$ACE, by = list(data$site), FUN = sd)
names( ACE_sd)<-  c('sample','ACE_sd' ) #修改列名
ACE_sd$ACE_sd <- sprintf("%0.2f", ACE_sd$ACE_sd)
#合并平均值和标准差数据
ACE <- merge(ACE_mean, ACE_sd, by = 'sample')
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

library(tidyr) 
#将添加列mean±sd
#ACE1 <- tidyr::unite(ACE,"library(tidyr) 
#将添加列mean±sd
#ACE1 <- tidyr::unite(ACE,"ACE_mean ± ACE_sd", ACE_mean, ACE_sd,sep = "±")#两种方法都可以,这个是好朋友振哥告诉我的
ACE <- unite(ACE, "ACE_mean ± ACE_sd", ACE_mean, ACE_sd, sep = "±", remove = FALSE)
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

到这里我们就成功的将ACEmean ± ACEsd给和在一起了,一般论文也是用这个。

2.4 计算其他多样性指数的标准差和平均值

其实方法和上面的一样,你可以不用看,我这里只是作为一个参考。

#计算chao1标准差和平均值
chao1_mean <- aggregate(data$Chao1, by = list(data$site), FUN = mean)
names(chao1_mean) <- c('sample','chao1_mean' ) #修改列名
chao1_mean$chao1_mean <- sprintf("%0.2f", chao1_mean$chao1_mean)

chao1_sd <- aggregate(data$Chao1, by = list(data$site), FUN = sd)
names( chao1_sd)<-  c('sample','chao1_sd' ) #修改列名
chao1_sd$chao1_sd <- sprintf("%0.2f", chao1_sd$chao1_sd)
#合并两个表格
chao1 <- merge(chao1_mean, chao1_sd, by = 'sample')
chao1 <- unite(chao1, "chao1_mean ± chao1_sd", chao1_mean, chao1_sd, sep = "±", remove = FALSE)

#计算Simpson标准差和平均值
Simpson_mean <- aggregate(data$Simpson, by = list(data$site), FUN = mean)
names(Simpson_mean) <- c('sample','Simpson_mean' ) #修改列名
Simpson_mean$Simpson_mean <- sprintf("%0.2f", Simpson_mean$Simpson_mean)

Simpson_sd <- aggregate(data$Simpson, by = list(data$site), FUN = sd)
names(Simpson_sd)<-  c('sample','Simpson_sd' ) #修改列名
Simpson_sd$Simpson_sd <- sprintf("%0.2f", Simpson_sd$Simpson_sd)
#合并两个表格
Simpson <- merge(Simpson_mean,Simpson_sd, by = 'sample')
Simpson <- unite(Simpson, "Simpson_mean ± Simpson_sd", Simpson_mean, Simpson_sd, sep = "±", remove = FALSE)

#计算Shannon标准差和平均值
Shannon_mean <- aggregate(data$Shannon, by = list(data$site), FUN = mean)
names(Shannon_mean) <- c('sample','Shannon_mean' ) #修改列名
Shannon_mean$Shannon_mean <- sprintf("%0.2f", Shannon_mean$Shannon_mean)

Shannon_sd <- aggregate(data$Shannon, by = list(data$site), FUN = sd)
names(Shannon_sd)<-  c('sample','Shannon_sd' ) #修改列名
Shannon_sd$Shannon_sd <- sprintf("%0.2f", Shannon_sd$Shannon_sd)

#合并两个表格
Shannon<- merge(Shannon_mean,Shannon_sd, by = 'sample')
Shannon <- unite(Shannon, "Shannon_mean ± Shannon_sd", Shannon_mean, Shannon_sd, sep = "±", remove = FALSE)
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

香农指数和辛普森指数的意义(辛普森指数和香农指数) 

2.5 将四个多样性指数合并并输出

Alpha1<- merge(ACE,chao1, by = 'sample')
Alpha1<- merge(Alpha1,Simpson, by = 'sample')
Alpha1<- merge(Alpha1,Shannon, by = 'sample')

write.table (Alpha1, file ="Alpha_1.csv",sep =",", quote =FALSE) #将数据导出
香农指数和辛普森指数的意义(辛普森指数和香农指数) 

抱歉了,由于时间问题,今天就先讲Alpha多样性指数求标准差及平均值,下一期带来使用单因素方差求Alpha多样性指数的显著性差异。

本文内容由互联网用户自发贡献,该文观点仅代表作者本人。爱坤分享仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 3252784850@qq.com 举报,一经查实,本站将于三个工作日内删除。文章链接:https://www.sfyjdwx.com/n/53316.html

(0)

相关推荐

  • 杜牧写江南春的背景故事(杜牧写江南春的背景资料)

    千里莺啼绿映红,水村山郭酒旗风。 南朝四百八十寺,多少楼台烟雨中。 这一首七绝《江南春》,是唐代大诗人杜牧的抒情诗代表作。此诗既有清新明快之调,又有沉郁顿挫之风,可谓身兼李、杜诗风之长。   全诗不发表议论,作者只是以丹青妙笔,自然点染出一幅江南春日图景,通过前后景物的反差,引发读者思考与回味。因为行文俊逸疏朗,表达虚实相生、凝练自然,历来为人称道…

    2022-11-25
    24
  • memtest86 教程(memtest86 v4.20怎么用)

    随着数字化的发展,现在市面上出现了很多内存检测工具,它们可以检测出内存的稳定度以及测试记忆的储存与检索资料的能力。国内使用最多的就是memtest。今天,我就将memtest的使用方法分享给你们 在计算机中,内存扮演着极其重要的角色,它影响着电脑的稳定性,因为所有程序的运行都是在内存中运行的。言归正传,下面,我就给大家介绍一下内存检测工具memtest的使用…

    2022-11-29
    30
  • 成化年制是什么年代(成化年制四字款真假)

    清代康熙、雍正、乾隆三朝,景德镇盛行仿成化御窑瓷器,御窑厂和民间窑场均有仿制。特别是由于这三朝皇帝对成化御窑瓷器均表现出浓厚兴趣,致使当时所仿成化御窑瓷器不但产量大,而且水平也高,有的仿品更达到使人真赝难辨的水平。 仿古三巨头——从左至右:乾隆、雍正、康熙 后仿成化御窑瓷器主要是指造型、纹饰、年款均模仿原作的一类仿品,如康熙朝仿成化斗彩鸡缸杯、雍正朝仿成化斗…

    2022-12-10
    41
  • 艾青作品有哪些(艾青作品风格)

    1.具有独特的意象和主题,其中心意象是太阳和土地,主题则是爱国主义。 艾青诗歌的中心意象是土地和太阳,而主题则是爱国主义。据统计,艾青首诗中,借土地来激发诗人情绪的诗,约占26%,全面直接抒写太阳及其边缘类的诗约占10%,且这类诗多是大篇幅的制作,如《大堰河——我的保姆》《我爱这土地》《雪落在中国的土地上》《向太阳》《吹号者》等。   土地的意象凝…

    2022-12-27
    9
  • 橙色怎么调出来的颜色(橙色怎么调出来的比例)

    如何使用不同的红色+黄色,混合得到鲜艳的橙色?   此练习其目的是加深对颜色混合、颜色偏见,和对颜料的局限性的理解。   橙色是什么颜色的?   橙色是二次色。 二次色即“间色”,是由红、黄、蓝三原色中任意两个原色调配成的色相。 由两种原色按不同比例,可调配多种二次色。   要混合橙色,需要混合黄色+红色。 在混合颜色…

    2023-01-26
    2
  • 关于描写秋雨的文章(名人描写秋雨的文章)

    秋 雨 作者:张爱玲 诵读:李进 来自四季风 雨,像银灰色粘湿的蛛丝,织成一片轻柔的网,网住了整个秋的世界。天也是暗沉沉的,像古老的住宅里缠满着蛛丝网的屋顶。 那堆在天上的灰白色的云片,就像屋顶上剥落的白粉。在这古旧的屋顶的笼罩下,一切都是异常的沉闷。   园子里绿翳翳的石榴,桑树,葡萄藤,都不过代表着过去盛夏的繁荣,现在已成了古罗马建筑的遗迹一样…

    2022-11-19
    11
  • 美白牙齿偏方有用吗(美白牙齿偏方有效果吗)

    夏天到了,小仙女们都裹上了防晒衣,带上了防晒帽,撑起了防晒伞,一切都是为了美白。但是爱美的小仙女们,对美白的追求已经从皮肤延伸到牙齿上啦。为了能拥有一口白牙,小仙女们各出奇招,有人买牙膏一定要买有美白效果的,也有人买了各种各样的美白牙贴,网上更是流传了很多牙齿美白的小妙招,那这些小妙招到底是好是坏呢?今天就来给大家揭晓一下。   1. 牙齿美白小妙…

    2022-11-21
    9
  • 怎么让安卓手机变砖(安卓手机变砖头怎么办)

    29日,复旦大学计算机学院教授、国家973首席科学家杨珉在微博发文,称自己和同事们于去年9月向谷歌提交了400多个漏洞,16个月后谷歌安全团队终于修复完毕。   随后他又补充说,原本谷歌安全团队自称能够2个月修复,实际却用了整整16个月才修复完成。 原文中,杨教授说这些漏洞可怕到能“让市面所有活着的安卓设备变砖头”,所有使用安卓代码的厂商都将受到这…

    2022-11-26
    16
  • QQ昵称个性(Qq昵称大全)

    晚初羡晚风捎信子雾尽风暖让风诠释爱与猫秋人间荒糖玫瑰典藏家未曾见山和海秋月星海 云边有只羊 晚颂海鲸 终将靠岸 下至人间白鹭柒巷子南边晨雾的鹿梦是彼方千岁合时九念于你云霭暮山向风和盛夏向晚意不适司月十渡绾挽卿 雾月的信 落日還在等可抵岁月漫长风季是过与你携手观海陪北看日落温酒故人归浪漫会有期风亭离晚故里长安下 写半截诗 妄想贩卖日落 屋外玫瑰遍地海和浪漫九月…

    2022-11-25
    12
  • 重庆有什么好玩的地方景点推荐(重庆有什么旅游的地方)

    重庆市旅游景点有:洪崖洞、长江索道、武隆喀斯特旅游区、大足石刻、重庆云阳龙缸国家地质公园、金佛山等。 重庆,地处中国内陆西南部,是西南部重要的城市,原来属于四川管辖,后来成为直辖市,目前是是中华人民共和国省级行政区、直辖市、国家中心城市、超大城市 ,国务院批复确定的国家重要的中心城市之一、长江上游地区经济中心、国家重要先进制造业中心、西部金融中心、西部国际综…

    2022-12-29
    2