[R] Matrix mulitplication

[Peter Dalgaard]

(Ted Harding) <Ted.Harding at nessie.mcc.ac.uk> writes:

> Sorry! Missed a trick here:
> 
>     At <- t(A)
>     Bt <- t(B)
>     E  <- B%*%A
>     Et <- t(E)
>     A%*%solve(Et%*%E + C)%*%Et%*%D
> 
> (saves 2 multiplications at the relatively cheap cost of 1 transpose)

Well, you might consider getting rid of the first two transposes since
you're not actually using them for anything....


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907