让我们假设以下向量:
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
欢迎任何提示!