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")
}