Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- for i=1:111
- knots=defaultknots(x,i);
- d=4 ;
- t=[zeros(1,d) knots repelem(max(x),d)];
- range=[min(x) max(x)];
- z=bspline_basismatrix(d,t,x);
- Z=[ones(1665,1) z];
- lme = fitlmematrix(X,y,Z,G,'CovariancePattern','Isotropic');
- AI=[];
- AI=lme.ModelCriterion.AIC
- AII=min(AI)
- i
- i=i+1;
- end
Add Comment
Please, Sign In to add comment