st: get predictions using mata

[adiallo5@worldbank.org]

Danielle ,
I  am not a mata specialist but for the second question,
I think you should use st_matrix.
Amadou.


-----------------------------------------------------------
Danielle wrote:


                                                                                 
                                                                                 
                                                                                 
 Dear Statalisters,                                                              
                                                                                 
 I am using -arima- with a panel dataset (possible if you specify the            
 "condition" option). But, the postestimation commands don't work, so I am       
 attempting to compute predictions and MSE myself, using mata, which I           
 am just beginning to learn.                                                     
                                                                                 
 I wrote the following little function to simulate getting predictions:          
                                                                                 
 --start--                                                                       
 real matrix test(real rowvector b, real matrix X)                               
 {                                                                               
              real matrix y                                                      
              for(i=1; i<=cols(X); i++) {                                        
                          real scalar yi                                         
                          yi = b * X[i, .]'                                      
              }                                                                  
              y = (y \ yi)                                                       
              return(y)                                                          
 }                                                                               
 --end--                                                                         
                                                                                 
 When I run:                                                                     
 b = (1, 2, 3, 4)                                                                
 X = (1, 2, 3, 4 \ 5, 6, 7, 8 \ 9, 10, 11, 12 \ 13, 14, 15, 16)                  
 test(b, X)                                                                      
                                                                                 
 I get the following error msg:                                                  
 : test(b, X)                                                                    
                   test():  3200  conformability error                           
                  <istmt>:     -  function returned error                        
 r(3200);                                                                        
                                                                                 
 I suspect that part of the problem is the "i" subscripting. In any              
 case, can someone please offer advice on how to fix this? Or, even better,      
 has anyone already written something to generate predictions?                   
                                                                                 
 Also, after I run -arima-, and get the coefficient matrix e(b), how do          
 I get this to be a mata matrix?                                                 
                                                                                 
 Thank you,                                                                      
 Danielle Ferry                                                                  
                                                                                 
 ------------------------------------------------------------------------------- 

*
*   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/