library(stats)
library(car)
## Warning: package 'car' was built under R version 4.5.3
## Loading required package: carData
## Warning: package 'carData' was built under R version 4.5.3
library(lmtest)
## Warning: package 'lmtest' was built under R version 4.5.3
## Loading required package: zoo
## Warning: package 'zoo' was built under R version 4.5.3
##
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric
library(zoo)
data <- read.table(file.choose(), header = TRUE)
data
## No TinggiBadan BeratBadan
## 1 1 150 45
## 2 2 153 47
## 3 3 156 48
## 4 4 158 50
## 5 5 160 51
## 6 6 163 53
## 7 7 166 55
## 8 8 168 56
## 9 9 171 59
## 10 10 174 61
summary(data)
## No TinggiBadan BeratBadan
## Min. : 1.00 Min. :150.0 Min. :45.00
## 1st Qu.: 3.25 1st Qu.:156.5 1st Qu.:48.50
## Median : 5.50 Median :161.5 Median :52.00
## Mean : 5.50 Mean :161.9 Mean :52.50
## 3rd Qu.: 7.75 3rd Qu.:167.5 3rd Qu.:55.75
## Max. :10.00 Max. :174.0 Max. :61.00
plot (data$TinggiBadan, data$BeratBadan, xlab = "Berat Badan", ylab = "Tinggi Badan", main = "Scatter Plot Berat Badan dan Tinggi Badan", pch = 19)
abline(lm(BeratBadan ~ TinggiBadan, data = data), col = "blue", lwd = 2)
shapiro.test(data$TinggiBadan)
##
## Shapiro-Wilk normality test
##
## data: data$TinggiBadan
## W = 0.97681, p-value = 0.9459
shapiro.test(data$BeratBadan)
##
## Shapiro-Wilk normality test
##
## data: data$BeratBadan
## W = 0.9713, p-value = 0.9026
cor.test(data$TinggiBadan, data$BeratBadan, method = "pearson")
##
## Pearson's product-moment correlation
##
## data: data$TinggiBadan and data$BeratBadan
## t = 35.534, df = 8, p-value = 4.307e-10
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.9862015 0.9992826
## sample estimates:
## cor
## 0.9968471