Dokumen R Markdown ini disusun sebagai teladan pembahasan analisis dengan tampilan tab berbeda seperti pada publikasi RPubs. Fokus utama adalah alur analisis (bukan hasil regresi), mulai dari pengumpulan data hingga evaluasi, sehingga mudah diikuti dan siap dipublish ke RPubs.
Pada tahap ini, data disiapkan untuk dianalisis. Data dapat berasal dari sumber sekunder seperti publikasi resmi atau data terbuka. Untuk keperluan teladan dan kemudahan replikasi, digunakan dataset contoh bawaan R. Tahap ini memastikan variabel respon dan prediktor tersedia serta konsisten.
data <- mtcars
head(data)
## 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
Eksplorasi data bertujuan memahami karakteristik awal data melalui statistik deskriptif dan visualisasi sederhana. Langkah ini membantu mengidentifikasi pola umum, rentang nilai, serta indikasi awal hubungan antar variabel.
summary(data)
## 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
plot(data$wt, data$mpg,
xlab = "Variabel Penjelas",
ylab = "Variabel Respon",
main = "Visualisasi Hubungan Awal")
Sebelum pemodelan, dilakukan pemeriksaan asumsi dasar regresi secara visual melalui grafik diagnostik. Tujuannya memastikan model yang dibangun berada pada kondisi yang layak untuk dianalisis lebih lanjut.
model_awal <- lm(mpg ~ wt, data = data)
par(mfrow = c(2,2))
plot(model_awal)
Tahap estimasi berfokus pada pembentukan model regresi linier sebagai representasi hubungan antara variabel. Pada teladan ini, estimasi ditampilkan seperlunya tanpa pembahasan angka secara mendalam.
model <- lm(mpg ~ wt, data = data)
summary(model)
##
## Call:
## lm(formula = mpg ~ wt, data = data)
##
## Residuals:
## Min 1Q Median 3Q Max
## -4.5432 -2.3647 -0.1252 1.4096 6.8727
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 37.2851 1.8776 19.858 < 2e-16 ***
## wt -5.3445 0.5591 -9.559 1.29e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.046 on 30 degrees of freedom
## Multiple R-squared: 0.7528, Adjusted R-squared: 0.7446
## F-statistic: 91.38 on 1 and 30 DF, p-value: 1.294e-10
Pengujian hipotesis digunakan untuk menilai signifikansi hubungan dalam model secara statistik. Tahap ini menegaskan kelayakan model sebelum dievaluasi secara menyeluruh.
anova(model)
## Analysis of Variance Table
##
## Response: mpg
## Df Sum Sq Mean Sq F value Pr(>F)
## wt 1 847.73 847.73 91.375 1.294e-10 ***
## Residuals 30 278.32 9.28
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Evaluasi dilakukan untuk menilai kinerja model secara umum melalui pola residual dan kecocokan model. Analisis ini membantu memastikan tidak terdapat pelanggaran asumsi yang mencolok.
par(mfrow = c(1,2))
hist(residuals(model), main = "Histogram Residual")
plot(fitted(model), residuals(model),
xlab = "Nilai Prediksi",
ylab = "Residual")
Secara keseluruhan, R Markdown ini menunjukkan bagaimana analisis disusun secara sistematis menggunakan tab terpisah untuk setiap tahapan. Pendekatan ini memudahkan pembaca mengikuti alur analisis dan sangat sesuai untuk dipublish ke RPubs sebagai laporan berbasis HTML.