question1
proc nlp data=t1;
max logf;
parms const pri sec man her10 size tang gea exch ind1 ind2 ind3 ind4 ind5 ind6 ind7 ind8 ind9 ind10 ind11 ind12 ind13 ind14 ind15 ind16 ind17 ind18 ind19 ind20 ind21 ind22 y1 y2 y3 y4 y5;
u=roa - (const + pri * prilarge + sec * second + man*manager+her10*h10+size*size_roa+tang*tangible+gea*gear+exch*exchange+ind1*B+ind2*C+ind3*D+ind4*E+ind5*F+ind6*G+ind7*H+ind8*I+ind9*J+ind10*K+ind11*L+ind12*M+ind13*C0+ind14*C1+ind15*C2+ind16*C3+ind17*C4+ind18*C5+ind19*C6+ind20*C7+ind21*C8+ind22*C99+y1*year99+y2*year00+y3*year01+y4*year02+y5*year03);
sigmau2=0.0604154**2;
sigmae2=0.0675433**2;
no=_n_;
f=-0.5*(((roa-const+pri*prilarge+sec*second+man*manager+her10*h10+size*size_roa+tang*tangible+gea*gear+exch*exchange+ind1*B+ind2*C+ind3*D+ind4*E+ind5*F+ind6*G+ind7*H+ind8*I+ind9*J+ind10*K+ind11*L+ind12*M+ind13*C0+ind14*C1+ind15*C2+ind16*C3+ind17*C4+ind18*C5+ind19*C6+ind20*C7+ind21*C8+ind22*C99+y1*year99+y2*year00+y3*year01+y4*year02+y5*year03)^2-sigmau2*u**2/(no*sigmau2+sigmae2))/sigmae2+log(no*sigmau2/sigmae2+1)+no*log(2π*sigmae2));
logf=log(f);
run;
询问,为什么我跑这个程序会出现提示:
proc nlp data=t1;
ERROR: 过程 NLP 没有找到。
我安装的SAS是9.1版的,是否与SAS的版本有关,
只在9.2中才能跑那个语句么?
还有,请帮忙看一下语句编写是否有误?
主要是想用极大似然面板估计法MLP法做,即通过
进而计算回归方程
,得出多个贝塔的估计值(X中有多个变量)。
question2
另,不知道MLP用上面的程序是否正确,我又编了以下的程序,
不知道哪一个才是正确的,烦请指正(如下述)
proc mixed data=group3.sampleroa method=ml covtest;
class stkcd;
model roa=exchange gear tangible size_roa manager h10 second prilarge year99 year00
year01 year02 year03 b c c0 c1 c2 c3 c4 c5 c6 c7 c8 c99 d e f g h i j k l m/solution;
random intercept year99 year00 year01 year02 year03 / sub=stkcd;
run;
quit;
question3
我用了上面的程序,但是结果
由于没有自由度DF,造成P值缺失,这又是怎么回事呢?
有哪些因素会造成自由度变为0啊???明明原始数据里有3000多个观测的,而我的因变量不超过20个。这个问题在网络上搜了好几次,一直搜不到解决的方法,希望论坛里的高手能帮忙解决一下,感激不尽。。
P.S如果哪个细节没有说清楚,请留言,我再详说。