Napisać funkcję statystyki(x), która dla zadanego wektora x będzie wyznaczać podstawowe statystyki opisowe z jego współrzędnych (przynajmniej średnią, odch. standardowe, skośność, kurtozę, mininum, maksimum).
Następnie wylosować po 100 liczb z rozkładu N(0,1) zapisując wyniki w zmiennych x i y oraz wykonać tę funkcję dla x, y, x+y, x-y, p(x)+p( y), p(x)-p( y), gdzie p(w) oznacza wektor zawierający uporządkowane współrzędne wektora w.

statystyki=function(x)
{n=length(x)
struktura = list(srednia=mean(x),odchylenie_standardowe=sd(x), minimum=min(x) , maksimum=max(x),
skosnosc=(n/(n-1)/(n-2)*sum(x)^3)/sd(x)^3 )
print(struktura)
}

x= c(3,5,6,7,8,9,6,5,4,3,34,56,3,4)
statystyki(x)
z= rnorm(100)
y= rnorm(100)
x=z+y
x=z-y

sort(x)







Napisać funkcję srednia_p(x, p) obliczającą średnią potęgową rzędu p elementów wektora x.
Następnie dla 1500 liczb pseudolosowych, wylosowanych z rozkładu jednostajnego na przedziale [1,10], obliczyć średnią arytmetyczną, średnią harmoniczną, średnią geometryczną oraz średnią kwadratową (potęgową rzędu 2).

-- rozwiązanie 1
srednia_p=function(x,p)
{
z=0
n=length(x)
for(i in 1:n)
{z=z+x[i]^p}
(z/n)^(1/p)
}
-- rozwiązanie 2
srednia_p=function(x,p)
{
z=0
n=length(x)
for(i in x)
{z=z+i^p}
(z/n)^(1/p)
}
-- rozwiązanie 3
srednia_p=function(x,p) (mean(x^p))^(1/p)

x=sample(2:10,8)+rnorm(8)
runif(1500)