Hi,
I am trying to use fitlm and rsquared ordinary functions and I can't figure out where the errors comes from.
Moreover, I want to plots 3 graphs but only the first one appears.
I would really appreciate your help.
here is my code :
function[]=reaction_order(A)
fit0=A;
fit1=log(A);
fit2=1./A;
time=[1:11]';
f0=fit(time,fit0','poly1');
figure
plot(f0,time,fit0);
xlabel('Time [s]');
ylabel('A [M]')
title('Zero order fit, A concentration as a function of time')
legend('A concentration', 'Fit curve')
mdl0=fitlm(time,f0);
rsquared0=mdl0.Rsquared.ordinary;
fprintf(rsquared0)
f2=fit(time,fit2','poly1');
figure;
plot(f2,time,fit2);
xlabel('Time [s]');
ylabel('1/A [1/M]')
title('Second order fit, 1/A as a function of time')
legend('1/A', 'Fit curve')
mdl2=fitlm(time,f2);
rsquared2=mdl2.Rsquared.ordinary;
The error is :
Error using classreg.regr.FitObject/assignData (line 134)
Predictor and response variables must have the same length.
Error in classreg.regr.TermsRegression/assignData (line 245)
model =
assignData@classreg.regr.ParametricRegression(model,X,y,w,asCat,varNames,excl);
Error in LinearModel.fit (line 1030)
model =
assignData(model,X,y,weights,asCatVar,dummyCoding,model.Formula.VariableNames,exclude);
Error in fitlm (line 121)
model = LinearModel.fit(X,varargin{:});
Error in chimieq2>reaction_order (line 50)
mdl0=fitlm(time,f0);
Error in chimieq2 (line 4)
reaction_order(A)
