我有兩個大的數據集,每個數據都有2000+數據並試圖找出每5行的協方差。如何找到R中每n行的協方差
x=c(1,2,3,4,5)
y=c(6,7,8,9,10)
df=data.frame(x,y)
group=rep(1:length(df),each=2,length=length(df))
什麼是我的下一個步驟,所以我可以找到像this`
cov(x[1:2,],y[1:2,])
和
cov(x[3:4,],y[3:4,])
我有兩個大的數據集,每個數據都有2000+數據並試圖找出每5行的協方差。如何找到R中每n行的協方差
x=c(1,2,3,4,5)
y=c(6,7,8,9,10)
df=data.frame(x,y)
group=rep(1:length(df),each=2,length=length(df))
什麼是我的下一個步驟,所以我可以找到像this`
cov(x[1:2,],y[1:2,])
和
cov(x[3:4,],y[3:4,])
library(zoo)
x = c(1,2,3,4,5)
y = c(6,7,8,9,10)
rows = 2
out = rollapply(data.frame(x,y), rows, function(x) cov(x[,1],x[,2]),
by.column=FALSE)
out
有在您的文章,@Ian了一些錯誤的協方差。 'a'沒有定義,'x'和'y'是開頭的向量,但後來他們不是,你有'}'而不是''''。另外,這個問題不是很清楚。 – Suren
對不起。我只是編輯它,是更好的? – Ian