Seol mar théacs é seo: Pairwise joint modeling of clustered and high-dimensional outcomes with covariate missingness in pediatric pneumonia care