%mc
fs = 1;
m = 1.3;
c = 4;
x = 1:30;
yunits = 'Pa';

ca = ao(plist('type', 'cdata', 'yvals', c.*ones(size(x)), 'yunits', yunits));
xa = ao(plist('type', 'cdata', 'yvals', x, 'yunits', yunits));
n = 1*randn(size(x));
sn = 2;
nats = sn.* ao(plist('tsfcn','randn(size(t))', 'fs', fs, 'nsecs', length(n), 'yunits', yunits));
ya = m.*xa + ca + nats;

% Fit with straightLineFit, pest output
yf_p = straightLineFit(ya, plist('output', 'pest'));
yf = yf_p.eval(xa, plist('type', 'tsdata'));

% Compare fitted data with original
iplot(ya,yf)

% Fit with straightLineFit, ao output
yf_a = straightLineFit(ya, plist('output', 'ao'));
