def zad5():
delikates = [
23.4, 30.9, 18.8, 23.0, 21.4, 1, 24.6, 23.8, 24.1, 18.7, 16.3, 20.3,
14.9, 35.4, 21.6, 21.2, 21.0, 15.0, 15.6, 24.0, 34.6, 40.9, 30.7, 24.5,
16.6, 1, 21.7, 1, 23.6, 1, 25.7, 19.3, 46.9, 23.3, 21.8, 33.3, 24.9,
24.4, 1, 19.8, 17.2, 21.5, 25.5, 23.3, 18.6, 22.0, 29.8, 33.3, 1, 21.3,
18.6, 26.8, 19.4, 21.1, 21.2, 20.5, 19.8, 26.3, 39.3, 21.4, 22.6, 1,
35.3, 7.0, 19.3, 21.3, 10.1, 20.2, 1, 36.2, 16.7, 21.1, 39.1, 19.9,
32.1, 23.1, 21.8, 30.4, 19.62, 15.5
]
reneta = [
16.5, 1, 22.6, 25.3, 23.7, 1, 23.3, 23.9, 16.2, 23.0, 21.6, 10.8, 12.2,
23.6, 10.1, 24.4, 16.4, 11.7, 17.7, 34.3, 24.3, 18.7, 27.5, 25.8, 22.5,
14.2, 21.7, 1, 31.2, 13.8, 29.7, 23.1, 26.1, 25.1, 23.4, 21.7, 24.4,
13.2, 22.1, 26.7, 22.7, 1, 18.2, 28.7, 29.1, 27.4, 22.3, 13.2, 22.5,
25.0, 1, 6.6, 23.7, 23.5, 17.3, 24.6, 27.8, 29.7, 25.3, 19.9, 18.2,
26.2, 20.4, 23.3, 26.7, 26.0, 1, 25.1, 33.1, 35.0, 25.3, 23.6, 23.2,
20.2, 24.7, 22.6, 39.1, 26.5, 22.7
]
# cdfplot(delikates)
# cdfplot(reneta)
x1, y1 = prepare_x_y_for_cdfplot(delikates)
x2, y2 = prepare_x_y_for_cdfplot(reneta)
plt.step(x1, y1)
plt.step(x2, y2)
plt.show()
eval_KS_test(delikates, "kstest delikates")
eval_KS_test(reneta, "kstest reneta")
eval_KS2_test(delikates, reneta, "kstest2 delikates vs reneta")
qqplot(delikates)
qqplot(reneta)
def zad9():
kondensatory = read_csv_file("kondensatory.csv")
pojemnosc = [float(x["pojemnosc"]) for x in kondensatory]
# eval_KS_test(pojemnosc, "KS test")
# eval_lilliefors(pojemnosc, "lilliefors test")
eval_SW_test(pojemnosc, "SW test")
qqplot(pojemnosc)
def zad8():
zarowki = read_csv_file("zarowki.csv")
times = [float(x["czas"]) for x in zarowki]
# eval_KS_test(times, "KS test")
# eval_lilliefors(times, "lilliefors test")
eval_SW_test(times, "SW test")
qqplot(times)
def zad7():
pacjenci = read_csv_file("pacjenci.csv")
sugar = [float(x["cukier"]) for x in pacjenci]
# eval_KS_test(sugar, "KS test")
# eval_lilliefors(sugar, "lilliefors test")
eval_SW_test(sugar, "SW test")
qqplot(sugar)
def zad2():
control_b = [
0.08, 0.10, 0.15, 0.17, 0.24, 0.34, 0.38, 0.42, 0.49, 0.50, 0.70, 0.94,
0.95, 1.26, 1.37, 1.55, 1.75, 3.20, 6.98, 50.57
]
cdfplot(control_b, xscale="log")
qqplot(control_b)
def zad6():
pacjenci = read_csv_file("pacjenci.csv")
men_heights = [int(x["wzrost"]) for x in pacjenci if x["plec"] == "M"]
women_heights = [int(x["wzrost"]) for x in pacjenci if x["plec"] == "K"]
eval_lilliefors(men_heights, "kstest men")
eval_lilliefors(women_heights, "kstest women")
qqplot(men_heights)
qqplot(women_heights)
def zad4():
filename = "pacjenci.csv"
pacjenci = read_csv_file(filename)
men_heights = [int(x["wzrost"]) for x in pacjenci if x["plec"] == "M"]
women_heights = [int(x["wzrost"]) for x in pacjenci if x["plec"] == "K"]
qqplot(men_heights)
qqplot(women_heights)
eval_KS_test(men_heights, "kstest men")
eval_KS_test(women_heights, "kstest women")
eval_KS2_test(men_heights, women_heights, "kstest2 men vs women")
def zad10():
absolwenci = read_csv_file("absolwenci.csv", delimiter=';')
rolnictwo = [
int(x["SALARY"]) for x in absolwenci if x["COLLEGE"] == "Rolnictwo"
]
pedagogika = [
int(x["SALARY"]) for x in absolwenci if x["COLLEGE"] == "Pedagogika"
]
qqplot(rolnictwo, title="Rolnictwo")
qqplot(pedagogika, title="Pedagogika")
eval_KS_test(rolnictwo, "Rolnictwo KS test")
eval_lilliefors(rolnictwo, "Rolnictwo lilliefors test")
eval_SW_test(rolnictwo, "Rolnictwo SW test")
print("\n\n\n")
eval_KS_test(pedagogika, "Pedagogika KS test")
eval_lilliefors(pedagogika, "Pedagogika lilliefors test")
eval_SW_test(pedagogika, "Pedagogika SW test")