[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

From |
"Gresch,Cornelia" <[email protected]> |

To |
<[email protected]> |

Subject |
st: extracting rownames of a matrix into a variable (preparatory work for MICE) |

Date |
Thu, 18 Sep 2008 20:52:58 +0200 |

Dear members, I want to extract a matrix into data-format in a way that the row-names are not lost but appear in a new variable as "observations" (contents of the first variable in string format). Does anybody of you have an idea how to get this running? Background of the question: I work on an imputation model with MICE (command ice). Since the dataset is very large I need to decide which variables to use as independent variables. Therefore I want to create a matrix including the R� (or LL) of bivariate regression models for all possible combination of the variables. The col-names of this matrix correspond to the dependent variables, the row-names of the matrix to the independent variables of the bivariate regression model. Because the final matrix will be very large I would like to convert it into variables (which can be easily used to define about 15-25 independent variables which again I want to use as predictors in the final imputation model).� Using�the command "svmat" results in the expected dataset but the row-names (which are necessary to identify the appropriate predictors) get lost. Below you can see the code and the resulting problem in more detail: (N.B.: variables are "e3his" "e220ex" "t_kfts" (which are metric); "tr10emp" "e121a" (which are categorial with corresponding dummies (tr10emp_1, tr10emp_2,...)); and "slspaet" "tr06sex" (which are binary)) matrix drop _all foreach av of var e3his e220ex t_kfts { � foreach katuv of any tr10emp e121a { ����� local `katuv' = "`katuv'_*" � } ����� matrix `av' = (1) �foreach uv of any e3his e220ex t_kfts /* metrische UVs */ /// ����� "`tr10emp'" "`e121a'"� /* polytome UVs */ /// ����� slspaet tr06sex /* Dummy-UVs */ { � ��� quietly reg `av' `uv' � ��� local r2 = "Fehler" ����� local r2 = e(r2) ����� matrix input R =(`r2') ����� matrix colnames `av' = `av' ����� matrix rownames R = `uv' ����� matrix `av' = (`av' \R) � }�� ����� if ("`av'" == "e3his") matrix ges_ols = (`av') ����� if ("`av'" ~= "e3his") matrix ges_ols = (ges_ols, `av') } matrix list ges_ols ges_ols[8,3] �������������� e3his���� e220ex���� t_kfts ������ r1��������� 1��������� 1��������� 1 ��� e3his��������� 1� .04320558� .05278807 �� e220ex� .04320558��������� 1� .01621444 �� t_kfts� .05278807� .01621444��������� 1 tr10emp_*� .16007611� .02387149� .16865035 ������ r1��������� 0��������� 0��������� 0 � slspaet� .01107509� .00863178� 4.931e-06 � tr06sex� 1.640e-06� .00104972� .00829559 If I convert this matrix to variables I get the following dataset: svmat ges_ols, names(col) �list ���� +--------------------------------+ ���� |��� e3his���� e220ex���� t_kfts | ���� |--------------------------------| � 1. |������� 1��������� 1��������� 1 | � 2. |������� 1�� .0432056�� .0527881 | � 3. | .0432056��������� 1�� .0162144 | � 4. | .0527881�� .0162144��������� 1 | � 5. | .1600761�� .0238715�� .1686504 | ���� |--------------------------------| � 6. |������� 0��������� 0��������� 0 | � 7. | .0110751�� .0086318�� 4.93e-06 | � 8. | 1.64e-06 ��.0010497�� .0082956 | ���� +--------------------------------+ But I also need a variable with the names of the variables in string format as I have it in the matrix (at least that's the only way I see to identify the 15-25 best fitting predictors). Does anybody of you have an idea how to get this running? Alternatively I tried to put the names as strings in the matrix, which also does not work and I also cannot transform the matrix (rows to columns and the other way around) and extract it that way since some of the independent variables are dummysets �(e.g. tr10emp_*) and therefore cannot be identified as variable name. Suggestions are more than welcome! Thanks and best regards from Europe, Germany, Berlin Cornelia ���������� * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:

- Prev by Date:
**Re: st: displaying only one legend with -graph combine-** - Next by Date:
**Re: st: test for structural break in panel** - Previous by thread:
**st: displaying only one legend with -graph combine-** - Next by thread:
**st: RE: extracting rownames of a matrix into a variable (preparatory work for MICE)** - Index(es):

© Copyright 1996–2024 StataCorp LLC | Terms of use | Privacy | Contact us | What's new | Site index |