R Markdown

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

summary(cars)
##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00
load("ProyectEstadistica.RData")

library(dplyr)     # Manipulación y transformación de datos
## 
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(ggplot2)   # Visualización gráfica de datos
library(psych)     # Estadísticos descriptivos avanzados
## 
## Adjuntando el paquete: 'psych'
## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha
library(corrplot)  # Visualización de matrices de correlación
## corrplot 0.95 loaded
library(readr)     # Importación eficiente de datos
library(GGally)    # Gráficos multivariados (ggpairs)
library(car)       # Diagnóstico de modelos de regresión
## Cargando paquete requerido: carData
## 
## Adjuntando el paquete: 'car'
## The following object is masked from 'package:psych':
## 
##     logit
## The following object is masked from 'package:dplyr':
## 
##     recode
library(MASS)      # Selección de modelos (stepAIC)
## 
## Adjuntando el paquete: 'MASS'
## The following object is masked from 'package:dplyr':
## 
##     select
head(Datos_Tasa_compensatoria_por_Caza_de_Fauna)
## # A tibble: 6 × 34
##   Año   `Autoridad Ambiental` `NOMBRE/RAZÓN SOCIAL` `PERSONA NATURAL O JURÍDICA`
##   <chr> <chr>                 <chr>                 <chr>                       
## 1 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## 2 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## 3 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## 4 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## 5 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## 6 2,021 CORTOLIMA             UNVERSIDAD DEL TOLIMA Persona jurídica            
## # ℹ 30 more variables: `REPRESENTANTE LEGAL` <chr>,
## #   `TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA` <chr>,
## #   `No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO` <chr>,
## #   `FECHA DE OTORGAMIENTO` <chr>, `VIGENCIA (meses)` <dbl>,
## #   `TIPO DE CAZA` <chr>, `VALOR TIPO DE CAZA` <dbl>,
## #   `ESPECIE O GRUPO TAXONÓMICO` <chr>, ORDEN <chr>, CLASE <chr>,
## #   `No. INDIVIDUOS y/o MUESTRAS APROBADOS` <chr>, DEPARTAMENTO <chr>, …
View(Datos_Tasa_compensatoria_por_Caza_de_Fauna)

#Exploración inicial
str(Datos_Tasa_compensatoria_por_Caza_de_Fauna)
## spc_tbl_ [1,526 × 34] (S3: spec_tbl_df/tbl_df/tbl/data.frame)
##  $ Año                                           : chr [1:1526] "2,021" "2,021" "2,021" "2,021" ...
##  $ Autoridad Ambiental                           : chr [1:1526] "CORTOLIMA" "CORTOLIMA" "CORTOLIMA" "CORTOLIMA" ...
##  $ NOMBRE/RAZÓN SOCIAL                           : chr [1:1526] "UNVERSIDAD DEL TOLIMA" "UNVERSIDAD DEL TOLIMA" "UNVERSIDAD DEL TOLIMA" "UNVERSIDAD DEL TOLIMA" ...
##  $ PERSONA NATURAL O JURÍDICA                    : chr [1:1526] "Persona jurídica" "Persona jurídica" "Persona jurídica" "Persona jurídica" ...
##  $ REPRESENTANTE LEGAL                           : chr [1:1526] "Omar A. Mejía Patiño" "Omar A. Mejía Patiño" "Omar A. Mejía Patiño" "Omar A. Mejía Patiño" ...
##  $ TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA     : chr [1:1526] "Permiso / Autoridad Ambiental diferente a ANLA" "Permiso / Autoridad Ambiental diferente a ANLA" "Permiso / Autoridad Ambiental diferente a ANLA" "Permiso / Autoridad Ambiental diferente a ANLA" ...
##  $ No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO: chr [1:1526] "Resolución 3758 de 2016" "Resolución 3758 de 2016" "Resolución 3758 de 2016" "Resolución 3758 de 2016" ...
##  $ FECHA DE OTORGAMIENTO                         : chr [1:1526] "2016-11-16T00:00:00.000" "2016-11-16T00:00:00.000" "2016-11-16T00:00:00.000" "2016-11-16T00:00:00.000" ...
##  $ VIGENCIA (meses)                              : num [1:1526] 120 120 120 120 120 120 120 6 120 120 ...
##  $ TIPO DE CAZA                                  : chr [1:1526] "Científica estudios ambientales" "Científica estudios ambientales" "Científica estudios ambientales" "Científica estudios ambientales" ...
##  $ VALOR TIPO DE CAZA                            : num [1:1526] 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 ...
##  $ ESPECIE O GRUPO TAXONÓMICO                    : chr [1:1526] "Panthiades bitias" "Actinote anteas" "Altinote dicaeus" "Catasticta tomyris tomyris" ...
##  $ ORDEN                                         : chr [1:1526] "Lepidoptera" "Lepidoptera" "Lepidoptera" "Lepidoptera" ...
##  $ CLASE                                         : chr [1:1526] "Insecta" "Insecta" "Insecta" "Insecta" ...
##  $ No. INDIVIDUOS y/o MUESTRAS APROBADOS         : chr [1:1526] "No registra" "No registra" "No registra" "No registra" ...
##  $ DEPARTAMENTO                                  : chr [1:1526] "Tolima" "Tolima" "Tolima" "Tolima" ...
##  $ MUNICIPIO                                     : chr [1:1526] "Dolores, Natagaima, Prado, Purificación, Cunday,Cunday, Venadillo, Alvarado, Piedras, Ortega, Espinal, Guamo, A"| __truncated__ "Dolores, Natagaima, Prado, Purificación, Cunday,Cunday, Venadillo, Alvarado, Piedras, Ortega, Espinal, Guamo, A"| __truncated__ "Dolores, Natagaima, Prado, Purificación, Cunday,Cunday, Venadillo, Alvarado, Piedras, Ortega, Espinal, Guamo, A"| __truncated__ "Dolores, Natagaima, Prado, Purificación, Cunday,Cunday, Venadillo, Alvarado, Piedras, Ortega, Espinal, Guamo, A"| __truncated__ ...
##  $ No. CERTIFICADO SiB                           : chr [1:1526] "1634DA1EAAC -16341C01996 -6321985BD7- 164F610CBDC- 16DBB9126D7- 16321A6F91C- 16321C92741- 16364C22A54- 1632181A"| __truncated__ "1634DA1EAAC -16341C01996 -6321985BD7- 164F610CBDC- 16DBB9126D7- 16321A6F91C- 16321C92741- 16364C22A54- 1632181A"| __truncated__ "1634DA1EAAC -16341C01996 -6321985BD7- 164F610CBDC- 16DBB9126D7- 16321A6F91C- 16321C92741- 16364C22A54- 1632181A"| __truncated__ "1634DA1EAAC -16341C01996 -6321985BD7- 164F610CBDC- 16DBB9126D7- 16321A6F91C- 16321C92741- 16364C22A54- 1632181A"| __truncated__ ...
##  $ ESTADO DE CONSERVACIÓN DE LA ESPECIE          : chr [1:1526] "Datos insuficientes (DD) o especies No evaluadas (NE)" "Datos insuficientes (DD) o especies No evaluadas (NE)" "Datos insuficientes (DD) o especies No evaluadas (NE)" "Datos insuficientes (DD) o especies No evaluadas (NE)" ...
##  $ ESTADO DE CONSERVACIÓN DEL HÁBITAT            : chr [1:1526] "Pobremente conservado" "Pobremente conservado" "Pobremente conservado" "Pobremente conservado" ...
##  $ PRESIÓN POR USO                               : chr [1:1526] NA NA NA NA ...
##  $ COEFICIENTE BIÓTICO (Cb)                      : num [1:1526] 1 1 1 1 1 1 1 1 1 1 ...
##  $ NACIONALIDAD (N)                              : num [1:1526] 0 0 0 0 0 0 0 0 0 0 ...
##  $ GRUPO TRÓFICO (Gt)                            : num [1:1526] 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.8 0.15 0.15 ...
##  $ COEFICIENTE DE VALORACIÓN 
## (V)               : num [1:1526] 1 1 1 1 1 1 1 1 1 1 ...
##  $ FACTOR REGIONAL (FR)                          : num [1:1526] 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.48 0.09 0.09 ...
##  $ No. ESPECÍMENES O MUESTRAS (Es)               : num [1:1526] 1 1 1 1 1 1 1 1 1 1 ...
##  $ UNIDAD DE Es                                  : chr [1:1526] "Individuo" "Individuo" "Individuo" "Individuo" ...
##  $ TARIFA MÍNIMA (TM)                            : num [1:1526] 10903 10903 10903 10903 10903 ...
##  $ COSTO DE IMPLEMENTACIÓN (CI)                  : num [1:1526] 29530 29530 29530 29530 29530 ...
##  $ MONTO A PAGAR (MP)                            : num [1:1526] 981 981 981 981 981 ...
##  $ RECLAMACIONES                                 : chr [1:1526] "No" "No" "No" "No" ...
##  $ MP DESPUÉS DE RECLAMACIÓN                     : num [1:1526] NA NA NA NA NA NA NA NA NA NA ...
##  $ RECAUDO                                       : num [1:1526] 981 981 981 981 981 ...
##  - attr(*, "spec")=
##   .. cols(
##   ..   Año = col_character(),
##   ..   `Autoridad Ambiental` = col_character(),
##   ..   `NOMBRE/RAZÓN SOCIAL` = col_character(),
##   ..   `PERSONA NATURAL O JURÍDICA` = col_character(),
##   ..   `REPRESENTANTE LEGAL` = col_character(),
##   ..   `TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA` = col_character(),
##   ..   `No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO` = col_character(),
##   ..   `FECHA DE OTORGAMIENTO` = col_character(),
##   ..   `VIGENCIA (meses)` = col_number(),
##   ..   `TIPO DE CAZA` = col_character(),
##   ..   `VALOR TIPO DE CAZA` = col_number(),
##   ..   `ESPECIE O GRUPO TAXONÓMICO` = col_character(),
##   ..   ORDEN = col_character(),
##   ..   CLASE = col_character(),
##   ..   `No. INDIVIDUOS y/o MUESTRAS APROBADOS` = col_character(),
##   ..   DEPARTAMENTO = col_character(),
##   ..   MUNICIPIO = col_character(),
##   ..   `No. CERTIFICADO SiB` = col_character(),
##   ..   `ESTADO DE CONSERVACIÓN DE LA ESPECIE` = col_character(),
##   ..   `ESTADO DE CONSERVACIÓN DEL HÁBITAT` = col_character(),
##   ..   `PRESIÓN POR USO` = col_character(),
##   ..   `COEFICIENTE BIÓTICO (Cb)` = col_number(),
##   ..   `NACIONALIDAD (N)` = col_number(),
##   ..   `GRUPO TRÓFICO (Gt)` = col_number(),
##   ..   `COEFICIENTE DE VALORACIÓN 
##   .. (V)` = col_number(),
##   ..   `FACTOR REGIONAL (FR)` = col_number(),
##   ..   `No. ESPECÍMENES O MUESTRAS (Es)` = col_number(),
##   ..   `UNIDAD DE Es` = col_character(),
##   ..   `TARIFA MÍNIMA (TM)` = col_number(),
##   ..   `COSTO DE IMPLEMENTACIÓN (CI)` = col_number(),
##   ..   `MONTO A PAGAR (MP)` = col_number(),
##   ..   RECLAMACIONES = col_character(),
##   ..   `MP DESPUÉS DE RECLAMACIÓN` = col_number(),
##   ..   RECAUDO = col_number()
##   .. )
##  - attr(*, "problems")=<externalptr>
dim(Datos_Tasa_compensatoria_por_Caza_de_Fauna)
## [1] 1526   34
names(Datos_Tasa_compensatoria_por_Caza_de_Fauna)
##  [1] "Año"                                           
##  [2] "Autoridad Ambiental"                           
##  [3] "NOMBRE/RAZÓN SOCIAL"                           
##  [4] "PERSONA NATURAL O JURÍDICA"                    
##  [5] "REPRESENTANTE LEGAL"                           
##  [6] "TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA"     
##  [7] "No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO"
##  [8] "FECHA DE OTORGAMIENTO"                         
##  [9] "VIGENCIA (meses)"                              
## [10] "TIPO DE CAZA"                                  
## [11] "VALOR TIPO DE CAZA"                            
## [12] "ESPECIE O GRUPO TAXONÓMICO"                    
## [13] "ORDEN"                                         
## [14] "CLASE"                                         
## [15] "No. INDIVIDUOS y/o MUESTRAS APROBADOS"         
## [16] "DEPARTAMENTO"                                  
## [17] "MUNICIPIO"                                     
## [18] "No. CERTIFICADO SiB"                           
## [19] "ESTADO DE CONSERVACIÓN DE LA ESPECIE"          
## [20] "ESTADO DE CONSERVACIÓN DEL HÁBITAT"            
## [21] "PRESIÓN POR USO"                               
## [22] "COEFICIENTE BIÓTICO (Cb)"                      
## [23] "NACIONALIDAD (N)"                              
## [24] "GRUPO TRÓFICO (Gt)"                            
## [25] "COEFICIENTE DE VALORACIÓN \n(V)"               
## [26] "FACTOR REGIONAL (FR)"                          
## [27] "No. ESPECÍMENES O MUESTRAS (Es)"               
## [28] "UNIDAD DE Es"                                  
## [29] "TARIFA MÍNIMA (TM)"                            
## [30] "COSTO DE IMPLEMENTACIÓN (CI)"                  
## [31] "MONTO A PAGAR (MP)"                            
## [32] "RECLAMACIONES"                                 
## [33] "MP DESPUÉS DE RECLAMACIÓN"                     
## [34] "RECAUDO"
summary(Datos_Tasa_compensatoria_por_Caza_de_Fauna)
##      Año            Autoridad Ambiental NOMBRE/RAZÓN SOCIAL
##  Length:1526        Length:1526         Length:1526        
##  Class :character   Class :character    Class :character   
##  Mode  :character   Mode  :character    Mode  :character   
##                                                            
##                                                            
##                                                            
##                                                            
##  PERSONA NATURAL O JURÍDICA REPRESENTANTE LEGAL
##  Length:1526                Length:1526        
##  Class :character           Class :character   
##  Mode  :character           Mode  :character   
##                                                
##                                                
##                                                
##                                                
##  TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA
##  Length:1526                              
##  Class :character                         
##  Mode  :character                         
##                                           
##                                           
##                                           
##                                           
##  No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO FECHA DE OTORGAMIENTO
##  Length:1526                                    Length:1526          
##  Class :character                               Class :character     
##  Mode  :character                               Mode  :character     
##                                                                      
##                                                                      
##                                                                      
##                                                                      
##  VIGENCIA (meses) TIPO DE CAZA       VALOR TIPO DE CAZA
##  Min.   :  0.00   Length:1526        Min.   :0.1000    
##  1st Qu.:  6.00   Class :character   1st Qu.:0.6000    
##  Median : 24.00   Mode  :character   Median :0.6000    
##  Mean   : 36.48                      Mean   :0.5965    
##  3rd Qu.: 24.00                      3rd Qu.:0.6000    
##  Max.   :120.00                      Max.   :1.2000    
##                                                        
##  ESPECIE O GRUPO TAXONÓMICO    ORDEN              CLASE          
##  Length:1526                Length:1526        Length:1526       
##  Class :character           Class :character   Class :character  
##  Mode  :character           Mode  :character   Mode  :character  
##                                                                  
##                                                                  
##                                                                  
##                                                                  
##  No. INDIVIDUOS y/o MUESTRAS APROBADOS DEPARTAMENTO        MUNICIPIO        
##  Length:1526                           Length:1526        Length:1526       
##  Class :character                      Class :character   Class :character  
##  Mode  :character                      Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##  No. CERTIFICADO SiB ESTADO DE CONSERVACIÓN DE LA ESPECIE
##  Length:1526         Length:1526                         
##  Class :character    Class :character                    
##  Mode  :character    Mode  :character                    
##                                                          
##                                                          
##                                                          
##                                                          
##  ESTADO DE CONSERVACIÓN DEL HÁBITAT PRESIÓN POR USO    COEFICIENTE BIÓTICO (Cb)
##  Length:1526                        Length:1526        Min.   :0.000           
##  Class :character                   Class :character   1st Qu.:1.000           
##  Mode  :character                   Mode  :character   Median :1.000           
##                                                        Mean   :1.209           
##                                                        3rd Qu.:1.000           
##                                                        Max.   :5.000           
##                                                                                
##  NACIONALIDAD (N) GRUPO TRÓFICO (Gt) COEFICIENTE DE VALORACIÓN \n(V)
##  Min.   :0        Min.   :0.0000     Min.   : 0.0000                
##  1st Qu.:0        1st Qu.:0.8000     1st Qu.: 1.0000                
##  Median :0        Median :0.8000     Median : 1.0000                
##  Mean   :0        Mean   :0.7282     Mean   : 0.9784                
##  3rd Qu.:0        3rd Qu.:0.8000     3rd Qu.: 1.0000                
##  Max.   :0        Max.   :1.0000     Max.   :20.0000                
##                                                                     
##  FACTOR REGIONAL (FR) No. ESPECÍMENES O MUESTRAS (Es) UNIDAD DE Es      
##  Min.   : 0.0000      Min.   :   0.00                 Length:1526       
##  1st Qu.: 0.4800      1st Qu.:   1.00                 Class :character  
##  Median : 0.4800      Median :   2.00                 Mode  :character  
##  Mean   : 0.6805      Mean   :  14.92                                   
##  3rd Qu.: 0.4800      3rd Qu.:   5.00                                   
##  Max.   :76.8000      Max.   :7261.00                                   
##                                                                         
##  TARIFA MÍNIMA (TM) COSTO DE IMPLEMENTACIÓN (CI) MONTO A PAGAR (MP)
##  Min.   : 9600      Min.   :26000                Min.   :    966   
##  1st Qu.:10567      1st Qu.:28620                1st Qu.:   5433   
##  Median :10903      Median :29530                Median :  11278   
##  Mean   :10961      Mean   :29686                Mean   :  59518   
##  3rd Qu.:11318      3rd Qu.:30652                3rd Qu.:  40577   
##  Max.   :11748      Max.   :31817                Max.   :7334339   
##  NA's   :50         NA's   :50                   NA's   :385       
##  RECLAMACIONES      MP DESPUÉS DE RECLAMACIÓN    RECAUDO       
##  Length:1526        Min.   :      0           Min.   :    981  
##  Class :character   1st Qu.:   6791           1st Qu.:   5433  
##  Mode  :character   Median :  16298           Median :  10865  
##                     Mean   :  41668           Mean   :  37914  
##                     3rd Qu.:  38028           3rd Qu.:  31817  
##                     Max.   :1249507           Max.   :1962540  
##                     NA's   :1012              NA's   :565

Including Plots

You can also embed plots, for example:

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.