请教一个stata的矩阵计算方法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请教如何写程序生成一个新变量Y2

数据表：

obs	X	Y1	Y2
1	A	4
2	B	5
3	C	9
4	D	6
5	E	3
6	F	8

关联矩阵：

A	B	C	D	E
B	E	F
C	A	D	F
D	B	C

请教如何生成一个新变量Y2

Y2的取值为根据关联矩阵计算数据表中x所对应的Y1的均值

如关联矩阵中第一行A所对应的BCDE的Y1的均值，则为数据表中第一个 样本 的Y2的取值                          

第一个样本的Y2=(5+9+6+3)/4=5.75

根据关联矩阵中第二行B所对应的EF的Y1的均值，则为数据表中第二个 样本 的Y2的取值

第二个 样本 的 Y2=(3+8)/2=5.5

根据关联矩阵 中第三行 C所对应的ADF的Y1的均值，则为数据表中第三个 样本 的Y2的取值

第三个样本的Y2=(4+6+8)/3=6

以此类推

如何用程序实现

非常感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：Stata tata 计算方法 非常感谢 数据表 计算方法 数据表 程序 如何 样本

请高手指点一下，之前想到用矩阵方法算，但是因为关联矩阵都不是对齐的，估计不行。应该是编程做循环，但是感觉太复杂了，还请高手指点一二。万分感激！

这个可以有很多种方法：Mata、matrix矩阵运算，也可以使用普通的Stata数据运算方式。假设你的矩阵叫作mat，下面是是一种方式：
snapshot erase _all
snapshot save,label("main data")
tempfile tmp
drop _all
matrix mat=mat'
svmat mat
drop in 1
save `tmp',replace
qui des
local k=r(k)

forva i=1/`k' {
   use `tmp',clear
   sort mat`i'
   save,replace
   snapshot restore 1
   gen mat`i'=X
   sort mat`i'
   merge 1:1 mat`i' using `tmp',keep(match master)
   sum Y1 if _merge==3
   local mean=s(mean)
   replace Y2=`mean' in `i'
   drop _merge mat`i'
   snapshot erase _all
   snapshot save, label("main data")
}
当然肯定还有更简单、更直接的方法，这只是一个不算聪明但是比较直观的办法。

老树皮 发表于 2013-7-17 13:03
这个可以有很多种方法：Mata、matrix矩阵运算，也可以使用普通的Stata数据运算方式。假设你的矩阵叫作mat， ...

非常感谢高手指点，只是编程不是太能看懂，能大概说一下思路么？
有几个问题需要说明：
关联矩阵中第一列是和数据表中的x对应的，是不同国家的名称，而每一行后面对应的都是乱的，没有规律，我重新弄了个关联矩阵，如A对应的是CDFG,B对应的是AC,等等，是指第一列不同国家的邻国名称。

A	C	D	F	G
B	A	C
C	E	F	G
D	B	C	E
E	A
F	C	D	E
G	E

因此，我上面说的关联矩阵实际上应该不算是矩阵，因为ABCD等都是字符型的，而且不整齐，如第一行有4个，第二行有2个，第三行有3个……
local k=r(k)    请问K是什么？


请问更简单更直接的办法是什么呢？  
再次感谢！

老树皮 发表于 2013-7-17 13:03
这个可以有很多种方法：Mata、matrix矩阵运算，也可以使用普通的Stata数据运算方式。假设你的矩阵叫作mat， ...

请问更简单更直接的办法是什么呢？

Brdic 发表于 2013-7-17 21:22
请问更简单更直接的办法是什么呢？

使用mata,应该程序更短一些,运算也更快一些,但是对编程技术要求更高,我没时间推敲,所以只能提供一个大体方向...

Brdic 发表于 2013-7-17 21:21
非常感谢高手指点，只是编程不是太能看懂，能大概说一下思路么？
有几个问题需要说明：
关联矩阵中第一 ...

今天很忙,晚上有时间的话给你细说一下思路. 另外,k表示有多少个变量。

老树皮 发表于 2013-7-18 09:39
今天很忙,晚上有时间的话给你细说一下思路. 另外,k表示有多少个变量。

好的，非常感谢！我大概弄懂了思路了。还有几个疑问：
mat矩阵如何录入进去？因为都是字符的，而且不整齐。
snapshot有什么作用？
谢谢！

Brdic 发表于 2013-7-18 09:59
好的，非常感谢！我大概弄懂了思路了。还有几个疑问：
mat矩阵如何录入进去？因为都是字符的，而且不整齐 ...

最近很忙，所以现在才答复。不知道你是否已经解决该问题。Anyway，答复如下：

因为你一开始说是关联矩阵，所以我以为你已经有矩阵了。但实际上你看我的那个程序，还是把矩阵转化成数据序列的，只不过转置了一下。所以如果你开始就将该矩阵保存成数据序列，就不是问题了。

snapshot 的作用是将数据暂时存在内存中。当然你也可以每一步都保存在硬盘中，然后调用，但是比较慢。

老树皮 发表于 2013-7-20 14:09
最近很忙，所以现在才答复。不知道你是否已经解决该问题。Anyway，答复如下：

因为你一开始说是关联矩 ...

感谢回复！
我重新研究了一下程序，
merge 1:1 mat`i' using `tmp',keep(match master) 意思是要把数据文件和矩阵文件匹配起来，因此应该有个临时文件来放数据文件，前面没有写这句，所以gen mat`i'=X 这句也执行不了。
我重新改了下，但是最后
replace y2=`mean' in `i' 这一句命令没有把均值放在第i行，我单独运行时提示是0个值发生改变，不知怎么回事
另外，根据
save `tmp',replace
qui des /* quietly  describe*/
local k=r(k)
k是根据tmp这个文件来定义，是不是一直要打开tmp这个文件k才能循环，因为中间我又打开了数据文件，这样才能提取x，并合并，最后执行这个循环好像循环不了了。
还有，实际上我是面板数据，上面的程序只是完成了一个年份的，如果要把所有年份的都弄一遍，是不是还要弄个大循环？