Analytic codes for "The Association between Vitamin D and Incident Herpes Zoster: A UK Biobank Study"

These are the codes for plotting the forest plot in the manuscript.

The hazard ratio of the Weibull regression model was obtained by using survival package (ver 3.2-11), and the regression output was organized using tidyr() function from the Broom package (ver 0.7.9). A forest plot organized these outputs was plotted using ggplot2 package (ver 3.3.5).

The results of sensitivity analysis were summarized in Supplementary Figures.

Sensitivity analysis	Figure
Stop follow-up by 31 August 2013	Supplementary Figures 2 Supplementary Figures 3
Using different covariates definitions	Supplementary Figures 4
Use Cox-regression to analyse the association between exposure and outcomes	Supplementary Figures 5 Supplementary Figures 6
Use Poisson regression to analyse the association between exposure and outcomes	Supplementary Figures 7 Supplementary Figures 8

Files Index:

Main analysis:

• Figure 3: The association between vitamin D status and the risk of herpes zoster.

• Figure 4: The association between vitamin D intake and the risk of herpes zoster.

Supplementary Figures:

• Supplementary Figure 2. Vitamin D status and the risk of herpes zoster excluding records after September 2013.

• Supplementary Figure 3. Vitamin D intake and the risk of herpes zoster excluding records after September 2013.

• Supplementary Figure 4. Sensitivity analysis of using different definitions of clinical covariates.

• Supplementary Figure 5. Sensitivity analysis of using stratified Cox regression to assess the association between vitamin D status and the hazards of incident herpes zoster before and after the vaccination program initiated.

• Supplementary Figure 6 Sensitivity analysis of using Cox proportional-hazards model to examine the association between vitaminD intake and the risk of herpes zoster.

• Supplementary Figure 7. Sensitivity analysis of using Poission regression model to examine the association between vitamin D status and the risk of herpes zoster

• Supplementary Figure 8. Sensitivity analysis of using Poission regression model to examine the association between vitamin D supplementation and the risk of herpes zoster

Code-lists:

All diagnostic codes are in CTV3, SNOMED-CT, or read 2. Prescription codes are in BNF, DM+D and Read 2 codes

• Comorbidities: including asthma, COPD, Chronic kidney disease, rheumatoid arthritis, inflammatory bowel diseases, depression,

• Immunosuppression: including including organ transplantation, chemoradiotherapy, cell-mediated immunosuppression, HIV, blood cancers, chemotherapy (biological and non-biological agents), bone marrow transplantation and long-term oral steroid.

• Vitamin D prescriptions: vitamin D prescriptions codes

• Herpes zoster diagnosis: diagnostic codes for diagnosing herpes zoster cases