BERIKUT Mean, Varians, Standar Deviasi, Skewness dan Kurtosis dengan menggunakan Program R :

· UKURAN PEMUSATAN (MEAN), UKURAN PENYEBARAN (VARIANS+STANDAR DEVIASI), UKURAN KESIMETRISAN (SKEWNESS) dan UKURAN KELANCIPAN (KURTOSIS) dengan metode biasa

stadeskmetodebiasa<-function(x)

{

n=length(x)

rata=sum(x)/n

kuadrat=(x-rata)^2

juml=sum(kuadrat)

varians=1/(n-1)*juml

stadev=sqrt(varians)

skewness=n/((n-1)*(n-2))*sum(((x-rata)/stadev)^3)

kurtosis=(n*(n+1))/((n-1)*(n-2)*(n-3))*sum(((x-rata)/stadev)^4)-(3*(n-1)^2)/((n-2)*(n-3))

cat("mean.biasa = ",rata,"\n")

cat("varians.biasa = ",varians,"\n")

cat("skewness.biasa = ",skewness,"\n")

cat("kurtosis.biasa = ",kurtosis,"\n")

}

· UKURAN PEMUSATAN (MEAN), UKURAN PENYEBARAN (VARIANS+STANDAR DEVIASI), UKURAN KESIMETRISAN (SKEWNESS) dan UKURAN KELANCIPAN (KURTOSIS) dengan metode MATRIKS

stadeskmetodematrix<-function(x)

{

n=length(x)

a=rep(1/n)

matrix1=rep(a,n)

matrix2=cbind(x)

rata=matrix1%*%matrix2

kuadrat=(x-rata)^2

juml=sum(kuadrat)

varians=1/(n-1)*juml

stdev=sqrt(varians)

b=n/((n-1)*(n-2))

matrix3=rep(b,n)

matrix4=cbind(((x-rata)/stdev)^3)

skewness=matrix3%*%matrix4

c=(n*(n+1))/((n-1)*(n-2)*(n-3))

matrix5=rep(c,n)

matrix6=cbind(((x-rata)/stdev)^4)

d=(3*(n-1)^2)/((n-2)*(n-3))

kurtosis=matrix5%*%matrix6-d

cat("mean = ",rata,"\n")

cat("varians = ",varians,"\n")

cat("skewness = ",skewness,"\n")

cat("kurtosis = ",kurtosis,"\n")

}