Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: Column vector into variable, accounting for -marksample-


From   "Lacy,Michael" <Michael.Lacy@colostate.edu>
To   "statalist-digest@hsphsun2.harvard.edu" <statalist-digest@hsphsun2.harvard.edu>
Subject   st: Column vector into variable, accounting for -marksample-
Date   Sun, 7 Oct 2012 19:45:01 +0000

I have a Stata program that obtains a column vector, say r(X) by calling another program. I want to put the
 values of r(X) into a variable in the original dataset, which for me is complicated by -marksample-. 

The values of r(X) arrive in the same order as the original data would be if cases not selected by 
-marksample- were dropped. So, if were there no use of marksample, the row index of r(X) would 
correspond to _n, and there would be no difficulty; -svmat- will align the contents of r(X) as desired. 
But this is not true because the data are processed -if `touse'-, so that the first element of r(X) might 
correspond to the 4th observation of the original data set, the second element to the 7th observation,
and so forth.  

My thinking so far involves reconstructing the r(X) matrix into another matrix, putting in missing values
as indicated by -touse-, per my fragment below.  Is there a less lunky way to accomplish this? 
(My program might be called many times in a simulation context, so time is an issue.)

This must be an utterly commonly task.  (And yes, I could do the same in Mata, but it's still more 
or less the same thing.)

program whatever
syntax ...[if] ...
marksample touse 

....Other program is called with if `touse', which returns r(X)
// Make r(X) into a variable
X = return(X)
mkmat `touse', matrix(tu)
mat `temp' = J(`=_N',1,.) // working matrix for repacking X
local pos = 1
forval i = 1/`=_N' {
   if (tu[`i',1] == 1) {
      mat `temp'[`i',1] = `D'[`pos',1] 
   	local ++pos
   }	  
}
svmat `temp', name(Xvar)


Thanks,

Mike Lacy
Dept. of Sociology
Colorado State University
Fort Collins CO 80523-1784

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index