Workshop #06

Wesley Burr & Dave Riegert

27/02/2023

library(psych)

## Warning: package ’psych’ was built under R version 4.3.2

library(rstatix)

## Warning: package ’rstatix’ was built under R version 4.3.2

##
## Attaching package: ’rstatix’

## The following object is masked from ’package:stats’:

##
## filter

Loading the data

data1 <- read.csv("data1.csv", header = TRUE)

#str(data1)

Company <- c(data1$rowname)

Alt <- c(data1$Alt)

EBIT <- c(data1$EBIT)

NWC <- c(data1$NWC)

Sale <- c(data1$Sale)

equity <- c(data1$equity)

RE <- c(data1$RE)

class.df<- data.frame(Alt,EBIT,NWC,Sale,equity,RE)

1
View(class.df)

Desribe the data

summary(data1)

## rowname Alt EBIT NWC

## Length:640 Min. :-9.610 Min. :-0.3800 Min. :-0.3400
## Class :character 1st Qu.: 3.200 1st Qu.: 0.0600 1st Qu.: 0.0600
## Mode :character Median : 5.485 Median : 0.1100 Median : 0.1700
## Mean : 7.438 Mean : 0.1142 Mean : 0.1997
## 3rd Qu.: 8.602 3rd Qu.: 0.1600 3rd Qu.: 0.3225
## Max. :58.580 Max. : 0.5400 Max. : 1.2200
## Sale equity RE
## Min. :0.1400 Min. :0.0000 Min. :-10.4700
## 1st Qu.:0.4700 1st Qu.:0.0400 1st Qu.: 0.0000
## Median :0.6150 Median :0.1100 Median : 0.1750
## Mean :0.6965 Mean :0.1569 Mean : 0.4215
## 3rd Qu.:0.7900 3rd Qu.:0.2000 3rd Qu.: 0.4725
## Max. :2.2700 Max. :1.5300 Max. : 47.0600

Correlation test

corr.test(x=class.df, method = "pearson")

## Call:corr.test(x = class.df, method = "pearson")

## Correlation matrix
## Alt EBIT NWC Sale equity RE
## Alt 1.00 0.21 0.15 0.06 -0.40 0.08
## EBIT 0.21 1.00 0.19 0.21 -0.02 -0.03
## NWC 0.15 0.19 1.00 0.12 -0.06 -0.10
## Sale 0.06 0.21 0.12 1.00 0.08 -0.04
## equity -0.40 -0.02 -0.06 0.08 1.00 -0.05
## RE 0.08 -0.03 -0.10 -0.04 -0.05 1.00
## Sample Size
## [1] 640
## Probability values (Entries above the diagonal are adjusted for multiple tests.)
## Alt EBIT NWC Sale equity RE
## Alt 0.00 0.00 0.00 0.72 0.00 0.32
## EBIT 0.00 0.00 0.00 0.00 0.99 0.99
## NWC 0.00 0.00 0.00 0.03 0.68 0.08
## Sale 0.14 0.00 0.00 0.00 0.32 0.93
## equity 0.00 0.66 0.11 0.04 0.00 0.93
## RE 0.04 0.49 0.01 0.28 0.23 0.00
##
## To see confidence intervals of the correlations, print with the short=FALSE option

2
Correlation matrices

cor_mat(class.df)

## # A tibble: 6 x 7
## rowname Alt EBIT NWC Sale equity RE
## * <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Alt 1 0.21 0.15 0.058 -0.4 0.081
## 2 EBIT 0.21 1 0.19 0.21 -0.017 -0.027
## 3 NWC 0.15 0.19 1 0.12 -0.063 -0.1
## 4 Sale 0.058 0.21 0.12 1 0.08 -0.043
## 5 equity -0.4 -0.017 -0.063 0.08 1 -0.047
## 6 RE 0.081 -0.027 -0.1 -0.043 -0.047 1

#plot the data and see the correlation

base.matrix <- cor_mat(class.df)
tidy.matrix <- pull_lower_triangle(base.matrix)
cor_plot(tidy.matrix)
t
Al

EBIT
IT
EB

NWC
C
W
N

Sale
le
Sa

equity
ty
ui
eq

RE

−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1

Identify the outliner

3
identify_outliers(data1, Alt)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 AKAM UW Equity 21.03 0.24 0.21 0.49 0.05 -0.56 TRUE FALSE
## 2 EPAM UN Equity 27.72 0.23 0.27 1.23 0.00 0.22 TRUE TRUE
## 3 MPWR UW Equity 27.41 0.12 0.31 0.71 0.00 0.15 TRUE TRUE
## 4 SWKS UW Equity 22.08 0.02 0.13 0.77 0.00 0.13 TRUE FALSE
## 5 VRSN UW Equity -8.11 0.08 0.41 0.53 0.22 -10.47 TRUE FALSE
## 6 ZBRA UW Equity 19.63 0.14 -0.04 0.30 0.79 0.25 TRUE FALSE
## 7 EPAM UN Equity 21.67 0.27 0.28 1.17 0.01 0.31 TRUE FALSE
## 8 MPWR UW Equity 30.78 0.11 0.36 0.77 0.00 0.23 TRUE TRUE
## 9 SWKS UW Equity 19.43 0.01 0.09 0.88 0.00 0.21 TRUE FALSE
## 10 TYL UN Equity 19.07 0.12 0.29 0.44 0.01 0.18 TRUE FALSE
## 11 VRSN UW Equity -9.61 0.08 0.43 0.45 0.19 -8.09 TRUE FALSE
## 12 MPWR UW Equity 34.13 0.15 0.32 0.77 0.00 0.23 TRUE TRUE
## 13 SWKS UW Equity 30.73 0.01 0.09 0.85 0.00 0.38 TRUE TRUE
## 14 VRSN UW Equity -8.35 0.03 0.28 0.45 0.19 -8.06 TRUE FALSE
## 15 ANET UN Equity 18.84 0.09 0.27 0.65 0.01 0.15 TRUE FALSE
## 16 ANSS UW Equity 17.24 -0.05 0.59 0.35 0.00 0.64 TRUE FALSE
## 17 EPAM UN Equity 27.31 0.30 0.34 1.25 0.01 0.37 TRUE TRUE
## 18 MPWR UW Equity 37.24 0.20 0.34 0.76 0.00 0.20 TRUE TRUE
## 19 NVDA UW Equity 20.53 0.16 0.49 0.70 0.04 0.44 TRUE FALSE
## 20 SWKS UW Equity 24.62 0.02 -0.03 0.80 0.00 0.49 TRUE FALSE
## 21 TXN UW Equity 17.97 0.17 0.29 0.81 0.05 1.90 TRUE FALSE
## 22 VRSN UW Equity -7.58 0.07 0.36 0.49 0.24 -7.98 TRUE FALSE
## 23 EPAM UN Equity 17.05 0.39 0.30 1.14 0.00 0.28 TRUE FALSE
## 24 INTU UW Equity 17.30 0.09 0.35 1.17 0.01 1.42 TRUE FALSE
## 25 MPWR UW Equity 24.81 0.25 0.50 0.73 0.00 0.15 TRUE TRUE
## 26 NVDA UW Equity 16.96 0.26 0.50 0.86 0.01 0.54 TRUE FALSE
## 27 ANSS UW Equity 26.77 -0.01 0.24 0.31 0.03 0.48 TRUE TRUE
## 28 EPAM UN Equity 23.53 0.32 0.37 1.02 0.02 0.23 TRUE FALSE
## 29 INTU UW Equity 22.34 0.19 0.69 1.08 0.01 1.29 TRUE FALSE
## 30 MPWR UW Equity 32.58 0.18 0.30 0.66 0.00 0.15 TRUE TRUE
## 31 NVDA UW Equity 26.07 0.21 0.52 0.88 0.02 0.66 TRUE TRUE
## 32 SWKS UW Equity 21.40 0.15 -0.04 0.70 0.00 0.77 TRUE FALSE
## 33 TXN UW Equity 16.94 0.15 0.21 0.80 0.05 1.92 TRUE FALSE
## 34 TYL UN Equity 18.52 0.15 0.07 0.50 0.00 0.27 TRUE FALSE
## 35 VRSN UW Equity -6.49 0.14 0.20 0.66 0.08 -9.54 TRUE FALSE
## 36 ADBE UW Equity 17.78 0.11 0.12 0.53 0.02 0.07 TRUE FALSE
## 37 AMD UW Equity 20.95 0.47 0.13 1.09 0.01 0.20 TRUE FALSE
## 38 ANSS UW Equity 18.13 -0.02 0.27 0.28 0.03 -1.28 TRUE FALSE
## 39 CDNS UW Equity 20.99 0.18 0.36 0.68 0.01 0.30 TRUE FALSE
## 40 ENPH UQ Equity 31.78 0.08 0.15 0.65 0.02 1.41 TRUE TRUE
## 41 EPAM UN Equity 21.60 0.30 0.41 0.98 0.01 0.70 TRUE FALSE
## 42 MPWR UW Equity 56.24 0.20 0.40 0.70 0.00 0.00 TRUE TRUE
## 43 NOW UN Equity 17.10 0.03 0.11 0.52 0.02 -0.32 TRUE FALSE
## 44 NVDA UW Equity 41.28 0.18 0.53 0.63 0.02 2.90 TRUE TRUE
## 45 SWKS UW Equity 20.12 0.12 -0.11 0.66 0.01 0.15 TRUE FALSE
## 46 TXN UW Equity 17.25 0.06 0.21 0.75 0.05 -0.22 TRUE FALSE
## 47 TYL UN Equity 21.56 0.13 0.08 0.43 0.00 -1.34 TRUE FALSE
## 48 VRSN UW Equity -6.01 0.13 0.09 0.72 0.07 1.71 TRUE FALSE
## 49 ADBE UW Equity 18.56 0.09 0.19 0.58 0.02 0.72 TRUE FALSE
## 50 AMD UW Equity 34.78 0.44 0.10 1.32 0.00 -0.57 TRUE TRUE

4
## 51 ANET UN Equity 21.20 0.13 0.13 0.51 0.00 0.31 TRUE FALSE
## 52 CDNS UW Equity 24.27 0.21 0.29 0.68 0.01 0.40 TRUE FALSE
## 53 ENPH UQ Equity 22.25 0.10 0.14 0.66 0.04 -0.09 TRUE FALSE
## 54 EPAM UN Equity 33.88 0.36 0.45 1.07 0.01 0.29 TRUE TRUE
## 55 INTU UW Equity 25.08 0.16 0.65 0.62 0.02 0.70 TRUE TRUE
## 56 MPWR UW Equity 58.58 0.27 0.38 0.76 0.00 0.14 TRUE TRUE
## 57 NVDA UW Equity 40.40 0.28 0.51 0.58 0.02 0.52 TRUE TRUE
## 58 TXN UW Equity 17.22 0.11 0.23 0.74 0.05 1.62 TRUE FALSE
## 59 TYL UN Equity 24.01 0.14 0.11 0.34 0.06 0.19 TRUE FALSE
## 60 VRSN UW Equity -4.95 0.11 0.27 0.67 0.06 -8.31 TRUE FALSE
## 61 CDNS UW Equity 19.39 0.23 0.20 0.68 0.01 0.46 TRUE FALSE
## 62 MPWR UW Equity 31.51 0.29 0.32 0.76 0.00 0.47 TRUE TRUE
## 63 ADBE UW Equity 18.05 0.07 0.09 0.65 0.01 0.25 TRUE FALSE
## 64 ANET UN Equity 26.23 0.15 0.12 0.65 0.00 -0.45 TRUE TRUE
## 65 CDNS UW Equity 22.40 0.21 0.19 0.69 0.02 0.48 TRUE FALSE
## 66 FICO UN Equity 18.37 0.06 0.16 0.96 0.09 0.86 TRUE FALSE
## 67 MPWR UW Equity 49.34 0.25 0.38 0.87 0.00 0.78 TRUE TRUE
## 68 NVDA UW Equity 40.60 0.28 0.48 0.65 0.02 1.66 TRUE TRUE
## 69 VRSN UW Equity -5.36 0.10 0.12 0.82 0.08 2.47 TRUE FALSE

identify_outliers(data1, EBIT)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 ANSS UW Equity 11.52 -0.15 0.41 0.34 0.00 0.38 TRUE FALSE
## 2 GEN UW Equity 2.90 -0.11 0.28 0.31 0.15 -0.15 TRUE FALSE
## 3 TYL UN Equity 14.29 -0.10 0.44 0.87 0.00 0.29 TRUE FALSE
## 4 GLW UN Equity 3.08 -0.16 0.30 0.32 0.22 0.43 TRUE FALSE
## 5 TER UW Equity 0.00 -0.12 0.30 0.64 0.00 0.29 TRUE FALSE
## 6 GLW UN Equity 2.92 -0.16 0.30 0.32 0.22 0.46 TRUE FALSE
## 7 TER UW Equity 0.00 -0.12 0.30 0.64 0.00 0.24 TRUE FALSE
## 8 GEN UW Equity 1.24 -0.21 0.16 0.22 0.44 -0.04 TRUE FALSE
## 9 GLW UN Equity 3.45 -0.11 0.47 0.34 0.17 0.50 TRUE FALSE
## 10 IBM UN Equity 5.31 -0.10 -0.05 0.68 0.22 1.24 TRUE FALSE
## 11 ROP UW Equity 0.00 0.35 0.17 0.26 0.33 0.29 TRUE FALSE
## 12 TER UW Equity 0.00 -0.38 0.46 0.63 0.07 0.17 TRUE TRUE
## 13 AMD UW Equity 2.70 0.40 0.17 1.42 0.07 -1.71 TRUE FALSE
## 14 EPAM UN Equity 17.05 0.39 0.30 1.14 0.00 0.28 TRUE FALSE
## 15 IBM UN Equity 4.61 -0.12 0.07 0.65 0.39 1.24 TRUE FALSE
## 16 TXN UW Equity 14.93 0.35 0.16 0.92 0.06 1.93 TRUE FALSE
## 17 AMD UW Equity 9.18 0.43 0.17 1.12 0.01 -1.29 TRUE FALSE
## 18 EPAM UN Equity 23.53 0.32 0.37 1.02 0.02 0.23 TRUE FALSE
## 19 AMD UW Equity 20.95 0.47 0.13 1.09 0.01 0.20 TRUE TRUE
## 20 AMD UW Equity 34.78 0.44 0.10 1.32 0.00 -0.57 TRUE FALSE
## 21 EPAM UN Equity 33.88 0.36 0.45 1.07 0.01 0.29 TRUE FALSE
## 22 QRVO UW Equity 6.43 0.32 0.46 0.56 0.08 0.00 TRUE FALSE
## 23 ZBRA UW Equity 10.81 0.32 -0.01 0.91 0.04 0.36 TRUE FALSE
## 24 AMD UW Equity 15.34 0.44 0.10 1.32 0.00 0.18 TRUE FALSE
## 25 EPAM UN Equity 14.44 0.36 0.45 1.07 0.01 0.97 TRUE FALSE
## 26 HPQ UN Equity 2.68 0.34 -0.05 1.63 0.45 0.00 TRUE FALSE
## 27 ORCL UN Equity 1.92 0.32 0.18 0.39 0.41 -0.01 TRUE FALSE
## 28 QRVO UW Equity 4.71 0.32 0.46 0.62 0.16 -0.09 TRUE FALSE
## 29 ZBRA UW Equity 7.99 0.38 0.11 0.91 0.04 1.98 TRUE FALSE
## 30 AMD UW Equity 13.11 0.54 -0.05 0.35 0.03 0.05 TRUE TRUE
## 31 EPAM UN Equity 13.84 0.37 0.41 1.20 0.01 1.14 TRUE FALSE

5
## 32 HPQ UN Equity 2.48 0.32 0.00 1.45 0.41 0.00 TRUE FALSE
## 33 ZBRA UW Equity 5.97 0.41 0.12 0.77 0.17 1.80 TRUE FALSE

identify_outliers(data1, NWC)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 CDW UW Equity 5.59 0.20 0.82 2.01 0.39 -0.24 TRUE FALSE
## 2 FSLR UW Equity 3.85 0.13 -0.34 0.43 0.06 0.41 TRUE FALSE
## 3 INTU UW Equity 15.96 0.14 1.22 1.28 0.01 1.64 TRUE TRUE
## 4 KLAC UW Equity 4.79 0.11 0.98 0.63 0.20 0.05 TRUE FALSE
## 5 SNPS UW Equity 6.45 0.14 1.01 0.50 0.01 0.36 TRUE FALSE
## 6 ORCL UN Equity 3.13 0.02 1.00 0.29 0.32 0.20 TRUE FALSE

identify_outliers(data1, Sale)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 ACN UN Equity 7.49 0.04 0.06 1.78 0.00 0.56 TRUE TRUE
## 2 AMD UW Equity 0.00 0.30 0.03 1.46 1.08 -1.64 TRUE FALSE
## 3 CDW UW Equity 5.18 0.07 0.65 1.98 0.53 -0.34 TRUE TRUE
## 4 ENPH UQ Equity 0.00 0.12 0.23 2.26 0.00 -0.83 TRUE TRUE
## 5 JBL UN Equity 2.88 0.04 0.27 1.86 0.40 0.13 TRUE TRUE
## 6 STX UW Equity 4.27 0.11 0.40 1.45 0.21 -0.19 TRUE FALSE
## 7 ACN UN Equity 8.26 0.04 0.09 1.81 0.00 0.66 TRUE TRUE
## 8 AMD UW Equity -0.43 0.26 0.20 1.29 0.97 -2.04 TRUE FALSE
## 9 CDW UW Equity 5.62 0.11 0.67 1.92 0.46 -0.27 TRUE TRUE
## 10 ENPH UQ Equity 2.08 0.11 0.22 2.16 0.11 -0.98 TRUE TRUE
## 11 JBL UN Equity 2.85 0.16 0.25 1.87 0.45 0.14 TRUE TRUE
## 12 STX UW Equity 3.36 0.09 0.48 1.40 0.27 -0.27 TRUE FALSE
## 13 ACN UN Equity 7.65 0.12 -0.01 1.69 0.00 0.65 TRUE FALSE
## 14 AMD UW Equity -0.59 0.26 0.20 1.29 0.97 -2.15 TRUE FALSE
## 15 CDW UW Equity 5.66 0.10 0.50 1.92 0.46 -0.26 TRUE TRUE
## 16 ENPH UQ Equity 0.90 0.11 0.22 2.16 0.11 -0.97 TRUE TRUE
## 17 HPQ UN Equity 2.80 0.19 0.09 1.66 0.27 1.11 TRUE FALSE
## 18 JBL UN Equity 2.74 -0.01 0.08 1.78 0.53 0.14 TRUE TRUE
## 19 STX UW Equity 2.25 0.06 0.49 1.36 0.58 -0.33 TRUE FALSE
## 20 ACN UN Equity 8.37 0.10 -0.11 1.62 0.00 0.35 TRUE FALSE
## 21 AMD UW Equity -0.08 0.29 0.16 1.30 0.14 -2.20 TRUE FALSE
## 22 CDW UW Equity 5.59 0.20 0.82 2.01 0.39 -0.24 TRUE TRUE
## 23 ENPH UQ Equity -0.45 0.11 0.30 1.97 0.54 -1.12 TRUE TRUE
## 24 HPQ UN Equity 2.78 0.16 0.04 1.58 0.22 -0.11 TRUE FALSE
## 25 INTU UW Equity 15.96 0.14 1.22 1.28 0.01 1.64 TRUE FALSE
## 26 JBL UN Equity 2.53 0.04 0.12 1.72 0.37 0.15 TRUE FALSE
## 27 ACN UN Equity 8.08 0.11 0.00 1.68 0.00 0.29 TRUE FALSE
## 28 AMD UW Equity 2.70 0.40 0.17 1.42 0.07 -1.71 TRUE FALSE
## 29 CDW UW Equity 5.83 0.29 0.60 2.27 0.27 -0.23 TRUE TRUE
## 30 HPQ UN Equity 2.91 0.19 0.08 1.69 0.16 -0.07 TRUE FALSE
## 31 JBL UN Equity 2.63 0.10 0.27 1.83 0.52 0.14 TRUE TRUE
## 32 ACN UN Equity 9.06 0.12 -0.15 1.45 0.00 0.27 TRUE FALSE
## 33 CDW UW Equity 6.51 0.29 0.69 2.25 0.17 -0.23 TRUE TRUE
## 34 HPQ UN Equity 2.95 0.19 0.17 1.76 0.20 -0.01 TRUE TRUE
## 35 JBL UN Equity 2.84 0.05 0.04 1.95 0.56 0.14 TRUE TRUE
## 36 CDW UW Equity 5.72 0.16 0.69 1.98 0.22 0.00 TRUE TRUE
## 37 HPQ UN Equity 2.63 0.20 0.12 1.63 0.32 0.00 TRUE FALSE

6
## 38 JBL UN Equity 2.67 0.05 0.07 1.89 0.64 0.05 TRUE TRUE
## 39 AMD UW Equity 34.78 0.44 0.10 1.32 0.00 -0.57 TRUE FALSE
## 40 CDW UW Equity 5.61 0.16 0.42 1.58 0.26 0.00 TRUE FALSE
## 41 HPQ UN Equity 2.84 0.31 0.03 1.64 0.27 0.00 TRUE FALSE
## 42 JBL UN Equity 2.75 0.11 0.14 1.76 0.41 0.12 TRUE TRUE
## 43 ACN UN Equity 8.12 0.12 -0.17 1.30 0.02 0.07 TRUE FALSE
## 44 AMD UW Equity 15.34 0.44 0.10 1.32 0.00 0.18 TRUE FALSE
## 45 CDW UW Equity 4.44 0.23 0.55 1.58 0.26 0.00 TRUE FALSE
## 46 HPQ UN Equity 2.68 0.34 -0.05 1.63 0.45 0.00 TRUE FALSE
## 47 JBL UN Equity 2.62 0.19 0.16 1.70 0.46 0.07 TRUE FALSE
## 48 STX UW Equity 1.87 0.11 0.07 1.30 0.39 6.38 TRUE FALSE
## 49 CDW UW Equity 5.55 0.10 0.40 1.81 0.25 0.00 TRUE TRUE
## 50 HPQ UN Equity 2.48 0.32 0.00 1.45 0.41 0.00 TRUE FALSE
## 51 JBL UN Equity 3.18 0.22 0.18 1.79 0.24 0.09 TRUE TRUE

identify_outliers(data1, equity)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 AMD UW Equity 0.00 0.30 0.03 1.46 1.08 -1.64 TRUE TRUE
## 2 CDW UW Equity 5.18 0.07 0.65 1.98 0.53 -0.34 TRUE FALSE
## 3 ZBRA UW Equity 19.63 0.14 -0.04 0.30 0.79 0.25 TRUE TRUE
## 4 AMD UW Equity -0.43 0.26 0.20 1.29 0.97 -2.04 TRUE TRUE
## 5 CDW UW Equity 5.62 0.11 0.67 1.92 0.46 -0.27 TRUE FALSE
## 6 HPE UN Equity 1.64 0.08 0.22 0.39 0.62 0.00 TRUE FALSE
## 7 JBL UN Equity 2.85 0.16 0.25 1.87 0.45 0.14 TRUE FALSE
## 8 ZBRA UW Equity 2.91 0.11 0.03 0.72 0.83 0.31 TRUE TRUE
## 9 AMD UW Equity -0.59 0.26 0.20 1.29 0.97 -2.15 TRUE TRUE
## 10 CDW UW Equity 5.66 0.10 0.50 1.92 0.46 -0.26 TRUE FALSE
## 11 JBL UN Equity 2.74 -0.01 0.08 1.78 0.53 0.14 TRUE FALSE
## 12 MU UW Equity 1.90 0.10 0.33 0.45 0.54 0.20 TRUE FALSE
## 13 STX UW Equity 2.25 0.06 0.49 1.36 0.58 -0.33 TRUE FALSE
## 14 WDC UW Equity 2.42 0.01 0.33 0.40 1.29 0.28 TRUE TRUE
## 15 ZBRA UW Equity 4.23 0.14 0.02 0.72 0.83 0.30 TRUE TRUE
## 16 ENPH UQ Equity -0.45 0.11 0.30 1.97 0.54 -1.12 TRUE FALSE
## 17 HPE UN Equity 0.91 0.15 0.04 0.47 0.63 0.05 TRUE FALSE
## 18 ON UW Equity 2.12 0.05 0.05 0.56 0.68 -0.10 TRUE TRUE
## 19 WDC UW Equity 3.62 0.02 0.45 0.64 0.50 0.30 TRUE FALSE
## 20 ZBRA UW Equity 4.35 0.14 -0.06 0.77 0.58 0.30 TRUE FALSE
## 21 HPE UN Equity 1.09 0.15 0.20 0.56 0.56 -0.13 TRUE FALSE
## 22 JBL UN Equity 2.63 0.10 0.27 1.83 0.52 0.14 TRUE FALSE
## 23 WDC UW Equity 3.53 0.08 0.23 0.71 0.49 0.30 TRUE FALSE
## 24 HPE UN Equity 0.91 0.16 -0.08 0.56 0.65 -0.11 TRUE FALSE
## 25 IBM UN Equity 4.90 -0.01 -0.14 0.38 0.57 1.01 TRUE FALSE
## 26 JBL UN Equity 2.84 0.05 0.04 1.95 0.56 0.14 TRUE FALSE
## 27 MCHP UW Equity 3.55 0.05 0.14 0.29 0.52 0.08 TRUE FALSE
## 28 WDC UW Equity 2.94 0.09 0.34 0.63 0.75 0.33 TRUE TRUE
## 29 HPE UN Equity 0.57 0.17 0.11 0.50 1.53 0.30 TRUE TRUE
## 30 IBM UN Equity 4.82 0.06 -0.09 0.35 0.58 -0.01 TRUE FALSE
## 31 JBL UN Equity 2.67 0.05 0.07 1.89 0.64 0.05 TRUE FALSE
## 32 MCHP UW Equity 4.37 0.02 0.14 0.30 0.58 -0.01 TRUE FALSE
## 33 ORCL UN Equity 2.72 0.18 0.28 0.34 0.45 -0.02 TRUE FALSE
## 34 WDC UW Equity 2.87 0.13 0.23 0.65 0.77 0.06 TRUE TRUE
## 35 HPE UN Equity 0.89 0.21 0.06 0.48 0.77 -0.15 TRUE TRUE
## 36 IBM UN Equity 4.42 0.09 -0.09 0.43 0.46 1.23 TRUE FALSE

7
## 37 HPE UN Equity 0.85 0.22 0.03 0.50 0.74 0.28 TRUE TRUE
## 38 HPQ UN Equity 2.68 0.34 -0.05 1.63 0.45 0.00 TRUE FALSE
## 39 IBM UN Equity 5.20 0.07 -0.14 0.43 0.46 -0.01 TRUE FALSE
## 40 JBL UN Equity 2.62 0.19 0.16 1.70 0.46 0.07 TRUE FALSE
## 41 WDC UW Equity 3.20 0.11 0.22 0.72 0.54 0.08 TRUE FALSE
## 42 GEN UW Equity 2.40 0.03 0.05 0.21 0.89 0.55 TRUE TRUE
## 43 HPE UN Equity 1.07 0.22 0.10 0.51 0.68 0.29 TRUE TRUE
## 44 STX UW Equity 0.35 0.19 0.00 0.98 0.45 11.15 TRUE FALSE
## 45 WDC UW Equity 2.30 0.16 0.12 0.50 0.60 0.10 TRUE FALSE

identify_outliers(data1, RE)

## rowname Alt EBIT NWC Sale equity RE is.outlier is.extreme

## 1 AMD UW Equity 0.00 0.30 0.03 1.46 1.08 -1.64 TRUE TRUE
## 2 ENPH UQ Equity 0.00 0.12 0.23 2.26 0.00 -0.83 TRUE FALSE
## 3 NXPI UW Equity 2.43 0.06 0.20 0.82 0.23 -0.77 TRUE FALSE
## 4 TXN UW Equity 9.67 0.14 0.24 0.75 0.08 1.57 TRUE FALSE
## 5 VRSN UW Equity -8.11 0.08 0.41 0.53 0.22 -10.47 TRUE TRUE
## 6 AMD UW Equity -0.43 0.26 0.20 1.29 0.97 -2.04 TRUE TRUE
## 7 ENPH UQ Equity 2.08 0.11 0.22 2.16 0.11 -0.98 TRUE FALSE
## 8 IBM UN Equity 5.01 0.02 0.26 0.74 0.24 1.21 TRUE FALSE
## 9 TXN UW Equity 11.03 0.12 0.19 0.80 0.07 1.80 TRUE FALSE
## 10 VRSN UW Equity -9.61 0.08 0.43 0.45 0.19 -8.09 TRUE TRUE
## 11 AMD UW Equity -0.59 0.26 0.20 1.29 0.97 -2.15 TRUE TRUE
## 12 ENPH UQ Equity 0.90 0.11 0.22 2.16 0.11 -0.97 TRUE FALSE
## 13 IBM UN Equity 5.60 0.00 0.25 0.74 0.24 1.25 TRUE FALSE
## 14 INTU UW Equity 11.66 0.13 0.64 1.10 0.03 1.42 TRUE FALSE
## 15 TXN UW Equity 14.40 0.01 0.17 0.80 0.07 1.83 TRUE FALSE
## 16 VRSN UW Equity -8.35 0.03 0.28 0.45 0.19 -8.06 TRUE TRUE
## 17 AMD UW Equity -0.08 0.29 0.16 1.30 0.14 -2.20 TRUE TRUE
## 18 ENPH UQ Equity -0.45 0.11 0.30 1.97 0.54 -1.12 TRUE FALSE
## 19 IBM UN Equity 5.31 -0.10 -0.05 0.68 0.22 1.24 TRUE FALSE
## 20 INTU UW Equity 15.96 0.14 1.22 1.28 0.01 1.64 TRUE FALSE
## 21 TXN UW Equity 17.97 0.17 0.29 0.81 0.05 1.90 TRUE TRUE
## 22 VRSN UW Equity -7.58 0.07 0.36 0.49 0.24 -7.98 TRUE TRUE
## 23 AMD UW Equity 2.70 0.40 0.17 1.42 0.07 -1.71 TRUE TRUE
## 24 ENPH UQ Equity 0.42 0.07 0.17 0.93 0.22 -0.74 TRUE FALSE
## 25 FICO UN Equity 12.49 0.06 0.24 0.75 0.12 1.20 TRUE FALSE
## 26 IBM UN Equity 4.61 -0.12 0.07 0.65 0.39 1.24 TRUE FALSE
## 27 INTU UW Equity 17.30 0.09 0.35 1.17 0.01 1.42 TRUE FALSE
## 28 TXN UW Equity 14.93 0.35 0.16 0.92 0.06 1.93 TRUE TRUE
## 29 VRSN UW Equity -4.44 0.12 0.21 0.63 0.10 -9.50 TRUE TRUE
## 30 AMD UW Equity 9.18 0.43 0.17 1.12 0.01 -1.29 TRUE FALSE
## 31 FICO UN Equity 13.32 0.05 0.19 0.81 0.09 1.21 TRUE FALSE
## 32 INTU UW Equity 22.34 0.19 0.69 1.08 0.01 1.29 TRUE FALSE
## 33 TXN UW Equity 16.94 0.15 0.21 0.80 0.05 1.92 TRUE TRUE
## 34 VRSN UW Equity -6.49 0.14 0.20 0.66 0.08 -9.54 TRUE TRUE
## 35 ADSK UW Equity 6.52 0.13 0.26 0.53 0.06 -1.20 TRUE FALSE
## 36 AKAM UW Equity 4.33 0.18 0.14 0.41 0.16 1.48 TRUE FALSE
## 37 ANET UN Equity 12.92 0.14 0.19 0.49 0.00 9.69 TRUE TRUE
## 38 ANSS UW Equity 18.13 -0.02 0.27 0.28 0.03 -1.28 TRUE FALSE
## 39 ENPH UQ Equity 31.78 0.08 0.15 0.65 0.02 1.41 TRUE FALSE
## 40 FFIV UW Equity 5.35 0.01 0.18 0.50 0.10 2.69 TRUE TRUE
## 41 FSLR UW Equity 3.75 0.07 -0.09 0.38 0.05 3.43 TRUE TRUE

8
## 42 FTNT UW Equity 6.99 0.01 0.02 0.64 0.00 7.61 TRUE TRUE
## 43 IT UN Equity 3.58 0.12 0.30 0.56 0.20 5.18 TRUE TRUE
## 44 KLAC UW Equity 6.26 0.13 0.58 0.63 0.12 1.60 TRUE FALSE
## 45 NVDA UW Equity 41.28 0.18 0.53 0.63 0.02 2.90 TRUE TRUE
## 46 ON UW Equity 3.21 0.07 0.05 0.61 0.27 1.41 TRUE FALSE
## 47 PTC UW Equity 5.79 0.05 0.06 0.43 0.13 47.06 TRUE TRUE
## 48 QRVO UW Equity 6.94 0.25 0.47 0.49 0.18 -0.90 TRUE FALSE
## 49 STX UW Equity 2.60 0.04 0.03 1.18 0.35 2.70 TRUE TRUE
## 50 TYL UN Equity 21.56 0.13 0.08 0.43 0.00 -1.34 TRUE FALSE
## 51 VRSN UW Equity -6.01 0.13 0.09 0.72 0.07 1.71 TRUE FALSE
## 52 ZBRA UW Equity 9.60 0.18 0.07 0.83 0.07 1.79 TRUE FALSE
## 53 FICO UN Equity 12.74 0.04 0.12 0.84 0.12 1.40 TRUE FALSE
## 54 IBM UN Equity 4.42 0.09 -0.09 0.43 0.46 1.23 TRUE FALSE
## 55 TXN UW Equity 17.22 0.11 0.23 0.74 0.05 1.62 TRUE FALSE
## 56 VRSN UW Equity -4.95 0.11 0.27 0.67 0.06 -8.31 TRUE TRUE
## 57 AKAM UW Equity 3.82 0.20 0.09 0.43 0.15 1.31 TRUE FALSE
## 58 ANSS UW Equity 11.79 -0.02 -0.15 0.30 0.03 -0.88 TRUE FALSE
## 59 ENPH UQ Equity 14.63 0.10 0.14 0.66 0.04 1.32 TRUE FALSE
## 60 FFIV UW Equity 4.33 0.09 0.16 0.51 0.08 3.58 TRUE TRUE
## 61 FSLR UW Equity 8.26 0.12 -0.10 0.39 0.04 4.35 TRUE TRUE
## 62 FTNT UW Equity 5.60 0.01 0.01 0.56 0.02 6.60 TRUE TRUE
## 63 IT UN Equity 5.23 0.10 0.11 0.64 0.12 5.67 TRUE TRUE
## 64 KEYS UN Equity 7.76 0.04 0.03 0.67 0.06 1.21 TRUE FALSE
## 65 KLAC UW Equity 5.83 0.17 0.57 0.73 0.15 1.90 TRUE TRUE
## 66 NVDA UW Equity 15.26 0.31 0.45 0.61 0.02 1.27 TRUE FALSE
## 67 ON UW Equity 5.60 0.07 0.05 0.70 0.11 1.22 TRUE FALSE
## 68 PTC UW Equity 6.00 0.14 -0.12 0.41 0.13 34.71 TRUE TRUE
## 69 STX UW Equity 1.87 0.11 0.07 1.30 0.39 6.38 TRUE TRUE
## 70 TYL UN Equity 6.55 0.14 0.11 0.34 0.06 -4.25 TRUE TRUE
## 71 VRSN UW Equity -4.47 0.12 0.20 0.67 0.06 1.87 TRUE FALSE
## 72 ZBRA UW Equity 7.99 0.38 0.11 0.91 0.04 1.98 TRUE TRUE
## 73 AKAM UW Equity 4.66 0.17 0.10 0.44 0.24 1.44 TRUE FALSE
## 74 ANSS UW Equity 15.68 0.03 -0.12 0.31 0.04 -0.80 TRUE FALSE
## 75 FFIV UW Equity 5.73 -0.05 0.10 0.54 0.03 3.09 TRUE TRUE
## 76 FSLR UW Equity 5.82 0.15 -0.11 0.32 0.01 4.59 TRUE TRUE
## 77 FTNT UW Equity 5.20 0.01 0.01 0.71 0.03 7.59 TRUE TRUE
## 78 IT UN Equity 6.56 0.10 -0.02 0.75 0.12 6.29 TRUE TRUE
## 79 KEYS UN Equity 7.21 0.05 -0.01 0.63 0.10 2.05 TRUE TRUE
## 80 KLAC UW Equity 7.05 0.26 0.56 0.75 0.09 2.01 TRUE TRUE
## 81 NVDA UW Equity 40.60 0.28 0.48 0.65 0.02 1.66 TRUE FALSE
## 82 PTC UW Equity 5.73 0.08 -0.10 0.33 0.11 24.52 TRUE TRUE
## 83 STX UW Equity 0.35 0.19 0.00 0.98 0.45 11.15 TRUE TRUE
## 84 TYL UN Equity 9.06 0.13 0.05 0.39 0.08 -6.69 TRUE TRUE
## 85 VRSN UW Equity -5.36 0.10 0.12 0.82 0.08 2.47 TRUE TRUE
## 86 ZBRA UW Equity 5.97 0.41 0.12 0.77 0.17 1.80 TRUE FALSE

Test Normality

shapiro_test(data1$Alt)

## # A tibble: 1 x 3
## variable statistic p.value

9
## <chr> <dbl> <dbl>
## 1 data1$Alt 0.776 1.08e-28

shapiro_test(data1$EBIT)

## # A tibble: 1 x 3
## variable statistic p.value
## <chr> <dbl> <dbl>
## 1 data1$EBIT 0.961 4.37e-12

shapiro_test(data1$NWC)

## # A tibble: 1 x 3
## variable statistic p.value
## <chr> <dbl> <dbl>
## 1 data1$NWC 0.960 2.88e-12

shapiro_test(data1$Sale)

## # A tibble: 1 x 3
## variable statistic p.value
## <chr> <dbl> <dbl>
## 1 data1$Sale 0.830 1.50e-25

shapiro_test(data1$equity)

## # A tibble: 1 x 3
## variable statistic p.value
## <chr> <dbl> <dbl>
## 1 data1$equity 0.752 7.22e-30

shapiro_test(data1$RE)

## # A tibble: 1 x 3
## variable statistic p.value
## <chr> <dbl> <dbl>
## 1 data1$RE 0.280 4.00e-44

Non of variable has the normal distribution.

Scatterplot and correlation

To make things easier, we will extract the two variables we are using into their own variable form.

# scatter plot
plot(data1$Alt, data1$EBIT)

10
 0.4
data1$EBIT

0.2
0.0
−0.2
−0.4

−10 0 10 20 30 40 50 60

data1$Alt

# scatter plot
plot(data1$Alt, data1$NWC)

11
 1.0
data1$NWC

0.5
0.0

−10 0 10 20 30 40 50 60

data1$Alt

# scatter plot
plot(data1$Alt, data1$Sale)

12
 2.0
1.5
data1$Sale

1.0
0.5

−10 0 10 20 30 40 50 60

data1$Alt

# scatter plot
plot(data1$Alt, data1$equity)

13
 1.5
1.0
data1$equity

0.5
0.0

−10 0 10 20 30 40 50 60

data1$Alt

# scatter plot
plot(data1$Alt, data1$RE)

14
 40
30
data1$RE

20
10
0
−10

−10 0 10 20 30 40 50 60

data1$Alt

# scatter plot
plot(data1$Alt, data1$EBIT)

15
 0.4
data1$EBIT

0.2
0.0
−0.2
−0.4

−10 0 10 20 30 40 50 60

data1$Alt

The scatterplot shows that there is a relationship between the Alt and explaintory variables, including EBIT,
NWC, Sale, Equity, RE. Let’s check the correlation:

Correlation test

# Correlation
cor.test(data1$Alt,data1$EBIT)

##
## Pearson’s product-moment correlation
##
## data: data1$Alt and data1$EBIT
## t = 5.5406, df = 638, p-value = 4.412e-08
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.1390699 0.2869966
## sample estimates:
## cor
## 0.2142615

# Correlation
cor.test(data1$Alt,data1$NWC)

##

16
## Pearson’s product-moment correlation
##
## data: data1$Alt and data1$NWC
## t = 3.9106, df = 638, p-value = 0.0001019
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.07640358 0.22779825
## sample estimates:
## cor
## 0.1529985

# Correlation
cor.test(data1$Alt,data1$Sale)

##
## Pearson’s product-moment correlation
##
## data: data1$Alt and data1$Sale
## t = 1.4655, df = 638, p-value = 0.1433
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.01966676 0.13481806
## sample estimates:
## cor
## 0.05792239

# Correlation
cor.test(data1$Alt,data1$equity)

##
## Pearson’s product-moment correlation
##
## data: data1$Alt and data1$equity
## t = -10.999, df = 638, p-value < 2.2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4624391 -0.3320195
## sample estimates:
## cor
## -0.3992471

# Correlation
cor.test(data1$Alt,data1$RE)

##
## Pearson’s product-moment correlation
##
## data: data1$Alt and data1$RE
## t = 2.0613, df = 638, p-value = 0.03968
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.003860737 0.157843147

17
## sample estimates:
## cor
## 0.08133727

Fitting the Line of Best Fit

The line of best fit is the statistical estimate (“best guess”) for the coefficients in the model

y = β0 + β1 x

which for us can be rewritten as

Alt = β0 + β1 · EBIT + β2 · NWC + β3 · Sale + β4 · Equity + β5 · RE

# fit the model and pull out the coefficients:

p_mod <- lm(Alt ~ EBIT+ NWC + Sale + equity + RE, data = data1)
p_mod$coefficients

## (Intercept) EBIT NWC Sale equity RE

## 6.9278115 14.2205870 3.7031355 0.8701934 -16.2305955 0.2071267

summary(p_mod)

##
## Call:
## lm(formula = Alt ~ EBIT + NWC + Sale + equity + RE, data = data1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -14.964 -3.372 -1.288 1.764 45.715
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.92781 0.66756 10.378 < 2e-16 ***
## EBIT 14.22059 2.84917 4.991 7.77e-07 ***
## NWC 3.70314 1.38313 2.677 0.00761 **
## Sale 0.87019 0.73486 1.184 0.23680
## equity -16.23060 1.47342 -11.016 < 2e-16 ***
## RE 0.20713 0.09188 2.254 0.02452 *
## ---
## Signif. codes: 0 ’***’ 0.001 ’**’ 0.01 ’*’ 0.05 ’.’ 0.1 ’ ’ 1
##
## Residual standard error: 6.67 on 634 degrees of freedom
## Multiple R-squared: 0.2184, Adjusted R-squared: 0.2122
## F-statistic: 35.42 on 5 and 634 DF, p-value: < 2.2e-16

So from the output of the model, we have the intercept (β0 ) equaling 149.7477, and the slope (β1 ) equalling
0.25924. We can round these - that much precision is silly - giving our model fit as

Alt = 6.92 + 14.2 · EBIT + 3.70 · NWC + 0.87 · Sale − 16.23 · Equity + 0.20 · RE

We can interpret this as follows:

18
 • RunTime coefficient: for every additional minute of run time, our factory will produce 2.8 more items
on average;
• Intercept (-367): It takes 367 items worth of production time to setup the factory to produce units.

plot(p_mod, which = 1:2)

Residuals vs Fitted
10 20 30 40 50

422486

614
Residuals

0
−20

−15 −10 −5 0 5 10 15

Fitted values
lm(Alt ~ EBIT + NWC + Sale + equity + RE)

19
 Q−Q Residuals

422 486
6
Standardized residuals

614
4
2
0
−2

−3 −2 −1 0 1 2 3

Theoretical Quantiles
lm(Alt ~ EBIT + NWC + Sale + equity + RE)

shapiro.test(p_mod$residuals)

##
## Shapiro-Wilk normality test
##
## data: p_mod$residuals
## W = 0.81856, p-value < 2.2e-16

hist(p_mod$residuals)

20
 Histogram of p_mod$residuals
300
200
Frequency

50 100
0

−10 0 10 20 30 40 50

p_mod$residuals

21