[R] Matrix multiplication

Steven Yen syen04 at gmail.com
Thu Jun 8 05:50:55 CEST 2017


I need to have all elements of a matrix multiplied by a weight before 
being post-multiplied by itself, as shown in the forst block of codes 
below. I can also multiply the matrix by the square root of the weight 
and then take the outer product.

Actually, what I need is this. Denote each row of the matrix by a row 
vector as xi and each element of the weighting vector as wi. Then, I 
need the sum of wi * t(xi) %*% xi over i.

Any suggestion for a compact approach would be appreciated.

set.seed(76543211)
w<-1:10; w
a<-matrix(rpois(20,2),nrow=10); a
b<-a
a<-w*a
t(a)%*%b

set.seed(76543211)
a<-matrix(rpois(20,2),nrow=10); a
a<-sqrt(w)*a; a
t(a)%*%a




On 1/4/2017 5:41 PM, Steven Yen wrote:
> I need help with gls{nlme}.
> Specifically, I am estimating an equation with AR(1) using 
> maximum-likelihood. I am not understanding the correlationoption 
> below. Help appreciated.
>
> ===
> library(nlme)
> eq1<-log(chnimp)~log(chempi)+log(gas)+log(rtwex)+befile6+
>                  affile6+afdec6
> reg1<-gls(eq1,data=mydata,correlation=corAR1(),method="ML",verbose=T)
>

	[[alternative HTML version deleted]]



More information about the R-help mailing list