| 5 | 1/1 | 返回列表 |
| 查看: 1348 | 回復(fù): 4 | |||
[求助]
請問哪位大神會用R軟件計(jì)算微生物多樣性指數(shù)? 已有1人參與
|
|
請問哪位大神會用R軟件計(jì)算微生物多樣性指數(shù)?我的一些高通量實(shí)驗(yàn)數(shù)據(jù)想做一些特殊的分析,為此需要學(xué)習(xí)一下用R軟件計(jì)算微生物多樣性指數(shù), 請會做的大神把詳細(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 |
金蟲 (著名寫手)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|