3 # (Slower, but trusted) R version of Moments_M3
4 .Moments_M3_R = function(X,Y)
12 M31=M3_1=M3_2=M3_3 = tensor(tensor(v,v),v)
14 M31 = M31 + (Y[i]*tensor(tensor(X[i,],X[i,]),X[i,]))
19 l1=l2=l3 = tensor(tensor(v,v),v)
22 l1 = l1 + Y[i]*tensor(tensor(e[,j],X[i,]),e[,j])
23 l2 = l2 + Y[i]*tensor(tensor(e[,j],e[,j]),X[i,])
24 l3 = l3 + Y[i]*tensor(tensor(X[i,],e[,j]),e[,j])
34 M3 = M31-(M3_1+M3_2+M3_3)
38 test_that("both versions of Moments_M3 agree on various inputs",
44 X = matrix( runif(n*d,min=-1,max=1), nrow=n )
45 Y = runif(n,min=-1,max=1)
47 M3_R = .Moments_M3_R(X,Y)
48 expect_equal(max(abs(M3 - M3_R)), 0)