library(titanic)
## Warning: package 'titanic' was built under R version 4.5.3
data(titanic_train)
total_selamat <- sum(titanic_train$Survived)
print(total_selamat)
## [1] 342
table(titanic_train$Survived)
##
## 0 1
## 549 342
data(faithful)
cor(faithful$eruptions, faithful$waiting)
## [1] 0.9008112
plot (mtcars$wt, mtcars$mpg)
cor(mtcars$mpg, mtcars$cyl)
## [1] -0.852162
library(ggplot2)
ggplot(iris, aes(x = Species, y =Sepal.Length, fill = Species)) +
geom_boxplot() +
labs(title = "Distribusi Sepal.Length antar Specier iris" ,
y = "Sepal Length (cm)") +
theme_minimal()
data(iris)
boxplot(Sepal.Length ~ Species, data = iris,
xlab = "Species",
ylab = "Sepal Length (cm)",
main = "Distribusi Sepal.Length antar Species",
col = c("red", "yellow", "blue"))
data(ChickWeight)
chicks_sample <- subset(ChickWeight, Chick %in% c(1,2,3,4,5))
plot(weight ~ Time, data = chicks_sample,
type = "n",
xlab = "Waktu (hari)",
ylab = "Berat (gram)",
main = "Tren Berat Anak Ayam dari Waktu")
for(i in unique(chicks_sample$Chick)) {
subset_data <- subset(chicks_sample, Chick == i)
lines(subset_data$Time, subset_data$weight, col = i, lwd = 1.5)
}
legend("topleft", legend = paste("Ayam", 1:5), col = 1:5, lwd = 1.5, cex =
0.8)
library(ggplot2)
ggplot(ChickWeight, aes(x = Time, y = weight, group = Chick, color = Diet)) +
geom_line(alpha = 0.5) +
geom_smooth(aes(group = 1), method = "loess", se = FALSE, color = "black",
size = 1.5) +
labs(title = "Tren Berat Anak Ayam dari Waktu",
x = "Waktu (hari)",
y = "Berat (gram)") +
theme_minimal()
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## `geom_smooth()` using formula = 'y ~ x'
data(airquality)
sum(is.na(airquality$Ozone))
## [1] 37
ozone_median <- median(airquality$Ozone, na.rm = TRUE)
airquality$Ozone_imputed <- ifelse(is.na(airquality$Ozone), ozone_median,
airquality$Ozone)
summary(airquality$Ozone_imputed)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.00 21.00 31.50 39.56 46.00 168.00
#soal 27
str(mtcars)
## 'data.frame': 32 obs. of 11 variables:
## $ mpg : num 21 21 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 ...
## $ cyl : num 6 6 4 6 8 6 8 4 4 6 ...
## $ disp: num 160 160 108 258 360 ...
## $ hp : num 110 110 93 110 175 105 245 62 95 123 ...
## $ drat: num 3.9 3.9 3.85 3.08 3.15 2.76 3.21 3.69 3.92 3.92 ...
## $ wt : num 2.62 2.88 2.32 3.21 3.44 ...
## $ qsec: num 16.5 17 18.6 19.4 17 ...
## $ vs : num 0 0 1 1 0 1 0 1 1 1 ...
## $ am : num 1 1 1 0 0 0 0 0 0 0 ...
## $ gear: num 4 4 4 3 3 3 3 4 4 4 ...
## $ carb: num 4 4 1 1 2 1 4 2 2 4 ...
summary(mtcars)
## mpg cyl disp hp
## Min. :10.40 Min. :4.000 Min. : 71.1 Min. : 52.0
## 1st Qu.:15.43 1st Qu.:4.000 1st Qu.:120.8 1st Qu.: 96.5
## Median :19.20 Median :6.000 Median :196.3 Median :123.0
## Mean :20.09 Mean :6.188 Mean :230.7 Mean :146.7
## 3rd Qu.:22.80 3rd Qu.:8.000 3rd Qu.:326.0 3rd Qu.:180.0
## Max. :33.90 Max. :8.000 Max. :472.0 Max. :335.0
## drat wt qsec vs
## Min. :2.760 Min. :1.513 Min. :14.50 Min. :0.0000
## 1st Qu.:3.080 1st Qu.:2.581 1st Qu.:16.89 1st Qu.:0.0000
## Median :3.695 Median :3.325 Median :17.71 Median :0.0000
## Mean :3.597 Mean :3.217 Mean :17.85 Mean :0.4375
## 3rd Qu.:3.920 3rd Qu.:3.610 3rd Qu.:18.90 3rd Qu.:1.0000
## Max. :4.930 Max. :5.424 Max. :22.90 Max. :1.0000
## am gear carb
## Min. :0.0000 Min. :3.000 Min. :1.000
## 1st Qu.:0.0000 1st Qu.:3.000 1st Qu.:2.000
## Median :0.0000 Median :4.000 Median :2.000
## Mean :0.4062 Mean :3.688 Mean :2.812
## 3rd Qu.:1.0000 3rd Qu.:4.000 3rd Qu.:4.000
## Max. :1.0000 Max. :5.000 Max. :8.000
head(mtcars)
## mpg cyl disp hp drat wt qsec vs am gear carb
## Mazda RX4 21.0 6 160 110 3.90 2.620 16.46 0 1 4 4
## Mazda RX4 Wag 21.0 6 160 110 3.90 2.875 17.02 0 1 4 4
## Datsun 710 22.8 4 108 93 3.85 2.320 18.61 1 1 4 1
## Hornet 4 Drive 21.4 6 258 110 3.08 3.215 19.44 1 0 3 1
## Hornet Sportabout 18.7 8 360 175 3.15 3.440 17.02 0 0 3 2
## Valiant 18.1 6 225 105 2.76 3.460 20.22 1 0 3 1
pairs(mtcars)
plot(mtcars$wt, mtcars$mpg)
#soal 30
library(ggplot2)
ggplot(diamonds, aes(x = cut, y = price)) +
geom_boxplot() +
labs(title = "Distribusi price berdasarkan cut",
x = "Cut",
y = "Price")
library(ggplot2)
data(diamonds)
class(diamonds$cut)
## [1] "ordered" "factor"
str(diamonds$cut)
## Ord.factor w/ 5 levels "Fair"<"Good"<..: 5 4 2 4 2 3 3 3 1 3 ...
unique(diamonds$cut)
## [1] Ideal Premium Good Very Good Fair
## Levels: Fair < Good < Very Good < Premium < Ideal