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 <feanor0@hotmail.com> 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]]

