Apuntes 96, 193-226
ISSN: 0252-1865
eISSN: 2223-1757
doi: 10.21678/apuntes.96.1841
© Creative Commons Attribution 3.0 Recibido el 15 de noviembre de 2022 Aceptado el 24 de mayo de 2023
Juan Walter Tudela-Mamani
Universidad Nacional del Altiplano-Facultad de Ingeniería Económica
Elias Cahui-Cahui
Universidad Nacional del Altiplano-Facultad de Ingeniería Económica
Javier Humpiri-Yucra
Universidad Nacional del Altiplano-Facultad de Ingeniería Económica
Resumen. El objetivo de la investigación es analizar la incidencia de la inversión pública sobre la recaudación del impuesto predial en gobiernos locales provinciales de la macrorregión sur del Perú durante el período 2010-2020. El tipo de estudio es aplicado, bajo el enfoque cuantitativo, diseño no experimental, tipo longitudinal de panel, de carácter descriptivo y explicativo. Los datos para cada municipalidad provincial fueron extraídos del portal informativo del Ministerio de Economía y Finanzas (MEF) segmentadas mediante el análisis de conglomerados. El método de estimación econométrica toma en cuenta la metodología para datos longitudinales de panel. Teniendo en cuenta el análisis de conglomerados y las estimaciones econométricas, los resultados muestran que un incremento del 1% en la inversión pública en transportes, saneamiento y cultura y deporte, provocó un aumento en la recaudación del impuesto predial de aproximadamente 0,3182%, 0,1208% y 0,0923% respectivamente (primer clúster); un incremento del 1% en la inversión pública en transportes, cultura y deporte, provocó un aumento en la recaudación del impuesto predial de aproximadamente 0,2588% y 0,1434% respectivamente (segundo clúster); y un incremento del 1% en la inversión pública en transportes y saneamiento, provocó un aumento en la recaudación del impuesto predial de aproximadamente 0,2923% y 0,1377% respectivamente (tercer clúster]).
Palabras clave. Impuesto predial, inversión pública, gobiernos locales, análisis de conglomerados, metodología para datos longitudinales de panel.
Public Investment and Property Tax Collection, an Application of Panel Data by Clusters in Provincial Local Governments of the Southern Macro-region of Peru, 2010-2020
Abstract. The objective of the research is to analyze the incidence of public investment on property tax collection in provincial local governments of the southern macro region of Peru during the period 2010 to 2020. The type of study is applied, under the quantitative approach, design non-expe-rimental, longitudinal panel type, descriptive and explanatory. The data for each provincial municipality were extracted from the information portal of the Ministry of Economy and Finance (MEF) segmented through cluster analysis. The econometric estimation method takes into account the metho-dology for longitudinal panel data. Taking into account the cluster analysis and econometric estimates, the results show that a 1% increase in public investment in transportation, sanitation, and culture and sports, caused an increase in property tax collection of approximately 0.3182%, 0.1208%. and 0.0923% respectively (first cluster), a 1% increase in public investment in transport, culture and sports, caused an increase in property tax collection of approximately 0.2588% and 0.1434% respectively (second cluster) and an increase in the 1% in public investment in transportation and sanitation, caused an increase in property tax collection of approximately 0.2923% and 0.1377% respectively (third cluster).
Keywords. Property tax, public investment, local governments, cluster analysis, longitudinal panel data methodology.
Introducción
El impuesto a la propiedad de bienes inmuebles, denominado «impuesto predial», constituye una herramienta esencial en el financiamiento del sector público, a favor del logro de objetivos de recaudación y planificación urbana, impulsando el desarrollo económico y el bienestar social (Torres & Terán, 2012). A nivel de América Latina, según las Naciones Unidas (2021), en 2018 la recaudación proveniente del impuesto predial representó en pro-medio un 3,8% del ingreso total, 1,7 puntos porcentuales por debajo de los países conformantes de la Organización para la Cooperación y el Desarrollo Económicos (OCDE). De igual manera, para esta organización (OECD, 2021), el impuesto a la propiedad de bienes inmuebles es de mayor aplica-ción en los países de América Latina, siendo administrado en mayoría de los casos por los gobiernos locales, cuya recaudación a este nivel oscila entre el 0,08% del PIB en Honduras y el 0,79% del PIB en Colombia.
En el Perú, bajo disposición del Decreto Legislativo N.º 776 y el Decreto Supremo N.º 156-2004-EF, se define el impuesto predial como aquel que grava el valor de los predios urbanos y rústicos (terrenos, edificaciones e instalaciones fijas y permanentes), siendo responsable de la recaudación, administración y fiscalización la municipalidad a la que pertenezca el predio. Cuantitativamente, la recaudación del impuesto predial en nuestro país en 2018 representó en promedio el 0,24% del PIB, ubicándose por encima de Ecuador, Guatemala, Honduras y México, y por debajo de Argentina, Brasil, Chile, Costa Rica, Colombia y Uruguay (OECD, 2021). Según Arias (2021), este porcentaje, relativamente bajo frente a sus similares, es explicado por tres razones fundamentales: (i) la base del impuesto predial no es el valor de mercado; (ii) la inexistencia de un catastro nacional de predios con un sistema de registro, que impide la identificación de estos; y (iii) un débil control de gran parte de las municipalidades en asuntos tributarios.
En nuestro país, con la promulgación de la Ley N.º 27783 en el año 2002, Ley de Descentralización, se establece la autonomía política y admi-nistrativa en los tres niveles de gobierno, mediante un equilibrado ejercicio de poderes. Con la Ley N.º 27972 en 2003, se establecen normas, meca-nismos y regímenes de carácter municipal, mediante lo cual se da inicio a un proceso normativo con la asignación de competencias, sean de carácter exclusivo o compartido, en busca del beneficio social de la población, siendo esta partícipe principal en el respeto de sus derechos y el cumplimiento de sus deberes.
El respeto de sus derechos en el sentido de ser atendidos de manera opor-tuna en la dotación de bienes y servicios públicos (proyectos de inversión) que permita mejorar su calidad de vida, y el cumplimiento de sus deberes
tiene que ver con las obligaciones establecidas por la ley, como son los tri-butos, entre ellos el pago del impuesto predial, el cual es un ingreso propio de la municipalidad. No obstante, según el BCRP (2019), en promedio, la recaudación por impuesto predial en las municipalidades alcanza el 0,2% del PIB, por lo que Yepes y De los Ríos (2017) sostienen que son cada vez más los municipios dependientes de las transferencias y asignaciones pre-supuestales, que, conforme al artículo 69 de la Ley N.º 27972, al igual que los tributos, son parte de las rentas municipales, con la diferencia de que son asignados por el gobierno central.
Estas transferencias han generado un limitante en el esfuerzo de los municipios en el recaudo de sus tributos locales, conocido como pereza fiscal (Muñoz, 2015), haciendo que las entidades locales antepongan la percep-ción de las transferencias a la puesta en marcha de medidas que mejoren su recaudación local; en razón de ello, el Ministerio de Economía y Finanzas, mediante la Ley N.º 29332 del año 2009, creó el Plan de Incentivos a la Mejora de la Gestión Municipal con el fin de estimular en los gobiernos locales la mejora de la recaudación y administración de los tributos muni-cipales, acompañado de la ejecución del gasto público de manera eficiente y efectiva para la dotación de servicios públicos de calidad y ejecución de inversiones a favor de la población; políticas de Estado que, para Mayoral y Uribe (2010), son esenciales en la disminución de los niveles de desigualdad en el largo plazo, y el desarrollo económico de los países.
Para analizar el efecto de la inversión pública sobre la recaudación del impuesto predial, el presente artículo parte de la hipótesis de que las inver-siones en infraestructura valorizan los predios, ayudan a la fiscalización municipal y, por ende, mejoran la cultura tributaria. El objetivo del presente artículo es analizar la incidencia de la inversión pública sobre la recaudación del impuesto predial en los gobiernos locales provinciales de la macrorregión sur del Perú en el período 2010-2020.
Entre los antecedentes, se cita a Ayala et al. (2017), quienes sostienen que los ingresos de la GAD municipal de Cantón Morona en Ecuador se dan principalmente por ingresos propios, transferencias del Estado y por las entidades financieras públicas; financiamiento que permite a la entidad la realización de proyectos, entre los principales, aquellos de infraestructura, mejora y prestación de servicios públicos, que favorecen en la calidad de vida de la población.
Asimismo, Camacho et al. (2017) concluyen que el recaudo del impuesto predial unificado en el distrito de Barranquilla, Colombia, se encuentra influenciado por los ingresos tributarios diferentes del predial, la tasa efectiva del recaudo, el avalúo catastral del predio y el tipo impositivo medio. En lo
referente a los determinantes económicos de la recaudación del impuesto predial, Unda y Moreno (2015) sostienen que la participación del Estado con la aportación de transferencias de recursos para fines específicos hacia los municipios determinó de manera negativa la recaudación del impuesto predial en México. Ello se ve respaldado por la investigación de Valenzuela y Hinojosa (2017), quienes concluyen que las transferencias federales presentan efectos negativos en los niveles de ingresos fiscales estatales con información de 31 entidades federativas mexicanas durante el período 2005-2014, lo que demuestra la necesidad de un ajuste estructural al sistema de transferencias. Por otra parte, Astudillo (2010) señala que las transferencias de recursos son la segunda fuente de ingreso entre los diferentes niveles de gobierno el sistema francés cuyo uso se da en favor de la igualdad; no obstante, evidencia que las transferencias generaron una baja recaudación fiscal de los contribuyentes.
Materiales y métodos
La investigación tiene un enfoque cuantitativo, diseño no experimental, tipo longitudinal transversal (panel data), de carácter descriptivo y expli-cativo (Hernández et al., 2014). El alcance de la presente investigación es la macrorregión sur del Perú, siendo las unidades de análisis 44 munici-palidades provinciales. La información de las variables para el análisis de conglomerados se extrajo de Ceplan (2017) e INEI (2017) respectivamente, mientras que para el análisis econométrico la información proviene del portal Transparencia Económica del Ministerio de Economía y Finanzas (MEF). De este modo, la base de datos en la estimación econométrica comprende información tipo panel que va desde el año 2010 hasta el año 2020.
Análisis de conglomerados
Yepes y De los Ríos (2017) sustentan la presencia de heterogeneidad en el ámbito municipal dadas las diferentes características y condiciones de cada uno de estos, entre ellas la ubicación geográfica, que tiene influencia en el grado de urbanidad y en las actividades económicas desempeñadas por la población. Frente a ello, el agrupamiento de cada unidad de análisis (conglomerado) debe realizarse en función de variables y/o características sobre las cuales las entidades municipales no posean intervención para poder modificarlas por lo menos en el corto plazo (Herrera & Francke, 2009).
El análisis de conglomerados es una técnica estadística multivariante que, a partir de casos y/o variables, busca agruparlos en grupos homogéneos denominados clústeres sugeridos por la propia particularidad que muestran los datos, de tal forma que los individuos con similares características sean
asignados a un mismo clúster (Pérez, 2004). Los métodos de agrupamiento en el análisis de conglomerados precisan de dos medidas de distancia: (i) entre las unidades, dentro de la cual destaca la medida de distancia euclidiana según la siguiente fórmula:
d(x, y) = �∑n (x ‒ y )2 (1)
i=1
i
i
donde d es la distancia entre las unidades de las muestras x e y, y n repre-senta el número de variables de agrupamiento; y (ii) entre los grupos que se van formando (asociativa), de lo que se obtiene la medida de asociación de Ward, cuya estimación, según Everitt, Landau y Leese (2001, citado en Herrera & Francke, 2009), es denotada como:
dk ij = � ni+nk � dki + � nj+nk � dkj + ( ‒nk ) dij (2)
( ) ni+nj+nk
ni+nj+nk
ni+nj+nk
donde dk(ij) representa la distancia entre los grupos k y por la unión de i y j los grupos y n respectivamente, y muestra el número de elementos del grupo i, j o k. Yepes y De los Ríos (2017) sostienen que Ward es factible con la medida de distancia euclidiana, cumpliendo las propiedades del agrupamiento de clústeres.
Metodología para datos longitudinales de panel
Según Greene (2018), el marco básico para la compresión de datos longi-tudinales es un modelo de regresión de la siguiente forma:
yit = x'itβ + z'i α + εit (3)
donde, i = 1, 2, …, N unidades sociales y t = 1, 2, …, T observaciones en el tiempo, y ε es el término de error que representa los efectos de las variables omitidas en el modelo. Asimismo, existen k regresores en x'it, sin incluir el término constante, la heterogeneidad (principal objetivo de aplicar data panel) está representada por , donde contiene un término constante y un conjunto de variables (observables y no observables).
Mínimos cuadrados ordinarios con errores estándar corregidos (PCSE)
Método propuesto por Beck y Katz (1995) como alternativa de solución ante la presencia de heteroscedasticidad y autocorrelación, así como mayor eficiencia frente a los estimadores de mínimos cuadrados generalizados factibles (MCGF). Matemáticamente, la fórmula idónea para el cálculo de la variabilidad de una estimación a través de mínimos cuadrados ordinarios (MCO) esta expresada por la raíz cuadrada de los términos diagonales de:
Cov(β� ) = (x'x)‒1[x'Ωx](x´x)‒1 (4)
Si se cumple el supuesto esférico en los errores, la ecuación (4) se reduce en un MCO común. Por otra parte, si los errores siguen la estructura de panel, la ecuación esta sujeta a proporcionar errores estándar incorrectos. Para solu-cionar esto, la ecuación (4) puede combinarse con la estructura de panel de errores, proporcionado así los errores estándar corregidos por el panel (PCSE). Para aquellos modelos de panel con errores contemporáneos correla-cionados y heteroscedásticos, Ω representa una matriz diagonal NT x NT con una matriz N x N de covarianzas contemporáneas, Σ, a lo largo de la diagonal. La estimación de la ecuación (4) requiere de la estimación de Σ, mediante MCO; asimismo, sea ei,t el residuo generado por MCO para la unidad i en el período t, se puede estimar un elemento típico de Σ denotado
por la siguiente ecuación:
ΣT e e
Σ�i,j =
t=1 i,t
T
j,t
(5)
Así, el cálculo de los errores estándar corregidos para panel se da tomando los elementos de la raíz cuadrada de la diagonal Σ, matemáticamente expre-sado como:
T
(x'x)‒1x' � (E'E) ∎ IT� x(x'x)‒1 (6)
donde E denota la matriz T x N de los residuos de MCO. La ecuación (6) permite obtener parámetros β� consistentes respecto a MCO tradicional.
Regresión de Prais-Winsten
Es una transformación de mínimos cuadrados generalizados factibles dise-ñada por Prais y Winsten en 1954, que permite la estimación de modelos consistentes frente a la presencia de heteroscedasticidad y autocorrelación en las perturbaciones aleatorias (Salvador, 2020), y que surge de:
yt = α + βxt + εt (7)
donde εt constituye la perturbación aleatoria y, bajo el supuesto de existencia de autocorrelación AR(1), se cumple:
εt = (ρεt‒1) + ct│ρ) < 1 (8)
donde ct representa la perturbación aleatoria ruido blanco. Para la obtención del modelo con ruido blanco, la técnica de Prais-Winsten asume un proceso AR(1) en el modelo, especificándose en la siguiente ecuación:
yt ‒ ρyt‒1 = α(1 ‒ ρ) + β(xt ‒ ρxt‒1) + ct (9)
donde ct es el ruido blanco.
Finalmente, con la introducción de raíces, la ecuación (9) sería reescrita como:
√1 ‒ ρ2y1 = α√1 ‒ ρ2 + �√1 ‒ ρ2x1�β + √1 ‒ ρ2c1 para t = 1 (10)
La ecuación (10) se estima mediante el método de mínimos cuadrados.
Resultados
- Análisis descriptivo
En la figura 1 se aprecia la orientación del comportamiento de la inversión pública según sectores de los gobiernos locales de la macroregión sur del Perú, donde se evidencia que gran parte de las municipalidades presentan una orientación activa de ejecución de proyectos en los sectores de transporte y saneamiento, seguidos por los proyectos en el sector educación, cultura y deporte.
Asimismo, se evidencia que el comportamiento del impuesto predial está altamente relacionado con la ejecución de proyectos enmarcados en el sector transporte, seguidos por los proyectos de saneamiento básico en la mayoría de las municipalidades de la macrorregión sur del Perú. Por lo que, dentro de la priorización de proyectos en los procesos de presupuesto participativo, es notoria la demanda de obras viales; sin embargo, estos tienen que ser corroborados con el cierre de brechas del nuevo Sistema de Programación Multianual y Gestión de Inversiones (Invierte.pe) y las políticas propias de los titulares de pliego. Un dato importante por resaltar es la caída de la recaudación del impuesto predial en el año 2020, resultado de la emergencia sanitaria de la COVID-19 que, según la ONU (2020), trajo como consecuencia una disminución en los ingresos de las familias y de las micro-, pequeñas y medianas empresas en todo el mundo. En el Perú esto ocasionó, de manera específica en la Programación Multianual de Inversiones (PMI), la poca celeridad de ejecución de proyectos en los sectores de transporte y saneamiento.
400
Gasto en inversiones (Millones)
350
300
250
200
150
100
50
Gastos de inversión y recaudación de impuestos, 2010-2020
Año 2020
Emergencia sanitaria por el Sars-Cov2
10.9 11.6
13.6
11.1
10.1
9.6
6.8
7.8
8.1
5.9
4.6
Recaudación del impuesto predial (Porcentaje)
16
14
12
10
8
6
4
2
- 0
2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020
Transporte Saneamiento Cultura y deporte Educación Impuesto predial
Fuente: Transparencia Económica del Ministerio de Economía y Finanzas (MEF). Elaboración: equipo de trabajo.
En la figura 2 literal (a), se observa que las inversiones de las munici-palidades de la región Moquegua están orientadas más a la ejecución de proyectos en el sector transporte, que representa alrededor del 60%, seguido del sector saneamiento con un 30% y, en menor proporción, por proyectos de educación; este último resultado puede deberse al cierre de brechas en dicho sector. Por otra parte, la ejecución de inversiones en las municipali-dades provinciales de la región Cusco se encuentran encaminadas al sector transporte, que representa el 45%, seguido por el sector saneamiento con el 20%. Por último, en las municipalidades de la región Puno, las inver-siones han venido siendo destinadas al sector transporte, con un promedio del 40%, y en un 30% al sector saneamiento. Estos resultados pueden ser factores determinantes de la recaudación del impuesto predial en dichos gobiernos locales.
En la figura 2 literal (b), se muestra que las inversiones de las munici-palidades de las regiones Arequipa, Madre de Dios y Tacna se encuentran orientadas en más del 50% a la ejecución de proyectos en el sector transporte y en menor cantidad a la ejecución de proyectos de saneamiento, que no superan el 30% de su inversión total. Estos resultados pueden deberse al grado de urbanidad propio de cada región y provincia, acompañado por los criterios de cierre de brechas de Invierte.pe.
Figura 2
Orientación de la inversión pública, 2010-2020 (millones de soles)
(a) (b)
70
60
50
40
30
20
10
0
Transporte Transporte
70
60
50
40
30
20
10
0
Educacion
Saneamiento Educación
Saneamiento
Cultura y deporte
Total Puno Cusco Moquegua
Cultura y deporte
Total Arequipa Madre de Dios Tacna
Fuente: Transparencia Económica del Ministerio de Economía y Finanzas (MEF). Elaboración: equipo de investigación.
En la figura 3, se evidencia el grado de participación que tienen los gobiernos locales en lo que corresponde a la recaudación del impuesto pre-dial por regiones. Las municipalidades provinciales de la región Arequipa alcanzan el 26% del total de la recaudación acumulada del impuesto predial total, seguidas por las municipalidades de la región Puno, que alcanzan el 23% de recaudación, y por las municipalidades de la región Tacna, que representan el 18%. Por su parte, las municipalidades de las regiones Cusco, Moquegua y Madre de Dios alcanzaron el 16%, 10% y 6% de la recaudación del impuesto predial total acumulado respectivamente. Según los datos, se observa que una mayor inversión en proyectos en el sector transporte por parte de los gobiernos locales permitirá aumentar la posibilidad de mejorar los niveles de recaudación del impuesto predial, situación que se observa en las municipalidades de las regiones Arequipa, Tacna, Puno y Cusco de la figura 2 literales (a) y (b). Por otro lado, el comportamiento de la recaudación acumulada del impuesto predial es influenciado por el tipo de urbanidad de cada localidad (provincia) y por el grado de responsabilidad ciudadana; es decir, la cultura tributaria en la población de cada provincia.
Figura 3
Recaudación del impuesto predial, 2010-2020 (millones de soles)
30
25
20
15
10
5
0
Puno
Tacna
Moquegua
Arequipa
Cusco
Madre de Dios
Fuente: Transparencia Económica del Ministerio de Economía y Finanzas (MEF). Elaboración: equipo de investigación.
- Análisis de conglomerados
El análisis de conglomerados fue realizado en función de las variables: (i) nivel de pobreza, dado que expresa el nivel de calidad de vida y las posibilidades de generación de ingresos que permitan la obligación de las personas con su comuna local (impuestos); y (ii) el nivel de urbanidad, que muestra el grado de concentración de la población en zonas con actividades económicas en pleno funcionamiento, generando mayor posibilidad de las autoridades competentes en la recaudación de impuestos. En el anexo 1, se ilustra la sintaxis de programación en el software R.
En la figura 4, se aprecia el número de clústeres óptimo y dendrograma, los cuales proponen la clasificación de los municipios provinciales en tres categorías o conglomerados, distribuidos como se presenta en la tabla 1.
Figura 4
10
Agrupación jerárquica, dendrograma
C1
C2
C3
Frecuencia entre todos los índices
6
8
7.5
Height
4
5.0
0 2
2.5
Moho La Union Paucartambo
Sandia Chumbivilcas Chucuito Huancané Paruro Condesuyos
Canas Tarata
San Antonio de Putina
Candarave Yunguyo El Collao Lampa
Quispicanchis
Melgar Azángaro Carabaya Acomayo Mariscal Nieto
Islay Tacna
San Román Camana Tahuamanu Jorge Basadre Arequipa
Ilo Cusco Tambopata
La Convención
Anta Calca Caylloma Puno Castilla
General Sanchez Cerro
Canchis Espinar Manu Carveli Urubamba
0.0
0 1 2 3 4 7 8 9
Número de clusters
Fuente: elaboración propia sobre la base de los resultados del software R.
Asimismo, en la tabla 1 se muestra que el primer clúster concentra municipalidades provinciales con menor grado de urbanidad y mayor nivel de pobreza, y que en la mayoría de los casos su población presenta una inadecuada provisión o carece de servicios básicos de infraestructura como comunicaciones, saneamiento, electricidad y transportes; donde destacan las provincias de Moho y Huancané en la región Puno; La Unión, en la región Arequipa; Acomayo, en la región Cusco; Tarata, en la región Tacna, entre otras. Por otro lado, en el segundo clúster se encuentran aquellas munici-palidades provinciales con un alto grado de urbanidad y menor nivel de pobreza, en las que destaca su grado de actividad económica, como son los casos de Ilo y Tacna, que son puntos clave en el transporte de mercancías dada su ubicación geográfica; Cusco, que brinda su potencial turístico como impulsor de su economía; o la provincia de San Román, considerada como zona de intercambio comercial. Por último, en el tercer clúster se encuentran municipalidades con un grado de urbanidad y nivel de pobreza por debajo del promedio general, entre ellas las provincias de Puno, Caravelí, Espinar y Canchis, cuyas economías tienen un grado medio de dinamismo.
Tabla 1
Categorías de las municipalidades provinciales
Categoría | Casos | Promedio | |
Urbanidad (%) | Pobreza (%) | ||
Clúster 1 | 21 | 37,72 | 48,29 |
Clúster 2 | 11 | 90,97 | 10,15 |
Clúster 3 | 12 | 54,09 | 23,38 |
Total | 44 | 55,50 | 31,96 |
Fuentes: Ceplan (2017) e INEI (2017). Elaboración: equipo de investigación.
- Análisis econométrico
A partir de la estructura de data panel denotada en la ecuación (3), se ha especificado el modelo econométrico de la siguiente manera:
LIPit = ai + β1LTit + β2LSit + β3LCDit + εit (11)
Donde:
LIPit: logaritmo natural del impuesto predial recaudada de la i-ésima muni-cipalidad provincial en el t-ésimo período.
LTit: logaritmo natural del gasto devengado en inversión en transportes de la i-ésima municipalidad provincial en el t-ésimo período.
LSit : logaritmo natural del gasto devengado en inversión en saneamiento de la i-ésima municipalidad provincial en el t-ésimo período.
LCDit: logaritmo natural del gasto devengado en inversión en cultura y deporte de la i-ésima municipalidad provincial en el t-ésimo período.
εit: término de error.
Para identificar el modelo econométrico apropiado, se han seguido los siguientes pasos: (i) la estimación del modelo por mínimos cuadrados ordinarios (MCO); (ii) la estimación por efectos aleatorios (EA), para luego aplicar la prueba de Breush Pagan (BP) que permitió elegir entre MCO y EA; con respecto a la prueba BP, se rechaza la hipótesis nula al 5%, lo que muestra evidencia estadística de heterogeneidad no observada del término de error, esto implica la necesidad de usar alguna de las metodologías que tengan en cuenta la existencia de efectos constantes en el tiempo (hacer uso de panel); (iii) se ha realizado la estimación por efectos fijos (EF); y (iv) se ha realizado la prueba de Hausman para la elección entre EA y EF; con respecto a la prueba de Hausman, se rechaza la hipótesis nula, lo que permite concluir que existe endogeneidad, que los parámetros de EF son eficientes e insesgados, y que prevalecen sobre los coeficientes de EA.
Los resultados de las regresiones econométricas por EF se aprecian en la tabla 2 (en el anexo 2, se ilustra la estimación econométrica completa que arroja el software Stata), donde se muestran las variables, sus coeficientes y estadísticos «z» para cada clúster. Asimismo, respecto a la significancia conjunta, al 1% de significancia se rechaza la hipótesis nula de que los β de las variables regresoras sean iguales a cero; es decir, conjuntamente las variables explicativas influyen en el modelo para cada clúster identificado. En lo que respecta al primer clúster, las variables transporte y saneamiento son significativas al 1%, y la variable cultura y deporte resultó significativa al 10%; en términos cuantitativos, un incremento del 1% en el gasto deven-gado en inversión en transportes supone un crecimiento del 0,3182% en la recaudación del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en transportes, la recaudación del impuesto predial aumentará en aproximadamente 0,014 soles). Un aumento del 1% en el gasto devengado en saneamiento implica un aumento en un 0,1208% en la recaudación del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en saneamiento, la recaudación del impuesto predial aumentará en aproximadamente 0,0048 soles). El aumento del 1% del gasto devengado en cultura y deporte supone un aumento del 0,0923% del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en cultura y deportes, la recaudación del impuesto predial
aumentará en aproximadamente 0,009 soles).
Las estimaciones del segundo clúster evidencian significancia individual en las variables transporte, y cultura y deporte, al 5%, así como signos posi-tivos en sus coeficientes. En términos cuantitativos, un incremento del 1% en el gasto devengado en transporte supone un crecimiento del 0,2588% de la recaudación del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en transporte, la recaudación del impuesto predial aumentará en aproximadamente 0,087 soles).
El aumento del 1% del gasto devengado en inversiones en cultura y deporte se traduce en un crecimiento del 0,1434% en la recaudación del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en cultura y deporte, la recaudación del impuesto predial aumentará en aproximadamente 0,3328 soles).
Tabla 2
Presentación de los resultados econométricos
Variables | Modelo logarítmico LIPit = ai + β1LTit + β2LSit + β3LCDit + εit | ||
Clúster 1 | Clúster 2 | Clúster 3 | |
Constante | 3,6048 (3,44)*** | 9,1427 (4,86)*** | 5,4245 (3,80)*** |
Transporte (LT) | 0,3182 (4,74)*** | 0,2588 (2,49)** | 0,2923 (3,44)*** |
Saneamiento (LS) | 0,1208 (3,20)*** | -0,0557 (-1,23) | 0,1377 (2,79)*** |
Cultura y deporte (LCD) | 0,0923 (1,93)* | 0,1434 (2,23)** | 0,0725 (1,52) |
R2 | 19,78 | 90,15 | 87,87 |
Estadístico Wald chi2 | 53,00*** | 15,38*** | 30,49*** |
Número de observaciones | 193 | 92 | 96 |
Número de grupos | 21 | 11 | 12 |
Notas. Las cifras entre paréntesis son los estadísticos z. *** Nivel de significancia al 1%, ** al 5% y * al 10%.
Fuente: elaboración propia sobre la base de los resultados del software Stata 16.
Por último, en el tercer clúster, las variables transporte y saneamiento resultaron significativas de manera individual a un nivel de significancia del 1% con signos positivos, lo que indica que un incremento del 1% del gasto devengado en transportes supone un crecimiento del 0,2923% en la recaudación del impuesto predial (el efecto marginal refleja que, por cada sol adicional en inversiones en transporte, la recaudación del impuesto predial aumentará en aproximadamente 0,038 soles).
Asimismo, un aumento del 1% en el gasto devengado en saneamiento redunda en una mejora en la recaudación del impuesto predial en 0,1377% (el efecto marginal refleja que, por cada sol adicional en inversiones en saneamiento, la recaudación del impuesto predial aumentará en aproxima-damente 0.064 soles).
No sobra enfatizar el hecho de que los modelos econométricos esti-mados en cada clúster han pasado por el análisis de heteroscedasticidad y autocorrelación. En efecto, para el clúster 1 se ha realizado el análisis de la heterocedasticidad a través de la prueba de Wald modificada para EF1, cuyo resultado evidenció la presencia de heteroscedasticidad. Por otro lado, el análisis de autocorrelación realizado mediante la prueba de Woldridge para EF evidenció ausencia de autocorrelación. En consecuencia, se ha proce-dido con la estimación de panel data para modelos con heteroscedasticidad mediante la técnica de mínimos cuadrados ordinarios con errores estándar corregido (PCSE), observándose en el modelo estimado significancia indi-vidual y conjunta.
Para los clústeres 1 y 2, se ha evidenciado presencia de autocorrelación y heterocedasticidad. Para garantizar estimaciones eficientes de los pará-metros, se ha recurrido a la estimación de panel data para modelos con heteroscedasticidad y autocorrelación. En este contexto, se ha utilizado el método de Prais-Winsten (PW), observándose en cada clúster significancia individual y conjunta.
- Análisis de política pública
En lo que concierne a las finanzas municipales dentro de los gobiernos loca-les, existen sectores que tienen mayor incidencia sobre la recaudación de los ingresos mediante impuestos municipales y recursos directamente recaudados. El hecho de la presencia del signo negativo en el sector saneamiento (a pesar de su no significancia), en el segundo clúster, supone que el reconocimiento de los proyectos de inversión en el mencionado sector aún es débil por parte de la ciudadanía, limitando la toma de conciencia que permita una cultura tributaria responsable en el interior de su jurisdicción. Asimismo, es sabido que, dentro de la secuencia de realización de obras de transporte, la ejecución de proyectos de saneamiento es condicionante para la puesta en marcha de dichos proyectos, que involucran en gran parte el asfalto y pavimento
1 El problema de heteroscedasticidad puede analizarse a través de la prueba del multiplicador de Lagrange de Breush y Pagan; no obstante, según Greene (2018), esta prueba es sensible al supuesto de normalidad de los errores; frente a ello, la prueba modificada de Wald para heteroscedasticidad es idónea aun cuando el mencionado supuesto es violado.
rígido, que son factores conducentes en su mayoría al cumplimiento de las obligaciones tributarias en la población (impuesto predial).
Una situación deseable es impulsar la inversión y ejecución de proyectos de saneamiento, así como de transportes (pistas y veredas), de manera que estos sean realizados en simultáneo; no obstante, prima la decisión y obli-gación en la dotación de bienes y servicios que permitan elevar la calidad de vida de las personas, y a su vez el cierre de brechas como parte del nuevo sistema de inversión pública, el cual involucra un diagnóstico eficiente y oportuno de las necesidades de inversión de la ciudadanía por parte de las autoridades competentes. Por otra parte, el sector cultura y deporte, presenta una influencia poco preponderante sobre la variable impuesto predial, pues centrarse en la ejecución de proyectos en el contexto de dicho sector hace que el impuesto predial tienda a crecer en menor proporción.
De todo el análisis realizado, se puede inferir que, en materia de recau-dación de impuestos municipales (el caso del impuesto predial), el sector transporte cumple un papel sumamente importante en el comportamiento del impuesto predial en la totalidad de los gobiernos locales de la macro-rregión sur del Perú.
De manera analítica, por clúster, se observa que, en el primero, los sectores transporte, saneamiento, y cultura y deporte, tienen una injeren-cia significativa sobre la variable recaudación del impuesto predial. Estos resultados nos llevan a manifestar que en las municipalidades provinciales situadas en este clúster aún existen brechas por cerrar o estas no vienen siendo atendidas en sus necesidades y proyectos identificados dentro de los procesos de presupuesto participativo; por ello, estas provincias tienen mayores niveles de pobreza y un menor grado de urbanidad. En ellas, las autoridades tienen que impulsar la ejecución de proyectos con la finalidad de que los ciudadanos cumplan con sus obligaciones tributarias locales.
En el segundo clúster, se puede observar que la variable saneamiento no es significativa para los gobiernos locales de este grupo porque se ha evidenciado dentro de estas municipalidades provinciales un alto grado de urbanidad y menores niveles de pobreza. En ellas, los gobiernos locales vienen ejecutando proyectos enfocados en los sectores transporte, y cultura y deporte, orientados al mejoramiento del ornato de la ciudad.
Por último, en los gobiernos locales del tercer clúster, y considerando los resultados econométricos, se observa que las variables de transporte y saneamiento tienen una injerencia significativa sobre la recaudación del impuesto predial a diferencia de la variable cultura y deporte. Lo anterior se justifica debido a que este grupo de municipalidades provinciales presentan un grado de urbanidad y nivel de pobreza por debajo del promedio nacional.
Asimismo, existe un grado medio de dinamismo en su economía local. En consecuencia, las autoridades locales deben priorizar proyectos dentro del sector transporte y saneamiento, que seguramente guardan relación con la cartera de inversiones y los resultados de los talleres de presupuesto partici-pativo, los cuales representan factores primordiales para que la ciudadanía pueda cumplir con sus obligaciones tributarias de pago del impuesto predial.
Conclusiones
Teniendo en cuenta el análisis de conglomerados y las estimaciones eco-nométricas, en primer lugar, se concluye que la variable con mayor relevancia en los modelos es la concerniente a la inversión en transporte, con una incidencia positiva y significativa sobre la recaudación del impuesto predial. En relación con el análisis de la inversión en saneamiento, se observa una incidencia positiva en el recaudo predial en el primer y el tercer clúster res-pectivamente. Por otro lado, se observa que la inversión en cultura y deporte incide positivamente sobre la recaudación del impuesto predial en los dos primeros clústeres, cuyas magnitudes varían entre cada clúster identificado según el nivel de pobreza y urbanidad.
El análisis público evidencia que la inversión en la macrorregión sur del Perú se orientó en mayor porcentaje hacia el sector transporte, por lo que mejorar el nivel de recaudo del impuesto predial dependerá de la identificación, priorización y ejecución de proyectos de inversión pública en este sector. Lo anterior implica una clara situación de análisis en la toma de decisiones en el marco de los procesos del presupuesto participativo y la normatividad de Invierte.pe. Con relación al sector saneamiento, la evidencia encontrada hace suponer que la ciudadanía aún no ve a este sector como un sector estratégico que los motive e impulse a cumplir con sus obligaciones tributarias, pese a ser condicionante para encaminar la ejecución de pro-yectos de pistas y veredas. Por último, los proyectos en el sector de cultura y deporte inciden en menor grado sobre el impuesto predial, debido a su relación con el nivel de urbanidad de cada espacio geográfico de los gobier-nos locales. Así, en aquellas localidades donde prevalece mayor urbanidad se pone énfasis en la ejecución de proyectos de transitabilidad y recreación. La recaudación de tributos, principalmente el impuesto predial, resulta importante dentro de los gobiernos locales dado que forma parte de sus ingre-sos propios. Una gestión administrativa orientada a una mayor recaudación posibilita el incremento de recursos necesarios para atender, por ejemplo, el financiamiento de proyectos de inversión sujetos a las necesidades de la población y, a través de la dotación de estos bienes o servicios, incentivar la
cultura tributaria en los ciudadanos.
Referencias
Arias, L. A. (2021). Política fiscal y tributaria frente a la pandemia global del coronavirus (L. Ráez, Ed.). 1.ª ed.. Proyecto Perú Debate 2021: propuestas hacia un mejor gobierno. https://www.cies.org.pe/sites/default/files/investigaciones/4._dp_politica_ fiscal_tributaria.pdf
Astudillo, M. (2010). Las transferencias intergubernamentales en el sistema fiscal francés. Dimensión Económicas, 2(4), 13-25. https://rde.iiec.unam.mx/revistas/4/articulos/ 2/21.php
Ayala, S., Zaruma, M., & Barragán, H. (2017). Origen y destino de recursos del presu-puesto de los GAD municipales: estudio del Cantón Morona. Killkana Sociales, 1(3), 29-36. https://doi.org/https://doi.org/10.26871/killkana_social.v1i3.60
BCRP. (2019). Reporte de inflación, marzo 2019. Panorama actual y proyecciones macro-económicas 2019-2020. https://www.bcrp.gob.pe/docs/Publicaciones/Reporte-Infla cion/2019/marzo/reporte-de-inflacion-marzo-2019.pdf
Beck, N., & Katz, J. (1995). What to do (and not to do) with time-series cross-section data. The American Political Science Review, 89(3), 634-647. https://doi.org/https:// doi.org/10.2307/2082979
Camacho, C., Hurtado, K., Navarro, E., Hurtado, J., & Nieves, S. (2017). Factores que inciden en el recaudo del impuesto predial unificado en el distrito de Barranquilla. Prospectiva, 15(1), 64-73. https://doi.org/https://doi.org/10.15665/rp.v15i1.591
Ceplan. (2017). Información departamental, provincial y distrital de población que requiere atención adicional y devengado per cápita, Perú. https://www.ceplan.gob.pe/wp-content/uploads/2017/08/Matriz-de-indicadores-nacionales-a-Julio-de-2017.pdf
D. L. N.º 776. (1993). Ley de Tributación Municipal. https://cdn.www.gob.pe/uploads/ document/file/251493/226857_file20181218-16260-14de8v2.pdf
D. S. N.º 156-2004-EF. (2004). Decreto que aprueba el Texto Único Ordenado de la Ley de Tributación Municipal. https://www.gob.pe/institucion/mef/normas-legales/ 226256-156-2004-ef
Greene, W. (2018). Econometric analysis (8.ª ed.). Pearson. https://www.pearson.com/ us/higher-education/program/Greene-Econometric-Analysis-8th-Edition/ PGM334862.html
Hernández, R., Fernández, C., & Baptista, P. (2014). Metodología de la investigación. (6.ª ed.), McGraw-Hill. https://www.uca.ac.cr/wp-content/uploads/2017/10/Investiga cion.pdf
Herrera, P., & Francke, P. (2009). Análisis de la eficiencia del gasto municipal y sus deter-minantes. Revista Economía, 32, 113-178. https://www.researchgate.net/publica tion/46532022_Analisis_de_la_eficiencia_del_gasto_municipal_y_de_sus_deter minantes/citation/download
INEI. (2017). Censos Nacionales 2017: XII de Población, VII de Vivienda y III de Comuni-dades Indígenas. https://censos2017.inei.gob.pe/redatam/
Ley N.º 27783. (2002). Ley de Bases de la Descentralización. https://www2.congreso.gob. pe/sicr/cendocbib/con4_uibd.nsf/8B3C23D0EB9643D605257FD400782856/
$FILE/27783.pdf
Ley N.º 27972. (2003). Ley Orgánica de Municipalidades. https://www.mef.gob.pe/conte nidos/presu_publ/capacita/programacion_formulacion_presupuestal2012/Anexos/ ley27972.pdf
Ley N.º 29332. (2009). Ley que crea el Plan de Incentivos a la Mejora de la Gestión Munici-pal. https://www.mef.gob.pe/en/por-instrumento/ley/9386-ley-n-29332/file
Mayoral, F., & Uribe, C. (2010). Determinantes económicos e institucionales del esfuerzo fiscal en América Latina. Investigación Económica, 69(273)(julio-septiembre), 85-113. http://www.scielo.org.mx/pdf/ineco/v69n273/v69n273a3.pdf
Muñoz, B. (2015). La pereza fiscal como problema relevante en la descentralización. Revista de Investigación en Administración, Contabilidad, Economía y Sociedad, 4(4), 141-
151. https://www.redalyc.org/journal/5518/551857280005/551857280005.pdf Naciones Unidas. (2021). Panorama fiscal de América Latina y el Caribe. Los desafíos de la
política fiscal en la recuperación transformadora pos-COVID-19.
OECD. (2021). Revenue statistics – Latin American countries: Comparative tables. OECD. Stat. https://stats.oecd.org/index.aspx?DataSetCode=RSLACT#
ONU. (2020). La pandemia de COVID-19 reduce gravemente los ingresos de los trabajadores. https://news.un.org/es/story/2020/09/1481042
Pérez, C. (2004). Técnica de análisis multivariante de datos. Pearson Educación, S. A. Salvador, J. (2020). Algunos determinantes del empleo en el sector público a través de los países.
Análisis econométricos en datos de panel 1995-2017. Cuadernos del INAP 44. Ins-tituto Nacional de la Administración Pública. https://www.argentina.gob.ar/sites/ default/files/cuinap_44.pdf
Torres, M., & Terán, H. (2012). El impuesto predial urbano. Análisis teórico y aproxima-ción inicial al caso de la ciudad de Caracas. Revista de Ciencias Sociales, 18(3), 416-
429. https://produccioncientificaluz.org/index.php/rcs/article/view/24988
Unda, M., & Moreno, C. (2015). La recaudación de impuesto predial en México: un aná-lisis de sus determinantes económicos en el período 1969-2010. Revista Mexicana de Ciencias Políticas y Sociales, 60(225), 45-78. http://dx.doi.org/10.1016/S0185-1918(15)30019-2
Valenzuela, R., & Hinojosa, A. (2017). Las transferencias federales, los contrapesos políti-cos y los ingresos fiscales estatales en México. Journal of Economic Literature, 14(42), 47-71. https://doi.org/10.1016/j.eunam.2017.09.002
Yepes, M., & De los Ríos, R. (2017). Determinantes de la recaudación de los impuestos muni-cipales. Cepal. https://www.cepal.org/sites/default/files/events/files/yepes_0.pdf
Anexo 1
Sintaxis y salidas del programa R-Studio (ML Conglomerados)
pkgs <- c(“factor extra”, “NoClust”, “cluster”, “fpc”, “dendextend”) install.packages(“pkgsearch”)
install.packages(“tidyverse”) library(factoextra);library(NbClust);library(fpc);library(cluster) library(readxl)
library(tidyverse)
Cluster<-read.csv(“C:/Users/LENOVO/Desktop/Investigación/R (BD y ST)/Data/Cluster1.csv”, row.names =“provincia”)
Filter
urbano pobreza
Azángaro
37.55
55.0
Carabaya 49.94 50.2
Chucuito 29.50 49.3
El Collao 39.32 42.3
Huancane 25.43 51.7
Lampa 44.26 43.6
Melgar 59.99 48.9
Moho 29.72 64.6
Puno 66.87 24.8
San Antonio de Putina 70.67 43.1
San Roman 91.31 18.9
View(Cluster)
ggplot(Cluster, aes(y=urbano))+geom_boxplot()
100
75
urbano
50
25
0
-0.4 -0.2 0.0 0.2 0.4
ggplot(Cluster,aes(y=pobreza))+geom_boxplot()
60
pobreza
40
20
-0.4 -0.2 0.0 0.2 0.4
df<-scale(Cluster) View(df)
head(df)
urbano | pobreza | |
Azangaro | -0.6535110 | l.3035976 |
Carabaya | -0.2024360 | l.0319987 |
Chucuito | -0.9465824 | 0.9810739 |
El Co11ao | -0.5890717 | o. 5849921 |
Huancane | -1.0947563 | 1.1168733 |
Lampa | -0.4092242 | 0.6585502 |
p1<-princomp(df) fviz_pca_ind(p1)
Individuals - PCA
Tarata | |||||||||
Candarave | |||||||||
Melgar | San Antonio | de Putina | |||||||
Mo | ho Az | Acoma angaro | Coarabaya | ||||||
Yun | guyo | Caylloma | San Rom | an Tacna Islay Are | quipa | ||||
Huancan | e Chucuito | Quispicanchis Lampa El Collao | Canchis | Puno | Camana | Mariscal Nieto | Ilo Cusco | ||
Paucartambo | La Union Chu | mbivilcas Sandia | Calca | Espinar Castilla General Sanchez | Cerro | Tam | bopata | ||
Carveli | Jorge Ba | Tahuamanu sadre | |||||||
Paruro | Canas | Condesuyos | Anta | Uruba | mba | ||||
La Co | nvencion | ||||||||
Manu |
1.0
y
Dim2 (13.8%)
0.5
0.0
-0.5
-1.0
-2 -1 0 1 2
Dim1 (86.2%)
k1<-kmeansruns(df, krange = 2, runs = 100) fviz_cluster(k1, data = df)
Cluster plot
Moho
Canas
Union Huancane
Chucuito Chumbivilcas
Azangaro
Acomayo
Carabaya
Paruro
Yunguyo
Sandia El Collao Lampa
Quispicanchis
San Antonio de Putina
Condesuyos
Candarave
Castilla Canchis
Espinar
Puno General Sanchez Cerro
San Roman ana
Tacna
Mariscal Nieto
Islay
Jorge Basadre
Arequipa
Manu
Tahuamanu
Tambopata
Cus Ilo
co
Cam
Carveli
Urubamba
on
La Convenci
Caylloma
Calca
Anta
Tarata
Melgar
Paucartambo
La
2
1
pobreza
cluster
a 1
0 a 2
-1
-2 -1
urbano0 1
k2<-kmeansruns(df, krange = 3, runs = 100)
fviz_cluster(k2, data = df)+labs(title=“”,x=“Nivel de urbanidad”,y=“Nivel de pobreza”)
Cluster plot
Moho
Canas
Union Huancane
Chucuito Chumbivilcas
Azangaro
Acomayo
Carabaya
Paruro
Yunguyo
Sandia El Collao Lampa
Quispicanchis
San Antonio de Putina
Candarave
Condesuyos
Anta
Castilla Canchis
Espinar
General Sanchez Cerro
Puno
La Convencion
Urubamba
Carveli
San Roman ana
Tacna
Jorge Basadre
Mariscal Nieto
Arequipa
Islay
Manu
Tahuamanu
Tambopata
Cus Ilo
co
Cam
Caylloma
Calca
Tarata
Melgar
Paucartambo
La
2
1
Nivel de pobreza
0
-1
cluster
a 1
a 2
a 3
-2 -1
Nivel de urba0nidad 1
k3<-kmeansruns(df, krange = 4, runs = 100) fviz_cluster(k3, data = df)
Cluster plot
Moho
Azangaro
Acomayo
Carabaya
Paruro
cuito
Chumbivilcas
Yunguyo
Canas
Quispicanchis Sandia El Collao Lampa
San Antonio de Putina
Candarave
Condesuyos
Anta
San Roman ana
Mariscal Nieto
Arequipa
Islay
Tacna
Jorge Basadre
Manu
Tahuamanu
Tambopata
Cus Ilo
co
Cam
Carveli
Urubamba
on
La Convenci
Puno
chez Cerro
lla Canchis
Espinar
Casti General San
Caylloma
Calca
Tarata
Melgar
Union Huancane
Chu
Paucartambo
La
2
1
pobreza
0
-1
cluster
a 1
a 2
a 3
a 4
-2 -1
urbano0 1
fviz_nbclust(df,kmeans, method=“wss”) ##No hay notoriedad del punto
de inflexión##
Optimal number of clusters
Total Within Sum of Square
75
50
25
1 2 3 4 5 6 7 8 9 10
Number of clusters k
fviz_nbclust(df,kmeans, method=“wss”)+
labs(subtitle = “Elbow method”)+geom_vline(xintercept = 2,linetype=2)
Optimal number of clusters
Elbow method
Total Within Sum of Square
75
50
25
1 2 3 4 5 6 7 8 9 10
Number of clusters k
fviz_nbclust(df,kmeans, method=“silhouette”)+
labs(subtitle = “Silhouette method”) ##la gráfica sugiere 2 cluster”
Optimal number of clusters
Silhouette method
0.5
Average silhouette width
0.4
0.3
0.2
0.1
0.0
1 2 3 4 5 6 7 8 9 10
Number of clusters k
set.seed(123)
fviz_nbclust(df, kmeans, nstart=25, method=“gap_stat”, nboot=50)+ labs(subtitle=“Gap statistic method”) ##sugiere 3 cluster
Optimal number of clusters
Gap statistic method
0.20
Gap statistic (k)
0.15
0.10
0.05
0.00
-0.05
1 2 3 4 5 6 7 8 9 10
Number of clusters k
res.nbclust<NbClust(df,distance=“euclidean”,min.nc=2,max.nc=9,metho-d=“complete”,index=“all”) ##sugiere 3 cluster##
************************************************************************
Among all indices:
5 proposed 2 as the best number of clusters
9 proposed 3 as the best number of clusters
4 proposed 4 as the best number of clusters
2 proposed 7 as the best number of clusters
2 proposed 8 as the best number of clusters
1 proposed 9 as the best number of clusters
***** conclusion *****
* According to the majority rule, the best number of clusters is 3
************************************************************************
Hubert Statistic values
0.013
0.015
0.017
0.019
Hubert statistic second differences
0.000
0.001
0.002
0.003
2 3 4 5 6 7 8 9 2 3 4 5 6 7 8 9
Number of clusters Number of clusters
factoextra::fviz_nbclust(res.nbclust)+theme_minimal()+ggtitle(“NbClut´s optimal number of cluster”)+labs(x=“Número de clusters”,y=“Frecuencia entre todos los índices”)+theme_classic()
NbClut’s optimal number of cluster
Frecuencia entre todos los índices
7.5
5.0
2.5
0.0
0 1 2 3 4 7 8 9
Número de clusters
cf2 <- clusterboot(df,B=100,bootmethod = c(“jitter”,“boot”),clustermethod
= kmeansCBI,krange=3,seed = 123) print(cf2)
Cluster %>%
ggplot(aes(x=factor(k2$cluster), y=pobreza, fill=factor(k2$cluster)))+-geom_boxplot()+geom_point()+xlab(“Cluster”)+labs(fill=“Cluster”)
60
pobreza
40
20
1 2 3
Cluster
Cluster
1
2
3
Cluster %>%
ggplot(aes(x=factor(k2$cluster), y=urbano, fill=factor(k2$cluster)))+geom_ boxplot()+geom_point()+xlab(“Cluster”)+labs(fill=“Cluster”)
100
75
urbano
50
25
0
1 2 3
Cluster
Cluster
1
2
3
d<-dist(df,method = “euclidean”) d
> d | |||||
Azangaro | Carabaya | Chucuito | El Collao | Huancane | |
Carabaya | 052,653,075 | ||||
Chucuito | 0.43578933 | 0.74588680 | |||
El Collao | 0.72148887 | 0.59101778 | 0.53356779 | ||
Huancane | 0.47912770 | 0.89634771 | 0.20099008 | 0.73390359 | |
Lampa | 0.68975516 | 0.42687837 | 0.62671787 | 0.19430880 | 0.82462982 |
Melgar | 0.88687963 | 0.37320497 | 1.11026115 | 0.84008919 | 1.26814006 |
Moho | 0.61345266 | 1.09808528 | 0.86575855 | 1.30931223 | 0.74644453 |
Puno | 2.01480724 | 1.56380135 | 1.94236106 | 1.40943427 | 2.14309600 |
San Antonio de Putina | 1.38104628 | 0.85496987 | 1.53935840 | 1.14223722 | 1.71740605 |
San Roman | 2.82897118 | 2.32487749 | 2.83241709 | 2.30990545 | 3.03266109 |
Sandia | 0.74734129 | 0.81811473 | 0.38822259 | 0.30602288 | 0.55761625 |
Yunguyo | 0.48617973 | 0.27458139 | 0.53088261 | 0.33111881 | 0.71098428 |
Arequipa | 3.43568399 | 2.94207444 | 3.40576251 | 2.87498427 | 3.60675254 |
Camana | 2.79368067 | 2.31265307 | 2.74797630 | 2.21614656 | 2.94892022 |
Carveli | 2.32969634 | 1.91947587 | 2.19373610 | 1.66568276 | 2.39028656 |
Castil1a | 1.77238540 | 1.39369552 | 1.62526824 | 1.09755410 | 1.82208456 |
Caylloma | 1.59346502 | 1.11782867 | 1.57906247 | 1.06206225 | 1.77898943 |
Condesuyos | 1.13966466 | 1.32522331 | 0.70458747 | 0.80788993 | 0.74741282 |
Islay | 3.12594462 | 2.62600704 | 3.11453818 | 2.58750649 | 3.31531069 |
La Union | 0.75178188 | 1.20946482 | 0.50597521 | 1.02155139 | 0.31753676 |
hc1<-hclust(d,method = “complete”) plot(hc1, cex=0.6, hang=-1)
Height
0 1
2
3
4
Paruro La Union Paucartambo Chumbivilcas
Chucuito Huancané Condesuyos
Canas Moho Sandia El Collao Lampa
Quispicanchis
Azángaro Yunguyo Carabaya Acomayo Arequipa
Ilo Cusco Tambopata Mariscal Nieto
Islay Tacna
San Román Camana Manu Tahuamanu Jorge Basadre
Melgar Tarata
San Antonio de Putina
Candarave
Puno Carveli Caylloma Castilla
General Sanchez Cerro
Canchis Espinar
La Convención Urubamba
Anta Calca
d
hclust (*, “complete”)
Moho La Union Paucartambo
Sandia Chumbivilcas
Chucuito Huancané Paruro Condesuyos
Canas Tarata
San Antonio de Putina
Candarave Yunguyo El Collao Lampa
d
hclust (*, “ward.D2”)
Quispicanchis
Melgar Azángaro Carabaya Acomayo Mariscal Nieto
Islay Tacna
San Román Camana Tahuamanu Jorge Basadre
Arequipa
Ilo Cusco Tambopata
La Convención
Anta Calca Caylloma
Puno Castilla
General Sanchez Cerro
Canchis Espinar Manu Carveli Urubamba
Height
Inversión pública y recaudación del impuesto predial, una aplicación de datos de panel por conglomerados en gobiernos locales provinciales de la macrorregión sur del Perú, 2010-2020
hc2<-hclust(d,method = “ward.D2”) hc2
plot(hc2, ces=0.6, hang=-1)
223
0 2 4 6 8 10
Anexo 2
Salidas econométricas del programa Stata
. xtpcse limpuesto ltransporte lsaneamiento lcultura, het Number of gaps in sample: 12
Linear regression, heteroskedastic panels corrected standard errors
Group variable: | departamento | Number of obs | = | 193 | ||
Time variable: | año | Number of groups | = | 21 | ||
Panels: | heteroskedastic (unbalanced) | Obs per group: | ||||
Autocorrelation: | no autocorrelation | min | = | 3 | ||
avg | = | 9.1904762 | ||||
max | = | 11 | ||||
Estimated covariances | = | 21 | R-squared | = | 0.1978 | |
Estimated autocorrelations | = | 0 | Wald chi2 (3) | = | 53.00 | |
Estimated coefficients | = | 4 | Prob > chi2 | = | 0.0000 | |
limpuesto | Coef. | Het-corrected Std. Err. | z | P> | z | | [95% Conf. Interval] | |
ltransporte | .3182316 | .0670773 | 4.74 | 0.000 | .1867625 | .4497006 |
lsaneamiento | .12081 | .0377421 | 3.20 | 0.001 | .0468369 | .194783 |
lcultura | .0923074 | .0479132 | l.93 | 0.054 | -.0016007 | .1862155 |
_cons | 3.604848 | 1.046442 | 3.44 | 0.001 | 1.553859 | 5.655837 |
. xtpcse limpuesto ltransporte lsaneamiento lcultura, het c(ar1)
Number of gaps in sample: 12
(note: computations for rho restarted at each gap)
(note: estimates for rho outside [‒1,1] bounded to be in the range [‒1,1]) Prais-Winsten regression, heteroskedastic panels corrected standard errors
Group variable: | departamento | Number of obs | = | 92 | ||
Time variable: | año | Number of groups | = | 11 | ||
Panels: | heteroskedastic (unbalanced) | Obs per group: | ||||
Autocorrelation: | common AR (1) | min | = | 4 | ||
avg | = | 8.3636364 | ||||
max | = | 11 | ||||
Estimated covariances | = | 11 | R-squared | = | 0.9015 | |
Estimated autocorrelations | = | 1 | Wald chi2 (3) | = | 15.38 | |
Estimated coefficients | = | 4 | Prob > chi2 | = | 0.0015 | |
limpuesto | Coef. | Het-corrected Std. Err. | z | P> | z | | [95% Conf. Interval] | |
ltransporte | .2588829 | .1041759 | 2.49 | 0.013 | .0547018 | .4630639 |
lsaneamiento | -.057518 | .0451628 | -1.23 | 0.217 | -.1442692 | .0327657 |
lcultura | .1434575 | .0642145 | 2.23 | 0.025 | .0175993 | .2693157 |
_cons | 9.142738 | 1.88045 | 4.86 | 0.000 | 5.457124 | 12.82835 |
rho | .716588 | |||||
. xtpcse limpuesto ltransporte lsaneamiento lcultura, het c(ar1)
Number of gaps in sample: 11
(note: computations for rho restarted at each gap)
(note: estimates for rho outside [‒1,1] bounded to be in the range [‒1,1]) Prais-Winsten regression, heteroskedastic panels corrected standard errors
Group variable: | departamento | Number of obs | = | 96 | ||
Time variable: | año | Number of groups | = | 12 | ||
Panels: | heteroskedastic (unbalanced) | Obs per group: | ||||
Autocorrelation: | common AR (1) | min | = | 1 | ||
avg | = | 8 | ||||
max | = | 11 | ||||
Estimated covariances | = | 12 | R-squared | = | 0.8787 | |
Estimated autocorrelations | = | 1 | Wald chi2 (3) | = | 30.49 | |
Estimated coefficients | = | 4 | Prob > chi2 | = | 0.0000 | |
limpuesto | Coef. | Het-corrected Std. Err. | z | P> | z | | [95% Conf. Interval] | |
ltransporte | .2923627 | .0849615 | 3.44 | 0.001 | .1258411 | .4588843 |
lsaneamiento | .13777 | .0493067 | 2.79 | 0.005 | .0411307 | .2344093 |
lcultura | .0725053 | .0476287 | 1.52 | 0.128 | -.0208452 | .1658558 |
_cons | 5.424502 | 1.428796 | 3.80 | 0.000 | 2.624113 | 8.22489 |
rho | .5650812 | |||||