Research article · Cross-sectional methodological study · STROBE-aligned

Socioeconomic and Housing Conditions in Rural Bolivia: A Privacy-First Reproducible Survey Analysis

Monica Cueto Tapia
Reproducible household survey analysis · Version 1.0 · 18 June 2026

Abstract

Background: Household surveys can inform housing and social policy, but public reproducibility may conflict with respondent confidentiality. Objective: To specify and evaluate a privacy-preserving, reproducible protocol for analysing housing adequacy and socioeconomic vulnerability. Methods: A cross-sectional synthetic population comprised 60 households and 271 members. Referential integrity and disclosure constraints were validated before analysis. Missingness was quantified, weighted descriptive estimands were accompanied by 2,000-resample bootstrap intervals, and conditional associations were estimated using parsimonious logistic regression. Results: Weighted adequate housing was 37.9% (95% bootstrap interval 26.7%-51.7%). Synthetic median monthly income was BOB 2210, with 15.0% missing. All exploratory regression intervals crossed the null value, indicating substantial statistical imprecision. Conclusions: A privacy-preserving analytical protocol can support computational reproducibility without releasing identifiable microdata. Valid empirical inference would additionally require authorised observations, a documented sampling design and validated measurement instruments.

Keywords: household survey; housing adequacy; reproducible research; synthetic data; missing data; logistic regression; research ethics

1. Introduction

Housing is multidimensional: construction materials, basic services, crowding and household resources jointly shape wellbeing. Survey analysis often compresses these dimensions into percentages without reporting denominators, uncertainty, missingness or the assumptions behind composite measures. These omissions weaken both scientific interpretation and operational decision-making.

A second challenge is disclosure. Household microdata may combine names, ages, occupations, income, health information and property characteristics. Removing a name column is insufficient when rare combinations can re-identify respondents. This project therefore treats privacy architecture as part of statistical quality rather than as an administrative appendix.

1.1 Research questions

  1. How can adequate housing be operationalised transparently from multiple components?
  2. Which socioeconomic characteristics are conditionally associated with adequate housing in a small synthetic cross-sectional sample?
  3. How should missingness and sampling uncertainty be quantified and communicated?
  4. How can survey analysis remain computationally reproducible without exposing confidential source records?

2. Methods

2.1 Study design and reporting framework

This simulation-based methodological study follows the Introduction-Methods-Results-Discussion structure and is mapped to the STROBE checklist for cross-sectional studies [1,2]. It is not a registered observational study and does not estimate population parameters for a real geographic area.

2.2 Data generation and units of analysis

A fixed pseudorandom seed generated 60 households and 271 nested members. Aggregate zones are fictional. No value was copied, perturbed or sampled from a real respondent. Household and person denominators are kept separate throughout.

2.3 Privacy and quality controls

Automated tests reject direct-identifier columns, duplicate identifiers, orphaned member records, invalid binary outcomes and impossible ages. Public file rules exclude SPSS, legacy Word and private raw-data directories from version control.

2.4 Estimands and uncertainty

For survey weights wi and outcome yi, the weighted mean is:

ȳw = Σi wiyi / Σiwi (1)

The 95% interval uses the 2.5th and 97.5th percentiles of 2,000 household bootstrap estimates:

CI.95(θ̂) = [ Q.025(θ̂*), Q.975(θ̂*) ] (2)

2.5 Housing and vulnerability measures

The housing score is the unweighted share of six documented components: electricity, improved water, improved sanitation, quality walls, quality floor and absence of severe crowding. Adequate housing requires a score of at least two thirds. The vulnerability index combines housing, head education and capped income:

Vi = 100[1 − (0.42Hi + 0.28Ei/16 + 0.30 min(Yi/6000,1))] + εi (3)

2.6 Exploratory model

A logistic generalised linear model estimates conditional associations with adequate housing:

log[pi/(1−pi)] = β0 + β1Ei + β2log(Yi+1) + β3Si + β4Gi (4)

Education, log income and household size are standardised. Missing education and income are median-imputed for computational completeness; empirical analysis should use design-aware multiple imputation and sensitivity analysis [5]. No coefficient is interpreted causally.

3. Results

Table 1. Synthetic sample characteristics
CharacteristicEstimate
Households, n60
Household members, n271
Household size, median (IQR)5.0 (3.0-6.0)
Age of household head, mean (SD)47.0 (11.1)
Woman-headed households, %36.7
Education of head, median years (IQR)7.0 (5.0-9.0)
Monthly income, median BOB (IQR)2210 (1690-3005)
Income missing, %15.0
Adequate housing, weighted % (95% CI)37.9 (26.7-51.7)
Vulnerability index, mean (SD)44.6 (12.0)
Missing data audit
Figure 1. Missingness by variable. Income has the largest deliberately introduced missing proportion. Missingness is displayed before outcome modelling.

3.1 Housing adequacy and services

The weighted adequate-housing estimate was 37.9% (95% bootstrap interval 26.7%-51.7%). Wide zone-specific intervals reflect the small number of households and discourage unstable rankings.

Housing adequacy by synthetic zone
Figure 2. Adequate housing by fictional aggregate zone. Points are estimates and bars are approximate 95% intervals.

3.2 Socioeconomic patterning

Median synthetic monthly household income was BOB 2210. The mean vulnerability index was 44.6 on a 0-100 scale. Correlations describe co-movement only and do not establish pathways of causation.

Income and vulnerability
Figure 3. Income and multidimensional vulnerability. Point size indicates household size and colour indicates housing classification.

3.3 Exploratory regression

Table 2. Exploratory logistic model for adequate housing
PredictorOdds ratio95% CIp-value
Education (+1 SD)0.930.51-1.690.811
Log income (+1 SD)0.870.46-1.630.655
Household size (+1 SD)0.920.53-1.580.759
Woman-headed household1.620.53-4.990.399
Interpretation: every confidence interval includes 1. The defensible conclusion is insufficient precision to distinguish these adjusted associations from the null, not evidence that the predictors have no relationship with housing.
Regression odds ratios
Figure 4. Adjusted odds ratios on a logarithmic scale. The dashed vertical line marks the null value.

4. Discussion

This analysis integrates disclosure controls, explicit denominators, uncertainty, documented measurement and automated reproduction. The wide intervals are analytically important: they show why complex machine learning and detailed subgroup ranking would be inappropriate for a 60-household study.

For a real survey, scientific interpretation would depend on information that cannot be reconstructed from the data alone: target population, sampling frame, selection probabilities, field dates, non-response process, questionnaire, interviewer procedures and ethics approval. These elements should be documented before substantive re-analysis.

5. Limitations

6. Conclusion

Strong research communication is not the accumulation of charts or equations. It is the alignment of a meaningful question, appropriate estimand, transparent method, uncertainty, ethical data handling and restrained conclusion. This workflow provides a publishable foundation once authorised, anonymised source data and a complete survey design are available.

Declarations

Ethics and privacy: No real participant record is processed by the public pipeline. Data availability: Synthetic CSV files and their generator are included. Code availability: The complete pipeline, tests and dependency lockfile accompany this paper. Funding: No external funding was received. Competing interests: None declared. Author contribution: Monica Cueto Tapia: conceptualisation, methodology, analysis design, interpretation and presentation.

References

  1. von Elm E, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement. Lancet. 2007;370:1453-1457. doi:10.1016/S0140-6736(07)61602-X.
  2. STROBE Initiative. Checklist for cross-sectional studies. Official checklist.
  3. Lohr SL. Sampling: Design and Analysis. 3rd ed. Chapman and Hall/CRC; 2021.
  4. Peng RD. Reproducible research in computational science. Science. 2011;334:1226-1227. doi:10.1126/science.1213847.
  5. Little RJA, Rubin DB. Statistical Analysis with Missing Data. 3rd ed. Wiley; 2019.
  6. Wickham H. Tidy data. Journal of Statistical Software. 2014;59:1-23. doi:10.18637/jss.v059.i10.