Vitamin D deficiency, Smoking, and Lung Function in the Normative Aging Study

Similar documents
Online supplement. Absolute Value of Lung Function (FEV 1 or FVC) Explains the Sex Difference in. Breathlessness in the General Population

Supplementary Appendix

Supplementary Appendix

Figure E1 Manhattan and QQ plots for FEV 1 meta-analyses across all Hispanic ancestry groups

Correlation and regression

COPD Phenotyping: Hyperpolarized Helium-3 Magnetic Resonance Functional Imaging

Biostatistics for biomedical profession. BIMM34 Karin Källen & Linda Hartman November-December 2015

Sociology 362 Data Exercise 6 Logistic Regression 2

The National Lung Health Education Program. Spirometric Reference Values for the 6-s FVC Maneuver*

Interplay Between S-adenosylmethionine, Folate, Cobalamin, and Arsenic Methylation in Bangladesh

Practical Biostatistics

An introduction to biostatistics: part 1

Suppose that we are concerned about the effects of smoking. How could we deal with this?

Statistics in medicine

Statistics in medicine

WORKSHOP 3 Measuring Association

Survey of Smoking Behavior. Survey of Smoking Behavior. Survey of Smoking Behavior

Introduction to Analysis of Genomic Data Using R Lecture 6: Review Statistics (Part II)

Analysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington

Disadvantages of using many pooled t procedures. The sampling distribution of the sample means. The variability between the sample means

10.4 Hypothesis Testing: Two Independent Samples Proportion

Nemours Biomedical Research Statistics Course. Li Xie Nemours Biostatistics Core October 14, 2014

The t-distribution. Patrick Breheny. October 13. z tests The χ 2 -distribution The t-distribution Summary

Einstein College of Medicine and Montefiore Medical Center, Bronx, New York

International Journal of Scientific & Engineering Research, Volume 5, Issue 4, April ISSN

The Influence Function of the Correlation Indexes in a Two-by-Two Table *

Lecture 9. Selected material from: Ch. 12 The analysis of categorical data and goodness of fit tests

Graphical models and causality: Directed acyclic graphs (DAGs) and conditional (in)dependence

Model and Working Correlation Structure Selection in GEE Analyses of Longitudinal Data

Readings Howitt & Cramer (2014) Overview

FORECASTING METHODS AND APPLICATIONS SPYROS MAKRIDAKIS STEVEN С WHEELWRIGHT. European Institute of Business Administration. Harvard Business School

Online Supplementary Material to: Perfluoroalkyl Chemicals, Menstrual Cycle Length and Fecundity: Findings from a. Prospective Pregnancy Study

Linear Regression Results Including Fish and Game Dioxin Congener Intake and Air Concentration Estimates

TriMatch: An R Package for Propensity Score Matching of Non-binary Treatments

Lecture 5: ANOVA and Correlation

3.1 Events, Sample Spaces, and Probability

Statistics in medicine

Inferences About Two Proportions

Analysis of Variance Bios 662

( ). Switch x and y and solve for y:

BIOL 51A - Biostatistics 1 1. Lecture 1: Intro to Biostatistics. Smoking: hazardous? FEV (l) Smoke

Achieving Optimal Covariate Balance Under General Treatment Regimes

Logistic Regression - problem 6.14

Chapter 14. Multiple Regression Models. Multiple Regression Models. Multiple Regression Models

Lecture 6 Multiple Linear Regression, cont.

1. Poisson distribution is widely used in statistics for modeling rare events.

STT 843 Key to Homework 1 Spring 2018

Describing Bivariate Relationships

Person-Time Data. Incidence. Cumulative Incidence: Example. Cumulative Incidence. Person-Time Data. Person-Time Data

Biostatistics 4: Trends and Differences

Beyond GLM and likelihood

Lecture 4 Multiple linear regression

Causal Inference with General Treatment Regimes: Generalizing the Propensity Score

Readings Howitt & Cramer (2014)

Nonparametric Mixture of Gaussian Processes with Constraints

Using Mixture Latent Markov Models for Analyzing Change in Longitudinal Data with the New Latent GOLD 5.0 GUI

Chapter 15 Sampling Distribution Models

Chapter 2: Describing Contingency Tables - I

5.1 Bivariate Relationships

For more information about how to cite these materials visit

CSC 411: Lecture 09: Naive Bayes

Correlation and Regression Bangkok, 14-18, Sept. 2015

10: Crosstabs & Independent Proportions

Epidemiology Principle of Biostatistics Chapter 14 - Dependent Samples and effect measures. John Koval

A spatial causal analysis of wildfire-contributed PM 2.5 using numerical model output. Brian Reich, NC State

Predict y from (possibly) many predictors x. Model Criticism Study the importance of columns

STA 291 Lecture 8. Probability. Probability Rules. Joint and Marginal Probability. STA Lecture 8 1

Survey of Smoking Behavior. Samples and Elements. Survey of Smoking Behavior. Samples and Elements

Conditional Probability

Binomial and Poisson Probability Distributions

Review. Number of variables. Standard Scores. Anecdotal / Clinical. Bivariate relationships. Ch. 3: Correlation & Linear Regression

Ph.D. Preliminary Examination Statistics June 2, 2014

You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?

Chapter 5: Data Transformation

= X+(z.95)SE(X) X-(z.95)SE(X)

Before and After Models in Observational Research Using Random Slopes and Intercepts

Introduction to logistic regression

CHAPTER6 LINEAR REGRESSION

Interpreting the coefficients

Use of Ultra-High Performance LC-MS for Challenging Assays. David Browne, Laboratory Manager (LC-MS-MS), Covance Laboratories, Harrogate

University of California at Berkeley Fall Introductory Applied Econometrics Final examination. Scores add up to 125 points

Foundations of Correlation and Regression

Chapter 2 Exercise Solutions. Chapter 2

Inference. ME104: Linear Regression Analysis Kenneth Benoit. August 15, August 15, 2012 Lecture 3 Multiple linear regression 1 1 / 58

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

Introduction to Structural Equation Modeling Dominique Zephyr Applied Statistics Lab

Time-Invariant Predictors in Longitudinal Models

Instructions. Information. Advice

Module 19: Simple Linear Regression

STA6938-Logistic Regression Model

Sample Size and Power Considerations for Longitudinal Studies

Statistical methods for comparing multiple groups. Lecture 7: ANOVA. ANOVA: Definition. ANOVA: Concepts

Department of Mathematics & Statistics Stat 2593 Final Examination 19 April 2001

Energy Flow in Technological Systems. December 01, 2014

Review of Regression Basics

Introductory Statistics with R: Linear models for continuous response (Chapters 6, 7, and 11)

Fling-Step vs Filtering: Effects on Response Spectra and Peak Motions. David M. Boore. Introduction

Data Analysis as a Decision Making Process

First & Second Term Booklet

Richard D Riley was supported by funding from a multivariate meta-analysis grant from

Transcription:

1 Vitamin D deficiency, Smoking, and Lung Function in the Normative Aging Study Nancy E. Lange MD, MPH, David Sparrow DSc, Pantel Vokonas MD, Augusto A. Litonjua, MD, MPH Online data supplement

2 For all tables, reference group for vitamin D deficiency (VDD) was the group of subjects without VDD, therefore negative effect estimates represent lower values in subjects with VDD compared to subjects without VDD. Table E1: Baseline characteristics of subjects included and not included in current analysis Baseline (Time 1: 1984 89) included not included (N=626) (N=798) Mean (SD) or Mean (SD) or p value* N(%) N(%) Age 59.9 (6.9) 65.3 <.1 BMI 27.2 (3.3) 26.9 (3.6).15 FEV 1 3.21 (.63) 2.93 (.66) <.1 FEV 1 %pred 89.5 (14.8) 86.8 (16.7).1 FVC 4.24 (.74) 3.89 (.74) <.1 FVC %pred 93.3 (12.8) 89.7 (14.3) <.1 FEV 1 /FVC ratio 75.5 (6.8) 75.1 (7.7).35 pack years 19.2 (22.3) 25.5 (27.) <.1 current smoking (N, % yes) 64 (1%) 12 (15%).8 GOLD 528 (84%) 634 (79%) 1 36 (6%) 44 (5.5%) 2 53 (8%) 9 (11.3%) 3 5 (1%) 16 (2%) 4 4 (.5%) 14 (1.8%) * t-test or chi-sq Table E2: Cross-sectional analysis, main effects of Vitamin D deficiency in relation to lung function, at each of three time points (Time 1: 1984-89, Time 2: 1992-1999, Time 3: 1995-23) Effect estimates (β) are for vitamin D deficiency. Multivariable analysis includes adjustment for age, height, pack-years of smoking, BMI and season of blood draw. Time 1 Time 2 Time 3 β p value β p value β p value FEV 1.53.21.12.84.73.15 FVC.29.53.1.11.45.42 ratio.811.14 1.66.3 1.34.83

3 Table E3: Bivariable and multivariable longitudinal models examining the main effect of vitamin D (both VDD and VDC) and lung function decline There were no significant relationships between vitamin D and lung function decline with the exception of an inverse relationship between VDC and FVC decline Vitamin D Deficiency Vitamin D Continuous Bivariate Multivariate Bivariate Multivariate β p value β p value β p value β p value FEV1.39.88.21.39.8.46.19.6 FVC.2.48.44.11.25.23.35.8 FEV1/FVC.57.23.46.37.31.9.216.25 Table E4: Cross-sectional analysis, vitamin D deficiency and smoking (packyears) in relation to lung function (Times 1-3) Effect estimates (β) are for the interaction between VDD and smoking as labeled at the top of the column, i.e. either VDD x pack-years or VDD x smoking status. Smoking status was analyzed as a dichotomous variable (current vs. ex- and never-smokers). Reference group is never & ex-smokers. Pack years Smoking Status β p value β p value Time 1 FEV 1.19.3.21.143 FVC.17.53.148.33 FEV 1 /FVC.225.35 2.23.21 Time 2 FEV 1.2.43.525.16 FVC.648.43.584.16 FEV 1 /FVC.21.53 3.524.23 Time 3 FEV 1.72.7.888.1 FVC.55.211.64.13 FEV 1 /FVC.938.39 14.491 <.1 Table E5: Cross-sectional analysis, vitamin D continuous and smoking (packyears) in relation to lung function Effect estimates (β) are for the interaction term between vitamin D as a continuous variable (VDC) and smoking (pack-years). Positive values imply incremental increase in lung function value for each 1ng/ml increase in vitamin D level.

4 β p value Time 1 FEV 1.14.19 FVC.38.74 FEV 1 /FVC.24.8 Time 2 FEV 1.8.32 FVC.1.3 FEV 1 /FVC.5.65 Time 3 FEV 1.16.1 FVC.15.14 FEV 1 /FVC.86.55 Table E6: Longitudinal Analysis showing effect of combination of vitamin D deficiency, smoking status, and change in lung function over time Effect estimates and p values are for the interaction term with vitamin D deficiency, time (in years) and smoking status. Negative effect estimates imply that in current smokers with vitamin D deficiency, lung function decline was more rapid as compared to current smokers who were vitamin D sufficient. Smoking Status β p value FEV1.26.52 FVC.9.54 FEV1/FVC.76.23 Table E7: Longitudinal Analysis showing effect of combination of vitamin D as continuous variable, pack-years of smoking, and change in lung function over time. Effect estimates and p values are for the interaction term with vitamin D as a continuous variable, time (in years), and pack-years of smoking. Positive effect estimates imply an incremental decrease in rate of decline (rate becomes less negative) per year per ng/ml increase in vitamin D level.

5 β p value FEV1 7.5E 6.7 FVC 8E 6.15 FEV1/FVC.2.85

6 Figure E1: Relationship between FEV 1 and pack-years of smoking There was a significant inverse relationship between FEV 1 and pack years of smoking 6 5 FEV 1 4 3 p=<.1 2 1 2 4 6 8 1 12 14 16 Pack years of Smoking

7 Figure E2: Relationship between lung function and vitamin D levels There was no significant relationship between vitamin D levels and A) FEV 1, B) FVC or C) FEV 1 /FVC. A) 6 5 4 p=.16 FEV 1 3 2 1 2 4 6 8 1 Vitamin D (ng/ml) B) 7 6 5 FVC 4 3 2 p=.15 1 2 4 6 8 1 Vitamin D (ng/ml) C)

8 FEV 1 /FVC 1 9 8 7 6 5 4 3 2 1 2 4 6 8 1 Vitamin D (ng/ml) p=.6 Figure E3: Relationship between vitamin D (continuous) and rates of change in lung function There was no significant relationship between vitamin D levels and rate of change of FEV 1 or FEV 1 /FVC, however there was a significant inverse relationship between vitamin D and rate of change of FVC (see Table E3) though the effect size was very small (β=-.35, p=.8 in adjusted model). Shown here are the bivariate relationships between vitamin D levels and lung function decline.

9 Rate of change of FEV1 (L/yr).1.5.5.1 FEV 1 1 2 3 4 5 6 7 8 9 1.15 Vitamin D level (ng/ml).15 FVC Rate of change of FVC (L/yr).1.5.5.1.15.2.25 1 2 3 4 5 6 7 8 9 1.3 Vitamin D level (ng/ml)

1 Rate of change of FEV1/FVC (%/yr) 14 12 1 8 6 4 2 FEV 1 /FVC 1 2 3 4 5 6 7 8 9 1 Vitamin D level (ng/ml)

11 Figure E4: Relationship between Lung Function, Vitamin D deficiency, and Smoking Status Shown here is the cross-sectional relationship between lung function (FEV 1 % predicted) and smoking status. In current smokers, the mean FEV 1 % predicted was significantly lower in those who were vitamin D deficient, compared to those who were not. In ex and non-smokers, there was no significant difference between subjects with and without vitamin D deficiency. p=<.1 NS NS 9 8 7 6 FEV 1 %pred 5 4 3 Vitamin D 2ng/ml Vitamin D >2ng/ml 2 1 Current smokers Ex smokers Never smokers Smoking Status