楼主: jyt1111
2058 2

[问答] 求助,有没有熟悉 ‘PVR’包的啊,有问题请教一下 [推广有奖]

  • 0关注
  • 1粉丝

博士生

57%

还不是VIP/贵宾

-

威望
0
论坛币
1282 个
通用积分
3.9600
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
52460 点
帖子
76
精华
0
在线时间
527 小时
注册时间
2010-3-7
最后登录
2025-7-28

楼主
jyt1111 发表于 2013-10-15 16:52:24 |AI写论文
20论坛币
各位大侠,有没有熟悉PVR(phylogenetic eigenvectors regression)包的啊?本人是R初学者,最近想做phylogenetic eigenvectors regression分析,对于程序包给出的应用实例看不明白,代码如下:library(ape)
library(PVR)
tree <- rcoal(10)
#Decomposing phylogenetic distance matrix derived from tree into a set of orthogonal vectors
x <- PVRdecomp(tree)
trait <- runif(10)
y <- PVR(x, trait = trait, method = "moran")
str(y)
其中我自己有物种之间的遗传距离数据,也就是不需要构建进化树进行分解,我不知道这种情况下,该怎么去进行数据分析,请知道的大侠不吝赐教,先谢谢了



关键词:有没有 Eigenvectors regression分析 eigenvector regression library 初学者 程序

沙发
DM小菜鸟 发表于 2015-1-7 18:29:39
一点一点来吧——

library(PVR)

tree <- rcoal(10)   generate trees by randomly clustering the tips

#Decomposing phylogenetic distance matrixderived from tree into a set of orthogonal vectors

x <- PVRdecomp(tree) 分解phylogenetic距离矩阵(计算基于phylogenetic)成一组正交特征向量

trait <- runif(10) generates random deviates

y <- PVR(x, trait = trait, method ="moran") The phylogenetic eigenvector regression (PVR) starts by performing aneigendecomposition of a pairwise double-centered phylogenetic distance matrixbetween species. The eigenvectors (representing the traits under analysis)estimated values express phylogenetic trends in data and residuals expressindependent evolution of each species.

X: n object of class PVR

trait :A vector, data frame or matrix that contains traits sets (for dataframes and matrices, each column must represent a trait set)

method Character string. A name for the eigenvectors selection method. It can be "moran", "stepwise", "psr" or "sequential".

   
str(y) 显示y的结构

藤椅
重于珠玉 发表于 2020-7-14 15:59:49 来自手机
jyt1111 发表于 2013-10-15 16:52
各位大侠,有没有熟悉PVR(phylogenetic eigenvectors regression)包的啊?本人是R初学者,最近想做phylog ...
您好,我想用pvr程序包做气候因子和系统发育对植物物候期的影响,请问您会用pvr进行操作吗,例子里的环境变量就是输入单纯的一个向量,可以输入矩阵吗?

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-20 17:21