data_new['e']=data_new['y']-model75.fittedvalues
sigma1=np.round((np.sum((data_new[data_new['region']==1]['e'])**2))/(len(data_new[data_new['region']==1]['e'])-1),2)
sigma2=np.round((np.sum((data_new[data_new['region']==2]['e'])**2))/(len(data_new[data_new['region']==2]['e'])-1),2)
sigma3=np.round((np.sum((data_new[data_new['region']==3]['e'])**2))/(len(data_new[data_new['region']==3]['e'])-1),2)
sigma4=np.round((np.sum((data_new[data_new['region']==4]
['e'])c1=np.round((sigma1/(((49)**-1)*(np.sum(data_new['e']**2))))**0.5,2)
c2=np.round((sigma2/(((49)**-1)*(np.sum(data_new['e']**2))))**0.5,2)
c3=np.round((sigma3/(((49)**-1)*(np.sum(data_new['e']**2))))**0.5,2)
c4=np.round((sigma4/(((49)**-1)*(