I think you can do the following to avoid using glsaccum and large number of
matrix.
Suppose m is a very big positive number, such that m+e will be all positive,
then do:
gen double _a1=a1*sqrt(m+e)
gen double _a2=a2*sqrt(m+e)
gen double _a3=a3*sqrt(m+e)
mat accum A = _a1 _a2 _a3, nocons
mat accum B = a1 a2 a3, nocons
mat C=A-m*B
will be the one you need.
-----Original Message-----
From: Alfonso Miranda [mailto:ecrgw@csv.warwick.ac.uk]
Sent: Tuesday, October 01, 2002 2:25 PM
To: statalist@hsphsun2.harvard.edu
Subject: st: matrix algebra
Dear statalisters,
I need to calculate a matrix product that has the following form
B = G'YG
where G is a n x k matrix and Y is a n x n diagonal matrix. My problem is:
(a) G should be formed with information for some variables, say, a1 a2 a3. I
have n observations for each variable.
(a) I have the information of the diagonal of Y in variable called e. So I
need to use e to create matrix Y.
(b) n is very large (more than 20,000 observations) so it is imposible to
create a nxn diagonal matrix with elements in the diagonal entries in e.
I guess I need to use something like:
matrix glsaccum B = a1 a2 a3, glsmat(Y)
But first I must define Y. Can anyone help me in this one?
Many Thanks
Alfonso Miranda.
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/