全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
10574 12
2009-07-20
悬赏 40 个论坛币 未解决
我用了一个小样本,数据你们可以下载试试
proc mixed可以用来估计混合模型,我试过可以用这个来估计固定效应,即估计结果与用proc tscsreg加fixone的一样。

程序一:
proc mixed data=test;
class stkcd ;
model detalever=stkcd tde_above tde_below/s;
run;

程序二:
proc tscsreg data=test;
model detalever=tde_above tde_below/fixone;
id stkcd Trdynt;
run;

即程序一和程序二估计出的结果是一样的。以下是问题:


问题一:在于如何用proc mixed 来估计随机效应模型呢?即如何设计程序,估计结果与下面的程序估计结果相同?
proc tscsreg data=test;
model detalever=tde_above tde_below/ranone;
id stkcd Trdynt;
run;

问题二:如果同时考虑固定效应和随机效应?程序怎么写?

test.rar

大小:38.58 KB

 马上下载

本附件包括:

  • test.sas7bdat

二维码

扫码加我 拉你入群

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

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

全部回复
2009-7-20 23:46:22
proc mixed有random这个statement,你查下手册
二维码

扫码加我 拉你入群

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

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

2009-7-20 23:52:09
我查了,可能是我基础比较差的原因,就是没试出来,另外,也怕自己的程序有问题
二维码

扫码加我 拉你入群

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

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

2009-7-21 04:56:51
从你的模型、数据和写的程序来看,本质上是分层随机抽样数据,同时包含随机固定和随机效应的混合效应模型,具体的程序如下,供参考(这里STKCD  TYDYNT我都当成随机效应;为修改任何变量属性):
proc mixed data=test METHOD=REML CL ALPHA=.05;           
class stkcd;   
model detalever = tde_above tde_below / HTYPE=3 DDFM=SATTERTH SOLUTION CL
ALPHA=.05;
random STKCD  TYDYNT/ SOLUTION CL ALPHA=.05;
run;
quit;

另外,你写的程序和提供的数据,我没有发现你对tydnt这个变量做了时间序列的定义(把tydnt转换成时间变量),当然,没有定义的情况下,前两个程序的参数是一样的,但proc tscsreg是要定义时间变量的。
大写的部分是补充的。钱一定要加,多少你看着办 [lol]
二维码

扫码加我 拉你入群

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

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

2009-7-21 10:43:24
4# 坐看云起时

TO:坐看云起时
呵,我只是附加最后执行的程序,另外,我这里没有考虑时间效应。

程序三:
proc tscsreg data=test;
model detalever=tde_above tde_below/ranone;
id stkcd Trdynt;
run;
结果如下:
                                                      Parameter Estimates
                                                                 Standard
                               Variable        DF    Estimate       Error    t Value    Pr > |t|
                               Intercept        1    -0.00609      0.0289      -0.21      0.8339
                               tde_above        1    1.223141      0.4427       2.76      0.0074
                               tde_below        1    0.375206      0.1822       2.06      0.0434
用你的程序执行后得到的结果为:
               The Mixed Procedure
                                                  Solution for Fixed Effects
                                       Standard
              Effect       Estimate       Error      DF    t Value    Pr > |t|     Alpha       Lower       Upper
              Intercept    0.002327     0.02614    17.8       0.09      0.9300      0.05    -0.05263     0.05729
              tde_above      1.2401      0.4410    66.3       2.81      0.0065      0.05      0.3597      2.1204
              tde_below      0.2910      0.1777    52.1       1.64      0.1076      0.05    -0.06565      0.6476
还是不一样啊
二维码

扫码加我 拉你入群

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

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

2009-7-21 12:39:49
要用proc tscsreg必须要把tydynt转换成时间变量属性(这是时间序列横截面回归,你现在的tydynt属性只是一般数量变量-best12.,不是日期属性-daye year monyy........;输入2001 2002.....机器不会自动识别这是时间变量,要从新定义)
,我给你的程序是假定你对tydynt转换成时间变量属性后再用(混合模型也要定义的)

你的结果是未定义tydynt时间属性的计算结果
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群