Proyecto análisis estadístico sobre la compensación por aprovechamiento de fauna silvestre en Colombia
María Paula Gómez Silva - Angie Maria Moreno Mantilla
Facultad de Ciencias Naturales e Ingeniería
San Gil, Santander - 2026
library(psych) # Estadísticos descriptivos avanzados
library(corrplot) # Visualización de matrices de correlación
library(readr) # Importación eficiente de datos
library(GGally) # Gráficos multivariados (ggpairs)
library(car) # Diagnóstico de modelos de regresión
library(MASS) # Selección de modelos (stepAIC)
library(ggplot2) # Visualización gráfica de datos
library(skimr) # Resumen más visual
library(dplyr) # Manipulación y transformación de datos1 Introducción al caso de estudio
A nivel nacional, se puede ejercer la caza de fauna silvestre nativa con diversos fines: de subsistencia, comerciales, de fomento, de control, deportiva y científica. Los usuarios que ejercen la caza deberán pagar la tasa compensatoria por caza de fauna silvestre nativa para coadyuvar a garantizar la renovabilidad de dicho recurso, reglamentada en el artículo 42 de la Ley 99 de 1993.
Nota: No se cobra a la caza de subsistencia. No aplica a las especies exóticas ni al recurso pesquero.
2 Carga de datos
Datos_Tasa_compensatoria_por_Caza_de_Fauna <- read_csv(
"Datos_de_la_Tasa_compensatoria_por_Caza_de_Fauna_Silvestre_20260509.csv",
col_types = cols(
Año = col_character(),
`FECHA DE OTORGAMIENTO` = col_character(),
`VIGENCIA (meses)` = col_number(),
`VALOR TIPO DE CAZA` = col_number(),
`COEFICIENTE BIÓTICO (Cb)` = col_number(),
`NACIONALIDAD (N)` = col_number(),
`GRUPO TRÓFICO (Gt)` = col_number(),
`COEFICIENTE DE VALORACIÓN \n(V)` = col_number(),
`FACTOR REGIONAL (FR)` = col_number(),
`No. ESPECÍMENES O MUESTRAS (Es)` = col_number(),
`TARIFA MÍNIMA (TM)` = col_number(),
`COSTO DE IMPLEMENTACIÓN (CI)` = col_number(),
`MONTO A PAGAR (MP)` = col_number(),
`MP DESPUÉS DE RECLAMACIÓN` = col_number(),
RECAUDO = col_number()
)
)3 Exploración inicial
3.1 Tipo de datos
## 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>3.3 Nombre de las variables
## [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"3.4 Resumen estadístico
## 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 :5654 Preparación de datos
# Asignamos nombre corto para facilitar el trabajo
datostasaCF <- Datos_Tasa_compensatoria_por_Caza_de_Fauna
summary(datostasaCF)## 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 :5655 Datos faltantes
## Año
## 0
## Autoridad Ambiental
## 0
## NOMBRE/RAZÓN SOCIAL
## 0
## PERSONA NATURAL O JURÍDICA
## 0
## REPRESENTANTE LEGAL
## 42
## TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA
## 0
## No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO
## 0
## FECHA DE OTORGAMIENTO
## 0
## VIGENCIA (meses)
## 0
## TIPO DE CAZA
## 0
## VALOR TIPO DE CAZA
## 0
## ESPECIE O GRUPO TAXONÓMICO
## 0
## ORDEN
## 4
## CLASE
## 4
## No. INDIVIDUOS y/o MUESTRAS APROBADOS
## 405
## DEPARTAMENTO
## 0
## MUNICIPIO
## 0
## No. CERTIFICADO SiB
## 0
## ESTADO DE CONSERVACIÓN DE LA ESPECIE
## 383
## ESTADO DE CONSERVACIÓN DEL HÁBITAT
## 352
## PRESIÓN POR USO
## 866
## COEFICIENTE BIÓTICO (Cb)
## 0
## NACIONALIDAD (N)
## 0
## GRUPO TRÓFICO (Gt)
## 0
## COEFICIENTE DE VALORACIÓN \n(V)
## 0
## FACTOR REGIONAL (FR)
## 0
## No. ESPECÍMENES O MUESTRAS (Es)
## 0
## UNIDAD DE Es
## 0
## TARIFA MÍNIMA (TM)
## 50
## COSTO DE IMPLEMENTACIÓN (CI)
## 50
## MONTO A PAGAR (MP)
## 385
## RECLAMACIONES
## 0
## MP DESPUÉS DE RECLAMACIÓN
## 1012
## RECAUDO
## 565# Porcentaje de NA por columna
porcentaje_na <- sapply(datostasaCF,function(x) mean(is.na(x))*100)
porcentaje_na## Año
## 0.0000000
## Autoridad Ambiental
## 0.0000000
## NOMBRE/RAZÓN SOCIAL
## 0.0000000
## PERSONA NATURAL O JURÍDICA
## 0.0000000
## REPRESENTANTE LEGAL
## 2.7522936
## TIPO DE AUTORIZACIÓN / ENTIDAD QUE OTORGA
## 0.0000000
## No. PERMISO/ LICENCIA AMBIENTAL/ SANCIONATORIO
## 0.0000000
## FECHA DE OTORGAMIENTO
## 0.0000000
## VIGENCIA (meses)
## 0.0000000
## TIPO DE CAZA
## 0.0000000
## VALOR TIPO DE CAZA
## 0.0000000
## ESPECIE O GRUPO TAXONÓMICO
## 0.0000000
## ORDEN
## 0.2621232
## CLASE
## 0.2621232
## No. INDIVIDUOS y/o MUESTRAS APROBADOS
## 26.5399738
## DEPARTAMENTO
## 0.0000000
## MUNICIPIO
## 0.0000000
## No. CERTIFICADO SiB
## 0.0000000
## ESTADO DE CONSERVACIÓN DE LA ESPECIE
## 25.0982962
## ESTADO DE CONSERVACIÓN DEL HÁBITAT
## 23.0668414
## PRESIÓN POR USO
## 56.7496723
## COEFICIENTE BIÓTICO (Cb)
## 0.0000000
## NACIONALIDAD (N)
## 0.0000000
## GRUPO TRÓFICO (Gt)
## 0.0000000
## COEFICIENTE DE VALORACIÓN \n(V)
## 0.0000000
## FACTOR REGIONAL (FR)
## 0.0000000
## No. ESPECÍMENES O MUESTRAS (Es)
## 0.0000000
## UNIDAD DE Es
## 0.0000000
## TARIFA MÍNIMA (TM)
## 3.2765400
## COSTO DE IMPLEMENTACIÓN (CI)
## 3.2765400
## MONTO A PAGAR (MP)
## 25.2293578
## RECLAMACIONES
## 0.0000000
## MP DESPUÉS DE RECLAMACIÓN
## 66.3171691
## RECAUDO
## 37.02490176 Selección de variables
tabla_analisis <- dplyr::select(datostasaCF,
"MONTO A PAGAR (MP)", # Variable depediente
"FACTOR REGIONAL (FR)",
"COEFICIENTE BIÓTICO (Cb)",
"GRUPO TRÓFICO (Gt)",
"COEFICIENTE DE VALORACIÓN \n(V)",
"No. ESPECÍMENES O MUESTRAS (Es)"
)
#Renombrar columnas
colnames(tabla_analisis) <- c("MP", "FR", "Cb", "Gt", "V", "Es")
#Vista previa
head(tabla_analisis)6.1 Resumen de variables seleccionadas
## MP FR Cb Gt
## Min. : 966 Min. : 0.0000 Min. :0.000 Min. :0.0000
## 1st Qu.: 5433 1st Qu.: 0.4800 1st Qu.:1.000 1st Qu.:0.8000
## Median : 11278 Median : 0.4800 Median :1.000 Median :0.8000
## Mean : 59518 Mean : 0.6805 Mean :1.209 Mean :0.7282
## 3rd Qu.: 40577 3rd Qu.: 0.4800 3rd Qu.:1.000 3rd Qu.:0.8000
## Max. :7334339 Max. :76.8000 Max. :5.000 Max. :1.0000
## NA's :385
## V Es
## Min. : 0.0000 Min. : 0.00
## 1st Qu.: 1.0000 1st Qu.: 1.00
## Median : 1.0000 Median : 2.00
## Mean : 0.9784 Mean : 14.92
## 3rd Qu.: 1.0000 3rd Qu.: 5.00
## Max. :20.0000 Max. :7261.00
## ## [1] 1526 67 Visualización de filas vacías
8 Análisis de correlación
8.1 Matriz de correlaciones
## MP FR Cb Gt V Es
## MP 1.00 0.18 0.10 -0.02 0.17 0.85
## FR 0.18 1.00 0.41 0.11 0.93 -0.02
## Cb 0.10 0.41 1.00 0.14 0.08 -0.02
## Gt -0.02 0.11 0.14 1.00 0.01 -0.12
## V 0.17 0.93 0.08 0.01 1.00 0.00
## Es 0.85 -0.02 -0.02 -0.12 0.00 1.00#Visualización de correlaciones
corrplot(matriz_cor,method = "color", addCoef.col = "black",main= "Matriz de correlaciones(Taza caza)", mar = c(0, 0, 1, 0))
