行列の係数と相関 pvalue のペアごとの組み合わせを列ごとに計算したいと考えています。
ここでは、次の 2 つの関数を使用します。
allCoef<- function(Y,X) { lm(Y~X+0)$coef }
allCorr.p<- function(Y,X) { cor.test(Y,X)$p.value }
たとえば、私は A の行列を持っています:
A= matrix(sample(1:100,16),4,4)
apply(Y=A,2,allCoef,X=A)
正常に動作します。
apply(Y=A,2,allCorr.p,X=A)
ただし、 にエラーが表示されcor.test.default(Y, X) : 'x' and 'y' must have the same lengthます。誰かが私がここで間違ったことを教えてもらえますか? 同じ行列を使用しているため、列の長さは同じでなければなりません。