计算得到Stata中矩阵的值
我在Stata数据集中有一个变量计算得到Stata中矩阵的值,stata,Stata,我在Stata数据集中有一个变量age，13个变量x1到x13，以及802个观察值年龄的值范围为1到9x1到x13的值范围为1到13 我想知道怎样数一的数字。。根据年龄的不同值，从x1到x13中的13个。例如，对于age1，在x1到x13中，计算1,2,3,4 我首先通过使用将x1更改为x13作为矩阵 mkmat x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13, matrix (a) 然后，我想使用以下循环进行计数： gen count = 0 quie
agex1x13x1x13x1x13agex1x13x1x13mkmat x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13, matrix (a)
gen count = 0
quietly forval i = 1/802 {
quietly forval j = 1/13 {
replace count = count + inrange(a[r'i', x'j'], 0, 1), if age==1
}
}
我失败了。我还是有点不确定你想要实现什么。但如果我理解正确的话，这里有一个方法 首先，一个简单的数据，其
agex1-x4clear
input age x1 x2 x3 x4
1 5 6 6 6
1 7 5 6 5
2 5 7 6 6
3 5 6 7 7
3 7 6 6 6
end
x1-x4forval i=5/7 {
egen n`i'=anycount(x1 x2 x3 x4),v(`i')
}
+----------------------------------------+
| age x1 x2 x3 x4 n5 n6 n7 |
|----------------------------------------|
1. | 1 5 6 6 6 1 3 0 |
2. | 1 7 5 6 5 2 1 1 |
3. | 2 5 7 6 6 1 2 1 |
4. | 3 5 6 7 7 1 1 2 |
5. | 3 7 6 6 6 0 3 1 |
+----------------------------------------+
n5x1-x4forval i=5/7 {
egen n`i'=anycount(x1 x2 x3 x4),v(`i')
}
+----------------------------------------+
| age x1 x2 x3 x4 n5 n6 n7 |
|----------------------------------------|
1. | 1 5 6 6 6 1 3 0 |
2. | 1 7 5 6 5 2 1 1 |
3. | 2 5 7 6 6 1 2 1 |
4. | 3 5 6 7 7 1 1 2 |
5. | 3 7 6 6 6 0 3 1 |
+----------------------------------------+
mat A=J(3,3,.) // age (1-3) and values (5-7)
mat rown A=age1 age2 age3
mat coln A=value5 value6 value7
forval i=5/7 {
forval j=1/3 {
qui su n`i' if age==`j'
loca k=`i'-4 // the first column for value5
mat A[`j',`k']=r(sum)
}
}
value5x1-x4A[3,3]
value5 value6 value7
age1 3 4 1
age2 1 2 1
age3 1 4 3
以Aspen为例，您可以这样做：
gen id = _n
reshape long x, i(id)
tab age x
countmatrix count = J(9, 13, 0)
forval i = 1/9 {
forval j = 1/13 {
forval J = 1/13 {
qui count if age == `i' & x`J' == `j'
matrix count[`i', `j'] = count[`i', `j'] + r(N)
}
}
}
我不太确定你努力取得什么结果。看起来您需要不止一个新的“count”变量。你能举一个简单数据的例子吗？对于当前值，删除“r”和“x”以引用矩阵中的元素。也可以检查
egenmat a=J（3,3，）mat list A