function R = stdevlines(X, lines)
Y = lines(:,1) * X + lines(:,2);
R = std(Y);