楼主: leany
594 1

[问答] 求助 [推广有奖]

  • 0关注
  • 0粉丝

初中生

52%

还不是VIP/贵宾

-

威望
0
论坛币
10 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
59 点
帖子
5
精华
0
在线时间
24 小时
注册时间
2010-1-20
最后登录
2014-10-22

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
大神可以帮忙看下吗,这里nz只有等于2时才可以运行,换成其他matlab就显示.*处错误,Matrix dimensions must agree.实在看不出哪里错了这里rawdata是一个957x5的矩阵。

data = rawdata(1:900,:);
nz =1;                           
T = rows(data)-nz;
R = data(1+nz:T+nz,1:2);
neq=2;
re=data(1+nz:T+nz,3);
r=abs(re);
p=data(1+nz:T+nz,4);
n=data(1+nz:T+nz,5);
Y = ones(T,neq);
X = R;
Z = ones(T,1);
for i = 1:nz
  Z = [Z rawdata(1+nz-i:T+nz-i,1:2)];
end

b=[1;10;1;1];
gmmopt.infoz.momt='leanygmm';
gmmopt.plot=0;
gmmopt.gmmit=10;
out=gmm(b,gmmopt,Y,X,Z);

function [m,e]=leanygmm(b,infoz,stat,Y,X,Z)
global r;
global p;
global n;
[T, neq] = size(Y);
k = cols(Z);
north = neq*k;
R=X;
imrs =b(1)/b(2)*p.*(r.^(b(1)-1))+b(3)*b(4)/b(2)*n.*(r.^(b(4)-1));
e = repmat(imrs,1,neq).*R - Y;
m = reshape(Z'*e/T,north,1);


已解决,谢谢
二维码

扫码加我 拉你入群

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

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

关键词:Dimensions dimension function reshape Global Matrix matlab

已有 1 人评分论坛币 收起 理由
Xaero + 60 自己动手,丰衣足食。

总评分: 论坛币 + 60   查看全部评分

沙发
Xaero 发表于 2013-11-24 10:29:34 |只看作者 |坛友微信交流群
done
十年一觉扬州梦。
智不足以Academy,才尚不够Industry,情无力于Life。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-2 09:20