[问答] 最大似然代码运行出错 [推广有奖]

3关注
0粉丝

博士生

90%

还不是VIP/贵宾

威望: 0 级
论坛币: 1256 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 16891 点
帖子: 237
精华: 0
在线时间: 486 小时
注册时间: 2008-4-1
最后登录: 2022-10-4

楼主

dtll666

发表于 2011-8-31 18:01:05 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

use gpe2;
data1=xlsreadm("research.xls","b2:ax50",12,0);
s=data1;
data2=xlsreadm("research.xls","b2:ac50",13,0);
x=data2;
y=s~x;
call reset;
_nplot=1;
_method=5;
_iter=100;
_tol=1.0e-5;
_vcov=1;
_b={1,2};
call estimate(&cesll,y);
end;

proc cesll(data,b);
   local s,j,u_j,sum_c;
   s=data[.,1:49];
   sum_c=0;
   j=1;
      do while j<29;
      u_j=sumc((data[.,49+j]-b[2].*s*data[.,49+j])^2);
      sum_c=u_j+sum_c;
      j=j+1;
      endo;
retp(-49*28*0.5*ln(b[1]^2)+28*ln(det(eye(49)-b[2].*s))-0.5/(b[1]^2)*sum_c);
endp;
运行时，老是出现Argument must be positive 的错误提示。请高手指点。