Cano-Ibáñez N, Martínez-Galiano JM, Luque-Fernández MA, Martín-Peláez S, Bueno-Cavanillas A, Delgado-Rodríguez M. Maternal Dietary Patterns during Pregnancy and Their Association with Gestational Weight Gain and Nutrient Adequacy [Internet]. International Journal of Environmental Research and Public Health 2020, Vol. 17, Page 7908 2020;17(21):7908. Publisher's VersionAbstract
Several epidemiologic studies have shown an association between Gestational Weight Gain (GWG) and offspring complications. The GWG is directly linked to maternal dietary intake and women’s nutritional status during pregnancy. The aim of this study was (1) to assess, in a sample of Spanish pregnant women, the association between maternal dietary patterns and GWG and (2) to assess maternal dietary patterns and nutrient adequate intake according to GWG. A retrospective study was conducted in a sample of 503 adult pregnant women in five hospitals in Eastern Andalusia (Spain). Data on demographic characteristics, anthropometric values, and dietary intake were collected from clinical records by trained midwives. Usual food intake was gathered through a validated Food Frequency Questionnaire (FFQ), and dietary patterns were obtained by principal component analysis. Nutrient adequacy was defined according to European dietary intake recommendations for pregnant women. Regression models adjusted by confounding factors were constructed to study the association between maternal dietary pattern and GWG, and maternal dietary patterns and nutritional adequacy. A negative association was found between GWG and the Mediterranean dietary pattern (crude β = −0.06, 95% CI: −0.11, −0.04). Independent of maternal dietary pattern, nutrient adequacy of dietary fiber, vitamin B9, D, E, and iodine was related to a Mediterranean dietary pattern (p \textless 0.05). A Mediterranean dietary pattern is related to lower GWG and better nutrient adequacy. The promotion of healthy dietary behavior consistent with the general advice promoted by the Mediterranean Diet (based on legumes, vegetables, nuts, olive oil, and whole cereals) will offer healthful, sustainable, and practical strategies to control GWG and ensure adequate nutrient intake during pregnancy.
Lee SF, Luque-Fernandez MA, Chen YH, Catalano PJ, Chiang CL, Wan EY-F, Wong IC-K, Chen MH, Ng AK. Doxorubicin and subsequent risk of cardiovascular diseases among survivors of diffuse large B-cell lymphoma in Hong Kong [Internet]. Blood Advances 2020;4(20):5107-5117. Publisher's VersionAbstract
Evidence regarding the dose-related impact of doxorubicin on subsequent cardiovascular diseases (CVDs) in Asian patients with diffuse large B-cell lymphoma (DLBCL) without preexisting CVDs is lacking. From a territory-wide electronic database in Hong Kong, we identified adults who were diagnosed with DLBCL and treated with chemotherapy between 2000 and 2018. We evaluated the patients for incident CVDs (including ischemic heart disease, heart failure, and cardiomyopathy). We evaluated the cause-specific cumulative incidence (csCI) of CVD with levels of doxorubicin exposure by using flexible parametric competing risk analysis and adjusting for demographics, comorbidities, therapeutic exposure, cardiovascular risk factors, and lifestyle factors. Controls were age- and sex-matched to DLBCL patients. We analyzed 2600 patients and 13 000 controls. The adjusted cause-specific hazard ratio (HR) for CVD in patients treated with ġt;500 mg doxorubicin compared with non-doxorubicin regimens was 2.65 (95\% confidence interval [CI], 1.23-5.74; P = .013). The 5-, 10-, and 15-year csCIs were 8.2\%, 11.3\%, and 12.8\% in patients vs 3.1\%, 4.4\%, and 5.2\% in controls, respectively. Hypertension (HR, 6.20; 95\% CI, 0.79-48.44; P = .082) and use of aspirin/angiotensin-converting enzyme inhibitor/beta-blocker at baseline (HR, 2.13-4.63; P \< .001 to .002) might confer a higher risk of subsequent CVDs. In this Hong Kong population-based study, doxorubicin exposure (absolute dose ġt;500 mg), together with hypertension or baseline use of medication for cardiovascular risk factors, was found to be associated with an increase in csCIs of CVDs. Tailoring therapeutic strategies to underlying CVD risk factors and risk-adapted monitoring and follow-up of susceptible DLBCL patients are advisable.
Nuñez O, Barranco MR, Fernández-Navarro P, Sanchez DR, Luque-Fernandez MA, Santamaría MP, Sánchez M-J. Deprivation gap in colorectal cancer survival attributable to stage at diagnosis: A population-based study in Spain [Internet]. Cancer Epidemiology 2020;68:101794. Publisher's VersionAbstract
Objective Socioeconomic inequalities in colorectal cancer (CRC) survival are a major concern of the Spanish public health system. If these inequalities were mainly due to differences in stage at diagnosis, population-based screening programs might reduce them substantially. We aimed to determine to what extent adverse stage distribution contributed to survival inequalities in a Spanish region before the implementation of a CRC screening program. Methods We analyzed data from a population-based cohort study that included all patients living in a region of southern Spain with CRC diagnosed between 2004 and 2013. The European Deprivation Index was used to assign each patient a socioeconomic level based on their area of residence. The role of tumor stage in survival disparities between socioeconomic groups was assessed using a causal mediation analysis. Results A total of 2802 men and 1957 women were included in the study. For men, the adjusted difference in deaths between the most deprived and the most affluent areas was 131 deaths per 1000 person-years by the first year after diagnosis. Of these deaths, 42 (per 1000 person-years) were attributable to differences in stage at diagnosis. No socioeconomic disparities in survival were detected among female patients. Conclusions In this study, we mainly detected socioeconomic disparities in short term survival of male patients. More than two thirds of these inequalities could not be attributed to differences in stage at diagnosis. Our results suggest that in addition to a screening program, other public health interventions are necessary to reduce the deprivation gap in survival.
Luque-Fernandez MA, Redondo-Sanchez D, RODRIGUEZ-BARRANCO MIGUEL, Yoe-Ling C-C, Salamanca-Fernandez E, Nunez O, Fernandez-Navarro P, Pollan M, Sanchez M-J. Socioeconomic Inequalities in Colorectal Cancer Survival in Southern Spain: A Multilevel Population-Based Cohort Study [Internet]. Clinical Epidemiology 2020;Volume 12:797–806. Publisher's VersionAbstract
Background: Colorectal cancer (CRC) is the most frequently diagnosed cancer in Spain. Socioeconomic inequalities in cancer survival are not documented in Spain. We aim to study the association of socioeconomic inequalities with overall mortality and survival among CRC patients in southern Spain.
Methods: We conducted a multilevel population-based cohort study, including CRC cases for the period 2011– 2013. The study time-to-event outcome was death, and the primary exposure was CRC patients’ socioeconomic status assessed by the Spanish deprivation index at the census tract level. We used a mixed-effects flexible hazard model, including census tract as a random intercept, to derive overall survival estimates by deprivation.
Results: Among 3589 CRC patients and 12,148 person-years at risk (pyr), 964 patients died before the end of the follow-up. Mortality by deprivation showed the highest mortality rate for the most deprived group (96.2 per 1000 pyr, 95% CI: 84.0– 110.2). After adjusting for sex, age, cancer stage, and the area of residence, the most deprived had a 60% higher excess mortality risk than the less deprived group (excess mortality risk ratio: 1.6, 95% CI: 1.1– 2.3).
Conclusions: We found a consistent association between deprivation and CRC excess mortality and survival. The reasons behind these inequalities need further investigation in order to improve equality cancer outcomes in all social groups.

Keywords: socioeconomic inequalities, colorectal cancer, survival, population-based epidemiology, epidemiological methods, multilevel
Luque-Fernández MÁ, Redondo-Sánchez D, Fernández P, Salamanca-Fernández E, Marcos-Gragera R, Guevara M, Carulla Aresté M, Jiménez R, Núñez O, Sabater C, López de Munain Marqués A, Chirlaque MD, Mateos A, Rodríguez-Barranco M, Espín Balbino J, Pollán M, Sánchez M-J. Protocolo del Estudio Poblacional Multinivel de las Desigualdades Socioeconómicas en la Distribución Geográfica de la Incidencia, la Mortalidad y la Supervivencia Neta del Cáncer en España: Estudio DESOCANES. Gaceta Sanitaria 2020;Abstract
La incidencia y la mortalidad brindan información sobre la carga de la morbilidad del cáncer y los años potenciales de vida perdidos debido al cáncer. Se ha desarrollado el Índice de Privación (IP) como una medida estandarizada para medir la privación socioeconómica en España a nivel de sección censal. Además, se puede combinar la información del IP con variables ecológicas poblacionales y los datos de los Estudios Europeos de Alta Resolución en Cáncer. El objetivo de este estudio es caracterizar las desigualdades socioeconómicas en la incidencia, el exceso de mortalidad, la mortalidad prematura y la supervivencia neta para tres de los cánceres más incidentes (pulmón, colon-recto y mama) en España usando el IP. Este estudio nacional multinivel de base poblacional evaluará el impacto de las desigualdades socioeconómicas. Se usarán el análisis espacial, la modelización multinivel, la supervivencia neta y la evaluación del impacto económico. Los resultados serán útiles para el apoyo a la toma de decisiones y la planificación y la gestión de intervenciones en salud pública destinadas a reducir el impacto de las desigualdades socioeconómicas en el diagnóstico y el pronóstico de los pacientes de cáncer en España. ER
Pazos M, Biarnes M, Blasco-Alberto A, Dyrda A, Luque-Fernandez MA, Gómez A, Mora C, Milla E, Muniesa JM, Anton A, Diaz-Aleman VT. SD-OCT peripapillary nerve fibre layer and ganglion cell complex parameters in glaucoma: principal component analysis [Internet]. British Journal of Ophthalmology 2020; Publisher's VersionAbstract
{Background/aims To identify objective glaucoma-related structural features based on peripapillary (p) and macular (m) spectral domain optical coherence tomography (SD-OCT) parameters and assess their discriminative ability between healthy and glaucoma patients.Methods Two hundred and sixty eyes (91 controls and 169 glaucoma) were included in this prospective study. After a complete examination, all participants underwent the posterior pole and the peripapillary retinal nerve fibre layer (pRNFL) protocols of the Spectralis SD-OCT. Principal component analysis (PCA), a data reduction method, was applied to identify and characterise the main information provided by the ganglion cell complex (GCC). The discriminative ability between healthy and glaucomatous eyes of the first principal components (PCs) was compared with that of conventional SD-OCT parameters (pRNFL, macular RNFL (mRNFL), macular ganglion cell layer (mGCL)and macular inner plexiform layer (mIPL)) using 10-fold cross-validated areas under the curve (AUC).Results The first PC explained 58% of the total information contained in the GCC and the pRNFL parameters and was the result of a general combination of almost all variables studied (diffuse distribution). Other PCs were driven mainly by pRNFL and mRNFL measurements. PCs and pRNFL had similar AUC (0.95 vs 0.96
Lee SF, Luk H, Wong A, Ng CK, Wong FCS, Luque-Fernandez MA. Prediction model for short-term mortality after palliative radiotherapy for patients having advanced cancer: a cohort study from routine electronic medical data [Internet]. Scientific Reports 2020;10(1):5779. Publisher's VersionAbstract
{We developed a predictive score system for 30-day mortality after palliative radiotherapy by using predictors from routine electronic medical record. Patients with metastatic cancer receiving first course palliative radiotherapy from 1 July, 2007 to 31 December, 2017 were identified. 30-day mortality odds ratios and probabilities of the death predictive score were obtained using multivariable logistic regression model. Overall, 5,795 patients participated. Median follow-up was 39.6 months (range, 24.5–69.3) for all surviving patients. 5,290 patients died over a median 110 days, of whom 995 (17.2%) died within 30 days of radiotherapy commencement. The most important mortality predictors were primary lung cancer (odds ratio: 1.73, 95% confidence interval: 1.47–2.04) and log peripheral blood neutrophil lymphocyte ratio (odds ratio: 1.71, 95% confidence interval: 1.52–1.92). The developed predictive scoring system had 10 predictor variables and 20 points. The cross-validated area under curve was 0.81 (95% confidence interval: 0.79–0.82). The calibration suggested a reasonably good fit for the model (likelihood-ratio statistic: 2.81
Luque-Fernandez MA, Gonçalves K, Salamanca-Fernández E, Redondo-Sanchez D, Lee SF, Rodríguez-Barranco M, Carmona-García MC, Marcos-Gragera R, Sánchez M-J. Multimorbidity and short-term overall mortality among colorectal cancer patients in Spain: A population-based cohort study [Internet]. European Journal of Cancer 2020;129:4 - 14. Publisher's VersionAbstract
Background Numerous studies have analysed the effect of comorbidity on cancer outcomes, but evidence on the association between multimorbidity and short-term mortality among colorectal cancer patients is limited. We aimed to assess this association and the most frequent patterns of multimorbidity associated with a higher short-term mortality risk among colorectal cancer patients in Spain. Methods Data were obtained from two Spanish population-based cancer registries and electronic health records. We estimated the unadjusted cumulative incidence of death by comorbidity status at 6 months and 1 year. We used a flexible parametric model to derive the excess mortality hazard ratios (HRs) for multimorbidity after adjusting for sex, age at diagnosis, cancer stage and treatment. We estimated the adjusted cumulative incidence of death by comorbidity status and identified multimorbidity patterns. Results Among the study participants, 1,048 cases were diagnosed with cancers of the colon and rectum, 2 cases with cancer of the anus with overlapping sites of the rectum and 11 cases with anal adenocarcinomas but treated as colorectal cancer patients. Among 1,061 colorectal cancer patients, 171 (16.2%) died before 6 months, 246 (23.3%) died before the 1-year follow-up, and 324 (30.5%) had multimorbidity. Patients with multimorbidity had two times higher mortality risk than those without comorbidities at 6 months (adjusted HR: 2.04; 95% confidence interval [CI]: 1.30–3.20, p = 0.002). The most frequent multimorbidity pattern was congestive heart failure + diabetes. However, patients with rheumatologic disease + diabetes had two times higher 1-year mortality risk than those without comorbidities (HR: 2.23; 95% CI: 1.23–4.07, p = 0.008). Conclusions Multimorbidity was a strong independent predictor of short-term mortality at 6 months and 1 year among the colorectal cancer patients in Spain. The identified multimorbidity pattern was consistent. Our findings might help identify patients at a higher risk for poor cancer and treatment outcomes.
Ma D, Luque-Fernandez MA, Bogdanet D, Desoye G, Dunne F, Halperin JA. Plasma Glycated CD59 Predicts Early Gestational Diabetes and Large for Gestational Age Newborns. The Journal of Clinical Endocrinology & Metabolism 2020;Abstract
Gestational diabetes mellitus (GDM) diagnosed in early pregnancy is a health care challenge because it increases the risk of adverse outcomes. Plasma glycated CD59 (pGCD59) is an emerging biomarker for diabetes and GDM. The aim of this study was to assess the performance of pGCD59 as a biomarker of early GDM and its association with delivering a large for gestational age (LGA) infant.Study: Blood levels of pGCD59Â were measured in samples from 693 obese women (BMI >29) undergoing a 75g, 2-hour oral glucose tolerance test (OGTT) at <20 weeks gestation in the DALI study: the main analyses included 486 subjects who had normal glucose tolerance throughout the pregnancy, 207 who met criteria for GDM at <20 weeks and 77 diagnosed with GDM at pregnancy week 24-28.Reference tests: 75g, 2-hour OGTT adjudicated based on IADPSG criteria.Index test: pGCD59 ELISA.Primary aim: To assess the performance of pGCD59 to identify women with GDM in early pregnancy (GDM<20).Secondary aim: To assess the association of pGCD59 with LGA and potentially others adverse neonatal outcomes linked to GDM. Mean pGCD59 levels were significantly higher (p-value <0.001) in women with GDM<20 (3.9 ±1.1 SPU) than in those without (2.7 ±0.7 SPU). pGCD59 accurately identified GDM in early pregnancy with an AUC ROC curves of 0.86 (95% CI: 0.83-0.90). One unit increase in maternal pGCD59 level was associated with 36% increased odds of delivering a large for gestational age (LGA) infant (odds ratio for LGA vs non-LGA infant: 1.4; 95%CI: 1.1-1.8; p-value: 0.016). Our results indicate that pGCD59 is a simple and accurate biomarker for detection of GDM in early pregnancy and risk assessment of LGA.
Fowler H, Belot A, Ellis L, Maringe C, Luque-Fernandez MA, Njagi EN, Navani N, Sarfati D, Rachet B. Comorbidity prevalence among cancer patients: a population-based cohort study of four cancers [Internet]. BMC Cancer 2020;20(1):2. Publisher's VersionAbstract
The presence of comorbidity affects the care of cancer patients, many of whom are living with multiple comorbidities. The prevalence of cancer comorbidity, beyond summary metrics, is not well known. This study aims to estimate the prevalence of comorbid conditions among cancer patients in England, and describe the association between cancer comorbidity and socio-economic position, using population-based electronic health records.
Luque-Fernandez MA, Redondo-Sanchez D, Lee SF, Rodriguez-Barranco M, Carmona-Garcia M, Marcos-Gragera R, Sanchez MJ. Multimorbidity by Patient and Tumor Factors and Time-to-Surgery Among Colorectal Cancer Patients in Spain : A Population-Based Study. Clinical Epidemiology 2020;2020:12:31-40.Abstract

Background: Cancer treatment and outcomes can be influenced by tumor characteristics, patient overall health status, and comorbidities. While previous studies have analyzed the influence of comorbidity on cancer outcomes, limited information is available regarding factors associated with the increased prevalence of comorbidities and multimorbidity among patients with colorectal cancer in Spain.

Patients and Methods: This cross-sectional study obtained data from all colorectal cancer cases diagnosed in two Spanish provinces in 2011 from two population-based cancer registries and electronic health records. We calculated the prevalence of comorbidities according to patient and tumor factors, identified factors associated with an increased prevalence of comorbidity and multimorbidity, analyzed the association between comorbidities and time-to-surgery, and developed an interactive web application (

Results: The most common comorbidities were diabetes (23.6%), chronic obstructive pulmonary disease (17.2%), and congestive heart failure (14.5%). Among all comorbidities, 52% of patients were diagnosed at more advanced stages (stage III/IV). Patients with advanced age, restricted performance status or who were disabled, obese, and smokers had a higher prevalence of multimorbidity. Patients with multimorbidity had a longer time-to-surgery than those without comorbidity (17 days, 95% confidence interval: 3– 29 days).

Conclusion: We identified a consistent pattern of factors associated with a higher prevalence of comorbidities and multimorbidity at diagnosis and an increased time-to-surgery among patients with colorectal cancer with multimorbidity in Spain. This pattern may provide insights for further etiological and preventive research and help to identify patients at a higher risk for poorer cancer outcomes and suboptimal treatment.

León-Gómez BB, Gotsens M, Mari-Dell'Olmo M, Domınguez-Berjón MF, Luque-Fernandez MÁ, Martin U, Rodrıguez-Sanz M, Pérez G. Bayesian smoothed small-areas analysis of urban inequalities in fertility across 1999–2013 [Internet]. Fertility Research and Practice 2019;5(1):17. Publisher's VersionAbstract
Since the 2008 economic crisis in Spain, overall fertility has continued to decrease, while urban inequalities have increased. There is a general lack of studies of fertility patterns in small-areas of Spanish cities. We explored the effects of the economic crisis on fertility during three time periods in urban settings in Spain.
Bellizzi S, Nivoli A, Salaris P, Ronzoni AR, Pichierri G, Palestra F, Wazwaz O, Luque-Fernandez MA. Sexual violence and eclampsia: analysis of data from Demographic and Health Surveys from seven low- and middle-income countries. Journal of Global Health 2019;9(2):020434.Abstract
BACKGROUND: Scientific literature has provided clear evidence of the profound impact of sexual violence on women's health, such as somatic disorders and mental adverse outcomes. However, consequences related to obstetric complications are not yet completely clarified. This study aimed to assess the association of lifetime exposure to intimate partner sexual violence with eclampsia. METHODS: We considered all the seven Demographic and Health Surveys (DHS) that included data on sexual violence and on signs and symptoms suggestive of eclampsia for women of reproductive age (15-49 years). We computed unadjusted and adjusted odds ratios (OR) to evaluate the risk of suggestive eclampsia by ever subjected to sexual violence. A sensitivity analysis was conducted restricting the study population to women who had their last live birth over the 12 months before the interview. RESULTS: Self-reported experience of sexual violence ranged from 3.7% in Mali to 9.2% in India while prevalence of women reporting signs and symptoms compatible with eclampsia ranged from 14.3% in Afghanistan to 0.7% in the Philippines. Reported sexual violence was associated with a 2-fold increased odd of signs and symptoms suggestive of eclampsia in the pooled analysis. The sensitivity analysis confirmed the strength of the association between sexual violence and eclampsia in Afghanistan and in India. CONCLUSIONS: Women and girls in low-and-middle-income countries are at high risk of sexual violence, which may represent a risk factor for hypertensive obstetric complication. Accurate counseling by health care providers during antenatal care consultations may represent an important opportunity to prevent adverse outcomes during pregnancy.
Hao M, Luque-Fernandez MA, Lopez D, Cote K, Newfield J, Connors M, Vaidya A. Benign Adrenocortical Tumors and the Detection of Nonadrenal Neoplasia [Internet]. International Journal of Endocrinology 2019;2019:9. Publisher's VersionAbstract
Context. Patients with adrenocortical tumors have been frequently observed to have nonadrenal neoplasia. Objective. To investigate whether patients with benign adrenocortical tumors have a higher likelihood of having nonadrenal neoplasia detected. Design and Participants. Case-control study of patients with benign adrenocortical tumors (cases; n = 400) and normal adrenal glands (controls; n = 400), who underwent repeated abdominal cross-sectional imaging. Main Outcomes. Primary analyses: association between case-control status and benign abdominal neoplasia detected via cross-sectional imaging. Secondary analyses: association between case-control status and tumors detected via other imaging modalities. Results. The mean interval of abdominal imaging was 4.7 (SD = 3.8) years for cases and 5.9 (4.8) years for controls. Cases were more likely to have detected intraductal papillary mucinous neoplasms (IPMNs) of the pancreas (8.5% vs. 4.5%, adjusted OR = 2.22, 95% CI (1.11, 4.63)) compared with controls. In secondary analyses, cases were more likely to have detected thyroid nodules (25.5% vs. 17.0%, adjusted OR = 1.77, 95% CI (1.15, 2.74)), hyperparathyroidism or parathyroid adenomas (3.5% vs. 1.3%, adjusted OR = 3.00, 95% CI (1.00, 11.64)), benign breast masses (6.0% vs. 3.3%, adjusted OR = 3.25, 95% CI (1.28, 8.78)), and benign prostatic hyperplasia (20.5% vs. 5.3%, adjusted OR = 3.20, 95% CI (1.14, 10.60)). Using a composite outcome, cases had higher odds of detection of the composite of IPMN, thyroid nodules, parathyroid tumors, benign breast masses, and prostate hyperplasia (adjusted OR = 2.36, 95% CI: 1.60, 3.50) when compared with controls. Conclusions. Patients with benign adrenocortical tumors had higher odds of detected pancreatic IPMN, as well as thyroid nodules, parathyroid tumors, benign breast masses, and prostate hyperplasia compared with patients with normal adrenal glands. These associations may have important implications for patient care and healthcare economics, regardless of whether they reflect incidental discoveries due to imaging detection or frequency bias, or a common risk for developing multiple neoplasia.
Puig-Barrachina V, Rodríguez-Sanz M, Domínguez-Berjón MF, Martín U, Luque MÁ, Ruiz M, Perez G. Decline in fertility induced by economic recession in Spain [Internet]. Gaceta Sanitaria 2019; Publisher's VersionAbstract
Objective To describe trends in fertility in Spain before (pre-recession; 1998-2008) and during (recession period; 2009-2013) the economic crisis of 2008, taking into account women's age and regional unemployment in 2010. Method The study consisted of a panel design including cross-sectional ecological data for the 17 regions of Spain. We describe fertility trends in Spain in two time periods, pre-recession (1998-2008) and recession (2009-2013). We used a cross-sectional, ecological study of Spanish-born women to calculate changes in fertility rates for each period using a linear regression model adjusted for year, period, and interaction between them. Results We found that compared to the pre-recession period, the fertility rate in Spain generally decreased during the economic recession. However, in some regions, such as the Canary Islands, this decrease began before the onset of the recession, while in other regions, such as the Basque country, the fertility rate continued to grow until 2011. The effects of the recession on the fertility rate are clearly observed in women aged 30-34 years. Conclusions The current economic recession has disrupted the positive trend in fertility that began at the start of this century. Since Spain already had very low fertility rates, the further decline caused by the economic recession could jeopardize the sustainability of welfare-state systems. Resumen Objetivo Describir las tendencias de la fecundidad en España en la época precrisis (1998-2008) y durante la crisis (2009-2013) económica, teniendo en cuenta la edad de las mujeres y el desempleo regional en 2010. Método Se utiliza un diseño panel que incluye datos ecológicos transversales para las 17 comunidades autónomas de España. Se describen las tendencias de fecundidad en los dos periodos. Para calcular los cambios en las tasas de fecundidad se utiliza un modelo de regresión lineal ajustado por año, periodo e interacción de ellas. Resultados En comparación con el periodo anterior, la tasa de fecundidad global en España disminuyó durante la crisis económica. Sin embargo, en algunas comunidades, como las Islas Canarias, esta disminución comenzó antes del inicio de la crisis, mientras que en otras, como el País Vasco, la tasa de fecundidad continuó creciendo hasta 2011. Los efectos de la crisis en la fecundidad se observan claramente en mujeres de 30 a 34 años. Conclusiones La crisis económica actual ha interrumpido la tendencia positiva en la fecundidad que comenzó a principios de este siglo. Dado que España ya tenía tasas de fecundidad muy bajas, el descenso causado por la crisis económica podría poner en peligro la sostenibilidad de los sistemas de bienestar social.
Luque-Fernandez MA, Redondo-Sánchez D, Maringe C. cvauroc: Command to compute cross-validated area under the curve for ROC analysis after predictive modeling for binary outcomes [Internet]. The Stata Journal 2019;19(3):615-625. Publisher's VersionAbstract
Receiver operating characteristic (ROC) analysis is used for comparing predictive models in both model selection and model evaluation. ROC analysis is often applied in clinical medicine and social science to assess the tradeoff between model sensitivity and specificity. After fitting a binary logistic or probit regression model with a set of independent variables, the predictive performance of this set of variables can be assessed by the area under the curve (AUC) from an ROC curve. An important aspect of predictive modeling (regardless of model type) is the ability of a model to generalize to new cases. Evaluating the predictive performance (AUC) of a set of independent variables using all cases from the original analysis sample often results in an overly optimistic estimate of predictive performance. One can use K-fold cross-validation to generate a more realistic estimate of predictive performance in situations with a small number of observations. AUC is estimated iteratively for k samples (the “test” samples) that are independent of the sample used to predict the dependent variable (the “training” sample). cvauroc implements k-fold cross-validation for the AUC for a binary outcome after fitting a logit or probit regression model, averaging the AUCs corresponding to each fold, and bootstrapping the cross-validated AUC to obtain statistical inference and 95% confidence intervals. Furthermore, cvauroc optionally provides the cross-validated fitted probabilities for the dependent variable or outcome, contained in a new variable named \_fit; the sensitivity and specificity for each of the levels of the predicted outcome, contained in two new variables named \_sen and \_spe; and the plot of the mean cross-validated AUC and k-fold ROC curves.
Schomaker M, Luque-Fernandez MA, Leroy V, Davies MA. Using longitudinal targeted maximum likelihood estimation in complex settings with dynamic interventions [Internet]. Statistics in Medicine 2019; Publisher's VersionAbstract
Longitudinal targeted maximum likelihood estimation (LTMLE) has very rarely been used to estimate dynamic treatment effects in the context of time-dependent confounding affected by prior treatment when faced with long follow-up times, multiple time-varying confounders, and complex associational relationships simultaneously. Reasons for this include the potential computational burden, technical challenges, restricted modeling options for long follow-up times, and limited practical guidance in the literature. However, LTMLE has desirable asymptotic properties, ie, it is doubly robust, and can yield valid inference when used in conjunction with machine learning. It also has the advantage of easy-to-calculate analytic standard errors in contrast to the g-formula, which requires bootstrapping. We use a topical and sophisticated question from HIV treatment research to show that LTMLE can be used successfully in complex realistic settings, and we compare results to competing estimators. Our example illustrates the following practical challenges common to many epidemiological studies: (1) long follow-up time (30 months); (2) gradually declining sample size; (3) limited support for some intervention rules of interest; (4) a high-dimensional set of potential adjustment variables, increasing both the need and the challenge of integrating appropriate machine learning methods; and (5) consideration of collider bias. Our analyses, as well as simulations, shed new light on the application of LTMLE in complex and realistic settings: We show that (1) LTMLE can yield stable and good estimates, even when confronted with small samples and limited modeling options; (2) machine learning utilized with a small set of simple learners (if more complex ones cannot be fitted) can outperform a single, complex model, which is tailored to incorporate prior clinical knowledge; and (3) performance can vary considerably depending on interventions and their support in the data, and therefore critical quality checks should accompany every LTMLE analysis. We provide guidance for the practical application of LTMLE.
Luque-Fernandez MA, Thomas A, Gelaye B, Racape J, Sanchez MJ, Williams MA. Secular trends in stillbirth by maternal socioeconomic status in Spain 2007–15: a population-based study of 4 million births [Internet]. European Journal of Public Health 2019; Publisher's VersionAbstract

Stillbirth, one of the urgent concerns of preventable perinatal deaths, has wide-reaching consequences for society. We studied secular stillbirth trends by maternal socioeconomic status (SES) in Spain.We developed a population-based observational study, including 4 083 919 births during 2007–15. We estimate stillbirth rates and secular trends by maternal SES. We also evaluated the joint effect of maternal educational attainment and the Human Development Index (HDI) of women’s country of origin on the risk of stillbirth. The data and statistical analysis can be accessed for reproducibility in a GitHub repository:

We found a consistent pattern of socioeconomic inequalities in the risk of delivering a stillborn, mainly characterized by a persistently higher risk, over time, among women with lower SES. Overall, women from countries with low HDIs and low educational attainments had approximately a four times higher risk of stillbirth (RR: 4.44; 95%CI: 3.71–5.32). Furthermore, we found a paradoxical reduction of the stillbirth gap over time between the highest and the lowest SESs, which is mostly due to the significant and increasing trend of stillbirth risk among highly educated women of advanced maternal age.

Our findings highlight no improvement in stillbirth rates among women of lower SES and an increasing trend among highly educated women of advanced maternal age over recent years. Public health policies developing preventive programmes to reduce stillbirth rates among women with lower SES are needed as well as the necessity of further study to understand the growing trend of age-related stillbirths among highly educated women in Spain.

Lee SF, Luque-Fernandez MA. Is cancer-related death associated with circadian rhythm? [Internet]. Cancer Communications 2019;39(1):27. Publisher's VersionAbstract

We investigated the temporal pattern of death in cancer patients using a large sample size and robust statistical methods to account for chronobiological periodicity.

We did not detect a circadian pattern of cancer death. The present study evaluated the temporal pattern of death among cancer patients using trigonometric functions and with time modeled in a circular scale.

To conclude, we found no evidence of a chronobiological circadian pattern in death among cancer patients by using robust statistical methods and data from a large population in a hospital setting. Increased understanding of the temporal pattern of deaths may yield important insights toward understanding external factors associated with cancer death.

Luque-Fernández MÁ, Calduch EN. Education in public health, epidemiology and biostatistics in Spain from a global and comparative perspective [Internet]. Gaceta Sanitaria 2019; Publisher's VersionAbstract
In the last ten years, there have been intense debates to boost the discipline, make it relevant to the genomic revolution, and place it at the forefront of the digital era. As a result, training in public health and epidemiology has been renewed, with marked statistical and methodological reinforcement, such as the current emphasis on causal inference, or the inclusion of master's in data science for health as a new academic degree in some schools of public health