Help with Matlab code-please

7 次查看（过去 30 天）

显示更早的评论

Ozmando 2019-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445582-help-with-matlab-code-please

关闭： MATLAB Answer Bot 2021-8-20

在 MATLAB Online 中打开

Good afternoon,

I am trying to run the following code in MATLAB but get an error.

Error using *

Inner matrix dimensions must agree.

Error in vcqr (line 27)

e = y - x*bhat; % Generate residuals

Any help to fix the problem please?

Thanks

clear
load fish1
Z = [ones(111,1) stormy mixed];
X = [ones(111,1) d];
pihat = inv(Z'*Z)*(Z'*d);
tau = ['15';'25';'35';'45';'50';'55';'65';'75';'85'];
for i = 1:size(tau,1),
    disp(i);
    hamster = str2num(strcat('.',tau(i,:)));
    boo = rq([ones(size(y)) Z*pihat],y,hamster);
    [b,vc] = vcqr(boo,y,[ones(size(y)) Z*pihat],hamster);
    [b,f,c,seq,conv] = mcmc_flat_fish1_mon2('iqrobj_fish',boo,vc,50000,y,[ones(size(y)),d],Z,hamster);
    save(strcat('fish',tau(i,:),'_storm_v3_mon2'),'b','f','c','seq','conv');
    [b,f,c,seq,conv] = mcmc_flat_fish2('iqrobj_fish',boo,vc,50000,y,[ones(size(y)),d],Z,hamster);
    save(strcat('fish',tau(i,:),'_storm_v3_2'),'b','f','c','seq','conv');
end
% Here is the portion of the function related to the error:
vc = zeros(k,k);
b = zeros(k,2);
S = (1/n)*x'*x; 
e = y - x*bhat; % Line 27