mean/median

# Среднее арифметическое (выборочное среднее)
def SampleMean(x):
    SampleMean = 0

    for i in range(len(x)):
        SampleMean += x[i]

    return SampleMean / len(x)

# Выборочная медиана
def SampleMedian(x):
    x_var = np.sort(x)

    if len(x_var) % 2 == 0:
        return (x_var[int(len(x_var) / 2)] + x_var[int(len(x_var) / 2 + 1)]) / 2
    else:
        return x_var[int(len(x_var) / 2 + 1)]

tetaAcc = 6 # teta
K = 50

# Оценки для распределений (выборочное среднее и медиана)
SampleMeanN01 = np.array([])
SampleMedianN01 = np.array([])
SampleMeanE = np.array([])
SampleMedianE = np.array([])
SampleMeanK3 = np.array([])
SampleMedianK3 = np.array([])

for i in range(K):
    K3 = KA(3, 50)
    yK = np.array([])

    for i in range(len(K3)):
        yK = np.append(yK, tetaAcc + K3[i])

    SampleMeanK3 = np.append(SampleMeanK3, SampleMean(yK))
    SampleMedianK3 = np.append(SampleMedianK3, SampleMedian(yK))

for i in range(K):
    N = NStandart(K)
    yN = np.array([])

    for i in range(len(N)):
        yN = np.append(yN, tetaAcc + N[i])

    SampleMeanN01 = np.append(SampleMeanN01, SampleMean(yN))
    SampleMedianN01 = np.append(SampleMedianN01, SampleMedian(yN))

for i in range(K):
    E = Exp(-1, 1, K)
    yR = np.array([])

    for i in range(len(E)):
        yR = np.append(yR, tetaAcc + E[i])

    SampleMeanE = np.append(SampleMeanE, SampleMean(yR))
    SampleMedianE = np.append(SampleMedianE, SampleMedian(yR))

# Среднеквадратическая погрешность
def RootMeanSquare(x):
    Mean = 0

    for i in range(len(x)):
        Mean += x[i]

    Mean /= len(x)
    s = 0

    for i in range(len(x)):
        s += (x[i] - Mean) ** 2

    s /= len(x) * (len(x) - 1)
    return np.sqrt(s)

# Среднеквадратические погрешности для оценок трех распределений
ErrorSampleMeanN01 = RootMeanSquare(SampleMeanN01)
ErrorSampleMedianN01 = RootMeanSquare(SampleMedianN01)
ErrorSampleMeanE = RootMeanSquare(SampleMeanE)
ErrorSampleMedianE = RootMeanSquare(SampleMedianE)
ErrorSampleMeanK3 = RootMeanSquare(SampleMeanK3)
ErrorSampleMedianK3 = RootMeanSquare(SampleMedianK3)

# Вывод результата
print("N01")
print("SampleMean")
print(SampleMeanN01)
print("SampleMedian")
print(SampleMedianN01)
print("E(-1, 1)")
print("SampleMean")
print(SampleMeanE)
print("SampleMedian")
print(SampleMedianE)
print("K(3)")
print("SampleMean")
print(SampleMeanK3)
print("SampleMedian")
print(SampleMedianK3)

print("N(0,1): SampleMean = %.5f , SampleMedian = %.5f" %(ErrorSampleMeanN01, ErrorSampleMedianN01))
print("E(-1, 1): SampleMean = %.5f , SampleMedian = %.5f" %(ErrorSampleMeanE, ErrorSampleMedianE))
print("K(3): SampleMean = %.5f , SampleMedian = %.5f" %(ErrorSampleMeanK3, ErrorSampleMedianK3))