library(readxl)
library(psych)
library(kableExtra)

base_pymes <- 
  read_excel("Base de datos depurada ELE IGS Análisis descriptivo.xlsx", 
  sheet = "BBDD ELE V")
## New names:
## * `` -> ...53
## * `Otras capacitaciones` -> `Otras capacitaciones...64`
## * `Otras capacitaciones` -> `Otras capacitaciones...70`
## * `` -> ...91
## * `` -> ...94
## * ...
#names(base_pymes)

####Ratio 1#########
####Tabla descriptiva 1##########
tabla_1_1<-
with(base_pymes, describe(`Activo C/ Pasivo C`))

kable(tabla_1_1, align = "c", 
      booktabs=TRUE, 
      linesep="", 
      caption = "tabla ratio 1")%>%
  kable_styling(font_size = 9, 
                full_width = FALSE, 
                latex_options = c("basic", "cambria", "hold", "hover"))%>%
  add_footnote("Fuente: Unidad de Análisis Institucional", 
               notation="alphabet")
tabla ratio 1
vars n mean sd median trimmed mad min max range skew kurtosis se
X1 1 2511 6.373253 77.25184 1.453381 1.794166 1.002968 0 3208.948 3208.948 33.12196 1255.896 1.541649
a Fuente: Unidad de Análisis Institucional
####Tabla descriptiva 2##########
tabla_1_2<-
with(base_pymes,describeBy(as.numeric(`Activo C/ Pasivo C`), 
           list(`Ïndice codificado`), mat = TRUE, digits = 2))

kable(tabla_1_2, align = "c", 
      booktabs=TRUE, 
      linesep="", 
      caption = "tabla ratio 1")%>%
  kable_styling(font_size = 9, 
                full_width = FALSE, 
                latex_options = c("basic", "cambria", "hold", "hover"))%>%
  add_footnote("Fuente: Unidad de Análisis Institucional", 
               notation="alphabet")
tabla ratio 1
item group1 vars n mean sd median trimmed mad min max range skew kurtosis se
X11 1 Alta probabilidad de quiebra 1 659 2.15 27.20 0.94 0.92 0.54 0 698.47 698.47 25.47 649.15 1.06
X12 2 Baja probabilidad de quiebra 1 1844 7.90 88.63 1.76 2.23 1.15 0 3208.95 3208.95 29.56 982.51 2.06
a Fuente: Unidad de Análisis Institucional
###Histograma 1#######
ratio_1<-
base_pymes$`Activo C/ Pasivo C`

ratio_1<-
subset(ratio_1, ratio_1<7)

hist(ratio_1,
     col = "rainbow"(20), breaks = 35)

####Ratio 2#########
####Tabla descriptiva 2##########
tabla_2_1<-
  with(base_pymes, describe(`Fondos propios/ PatN+ Pasivo`))

kable(tabla_2_1[,1:10], align = "c", 
      booktabs=TRUE, 
      linesep="", 
      caption = "tabla ratio 1")%>%
  kable_styling(font_size = 12, 
                full_width = FALSE, 
                latex_options = c("basic", "cambria", "hold", "hover"))%>%
  add_footnote("Fuente: Unidad de Análisis Institucional", 
               notation="alphabet")
tabla ratio 1
vars n mean sd median trimmed mad min max range
X1 1 2517 0.3831313 0.5853572 0.35044 0.3646841 0.283316 -2.062886 23.34344 25.40632
a Fuente: Unidad de Análisis Institucional
####Tabla descriptiva 2##########
tabla_2_2<-
  with(base_pymes,describeBy(as.numeric(`Fondos propios/ PatN+ Pasivo`), 
                             list(`Ïndice codificado`), mat = TRUE, digits = 2))

kable(tabla_2_2[,1:10], align = "c", 
      booktabs=TRUE, 
      linesep="", 
      caption = "tabla ratio 1")%>%
  kable_styling(font_size = 12, 
                full_width = FALSE, 
                latex_options = c("basic", "cambria", "hold", "hover"))%>%
  add_footnote("Fuente: Unidad de Análisis Institucional", 
               notation="alphabet")
tabla ratio 1
item group1 vars n mean sd median trimmed mad min
X11 1 Alta probabilidad de quiebra 1 659 0.18 0.30 0.18 0.19 0.17 -2.06
X12 2 Baja probabilidad de quiebra 1 1847 0.45 0.64 0.43 0.43 0.26 -1.06
a Fuente: Unidad de Análisis Institucional
###Histograma 2#######

ratio_2<-
base_pymes$`Fondos propios/ PatN+ Pasivo`

ratio_2<-
subset(ratio_2, ratio_2<1 & 
                ratio_2>0)

hist(ratio_2,
     col = "rainbow"(20), breaks = 35)