%estimation of standard deviation
d_res=d_data-d_model;

n=length(d_res); %number of observations

sigma2=sum(d_res.^2)./(n-2); %variance

sigma=sqrt(sigma2); %standard deviation