Calculating standard deviation of multivarible regression analysis

So I have been working on a multivariable regression analysis with 13 observations and 4 variables. I am calculating coefficients for each hour of the day by the way. I calculated the coefficients with b = regress(y,X) function. In addition to the coefficients, I would like to get standard deviations of the residuals for each variable. I could get for example confidence intervals, but the regress() only gives me one value for each loop/hour. There are functions like fitlm which gives me more information, but they are more difficult to use with loops.

So do you have any idea what is the best way to calculate 24 hourly coefficints for 4 variables, as well as the standard deviation of errors/standard deviation of coefficints for each variable?