R计算 PD Calculate Faith’s Phylogenetic Diversity

网址: https://cran.r-project.org/web/packages/picante/index.html

Description: Calculate the sum of the total phylogenetic branch length for one or multiple samples.

Usage: pd(samp, tree, include.root=TRUE)

Arguments

samp: Community data matrix

tree:A phylo tree object

include.root: Should the root node be included in all PD calculations (default = TRUE)

#加载包
library(picante)

library(ape)

#测试数据
data(phylocom)

pd(phylocomsample,phylocomsample, phylocomsample,phylocomphylo)

  1. 读取树文件:
    tree= read.tree(“D:/timberproject2021/sp_tree20210707.tre”)

  2. 读取samp数据矩阵
    b = read.csv(“D:/timberproject2021/GridsID-species2.csv”, header = T, row.names=1)

  3. 计算pd
    c = pd(b, tree, include.root=TRUE)

  4. 导出计算结果
    write.csv(c,file = “D:/aaa.csv”)

报错:修改树:
If youwant the root ofyour tree to correspond to the most recent ancestor of the taxaactually present in your sample, you should prune the tree before running pd:

prunedTree<-prune.sample(b,tree)

参考资料:

https://daijiang.name/en/2014/05/04/notes-func-phylo-book-1/

原文网址:https://www.jianshu.com/p/789a93c2f62a

R—计算系统发育多样性PD (Calculate Faith’s Phylogenetic Diversity)相关推荐

  1. R计算功能多样性— functional diversity (FD)

    一般植物功能特征被划分为3类:一是植物形态特征, 包括生长型.生活型.植株高度等:二是植物生殖特征,包括传粉方式.扩散方式.种子重量等:三是植物生理特征, 如植物固氮能力等[30].--参考 为研究中 ...

  2. C语言学习之购房从银行贷了一笔款d,准备每月还款额为p,月利率为r,计算多少月能还清。

    购房从银行贷了一笔款d,准备每月还款额为p,月利率为r,计算多少月能还清. 设d为300 000元,p为6000元,r为1%.对求得的月份取小数点后一位,对第2位按四舍五入处理. #include & ...

  3. R 计算均方差MSE(mean squared error)

    本文介绍MSE(均方差),并使用两种R方法实现. MSE(均方差) 判定预测模型的准确度的常用方法是均方差MSE( mean squared error).计算公示为: MSE = (1/n) * Σ ...

  4. R计算贝叶斯信息指标(Bayesian Information Criterion,BIC)实战

    R计算贝叶斯信息指标(Bayesian Information Criterion,BIC)实战 目录 R计算贝叶斯信息指标(Bayesian Information Criterion,BIC)实战

  5. R计算获取决策曲线数据(Decision Curve Analysis,DCA)并使用python进行可视化

    R计算获取决策曲线数据(Decision Curve Analysis,DCA)并使用python进行可视化 目录 R计算获取决策曲线数据(Decision Curve Analysis,DCA)并使 ...

  6. R计算数组的累计加和(Cumulative Sums)

    R计算数组的累计加和(Cumulative Sums) 目录 R计算数组的累计加和(Cumulative Sums) R中的累加和 可视化R中的累加和

  7. R计算混淆矩阵(Confusion Matrix)

    R计算混淆矩阵(Confusion Matrix) 目录 R计算混淆矩阵(Confusion Matrix) 拟合Logistic回归模型 创建混淆矩阵

  8. R计算赤信息指标(Akaike information criterion,AIC)

    R计算赤信息指标(Akaike information criterion,AIC) 目录 R计算赤信息指标(Akaike information criterion,AIC) 计算AIC AIC结果 ...

  9. R计算回归模型Mallows’ Cp指标

    R计算回归模型Mallows' Cp指标 目录 R计算回归模型Mallows' Cp指标 R中Mallows'Cp的计算 Mallow's Cp指标要点

最新文章

  1. 清华孵化,打造智能驾驶感控平台,超星未来获千万A轮融资
  2. Socket,SocketImpl与SocketImplFactory的关系
  3. Boost TCP serverclient 有回调无发送
  4. Java批量完成对文件夹下全部的css与js压缩,利用yuicompressor
  5. java 实体类 临时注解_JPA:Java持久层API--配置流程
  6. css案例学习之id要唯一
  7. 大型综合体弱电智能化解决方案标书
  8. Pycharm安装包(类库)的方法总结及解决包下载慢的问题
  9. 前端使用canvas拼接多张图片
  10. Selenium WebDriver使用IE浏览器
  11. 【zznu-夏季队内积分赛3-G】2333
  12. 终端模拟器常用快捷键
  13. eclipse中转换项目编码时没有GBK的问题
  14. python账号怎么注销_登录后如何注销
  15. HTML怎么跟随页面缩放,如何让网页跟着 浏览器全比例缩小(示例代码)
  16. ISP概述、工作原理及架构
  17. Texlive配置华文新魏字体
  18. https自签证书tls握手时错误或go系统错误处理
  19. Google advertiser api开发概述——最佳做法建议
  20. JVM中对象的分代年龄

热门文章

  1. log4j 源码解析_log4j1.x设置自动加载log4j.xml
  2. linux新增ssh端口80,Linux(CentOS 7) 新增或修改 SSH默认端口
  3. Ubuntu下如何用Empathy上QQ
  4. 自编码clone图解
  5. Hey, 看看小程序的page-frame.html把~
  6. Java基础-等待子线程结束
  7. c++,不能声明为虚函数的函数
  8. 用dmidecode - 查看硬件信息
  9. (转)oracle 11g安装后用户名忘记怎么办
  10. JavaScript求当月天数