Algorithmic Correlations.R

Sergio Mora — Feb 11, 2018, 10:06 PM

library(readxl)
library(ggplot2)
library(dplyr)

Attaching package: 'dplyr'

The following objects are masked from 'package:stats':

    filter, lag

The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union

library(writexl)
Acciones <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                       sheet = "COLCAP")
View(Acciones)
BCOLOMBIA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                     sheet = "BCOLOMBIA")
BOGOTA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                       sheet = "BOGOTA")
CELSIA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "CELSIA")
CEMARGOS <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "CEMARGOS")
CLH <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "CLH")
CNEC <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "CNEC")
CONCONCRET <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "CONCONCRET")
ECOPETROL <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "ECOPETROL")
EEB <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "EEB")
EXITO <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "EXITO")
GRUPOARGOS <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "GRUPOARGOS")
GRUPOAVAL <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "GRUPOAVAL")
GRUPOSURA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "GRUPOSURA")
ISA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "ISA")
NUTRESA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "NUTRESA")
PFAVAL <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFAVAL")
PFAVH <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFAVH")
PFBCOLOM <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFBCOLOMBIA")
PFCEMARGOS <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFCEMARGOS")
PFDAVVNDA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFDAVVNDA")
PFGRUPOARG <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFGRUPOARG")
PFGRUPSURA <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                        sheet = "PFGRUPSURA")
PROMIGAS <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                         sheet = "PROMIGAS")
CORFICOLCF <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                       sheet = "CORFICOLCF")
ETB <- read_excel("C:/Users/Sergio Mora/OneDrive/Estudios/Mercado de Valores 002/ACCIONES/e-bvc/Acciones.xlsx", 
                       sheet = "ETB")

BCOLOMBIA$Accion <- c("BCOLOMBIA")
BOGOTA$Accion <-  c("BOGOTA")
CELSIA$Accion <-  c("CELSIA")
CEMARGOS$Accion <-  c("CEMARGOS")
CLH$Accion <-  c("CLH")
CNEC$Accion <-  c("CNEC")
CONCONCRET$Accion <-  c("CONCONCRET")
CORFICOLCF$Accion <-  c("CORFICOLCF")
ECOPETROL$Accion <-  c("ECOPETROL")
EEB$Accion <-  c("EEB")
ETB$Accion <-  c("ETB")
EXITO$Accion <-  c("EXITO")
GRUPOARGOS$Accion <-  c("GRUPOARGOS")
GRUPOAVAL$Accion <-  c("GRUPOAVAL")
GRUPOSURA$Accion <-  c("GRUPOSURA")
ISA$Accion <-  c("ISA")
NUTRESA$Accion <-  c("NUTRESA")
PFAVAL$Accion <-  c("PFAVAL")
PFAVH$Accion <-  c("PFAVH")
PFBCOLOM$Accion <-  c("PFBCOLOM")
PFCEMARGOS$Accion <-  c("PFCEMARGOS")
PFDAVVNDA$Accion <-  c("PFDAVVNDA")
PFGRUPOARG$Accion <-  c("PFGRUPOARG")
PFGRUPSURA$Accion <-  c("PFGRUPSURA")
PROMIGAS$Accion <-  c("PROMIGAS")



COLCAP <- data.frame(rbind(BCOLOMBIA, 
                           BOGOTA,
                           CELSIA,
                           CEMARGOS,
                           CLH,
                           CNEC,
                           CONCONCRET,
                           CORFICOLCF,
                           ECOPETROL,
                           EEB,
                           ETB,
                           EXITO,
                           GRUPOARGOS,
                           GRUPOAVAL,
                           GRUPOSURA,
                           ISA,
                           NUTRESA,
                           PFAVAL,
                           PFAVH,
                           PFBCOLOM,
                           PFCEMARGOS,
                           PFDAVVNDA,
                           PFGRUPSURA,
                           PROMIGAS
                           )
                     )

qplot(COLCAP$Fecha,
      COLCAP$Apertura,
    data = COLCAP, color=COLCAP$Accion)+geom_smooth()

`geom_smooth()` using method = 'loess'

NROW(BCOLOMBIA)

[1] 242

NROW(BOGOTA) # 241 rows ok

[1] 242

NROW(CELSIA)

[1] 242

NROW(CEMARGOS)

[1] 242

NROW(CLH)

[1] 242

NROW(CNEC)

[1] 242

NROW(CONCONCRET) # 236 rows

[1] 236

NROW(CORFICOLCF)

[1] 242

NROW(ECOPETROL)

[1] 242

NROW(EEB)

[1] 242

NROW(ETB)

[1] 242

NROW(EXITO)

[1] 242

NROW(GRUPOARGOS)

[1] 242

NROW(GRUPOAVAL) # 240 rows

[1] 240

NROW(GRUPOSURA)

[1] 242

NROW(ISA)

[1] 242

NROW(NUTRESA)

[1] 242

NROW(PFAVAL)

[1] 242

NROW(PFAVH)

[1] 242

NROW(PFBCOLOM)

[1] 242

NROW(PFCEMARGOS) # 241 rows

[1] 241

NROW(PFDAVVNDA)

[1] 242

NROW(PFGRUPSURA)

[1] 242

NROW(PROMIGAS) # 169 rows

[1] 169

Open_Price <- data.frame(BCOLOMBIA$Apertura, 
    BOGOTA$Apertura,
    CELSIA$Apertura,
    CEMARGOS$Apertura,
    CLH$Apertura,
    CNEC$Apertura,
    CORFICOLCF$Apertura,
    ECOPETROL$Apertura,
    EEB$Apertura,
    ETB$Apertura,
    EXITO$Apertura,
    GRUPOARGOS$Apertura,
    GRUPOSURA$Apertura,
    ISA$Apertura,
    NUTRESA$Apertura,
    PFAVAL$Apertura,
    PFAVH$Apertura,
    PFBCOLOM$Apertura,
    PFDAVVNDA$Apertura,
    PFGRUPSURA$Apertura
    )

Correlaciones <- data.frame(cor(Open_Price))
head(Correlaciones)

                   BCOLOMBIA.Apertura BOGOTA.Apertura CELSIA.Apertura
BCOLOMBIA.Apertura         1.00000000       0.7891692       0.7735117
BOGOTA.Apertura            0.78916918       1.0000000       0.8386523
CELSIA.Apertura            0.77351168       0.8386523       1.0000000
CEMARGOS.Apertura          0.05677931      -0.3615567      -0.2518413
CLH.Apertura               0.64979303       0.3042414       0.3275128
CNEC.Apertura              0.68520757       0.8119471       0.7143185
                   CEMARGOS.Apertura CLH.Apertura CNEC.Apertura
BCOLOMBIA.Apertura        0.05677931    0.6497930     0.6852076
BOGOTA.Apertura          -0.36155674    0.3042414     0.8119471
CELSIA.Apertura          -0.25184126    0.3275128     0.7143185
CEMARGOS.Apertura         1.00000000    0.2265111    -0.3458689
CLH.Apertura              0.22651106    1.0000000     0.2006345
CNEC.Apertura            -0.34586889    0.2006345     1.0000000
                   CORFICOLCF.Apertura ECOPETROL.Apertura EEB.Apertura
BCOLOMBIA.Apertura         -0.11647491          0.2040670    0.6826023
BOGOTA.Apertura             0.16860399          0.4042231    0.8533371
CELSIA.Apertura             0.01482925          0.4173286    0.8132259
CEMARGOS.Apertura          -0.10814527         -0.4460179   -0.5520147
CLH.Apertura               -0.05490977         -0.1778202    0.3370601
CNEC.Apertura               0.04024952          0.4684309    0.7836001
                   ETB.Apertura EXITO.Apertura GRUPOARGOS.Apertura
BCOLOMBIA.Apertura  -0.54200160      0.1220477           0.7946032
BOGOTA.Apertura     -0.84162573      0.3496684           0.5723886
CELSIA.Apertura     -0.71436730      0.3575046           0.6564179
CEMARGOS.Apertura    0.59280989     -0.5181229           0.2691384
CLH.Apertura        -0.06807208     -0.1025821           0.4379229
CNEC.Apertura       -0.74343679      0.3167672           0.5223298
                   GRUPOSURA.Apertura ISA.Apertura NUTRESA.Apertura
BCOLOMBIA.Apertura          0.8208942    0.8208615        0.7757505
BOGOTA.Apertura             0.7532662    0.8405057        0.9282651
CELSIA.Apertura             0.6477262    0.8481945        0.8823657
CEMARGOS.Apertura           0.1454458   -0.3031854       -0.4368978
CLH.Apertura                0.3689217    0.3948379        0.3391512
CNEC.Apertura               0.6346022    0.7825066        0.8137809
                   PFAVAL.Apertura PFAVH.Apertura PFBCOLOM.Apertura
BCOLOMBIA.Apertura       0.8555142      0.1858221         0.9331207
BOGOTA.Apertura          0.9452896      0.4848654         0.5887050
CELSIA.Apertura          0.7912607      0.3656806         0.6054434
CEMARGOS.Apertura       -0.1991119     -0.4756692         0.3059038
CLH.Apertura             0.4188318     -0.1612028         0.7187572
CNEC.Apertura            0.8187897      0.3633069         0.5199950
                   PFDAVVNDA.Apertura PFGRUPSURA.Apertura
BCOLOMBIA.Apertura          0.7502882           0.7725101
BOGOTA.Apertura             0.3984027           0.6739623
CELSIA.Apertura             0.4232729           0.5690196
CEMARGOS.Apertura           0.4457623           0.2402639
CLH.Apertura                0.6556941           0.3880688
CNEC.Apertura               0.3470615           0.5612687

write_xlsx(Correlaciones, 
           path = "Correlaciones.xlsx", 
           col_names = TRUE)

[1] "C:\\Users\\Sergio Mora\\OneDrive\\Estudios\\Mercado de Valores 002\\ACCIONES\\e-bvc\\Correlaciones.xlsx"