Load Package

library(readxl)

antar=read_excel("E:\\Praktikum-5.xlsx",sheet="Sheet1")
eeo=read_excel("E:\\praktikum 7 EEO.xlsx",sheet="data")

Data 1

#===============data antar===================
X_1=cbind(antar$x1,antar$x2);str(X_1)
##  num [1:25, 1:2] 7 3 3 4 6 7 2 7 30 5 ...
y_1=antar$y

Nilai eigen

mat_xx1=t(X_1) %*% X_1
ac_antar=eigen(mat_xx1)
ac_antar
## eigen() decomposition
## $values
## [1] 6728353.8952     389.1048
## 
## $vectors
##           [,1]       [,2]
## [1,] 0.0199058 -0.9998019
## [2,] 0.9998019  0.0199058

Condition Number of X'X

#Condition Number of X'X:
k_antar=max(ac_antar$values)/min(ac_antar$values)
k_antar
## [1] 17291.88

Dari hasil Condition Number terlihat bahwa nilainya lebih besar 1000 sehingga diindikasikan terjadi multikolinier kuat ##Condition indices of X'X

#Condition Number of X'X:
k_1=max(ac_antar$values)/min(ac_antar$values[1])
k_1
## [1] 1
k_2=max(ac_antar$values)/min(ac_antar$values[2])
k_2
## [1] 17291.88

multikolinier terjadi pada peubah penjelas ke 2

Data 2

#==================data eeo==================
X_2=cbind(eeo$FAM,eeo$PEER,eeo$SCHOOL)
y_2=eeo$ACHV

Nilai eigen

mat_xx2=t(X_2) %*% X_2
ac_eeo=eigen(mat_xx2)
ac_eeo
## eigen() decomposition
## $values
## [1] 209.3295610   2.7303093   0.5833274
## 
## $vectors
##           [,1]       [,2]       [,3]
## [1,] 0.6174223  0.6698093 -0.4124865
## [2,] 0.5243779 -0.7413256 -0.4188844
## [3,] 0.5863595 -0.0423298  0.8089442

Condition Number of X'X:

k_eeo=max(ac_eeo$values)/min(ac_eeo$values)
k_eeo
## [1] 358.8543

Dari hasil Condition Number terlihat bahwa nilainya berada diantara 100 dan 1000 sehingga diindikasikan terjadi multikolinier lemah

Condition Indices of X'X:

k1_eeo=max(ac_eeo$values)/ac_eeo$values[1]
k1_eeo
## [1] 1
k2_eeo=max(ac_eeo$values)/ac_eeo$values[2]
k2_eeo
## [1] 76.66881
k3_eeo=max(ac_eeo$values)/ac_eeo$values[3]
k3_eeo
## [1] 358.8543

multikolinier terjadi pada peubah penjelas ke 3