雙總和我想計算這些量到編碼中的R
a12=sum_(i from 1 to m)sum_(j1<j2)(I(X[i]>Y[j1] and X[i]>Y[j2]))
a13=sum_(j from 1 to n)sum_(i1<i2)(I(X[i1]>Y[j] and X[i2]>Y[j]))
其中I是指標函數。
所以我想出了這個R代碼裏面
a12=0; a13=0
for (l in 1:(length(Z1)-1)){
for (m in 1:(length(Z2)-1)){
a12<-a12+(Z1[l]<Z2[m])*(Z1[l+1]<Z2[m])*1
a13<-a13+(Z1[l]<Z2[m])*(Z1[l]<Z2[m+1])*1
} # closing m
} # closing l
a12=a12+sum((Z1[-length(Z1)]<Z2[length(Z2)])*(Z1[-1]<Z2[length(Z2)])*1)
a13=a13+sum((Z1[length(Z1)]<Z2[-length(Z2)])*(Z1[length(Z1)]<Z2[-1])*1)
a12;
a13
不幸的是,這不僅是非常緩慢的,但我沒有得到什麼,我應該得到的。
請問你能幫助我嗎?
感謝,
羅蘭
例子和'Z2'服用
z
下部三角形的總和得到最終結果和你所期望的結果將有助於。 – 2011-03-28 21:26:00也解釋你的雙倍更好。 j1和j2只有在你在a12中使用它們之後才被定義**。你到底想要做什麼?鏈接到一張紙或更好的公式也會有所幫助。 – 2011-03-28 21:45:30
@Joris我認爲'Sum_(j1
2011-03-29 02:22:38