Tugas_4 RLB_Dataset1

#Membuat Scatter Plot X1,X2,X3
#Scatter Plot X1
library(ggplot2)
ggplot(data, aes(x = X1, y = Y)) + 
  geom_point(size = 3, color = "black") + 
  geom_smooth(
    method = "lm",
    se = FALSE,
    color = "darkblue",
    linetype = "solid",
    linewidth = 1
  ) + 
  labs(
    title = "Scatter plot Y terhadap X1",
    x = "X1",
    y = "y"
  ) +
  theme_minimal()

## `geom_smooth()` using formula = 'y ~ x'

#Scatter Plot X2
library(ggplot2)
ggplot(data, aes(x = X2, y = Y)) + 
  geom_point(size = 3, color = "black") + 
  geom_smooth(
    method = "lm",
    se = FALSE,
    color = "darkviolet",
    linetype = "solid",
    linewidth = 1
  ) + 
  labs(
    title = "Scatter plot Y terhadap X2",
    x = "X2",
    y = "y"
  ) +
  theme_minimal()

## `geom_smooth()` using formula = 'y ~ x'

#Scatter Plot X3
library(ggplot2)
ggplot(data, aes(x = X3, y = Y)) + 
  geom_point(size = 3, color = "black") + 
  geom_smooth(
    method = "lm",
    se = FALSE,
    color = "darkgreen",
    linetype = "solid",
    linewidth = 1
  ) + 
  labs(
    title = "Scatter plot Y terhadap X3",
    x = "X3",
    y = "y"
  ) +
  theme_minimal()

## `geom_smooth()` using formula = 'y ~ x'

#Analisis Model Regresi
#Estimasi Parameter
rlb=lm(Y~X1+X2+X3, data=data)
summary(rlb)

## 
## Call:
## lm(formula = Y ~ X1 + X2 + X3, data = data)
## 
## Residuals:
##       1       2       3       4       5       6       7       8       9 
## -0.2981 -0.2142  0.4714 -0.4746 -0.9524  0.4761  0.9274  0.7466 -0.6823 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 21.87353    4.07389   5.369  0.00302 ** 
## X1           0.41277    0.03067  13.460 4.05e-05 ***
## X2           2.20267    0.47199   4.667  0.00550 ** 
## X3          -0.07895    0.15551  -0.508  0.63330    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8484 on 5 degrees of freedom
## Multiple R-squared:  0.9888, Adjusted R-squared:  0.9821 
## F-statistic: 147.1 on 3 and 5 DF,  p-value: 2.696e-05

#Pendekatan Multikolinieritas
library(car)

## Loading required package: carData

vif(rlb)

##       X1       X2       X3 
## 1.403219 3.093549 2.607991

#Uji Asumsi
#Uji Non-Heterokedastisitas
#Breusch-Pagan Test
library(lmtest)

## Loading required package: zoo

## 
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric

bptest(rlb)

## 
##  studentized Breusch-Pagan test
## 
## data:  rlb
## BP = 0.52402, df = 3, p-value = 0.9136

#Uji Glejser
resrlb = resid(rlb)
resrlb

##          1          2          3          4          5          6          7 
## -0.2980857 -0.2141744  0.4714083 -0.4745650 -0.9524106  0.4760829  0.9274363 
##          8          9 
##  0.7466056 -0.6822975

aresrlb1 = abs(resrlb)
aresrlb1

##         1         2         3         4         5         6         7         8 
## 0.2980857 0.2141744 0.4714083 0.4745650 0.9524106 0.4760829 0.9274363 0.7466056 
##         9 
## 0.6822975

aresrlb2 = lm(aresrlb1~X1+X2+X3, data = data)
summary(aresrlb2)

## 
## Call:
## lm(formula = aresrlb1 ~ X1 + X2 + X3, data = data)
## 
## Residuals:
##        1        2        3        4        5        6        7        8 
## -0.24968 -0.37078 -0.06495 -0.01987  0.35744 -0.17650  0.37477  0.12608 
##        9 
##  0.02347 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)
## (Intercept)  1.000238   1.548904   0.646    0.547
## X1           0.002391   0.011660   0.205    0.846
## X2           0.035130   0.179452   0.196    0.853
## X3          -0.024936   0.059127  -0.422    0.691
## 
## Residual standard error: 0.3226 on 5 degrees of freedom
## Multiple R-squared:  0.04467,    Adjusted R-squared:  -0.5285 
## F-statistic: 0.07793 on 3 and 5 DF,  p-value: 0.9692

#Uji Non-Autokorelasi
#Durbin-Watson Test
library(lmtest)
dwtest(rlb)

## 
##  Durbin-Watson test
## 
## data:  rlb
## DW = 1.6446, p-value = 0.1165
## alternative hypothesis: true autocorrelation is greater than 0

#Bruesch-Godfrey Test
bgtest(rlb)

## 
##  Breusch-Godfrey test for serial correlation of order up to 1
## 
## data:  rlb
## LM test = 0.51249, df = 1, p-value = 0.4741

#Uji Normalitas Residual
shapiro.test(resid(rlb))

## 
##  Shapiro-Wilk normality test
## 
## data:  resid(rlb)
## W = 0.93208, p-value = 0.5014