It looks like code is supposed to do exactly what I am after, unfortunately I am having problems with your York_fit.m code:
at 40 ==> tmp=Y/[X; ones(1,N)]; shouldn't it be tmp=[Y/X; ones(1,N)]; ?
but even if I change this I am still encountering more errors:
??? Error using ==> times
Matrix dimensions must agree.
Error in ==> york_fit at 58
Can I use multiple regression and
can I omit the offset in the model?
E.g. Y is explained by two betas corresponding with two different x (of same type of measurement, i.e. same error), but no offset is fitted?
Y = X*beta + error
where X = [x1; x2]