| 5 | 1/1 | 返回列表 |
| 查看: 1353 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
請問哪位大神會(huì)用R軟件計(jì)算微生物多樣性指數(shù)? 已有1人參與
|
|||
|
請問哪位大神會(huì)用R軟件計(jì)算微生物多樣性指數(shù)?我的一些高通量實(shí)驗(yàn)數(shù)據(jù)想做一些特殊的分析,為此需要學(xué)習(xí)一下用R軟件計(jì)算微生物多樣性指數(shù), 請會(huì)做的大神把詳細(xì)的過程寫給我好嗎?非常感謝! |
高通量測序 |
|
library(vegan) otu=read.table("otu.txt",header=T) otu=otu.matrix(otu)#將數(shù)據(jù)類型轉(zhuǎn)化為矩陣 Shannon.Wiener <- diversity(otu,index="shannon" #計(jì)算香農(nóng)指數(shù)Simpson <- diversity(otu,index="simpson" ![]() Inverse.Simpson <- diversity(otu,index="inv" ![]() S <- specnumber(otu)#計(jì)算物種累計(jì)數(shù) plot(S) J <- Shannon.Wiener/log(S) #計(jì)算Pielou均勻度指數(shù) J |
|
otu=read.table("otu.txt",header=T) otu=otu.matrix(otu) #將數(shù)據(jù)類型轉(zhuǎn)化為矩陣 Shannon.Wiener <- diversity(otu,index="shannon" #計(jì)算香農(nóng)指數(shù)Simpson <- diversity(otu,index="simpson" ![]() Inverse.Simpson <- diversity(otu,index="inv" ![]() S <- specnumber(otu) #計(jì)算物種累計(jì)數(shù) plot(S) J <- Shannon.Wiener/log(S) #計(jì)算Pielou均勻度指數(shù) J |
金蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +6 | 一定上岸aaa 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +10 | 路癡小琪 2026-03-27 | 10/500 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-27 | 3/150 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|