Memanggil Data Iris

z = read.csv("D:/StatistikaSemester3/iris",header =TRUE )
head(z)
##   Id SepalLengthCm SepalWidthCm PetalLengthCm PetalWidthCm     Species
## 1  1           5.1          3.5           1.4          0.2 Iris-setosa
## 2  2           4.9          3.0           1.4          0.2 Iris-setosa
## 3  3           4.7          3.2           1.3          0.2 Iris-setosa
## 4  4           4.6          3.1           1.5          0.2 Iris-setosa
## 5  5           5.0          3.6           1.4          0.2 Iris-setosa
## 6  6           5.4          3.9           1.7          0.4 Iris-setosa

Menentukan Summary

summary(z)
##        Id         SepalLengthCm    SepalWidthCm   PetalLengthCm  
##  Min.   :  1.00   Min.   :4.300   Min.   :2.000   Min.   :1.000  
##  1st Qu.: 38.25   1st Qu.:5.100   1st Qu.:2.800   1st Qu.:1.600  
##  Median : 75.50   Median :5.800   Median :3.000   Median :4.350  
##  Mean   : 75.50   Mean   :5.843   Mean   :3.054   Mean   :3.759  
##  3rd Qu.:112.75   3rd Qu.:6.400   3rd Qu.:3.300   3rd Qu.:5.100  
##  Max.   :150.00   Max.   :7.900   Max.   :4.400   Max.   :6.900  
##   PetalWidthCm     Species         
##  Min.   :0.100   Length:150        
##  1st Qu.:0.300   Class :character  
##  Median :1.300   Mode  :character  
##  Mean   :1.199                     
##  3rd Qu.:1.800                     
##  Max.   :2.500

Menentukan Nilai Modus

mod <- function(a){uniqx <- unique (a)
                   uniqx[which.max(tabulate(match(a, uniqx)))]}
mod(z$SepalLengthCm)
## [1] 5
mod(z$SepalWidthCm)
## [1] 3
mod(z$PetalLengthCm)
## [1] 1.5
mod(z$PetalWidthCm)
## [1] 0.2

Menghitung Nilai Rata-Rata Geometri

meangeo <- function(a){
  exp(mean(log(a)))
}
meangeo(z$SepalLengthCm)
## [1] 5.78572
meangeo(z$SepalWidthCm)
## [1] 3.023582
meangeo(z$PetalLengthCm)
## [1] 3.239757
meangeo(z$PetalWidthCm)
## [1] 0.837827

Menghitung Nilai Rata-Rata Harmoni

meanharm <- function(a){1/(mean(1/a))}
meanharm(z$SepalLengthCm)
## [1] 5.728905
meanharm(z$SepalWidthCm)
## [1] 2.993137
meanharm(z$PetalLengthCm)
## [1] 2.696472
meanharm(z$PetalWidthCm)
## [1] 0.4866465

Menghitung Nilai Skewness dan Kurtosis

library(e1071)
skewness(z$SepalLengthCm)
## [1] 0.3086407
skewness(z$SepalWidthCm)
## [1] 0.3274013
skewness(z$PetalLengthCm)
## [1] -0.2689994
skewness(z$PetalWidthCm)
## [1] -0.102906
kurtosis(z$SepalLengthCm)
## [1] -0.6058125
kurtosis(z$SepalWidthCm)
## [1] 0.1983681
kurtosis(z$PetalLengthCm)
## [1] -1.416683
kurtosis(z$PetalWidthCm)
## [1] -1.357368

Menampilkan Suatu Grafik

Menampilkan Boxplot

boxplot(z$SepalLengthCm, main = "SepalLengthCm", xlab = "SepalLength(Cm)", col = "coral")

boxplot(z$SepalWidthCm, main = "SepalWidthCm", xlab = "SepalWidth(Cm)", col = "Cadet Blue")

boxplot(z$PetalLengthCm, main = "PetalLengthCm", xlab = "PetalLength(Cm)",col = "Dark Sea Green")

boxplot(z$PetalWidthCm, main = "PetalWidthCm", xlab= "PetalWidth(Cm)", col = "Khaki")

Menampilkan Histogram

hist(z$SepalLengthCm, main = "SepalLengthCm", xlab= "SepalLength(Cm)", col = "Light Steel Blue")

hist(z$SepalWidthCm, main = "SepalWidthCm", xlab = "SepalWidth(Cm)", col = "Medium Slate Blue")

hist(z$PetalLengthCm, main = "PetalLengthCm", xlab = "PetalLength(Cm)", col = "Moccasin")

hist(z$PetalWidthCm, main = "PetalWidthCm", xlab = "PetalWidth(Cm)", col = "Pale Green")

## Menampilkan Scatter Plot

par(mfrow=c(2,2))
plot(z$SepalLengthCm, main = "SepalLengthCm", xlab = "SepalLength(Cm)", col = "Blue" )

plot(z$SepalWidthCm, main = "SepalWidthCm", xlab = "WidthLength(Cm)", col = "Red")

plot(z$PetalLengthCm, main = "PetalLengthCm", xlab = "PetalLength(Cm)", col = "Yellow")

plot(z$PetalWidthCm, main = "PetalWidthCm", xlab = "PetalWidth(Cm)", col = "Orange")

## Menampikan Barplot

barplot(z$SepalLengthCm, main = "SepalLengthCm", xlab = "SepalLength(Cm)", col = "Pale Green")

barplot(z$SepalWidthCm, main = "SepalWidthCm", xlab = "SepalWidth(Cm)", col = "Pale Turquoise")

barplot(z$PetalLengthCm, main = "PetalLengthCm", xlab = "PetalLength(Cm", col="Peach Puff")

barplot(z$PetalWidthCm, main = "PetalWidthCm", xlab = "PetalWidth(Cm)", col="Violet")