Hi Murali, Just an idea, probably not the best :
x<-1:4 y<-1:6 z<-matrix(1:(length(x)*length(y)),nrow=length(x)) I <- matrix(1,nrow=length(x),ncol=length(y)) I[row(I)==col(I)] <- 0 sum (outer (x, y, '*') * I) sum (outer (x, y, '*') * z * I) Hope this helps, Nael On Wed, Jul 2, 2008 at 6:30 PM, Murali Menon <[EMAIL PROTECTED]> wrote: > > folks, > > is there a clever way to compute the sum of the product of two vectors such > that the common indices are not multiplied together? > > i.e. if i have vectors X, Y, how can i compute > > Sum (X[i] * Y[j]) > i != j > > where i != j > > also, what if i wanted > > Sum (X[i] * Y[j] * R[i, j]) > i != j > > where R is a matrix? > > thanks, > > murali > > > _________________________________________________________________ > Enter the Zune-A-Day Giveaway for your chance to win day after day after > day > > M_Mobile_Zune_V1 > [[alternative HTML version deleted]] > > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.