-4

让我们假设以下向量:

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)。我寻找使用用户定义的函数,但它应该返回要拟合的参数......</p>

作为一个具体的例子,让我们从进化树中取出三个物种:

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

欢迎任何提示!

4

1 回答 1

0

答案可以在 CERN 论文中找到:

ftp://ftp.desy.de/pub/preprints/cern/ppe/ppe94-185.ps.gz

该过程是广义最小二乘回归。

有关结果,请参见第 (1) 页的公式 (2)。

于 2015-01-06T16:33:32.017 回答