Correlation Analysis for Total motor Vehicle Sales

Does exchange rate influence motor vehicle sales?

Vijeta Tulsiyan - s3398979

Last updated: 21 October, 2017

Introduction

Total motor vehicle sales in Australia for the month of August 2017 alone was 100,809 which increased by 2.7% compared to last year. The sale of Vehicles, both passenger cars and sports utility, was robust in last few years.

Despite news in media that local automobile companies are struggling to cope competitive prices from import of cheaper car models and looking for government support to sustain, wondering what is giving momentum to total sales?

Introduction Cont.

Problem Statement

Data

Data Cont.

Descriptive Statistics and Visualisation

#Check for any Linear relationship
plot(Vehicle_Sales ~ AUD_USD, data = Vehicle)

plot(log(Vehicle_Sales) ~ log(AUD_USD), data = Vehicle)

Data Visualisation continued…

par(mfrow=c(2,2))
Vehicle$AUD_USD %>%  hist(main = "AUD_USD")
log(Vehicle$AUD_USD) %>%  hist(main = "log(AUD_USD)")
Vehicle$Vehicle_Sales %>%  hist(main = "Vehicle_Sales")
log(Vehicle$Vehicle_Sales) %>%  hist(main = "log(Vehicle_Sales)")

Data Visualisation continued…

#Linear Model FIT

model <- lm(Vehicle_Sales ~ AUD_USD, data = Vehicle)
model %>% summary()
## 
## Call:
## lm(formula = Vehicle_Sales ~ AUD_USD, data = Vehicle)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -11920  -5350  -1457   2611  17201 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 52619.94    2875.59   18.30   <2e-16 ***
## AUD_USD       401.55      35.33   11.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7646 on 198 degrees of freedom
## Multiple R-squared:  0.3949, Adjusted R-squared:  0.3918 
## F-statistic: 129.2 on 1 and 198 DF,  p-value: < 2.2e-16

Data Visualisation continued…

plot(Vehicle_Sales ~ AUD_USD, data = Vehicle, xlab = "AUD_USD", ylab = "Total Vehicle Sales in Australia")
abline(model, col = "red")

Hypothesis Testing

#Hypothesis Testing
library(Hmisc)
bivariate<-as.matrix(dplyr::select(Vehicle, Vehicle_Sales,AUD_USD)) #Create a matrix of the variables to be correlated
rcorr(bivariate, type = "pearson")
##               Vehicle_Sales AUD_USD
## Vehicle_Sales          1.00    0.63
## AUD_USD                0.63    1.00
## 
## n= 200 
## 
## 
## P
##               Vehicle_Sales AUD_USD
## Vehicle_Sales                0     
## AUD_USD        0
# COnfidence Interval
library(psychometric)
r=cor(Vehicle$Vehicle_Sales, Vehicle$AUD_USD)
CIr(r = r, n = 200, level = .95)
## [1] 0.5364271 0.7056214
#Correlation
cor(Vehicle$Vehicle_Sales, Vehicle$AUD_USD)
## [1] 0.6283998

Discussion

References

http://dfat.gov.au

https://www.aph.gov.au/About_Parliament/Parliamentary_Departments/Parliamentary_Library/pubs/MSB

https://www.eia.gov/dnav/pet/hist/LeafHandler.ashx?n=pet&s=rwtc&f=m