Skip to content

matrix multiplication

1 message · arun

#
Hi,
Try this:
#mat1 is the data

res<-do.call(cbind,lapply(seq_len(nrow(mat1)),function(i) {new1<-do.call(rbind,lapply(seq_len(nrow(mat1[-i,])),function(j) {x1<-rbind(mat1[i,],mat1[j,]); x2<-(abs(x1[1,1]-x1[2,1])*abs(x1[1,5]-x1[2,5]))+(abs(x1[1,2]-x1[2,2])*abs(x1[1,6]-x1[2,6]))+(abs(x1[1,3]-x1[2,3])*abs(x1[1,7]-x1[2,7]))+(abs(x1[1,4]-x1[2,4])*abs(x1[1,8]-x1[2,8]))}));new1}))
head(res,3)
#??? [,1]?? [,2]?? [,3]?? [,4]?? [,5]?? [,6]?? [,7]?? [,8]?? [,9]? [,10]
#[1,]????? 0? 79745 146483? 40874 206818 229688 149126 230028? 90188 157418
#[2,]? 79745????? 0 117950 179977 187223 343090 119664 137120 134010? 74077
#[3,] 146483 117950????? 0 123746? 54717 124372 132140? 36383? 70878 112656
-----------------------------------------------------------------------------------
A.K.