mata可不可以处理这样的式子？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

    做论文遇到瓶颈，看过了mata的手册，但是还是不会做具体的操作，所以希望论坛里的高手不要潜水，出来支几招，帮我渡过难关，先谢谢了！！
    公式如下：

      
    其中b、Φ、∑是已知数（第一个∑是求和符号，第二个∑是个已知数），求的是δ。b是1*4矩阵，b的上标n代表的是期数，不参与b的计算，只在等式右边用到。e=（0,0,1,0)'，Φ、∑是4*4矩阵，δ也是4*4矩阵，j是指数。
    因为未知数δ在等式右侧，公式又比较复杂，我不知道该怎么处理，是简单的现成的程序就可以还是要有复杂的编程才能算出来？这里求大家帮帮忙，谢谢了。

这里也把公式的推导过程贴出来供参考。



上式中n表示利率的期数，r代表利率。b表示长期利率对经济增长、通货膨胀、基准利率和基准利差等变量的敏感程度，x表示以上几个宏观变量，是已知数。
我参考的文献是《中国主要宏观经济变量与利率期限结构的关系:基于VAR-ATSM模型的分析》

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATA Mat 利率期限结构 宏观变量 宏观经济 求和符号 未知数 我不知道 程序 论文

这个，这个，我以前没有接触过mata，用stata处理矩阵问题也是第一次，所以实在是很难琢磨出解决的方法，论文又有时间要求，现在本人很是着急，请大家能帮忙就帮帮忙啊。谢谢啦。

虽然我不懂，但我觉得
1、你最好把公式的出处都贴出来，全面看一个问题可能会更容易理解，知道你到底需要做什么
2、stata不是万能的，有肯解决不了
3、看你文献中是如何解决的

蓝色 发表于 2013-6-2 22:48
虽然我不懂，但我觉得
1、你最好把公式的出处都贴出来，全面看一个问题可能会更容易理解，知道你到底需要做 ...

哦，谢谢你的回复，好歹是有人看了我的问题。
1，我不太理解第一条，这个公式是一个模型中几个公式推导出来的。你是指应该把公式的推导过程也贴出来么？
2，这个问题应该是解决的了的，因为文献里用同样的公式计算参数，我想应该是可以处理的。
3，额，英文文献里根本就没有提及过软件的问题，只给出结果。中文文献里只说用过gauss6.0软件，也没有具体的计算过程，直接给结果。我查过资料，guass的功能stata是可以实现的，而且我之前也用过stata，另一个软件我没接触过，所以想用stata来算。

你吧公式推导和应用的文献贴出来
或许根据文章中其它线索可以解决

蓝色 发表于 2013-6-3 22:42
你吧公式推导和应用的文献贴出来
或许根据文章中其它线索可以解决

哦哦，谢谢啊，我试试。

难是难在Affine Term Structure Models【ATSM】的部份，
毕竟ATSM都有专章在 Handbook of Financial Econometrics 罗！

目前Stata应当没有专门的程序啦！ 编写有点难度…
虽然有人在问，请参见
http://www.stata.com/statalist/archive/2013-02/msg00537.html

看起来Matlab那边是有相关的程序 【Matlab不愧是金融计量的好东西】
请参见
http://www.stanford.edu/~piazzesi/

个人感觉，如果各矩阵的阶数不高，直接在stata中操作矩阵即可。这样在循环中可以用到宏代换。

mata中的循环难以用上类似宏代换的功能，这使得变量序列的循环处理比较麻烦。

楼主未贴出数据。只能猜楼主想要这样的结果(其中几个方阵都是4阶的)：

mat e=(0,0,1,0)
forv n=1/16{
mat b`n'=J(4,4,0)
forv j=0/`=`n'-1'{
mat b`n'=b`n'+(phi`j')-(sigma`j')*(delta`j')
}
mat b`n'=e*b`n'/`n'
}

mat l b1
mat l b2
……
mat l b16

sungmoo 发表于 2013-6-5 11:29
个人感觉，如果各矩阵的阶数不高，直接在stata中操作矩阵即可。这样在循环中可以用到宏代换。

mata中的循 ...

直接操作？可以计算出这样的式子？
式子里n的最大值是16，但是这个式子涉及到了循环么？我编程现在基本是一片空白，麻烦你了。

本式不是求和吗？可用循环求和。

另外，n并不是矩阵的阶数。