Studentized residuals python

def internally_studentized_residual(X,Y):
    X = np.array(X, dtype=float)
    Y = np.array(Y, dtype=float)
    mean_X = np.mean(X)
    mean_Y = np.mean(Y)
    n = len(X)
    diff_mean_sqr = np.dot((X - mean_X), (X - mean_X))
    beta1 = np.dot((X - mean_X), (Y - mean_Y)) / diff_mean_sqr
    beta0 = mean_Y - beta1 * mean_X
    y_hat = beta0 + beta1 * X
    residuals = Y - y_hat
    h_ii = (X - mean_X) ** 2 / diff_mean_sqr + (1 / n)
    Var_e = math.sqrt(sum((Y - y_hat) ** 2)/(n-2))
    SE_regression = Var_e*((1-h_ii) ** 0.5)
    studentized_residuals = residuals/SE_regression
    return studentized_residuals

def deleted_studentized_residual(X,Y):
    #formula from https://newonlinecourses.science.psu.edu/stat501/node/401/
    r = internally_studentized_residual(X,Y)
    n = len(r)
    return [r_i*math.sqrt((n-2-1)/(n-2-r_i**2)) for r_i in r]

Source: stackoverflow.com

Similar codes

1 vote124 views1 answer
How to make a crosshair in python
Python
0 vote123 views1 answer
ValueError: query data dimension must match training data dimension
Python
0 vote113 views1 answer
Slicing time series
Python
4 votes137 views1 answer
Get all count rows pandas
Python
4 votes127 views1 answer
Print type of exception python
Python
0 vote126 views1 answer
Nptel swayam
Python
0 vote146 views1 answer
Django amzon like app
Python Django
0 vote107 views1 answer
How to get the access of python on cmd
Python
0 vote116 views1 answer
Creating a record in python
Python
0 vote132 views1 answer
Opencv loop video
Python