# -*- coding: utf-8 -*-
"""
Created on Sun May 27 10:36:11 2018

@author: Kinga
"""

import numpy as np
from statistics import mode

x = []

x.extend(np.full((10), 0))
x.extend(np.full((18), 1))
x.extend(np.full((26), 2))
x.extend(np.full((29), 3))
x.extend(np.full((11), 4))
x.extend(np.full((3), 5))
x.extend(np.full((0), 6))
x.extend(np.full((1), 7))
x.extend(np.full((1), 8))


x=np.sort(x)
xsr = np.mean(x) #srednia arytmetyczna
print('Wartosc sredniej arytmetycznej wynosi:',xsr)

odchx = np.std(x) # odchylenie standardowe
print('Wartosc odchylenia stanardowego wynosi:',odchx)

modax = mode(x) #modalna
print('Wartosc modalnej wynosi:',modax)

medx = np.median(x) #mediana
print('Wartosc mediany wynosi:',medx)

#kwartyle

n1 = (len(x)+3)/4
n1c = int(np.floor(n1)) #zaokraglenie
n1u=n1-n1c
#print('Czesc calkowita wynosi:',"\n",n1c,"\n") #czesc calkowita
#print('Czesc ulamkowa wynosi:',"\n",n1u,"\n") #czesc ulamkowa

Q1 = x[n1c-1]+(x[n1c]-x[n1c-1])*n1u
print('Pierwszy kwartyl:',Q1)

print('Drugi kwartyl:',medx)#drugi kwartyl to mediana

n1 = (3*len(x)+1)/4
n1c = int(np.floor(n1)) #zaokraglenie w dół
n1u=n1-n1c

Q3 = x[n1c-1]+(x[n1c]-x[n1c-1])*n1u
print('Trzeci kwartyl:',Q3)

max = np.max(x)
min = np.min(x)
R = max - min #rozstep
print('Rozstęp wynosi:',R)

bwz = np.abs(x-xsr)
opos = np.mean(bwz)
print('Odchylenie przeciętne od sredniej:',opos)

u = np.abs(x-medx)
opom = np.mean(u)
print('Odchylenie przeciętne od mediany:',opom)

RQ = Q3 - Q1
print('Rozstep cwiartkowy:',RQ)

DQR = RQ/2
print('Odchylenie cwiartkowe:',DQR)

wzws = odchx/xsr*100 #wspolczynnik zmiennosci wzgledem sredniej
print('wspolczynnik zmiennosci wzgledem sredniej:',wzws)

wzwm = DQR/medx*100 #wspolczynnik zmiennosci wzgledem mediany
print('wspolczynnik zmiennosci wzgledem mediany:',wzwm)


#wspolczynnik asymetrii
wa=np.sum((x-xsr)**3)/len(x)/odchx**3
print('Wspólczynnik asymetrii',wa)

#wspolczynnik skosnosci Pearsona

Ws=(xsr-medx)/odchx
print('wspolczynnik skosnosci Pearsona:', Ws)

#pozycyjny wspolczynnik skosnosci
Wp=(Q1+Q3-2*medx)/(Q3-Q1)
print('pozycyjny wspolczynnik skosnosci:', Wp)

if wa ==0:
    print('Rozkład symetryczny')
elif wa>0:
    print('Asymetria prawostronna')
else:
    print('asymetria lewostronna')

#tworzenie funkcji
def mz(dane, rzad):
    momz=np.mean(dane**rzad)
    return momz

print(xsr,"=",mz(x,1))