Introducción

La gestión y planificación integral de espacios urbanos como calles, banquetas, plazas, parques urbanos (PU) y áreas verdes (AV) es fundamental para promover una mejor calidad de vida (CV) en ciudades (^{Almeida, Mariano, Agostinho, Liu y Giannetti, 2018}; ^{Bertram y Rehdanz, 2015}; ^{González y Sánchez, 2014}; ^{Harnik y Crompton, 2014}; ^{Kim y Jin, 2018}; ^{Obando y Salcedo, 2015}). La importancia económica, social y ambiental de los parques y áreas verdes, como la describen ^{Holguín y Campos (2017)} y ^{Morales-Cerdas, Piedra y Romero (2018)}, radica en que proveen servicios ambientales y socioeconómicos a la población como oxígeno, captura de carbono, biodiversidad, recreación, salud pública (física y psicológica) y turismo, entre otros, que inciden en la calidad de vida citadina (^{Obando y Salcedo, 2015}; ^{Rivera, 2014}; ^{White et al., 2013}; ^{Zhang y Zhou, 2018}). Puesto que uno de los objetivos de la planificación urbana es proveer calidad de vida, entonces es necesario conocerla como parte de la evaluación de programas de provisión, conservación y/o mejoramiento tanto de parques urbanos como de árboles y plantas de ornato (APO) en banquetas, jardines, plazas ajardinadas y camellones. De acuerdo con ^{Urzúa y Caqueo-Urízar (2012)}, el concepto de calidad de vida considera perspectivas científicas como la medicina, la economía y otras, y cada una expresa un enfoque diferente; también señalan que, si bien no existe una definición única sobre CV, hay categorías generales que permiten agrupar definiciones en un marco teórico que la considera como sinónimo de satisfacción personal o satisfacción con la vida, cuyo foco principal de interés es el bienestar subjetivo (BS), entendido como lo que piensan y sienten las personas sobre su vida, cuando se toman en cuenta conclusiones en los ámbitos cognitivo y afectivo (^{Aparicio, 2011}). Skevington y Böhnke, con base en ^{Diener y Diener (1995)}, entienden el BS como “aquel que incluye evaluaciones cognitivas o una evaluación de la satisfacción con la vida como un todo, además de reacciones emocionales a los acontecimientos de la vida” (2018, 23). ^{Ferrans (1990)} define calidad de vida como una sensación de bienestar que se deriva de la satisfacción o insatisfacción con diferentes áreas de la existencia, importantes para la persona. Hay diversas perspectivas teóricas en las ciencias sociales para conocer, analizar y estimar la CV y el BS (^{Denegri, García y González, 2015}; ^{Easterlin, 2015}; ^{Gómez-Álvarez y Ortiz, 2016}). Ambos conceptos se aplican en la valoración económica de bienes y servicios ambientales, que es el tema de este artículo. La economía está investigando la estimación de valores monetarios de bienes/servicios ambientales considerando una medida de CV como variable proxy de la utilidad económica,¹ que se obtiene al preguntarle a la persona (^{Ferreira y Moro, 2010}; ^{Frey, Luechinger y Stutzer, 2010}; ^{Powdthavee, 2009}; ^{Skevington y Böhnke, 2018}). La estimación de dichos valores consiste en obtener estadísticamente un modelo de calidad de vida o de satisfacción personal en función de bienes/servicios (que incluye los ambientales), del ingreso y de variables sociodemográficas y de percepción controladas. Hay trabajos que emplean el enfoque de CV y de satisfacción y utilidad económica, que valoran la contaminación atmosférica, los parques urbanos, las áreas verdes, el ruido citadino, la salud pública, la belleza escénica y la biodiversidad (^{Ambrey y Fleming, 2011}; ^{Du, Shin y Managi, 2018}; ^{Jones, 2017}; ^{Kopmann y Rehdanz, 2013}; ^{Levinson, 2012}; ^{Luechinger y Raschky, 2009}; ^{MacKerron y Mourato, 2009}). El objetivo de esta investigación es modelar la percepción de calidad de vida de habitantes de Ciudad de México en función de su percepción sobre la presencia de parques, árboles y plantas de ornato en su entorno, considerando también variables sociodemográficas, provisión de servicios públicos y de percepción de seguridad. Y, con los resultados, obtener y analizar el valor monetario que tiene la presencia de parques, árboles y plantas de ornato en el entorno donde residen los habitantes, tomando en cuenta deciles altos, bajos y promedio de ingresos mensuales del hogar.

Antecedentes

El valor de los bienes y servicios ambientales sin precio de mercado se estima a partir de la perspectiva metodológica de la economía ambiental. Para los servicios que prestan los parques urbanos y las áreas verdes se citan los métodos de valoración contingente y de precios hedónicos (^{Chen y Qi, 2018}; ^{Haro-Martínez y Taddei-Bringas, 2014}; ^{Kim, Lee, Lee y Choi 2019}; ^{Latinopoulos, Mallios y Latinopoulos, 2016}; ^{Raffo y Mayta, 2015}; ^{Romero y Vargas, 2016}; ^{Verbič, Slabe y Klun, 2016}). Sin embargo, la investigación económica y sus resultados sobre calidad ambiental, uso/consumo de bienes/servicios ambientales, bienestar subjetivo y calidad de vida de las personas es un área emergente para valorar monetariamente los bienes y servicios ambientales (^{Ambrey y Fleming, 2013}; ^{Ferreira y Moro, 2010}; ^{Kopmann y Rehdanz, 2013}; ^{Levinson, 2012}).

Los trabajos de ^{Smyth, Mishra y Qian (2008)}; ^{Ambrey y Fleming (2013)}; ^{Bertram y Rehdanz (2015)} y ^{Vargas y Roldán (2018)} abordan el BS, la CV y la valoración económica de los servicios que prestan los parques urbanos y en general las áreas verdes. De acuerdo con ^{Vargas y Roldán (2018)}, existen argumentos que sugieren una relación positiva entre los PU, y los APO con la CV de las personas, sin embargo, es necesario estimar su importancia. Según ^{Bertram y Rehdanz (2015)}; ^{Kim y Jin (2018)} y ^{Vargas y Roldán (2018)}, esta cuestión se puede abordar a partir de la valoración económica de atributos ambientales con datos de CV o de BS, a través de una medida de satisfacción con la vida y/o con los servicios. De esta forma, los estudios que consideran la satisfacción personal en función de PU y AV toman en cuenta una función de utilidad indirecta (de máxima satisfacción posible por el consumo de bienes y servicios, cuando se sopesa el ingreso de las personas), estimada con una variable proxy como la CV del consumidor derivada del bienestar que pueden proporcionar los servicios ambientales y sociales de parques urbanos y áreas verdes (^{Ambrey y Fleming, 2011}; ^{Kim y Jin, 2018}; ^{Vargas y Roldán, 2018}). Si se observa lo anterior, el modelo general para esta investigación es:

En este modelo, la utilidad económica se mide a través de una variable proxy, en forma de percepción de calidad de vida (CV_i), la cual se le pregunta a la persona; Y es el ingreso del consumidor trasformado con logaritmos naturales (LN), ya que se asume utilidad marginal decreciente del ingreso² (^{Kim y Jin, 2018}; ^{Vargas y Roldán, 2018}); BA_i es el bien de importancia ambiental en el entorno como variable por valorar (que en esta investigación serán los bienes PU y APO); X _i es un conjunto de variables sociodemográficas (sexo, edad, educación, residencia, estado civil, ocupación) y Z _i son variables de la calle donde reside, como la calidad del recubrimiento, rampas en las banquetas y ambulantaje en el espacio público. Además de la presencia del mercado público y el centro de salud en su entorno y la percepción de seguridad de la calle, cada variable se introduciría en el modelo de manera individual. Para esta investigación, si se sigue a ^{Levinson (2012)} y ^{Bertram y Rehdanz (2015)}, las variables a valorar parque urbano y árboles y plantas de ornato ubicados en el entorno donde reside el habitante de la ciudad se modelan por separado. Las variables PU y APO, al igual que el resto inciden en la utilidad económica en forma de calidad de vida del consumidor, y εi es el error aleatorio. De acuerdo con ^{Ambrey y Fleming (2013)} y ^{Bertram y Rehdanz (2015)}, la ecuación 1 se puede estimar a partir de regresión lineal si la variable dependiente (CV) es cuantitativa o, a través de regresiones del tipo probit ordinal, si la variable dependiente se considera en escala de medición ordinal. ^{Welsch y Ferreira (2014)} y ^{Vargas y Roldán (2018)} fundamentan la elección de la regresión lineal, si se asume que la comparación entre niveles diferentes de CV_i no es ordinal sino cardinal, lo anterior significa que la distancia entre dos niveles de calidad de vida es proporcional a la distancia que hay entre dos de CV_i del consumidor, en consecuencia la estimación se puede realizar por mínimos cuadrados ordinarios. Si se considera a ^{Ambrey y Fleming (2013)}, ^{Kopmann y Rehdanz (2013)} y ^{Bertram y Rehdanz (2015)}, los resultados, al emplear la regresión lineal y los modelos ordinales, tenderán a ser los mismos. Por ello sugieren el modelo de regresión lineal, por la interpretación directa del efecto de los parámetros obtenidos sobre la variable dependiente. Por lo anterior, en este trabajo se emplean estimaciones de regresión lineal, sin embargo, se reportan también los resumidos del probit ordinal. No se debe perder de vista que con la ecuación 1 y de acuerdo con ^{Freeman, Herriges y Kling (2014)} y ^{Vargas y Roldán (2018)}, el valor de un cambio en la presencia/calidad/cantidad del bien ambiental corresponde al cambio que debe ocurrir en el ingreso para que la CV sea constante. Ante la presencia o ausencia de un bien ambiental correspondería un cambio en el ingreso, que deja al individuo en el mismo nivel de utilidad o calidad de vida. Por lo que según este paradigma de utilidad/CV, ^{Frey et al. (2010)}, ^{Ferreira y Moro (2010)}, ^{Levinson (2012)}, ^{Ambrey y Fleming (2013)} y ^{Vargas y Roldán (2018)} señalan que es posible estimar el valor monetario de un cambio en el bien ambiental y obtener una tasa marginal de sustitución entre el ingreso y la presencia/calidad/cantidad del bien/servicio que se quiere valorar monetariamente. En la literatura, a dicho valor monetario se le denomina disposición a pagar (DAP), sin que ésta se entienda necesariamente en el sentido de los modelos de valoración contingente³ (^{Bertram y Rehdanz, 2015}; ^{Ferreira y Moro, 2010}; ^{Kim y Jin, 2018}; ^{Kopmann y Rehdanz, 2013}). ^{Welsch y Ferreira (2014)} se refieren al enfoque de valoración basado en CV de tipo subjetivo como “método de preferencia experimentado”, y señalan también que, a diferencia de la valoración contingente, el método empleado en este trabajo no se basa en opciones sino en asociación estadística entre bienestar o CV subjetiva con indicadores de calidad ambiental. Por consiguiente, la DAP a la que se hace referencia aquí se obtendría a partir de la razón de las derivadas parciales de la calidad de vida considerando un cambio en el bien ambiental (numerador), con respecto a uno en el logaritmo del ingreso en el denominador, con el resto de variables constantes. De esta forma, al relacionar una medida de cambio de un bien con respecto a una en el ingreso (la tasa marginal de sustitución), se obtiene una medida del valor monetario del bien (^{Freeman et al., 2014}; ^{Riera et al., 2016}). Entonces, la DAP se obtendría a partir de la división de los parámetros estimados para el BA y el logaritmo del ingreso tomando en cuenta el valor promedio del ingreso (Y-) o algún cuantil de éste; por lo tanto, su forma básica es:

donde Y- representa el promedio de la variable ingreso, ρ es el parámetro estimado por la presencia, cantidad o calidad de un bien o servicio ambiental y β es el parámetro correspondiente para el logaritmo natural del ingreso. También ^{Levinson (2012)}, ^{Ambrey y Fleming (2013)} y ^{Vargas y Roldán (2018)} estiman este tipo de modelos considerando interacciones entre el ingreso y el bien para determinar si este último se comporta como uno normal o inferior.⁴ En consecuencia, la fórmula 2 cambiaría en el número de términos del numerador y del denominador, en caso de que exista una interacción entre el bien ambiental y el ingreso según el debido proceso de derivación matemática, de acuerdo con lo señalado por ^{Ambrey y Fleming (2013)} y ^{Bertram y Rehdanz (2015)}. Además, para ^{Ferreira y Moro (2010)}, ^{Kopmann y Rehdanz (2013)}, ^{Bertram y Rehdanz (2015)} y ^{Vargas y Roldán (2018)}, en este tipo de modelos puede haber percepciones y características de personalidad no observables, correlacionadas con las variables incluidas en la regresión. Lo anterior puede causar problemas de endogeneidad relacionados con la variable ingreso, como se explicará después. Por lo tanto, para leer los resultados de este trabajo se deben considerar tales limitaciones.

Área de estudio

De acuerdo con ^{Meza y Moncada (2010)} y ^{Fernández-Álvarez (2017)}, Ciudad de México es una megaciudad del siglo XXI, densamente poblada y en donde es necesario que sus habitantes tengan acceso a AV en buenas condiciones, para que mantengan o mejoren su CV (véase Figura 1).

Fuente: Instituto Nacional de Estadística y Geografía (^{INEGI, 2018}).

Figura 1 Mapa de Ciudad de México y sus alcaldías 

De acuerdo con ^{Checa-Artasu (2016)}, ^{Martínez-Soto et al. (2016)} y ^{Morales-Cerdas et al. (2018)}, las AV son un patrimonio natural importante para la sustentabilidad urbana, debido a los servicios ambientales que proporcionan, como captura de carbono y mitigación del efecto “isla de calor”, además de ofrecer beneficios sociales, económicos y médicos. Aunque no es objetivo de esta investigación discutir y analizar aspectos de planificación y problemática de las AV en Ciudad de México, ^{Fernández-Álvarez (2017)} señala que su distribución y accesibilidad se correlacionan con las características demográficas de los sectores socioeconómicos. En la actualidad es muy complicado, por cuestiones teórico-prácticas, señalar un óptimo de superficie verde per cápita para todas las ciudades, debido a la naturaleza multiobjetivo de las AV urbanas y de los PU (servicios recreativos diversos y ambientales múltiples) y también a las diferencias económicas, sociales, geográficas y ambientales (^{Garvin, 2013}; ^{Harnik, 2010}; ^{World Health Organization, 2010}). Sin embargo, es necesario especificar cantidades de superficie verde que guíen al especialista en planeación urbana y/o de AV en temas de acceso, justicia ambiental y de CV, como lo hacen ^{Checa-Artasu (2016)}, ^{Martínez-Soto et al. (2016)} y ^{Fernández-Álvarez (2017)}. Las estadísticas oficiales por alcaldía sobre superficie verde del gobierno de Ciudad de México son de 2009 (véase Tabla 1) (Procuraduría Ambiental y del Ordenamiento Territorial de la Ciudad de México ^{[PAOT], 2018}), por lo que se deben manejar con reserva, y esperar actualizaciones para la segunda década del siglo XXI. Además, tales cifras separan las superficies verdes en suelo urbano de las que tienen estatus de conservación, pero no especifica los criterios para hacerlo, y tampoco el estatus de accesibilidad a ellas (si es libre o restringido), ni hay claridad sobre los aspectos técnicos y legales de lo que es un área verde (^{Checa-Artasu, 2016}).

Aquí se optó por abordar la relación entre CV y AV, a través de la percepción subjetiva de los beneficiarios en relación con la presencia de PU y APO en el entorno donde residen los citadinos, en lugar de manejar datos de superficie verde total, superficie verde por persona o bien de distancia de los hogares a parques urbanos más cercanos. Aunque esto no excluye la posibilidad de que estudios futuros aborden tales perspectivas con datos de superficie y/o distancia de los PU/AV, como lo hacen ^{Ambrey y Fleming (2011}, ²⁰¹³⁾; ^{Bertram y Rehdanz (2015)}; ^{Aoshima, Uchida y Ushimaru (2018)} y ^{Vargas y Roldán (2018)}.

Metodología

Precisiones sobre las variables

Estimación de modelos de regresión

Los modelos se estimaron en pasos, para ajustar por separado las variables PU y APO, acompañadas siempre del LN del ingreso total mensual del hogar, como lo sugieren ^{Ambrey y Fleming (2013)}, ^{Bertram y Rehdanz (2015)}, ^{Aoshima et al. (2018)}, ^{Vargas y Roldán (2018)} y ^{Du et al. (2018)}. Dicho ingreso comprendió aportaciones (sueldos/salarios) de los integrantes del hogar, jubilaciones, becas, apoyos y trasferencias gubernamentales de ayuda social y de otras fuentes, como actividades comerciales y trasferencias/contribuciones de familiares que viven en otro país. En la primera etapa se consideraron solo variables sociodemográficas, en la segunda se incorporaron las del entorno urbano (conjunto Z), y en la tercera se incluyeron interacciones entre el LN del ingreso y las de PU y APO, para determinar la naturaleza normal o inferior de dicho bien. Sin embargo, las interacciones no resultaron significativas, por lo que no se consideraron en los modelos finales. También se probaron variables independientes, como el LN de habitantes por hogar y la presencia de servicios públicos y características del entorno: letrero en calles, puestos de venta semifijos, alumbrado público, botes de basura, de alcantarillado, biblioteca, deportivo, módulo de seguridad y de limpieza. Sin embargo, no fueron significativas ni aportaban a la bondad de ajuste o a la parsimonia del modelo, por lo que no estuvieron en los modelos finales estimados. Lo anterior se realizó con el criterio de información de Aikake (AIC) y el coeficiente de determinación R² para la selección de variables. Todas las estimaciones consideraron errores estándar robustos de clúster (supuso la autocorrelación de las observaciones dentro de las alcaldías), como lo sugieren ^{Kopmann y Rehdanz (2013)}, ^{Scott y Freese (2014)}, ^{Bertram y Rehdanz (2015)} y ^{Vargas y Roldán (2018)}. Lo anterior, se realizó conforme al procedimiento adoptado en investigaciones similares de valoración económica de bienes y servicios ambientales a partir de modelos de CV, BS y satisfacción (^{Ambrey y Fleming, 2013}; ^{Bertram y Rehdanz, 2015}; ^{Kopmann y Rehdanz, 2013}; ^{Vargas y Roldán, 2018}).

Estimación de la disposición a pagar

La DAP se estimó, por separado, para la presencia de PU y APO en el entorno. Para ello se emplearon los resultados obtenidos a partir de regresión lineal, además del probit ordinal. De forma semejante a como lo realizan ^{Smyth et al. (2008}), ^{Ambrey y Fleming (2013)}, ^{Bertram y Rehdanz (2015)}, ^{Vargas y Roldán (2018)} y ^{Aoshima et al. (2018)}, se estimaron las DAP para PU y APO, mediante las fórmulas 3 y 4 respectivamente:

En estas fórmulas, en el numerador se introduce el parámetro del bien por valorar (PU y APO), y en el denominador el parámetro obtenido para la variable LN del ingreso mensual del hogar. De acuerdo con ^{Ambrey y Fleming (2011}, ²⁰¹³⁾, los valores de DAP estimados se dividieron entre el número promedio de habitantes por hogar a partir de los resultados obtenidos en la misma encuesta. Para las estimaciones se consideró el promedio del ingreso mensual del hogar, y éstas se extendieron a los deciles primero, tercero, séptimo y noveno del ingreso, para determinar su comportamiento ante altos y bajos ingresos mensuales.

En los análisis de regresión y estimaciones se utilizó Stata, versión 15, de acuerdo con lo recomendado por ^{Cameron y Trivedi (2013)}, ^{Hilbe (2014)} y ^{Scott y Freese (2014)}.

Resultados y discusión

Resultados de modelos estimados

En los modelos estimados se introdujo la edad junto con su forma cuadrática, con base en los argumentos de ^{Blanchflower y Oswald (2008)}, ^{Dolan y White (2007)}, ^{Carstensen et al. (2011)}, ^{Ambrey y Fleming (2013)} y ^{Graham (2016)}, quienes señalan que las personas empiezan a incrementar su satisfacción por la vida cuando tienen entre 40 y 50 años, en comparación con la caída que experimentaban antes. Esto se describe como un comportamiento gráfico en U de la CV respecto a la edad (^{Ambrey y Fleming, 2013}), que ha ocurrido en Estados Unidos, Canadá y Reino Unido (^{Deaton, 2008}). El ^{INEGI (2012)} reporta un comportamiento en U en la relación de satisfacción con la vida y la edad; encontró el nivel más bajo en el grupo de 45 a 59 años del mexicano promedio. Con los resultados de este trabajo, los parámetros y signos obtenidos para la edad y edad al cuadrado en los modelos estimados de la población estudiada, el nivel más bajo de CV se presenta a los 47 años, lo que coincide con el comportamiento descrito por el ^{INEGI (2012)}, ^{Deaton (2008)} y ^{Graham (2016)}. Se optó por introducir a la educación como variable nominal, debido a que los años cursados en cada nivel son diferentes: seis de primaria, tres de secundaria, tres de bachillerato y cuatro o más en el caso de licenciatura, maestría, doctorado y posdoctorado. En las Tablas 5 y 6 se observa una significancia estadística positiva entre los grados de secundaria en adelante, con respecto al grupo de la población que no tiene estudios, en coincidencia con los resultados encontrados por ^{Kopmann y Rehdanz (2013)}, ^{Bertram y Rehdanz (2015)} y ^{Vargas y Roldán (2018)}. Ninguna de las variables categóricas resultó significativa en el rubro de ocupación, no obstante, el signo positivo correspondiente a la población trabajadora coincide con el comportamiento encontrado por ^{Vargas y Roldán (2018)} y ^{Kim y Jin (2018}) respecto a que las personas que trabajan se sentirán más satisfechas o tendrán mejor CV que las desempleadas y/o que buscan trabajo. La variable LN del ingreso en los modelos ajustados tuvo un signo positivo, como se reportó en estudios de CV y valoración de AV y de otros servicios ambientales (^{Ambrey y Fleming, 2013}; ^{Bertram y Rehdanz, 2015}; ^{Kim y Jin, 2018}; ^{Smyth et al., 2008}; ^{Vargas y Roldán, 2018}), y fue estadísticamente significativa en todos los modelos, al menos con un α ≤ 0.01. En cuanto a la de estado civil, los resultados estadísticamente significativos en ambos modelos (regresión lineal y probit ordinal), señalan que las personas casadas/unión libre (al menos con un α de 0.05) y solteras (al menos con un α de 0.01) tendieron a reportar mayor CV que las divorciadas o separadas. En el caso de viudez, el parámetro no fue significativo.

Sobre el efecto de un PU en el entorno, los resultados en ambas regresiones (véanse Tablas 5 y 6) sustentan, al menos con α=0.05, la hipótesis de que la existencia de uno incide positivamente en la CV de los habitantes de Ciudad de México. La presencia de APO también tiene un efecto positivo en la calidad de vida, al menos con α=0.05, en los modelos de regresión (véanse Tablas 5 y 6). Además, una prueba de ji-cuadrada para analizar la dependencia entre las variables PU y presencia de APO sustenta la hipótesis de que la existencia de la primera en el entorno está asociada con la presencia de la segunda en la calle y el entorno (Ji²= 73.8468, al menos con un α≤0.01).

Provisión de otros servicios urbanos en el entorno

Los resultados de las regresiones lineal y probit ordinal (véanse Tablas 5 y 6) establecen que la presencia de rampas para sillas de ruedas, adoquín y pavimento en las calles (recubrimiento) repercute de manera positiva y significativa en la CV de las personas, con respecto a quienes viven en calles de tierra. La calidad de vida se incrementa cuando mejora el recubrimiento (de muy deficiente a excelente).

El modelo de regresión lineal (véase Tabla 5) permite observar un efecto estadísticamente significativo y positivo con respecto a la categoría base (calle muy insegura) en la clasificación de segura y altamente segura (al menos con un α≤0.1). En los resultados del modelo probit (véase Tabla 6), el efecto positivo se extiende a la categoría de calle más o menos segura, con respecto a la de base.

Los habitantes de Coyoacán, Iztapalapa, Magdalena Contreras, Álvaro Obregón, Tláhuac, Xochimilco, Benito Juárez y Miguel Hidalgo tienden a reportar más CV en términos estadísticamente significativos (al menos un α= 0.05), con respecto a Milpa Alta (base) (véase Tabla 5). Si bien el objetivo de este trabajo no es analizar la calidad de vida por alcaldía, este comportamiento coincide de manera parcial con los resultados del índice de bienestar social del gobierno de Ciudad de México (^{CEDSCM, 2016}), que es mayor en las ocho alcaldías mencionadas, en comparación con el de Milpa Alta que ocupa el último lugar en dicho índice, según el ^{CEDSCM (2016)}. Caso contrario al de Iztacalco y Cuauhtémoc, cuyos habitantes tienden a manifestar CV inferior a la población de Milpa Alta, al menos un α= 0.05. El resto de las variables dicotómicas no fueron significativas, al menos con α= 0.1. Sin embargo, la comparación de resultados de esta investigación con el índice de bienestar del ^{CEDSCM (2016)} se debe acotar, y considerar que comprende variables indicadoras tanto de calidad del entorno (AV, banquetas y alumbrado) como de tiempo libre (acceso a servicios de salud, internet, a alimentación, calidad de vivienda y felicidad) y no incluye otras, como seguridad pública. En esta investigación, los resultados no deben interpretarse en términos de superficie de área verde por alcaldía, ya que esta variable no formó parte de los modelos estimados; tampoco se contempló la de accesibilidad a AV a partir de los hogares. Sin embargo, Milpa Alta tiene la mayor superficie total de áreas verdes en Ciudad de México (véase Tabla 1), pero parte de ella está sujeta a conservación. Los resultados de las variables dicotómicas por alcaldía son semejantes entre modelos de regresión lineal y probit estimados para PU y APO aunque, por razones de espacio, los arrojados por el probit no se reportan. En estudios futuros sobre el tema, en los modelos de CV se podrían incorporar datos de superficie verde e incluso considerar medidas de distancia a los PU y alguna variable de accesibilidad a las AV.

Obtención de la disposición a pagar

En la Tabla 7 aparecen los resultados de la DAP, estimados con los resultados de regresión lineal y probit ordinal. Los cálculos se dividieron por el promedio de habitantes en hogar, el cual fue de 4, y se encontró lo siguiente:

Los resultados de la DAP se deben acotar por el posible efecto de endogeneidad sobre la variable ingreso, por lo que el valor del parámetro ingreso pudiera estar subestimado y, por tanto, producir una sobrevaloración de la DAP, como lo señalan ^{Luechinger (2009)}; ^{Ambrey y Fleming (2014)} y ^{Vargas y Roldán (2018)}. Sobre este tema, según ^{Luechinger (2009)}, una endogeneidad grave y problemas de variables omitidas pueden incidir en la estimación del efecto del ingreso en la satisfacción con la vida. De acuerdo con ^{Vargas y Roldán (2018)}, puede ocurrir que las personas más productivas reciban más ingreso y, por lo tanto, que estén más satisfechas. Según ^{Clark, Frijters y Shields (2008)} y ^{Luechinger (2009)}, la evidencia apunta a que las personas más felices ganan más, y que hay variables difíciles de controlar en el análisis estadístico, que influyen en la felicidad, la satisfacción y la CV, entre ellas están los factores laborales, como horas de trabajo, y los psicológicos, como estrés y salud física. Los resultados de esta investigación deben leerse tomando en cuenta esta restricción (las variables difíciles de controlar), y los valores estimados de la DAP se deben considerar solo como indicativos del monetario de los PU y los APO en Ciudad de México, como sugieren ^{Vargas y Roldán (2018)} para los valores estimados en parques urbanos de Colombia. Además, estos valores pudieran ser complementarios a los monetarios, que se estiman a través de precios hedónicos de viviendas.

De acuerdo con ^{Ambrey y Fleming (2013)} y ^{Vargas y Roldán (2018)}, para valorar una AV y otros servicios públicos citadinos, los resultados provenientes del bienestar subjetivo y la calidad de vida quizá sean de naturaleza hedónica. Por ello, esos resultados pueden ser estimaciones complementarias a las de los precios hedónicos del mercado de viviendas (^{Bertram y Rehdanz, 2015}; ^{Kim y Jin, 2018}, ^{Kim et al., 2019}). Según ^{Vargas y Roldán (2018)}, en el modelo de precios hedónicos, la utilidad/satisfacción generada por la presencia de PU y APO constituiría parte del precio de la vivienda. En las regresiones estimadas aquí, al no incluir dicho precio, podrían indicar que los valores monetarios calculados son el de PU y de APO, que no es capturado por el precio de la vivienda. En consecuencia, la complementariedad de los resultados, de acuerdo con ^{Bertram y Rehdanz (2015)} y ^{Kim et al. (2019)}, dependerá de la existencia de equilibrio (el punto donde la oferta se iguala a la demanda)⁵ en el mercado de vivienda. En caso de no existir un equilibrio de mercado, las estimaciones obtenidas a través de funciones de BS/CV podrían representar el valor del bien o servicio ambiental que se está evaluando. Por tanto, es necesario abordar la valoración económica de servicios ambientales con el método de precios hedónicos, y analizar el equilibrio de mercado de viviendas, que permitirá conocer cómo y en cuánto el valor monetario de los PU o de los APO constituye parte del precio de una vivienda. Este es un tema pendiente, empero, los resultados de esta investigación permiten señalar que la presencia de parques urbanos y árboles y plantas de ornato en el entorno incide positivamente en la CV de los citadinos y que los valores monetarios de la DAP estimados son un indicador posible del valor económico de dichos bienes ambientales. Sin embargo, los resultados de este trabajo están acotados por los posibles problemas de endogeneidad descritos. Además, se debe considerar que las AV también causan males sociales como congestión del espacio, inseguridad y delincuencia, además de la falta de mantenimiento, como lo señalan ^{Kim y Jin (2018)} y ^{Vargas y Roldán (2018)}. En investigaciones futuras se podrían abordar cuestiones de valoración de PU y AV a través de precios hedónicos y análisis de mercados de vivienda, y estimar valores monetarios (DAP) tomando en cuenta distancias de los parques urbanos a los hogares, superficie de áreas verdes y accesibilidad a ellas e incluso densidades poblacionales de los lugares de residencia. Estos elementos podrían servir para tratar problemas de planificación urbana relacionados con dotación de AV y de congestión. En este aspecto, los métodos de BS/CV pueden ser útiles porque complementarían los hedónicos de valoración (^{Ambrey y Fleming, 2011}; ^{Aoshima et al., 2018}; ^{Kim y Jin, 2018}). No obstante, los valores obtenidos aquí son un indicador de la importancia de las áreas verdes (parques urbanos y árboles y plantas de ornato) en la vida de los citadinos.

Conclusiones

La estimación de los modelos de regresión lineal y probit ordinal de CV en función de variables sociodemográficas, servicios públicos y percepción de inseguridad permite concluir que los servicios de tipo ambiental, presencia de PU y de APO en el entorno inciden significativa y positivamente en la calidad de vida de los citadinos.

El signo de los parámetros estimados para la edad coincide con el comportamiento en U, encontrado en investigaciones sobre BS y CV para otros países. El valor mínimo de CV, ubicado en los 47 años, coincide con lo reportado sobre BS por el ^{INEGI (2012)}. El signo y la significancia estadística del ingreso permiten establecer que a mayores ingresos habrá mejor calidad de vida.

Para la presencia de parques en el entorno de los lugares de residencia de los citadinos, la DAP estimada con los parámetros por regresión lineal y probit ordinal, respectivamente, fluctuó de 277.9 a 284.8 pesos, en ingresos del primer decil, de 719.1 a 736.8, en media del ingreso, y en el noveno decil de 1 945.5 a 1 993.5. Este cálculo se hizo por habitante del hogar, como en investigaciones similares. En este caso, los valores mayores de la disposición a pagar se obtienen a partir de los parámetros del probit ordinal.

En cuanto a la presencia de APO en el entorno, la DAP estimada por probit ordinal y regresión lineal, respectivamente, para la presencia de APO en el entorno, fluctuó de 648.3 a 666.1 pesos para el primer decil; para la media, de 1 677.4 de 1 723.3 y, en el noveno decil, de 4 538.2 a 4 662.5. Aquí los valores mayores de la DAP se obtuvieron a partir de los parámetros de la regresión lineal.

Cuando existen árboles y plantas de ornato en el entorno, la DAP tiende a ser mayor que para la de parques urbanos.

Los valores de la DAP estimados de bienes ambientales del entorno donde residen los citadinos pueden ser complementarios a los de PU y APO, calculados por el método de precios hedónicos, al formar parte o capitalizarse en el precio de las viviendas. Sin embargo, hacen falta investigaciones nuevas sobre el tema, para determinar la existencia de equilibrios en mercados de vivienda, de no existir dicho equilibrio los valores estimados por los métodos de CV/BS podrían representar el valor de los bienes ambientales.

Los resultados de esta investigación están delimitados por la posible endogeneidad presente en el modelo, debido a que las personas con mayor CV tienden a ser las más productivas y, por lo tanto, a ganar más. Lo que provocaría una sobreestimación de los valores de la DAP. En estudios futuros se tendrán que incorporar variables de tipo laboral (horas de trabajo) y factores psicológicos, como estrés.

Los resultados obtenidos también están limitados a considerar solo la presencia de servicios ambientales y no la superficie, distancia, ni la accesibilidad a dichos servicios. Lo anterior es fundamental para abordar problemas relacionados con presencia de males ambientales derivados de los PU y las AV como congestión, inseguridad y falta de mantenimiento, lo que repercutiría en los valores de la DAP estimados. No obstante las limitantes del estudio, los resultados y los valores monetarios calculados aquí son un indicador de la importancia que tienen los parques urbanos y los árboles y plantas de ornato en la calidad de vida de los habitantes de la ciudad.