让我们假设以下向量:
x = c( 0.5, 0.4, 0.8 )
其中x [1]和x [2]值与相关矩阵相关:
x[1] x[2] x[3]
x[1] 1 0.8 0
x[2] 0.8 1 0
x[3] 0 0 1
我想计算x的平均值,但要考虑相关性。
我尝试使用lm()的广义最小二乘法,但这意味着使用水平函数,而lm()不喜欢poly(x,0)。我寻找使用用户定义的函数,但它应该返回要拟合的参数...
作为一个具体的例子,让我们从进化树中取三个物种:
library(ape)
## The evolution tree
t=rtree(3)
## Plot it, you notice that two are closer to each other than the 3rd one
plot(t)
## Correlation matrix
vcv.phylo(t,corr=T)
t1 t3 t2
t1 1.0000000 0.4019544 0
t3 0.4019544 1.0000000 0
t2 0.0000000 0.0000000 1
欢迎提示!
答案 0 :(得分:0)