treatment effectiveness and healthcare resource use in daily practice may be different from effectiveness and resource use in clinical trials . this is increasingly recognized by researchers , policy makers and decision makers responsible for pricing and reimbursement of healthcare interventions . hence , the importance of data generated from sources that reflect the use and the associated outcomes in routine practice settings is growing . suitable sources of such real - life data could be health records kept by physicians ( eg , their routine records or specifically established databases ; paper - based or electronic ) , patient registries that enroll patients with specific diseases or other characteristics of interest ( eg , cancer registries ) , or administrative claims databases of healthcare insurers and provider organizations set up for the purpose of reimbursement of providers for their expenses . the inherent problem of analyses conducted on such databases relates to the fact that the data are often collected for other purposes . when treatment is not assigned to patients at random , disease severity and prognosis of patients the currently most widely used severity classification is the spirometric classification of the global initiative for chronic obstructive lung disease ( gold),1 which is based on lung function in terms of the forced expiratory volume in one second as a percentage of the value predicted for sex , age and height ( fev1-%predicted ) . it is used to diagnose copd , to monitor disease progression and to aid in treatment decisions . lung function is an important indicator of copd severity because it is associated with mortality,2 exacerbations,3,4 health care utilization and costs.5 there is also an association between lung function and quality of life6 although this relationship is not as strong.7,8 consequently , the same associations have been shown for the gold classification as it is based on fev1.913 it is widely recognized , also by gold,1 that the impact of copd does not only depend on lung function and further - more that is not only a lung disease but a multi - dimensional disease with many systemic , extra - pulmonary effects . therefore , composite measures of copd severity have been proposed . however , almost all of these still accept the importance of airflow limitation and they are partly based on the fev1-%predicted.1418 however , lung function measures obtained through spirometry are not routinely available in databases and registries . most routine databases also lack information on the other parameters that constitute the more recent composite measures of copd severity . previous retrospective database research has approximated copd severity using demographic ( eg , age and/or smoking status ) and resource utilization data ( eg , medications used and/or hospital admission).1924 in this study , we attempted to develop a multivariable predictive algorithm to derive the severity of copd as classified by gold , using variables that are commonly available in routine databases . the analysis was performed on the baseline data of the  understanding potential long - term impacts on function with tiotropium  ( uplift ) trial , a large , 4-year trial in copd patients with rate of decline in fev1 as primary endpoint.25 this trial was especially suitable for this purpose because of its size ( with almost 6000 patients randomized ) , and because among the baseline data collected are many data commonly available in routine databases ( eg , demographics , medications used in the past , hospital admissions in the year preceding enrollment ) . moreover , fev1 was collected with high quality , thus providing confidence in the gold severity assignment of trial patients . the uplift trial was a randomized , double - blind , placebo - controlled 4-year trial , investigating the effect of tiotropium on the yearly rate of decline in fev1 in patients with moderate to very severe copd according to gold ( stages ii to iv ) . 25 main inclusion criteria besides a diagnosis of copd included age of 40 years or over , a smoking history of at least 10 pack - years , a post - bronchodilator fev1-%predicted of 70% or less , and an fev1 of 70% or less of the forced vital capacity . key exclusion criteria were a history of asthma , a copd exacerbation or respiratory infection within 4 weeks before screening , use of supplemental oxygen for more than 12 hours per day . uplift collected information on patients demographics , comorbidities and co - medications as well as copd characteristics and smoking history , along with contacts with health care providers in various settings in the year preceding enrollment . in order to enhance the homogeneity of treatment patterns , the sample was limited to 3698 patients from western europe , the united states , australia and new zealand with copd severity stage ii , iii or iv . the sample was randomly split into two sets : one for developing the algorithm ( 2/3 of patients ) and the other for validation ( 1/3 of patients ) . potential variables for inclusion in the algorithm to approximate disease severity were those used in previous database research and those used in diagnostic and severity classifications.1417,2629 variables were then selected if available in uplift and likely to be available in routine databases . the final list of variables to be tested consisted of age ( continuously in years , squared , and in categories : 4049 , 5059 , 6069 , 7079 , 8090 , 90 + ) . bmi ( continuously and categorized : < 18.5 , 18.521.0 , 21.025.0 , 25.030.0 , > 30.0 ) , smoking status ( current / former smoker ) , pack - years of smoking ( continuously and categorized : < 20 , 2060 , > 60 ) , sex , type of respiratory maintenance medications ( short - acting and long - acting bronchodilators , inhaled corticosteroids ( ics ) , mucolytics , leukotriene modifiers , xanthine ) , number of medication types ( 0 to 5 ) , treatment of exacerbations ( number of courses of antibiotics or oral steroids , categorized : 0 , 1 , 2 , 3 in one year ) , resource use in one year ( number of emergency room visits with / without hospital admission , categorized ( 0 , 1 , 2 ) , number of scheduled and unscheduled general practitioner ( gp ) visits ( categorized : 0 , 1 , 2 , 34 , 5 and 0 , 1 , 2 , 3 , respectively ) , hospital admissions ( yes / no ) , use of oxygen at home ( yes / no ) , statin use ( yes / no ) , use of other cardiovascular medication ( yes / no ) , number of co - morbidities ( categorized : 0 , 1 , 2 , 35 , 69 , 10 ) , charlson comorbidity index30 ( categorized : 1 , 2 , 3 , 4 ) , presence of selected co - morbidities ( arrhythmia , coronary heart disease , vascular disease , hypertension , disorders of nervous system , stroke , diabetes , depression , anemia , and platelet disorders , osteoporosis and cataract ) , as well as interaction terms ( age  gender , age  bmi , bmi  gender ) . several partial proportional odds ( ppo ) ordered logit models were used to estimate the probabilities of being in gold stages ii ( fev1%-predicted > 50% ) , iii ( fev1%-predicted fev1%-predicted 30%50% ) , and iv ( fev1%-predicted < 30% ) , as copd severity categories are ordered . the ppo resembles the standard ordered logit model , which is the best - known ordered regression model . in our case , with three possible outcomes y , it estimates the probability that y > stage ii and the probability that y > stage iii . ordered regression models assume that the observed ordered outcome y is a function of a continuous and unobservable variable y*. two thresholds j are assumed to determine in which stage a patient is classified : y > yj if y * > j , where j denotes a specific outcome and tj is the upper limit of y * for this outcome . y * is related to the explanatory variables x : yi*=xij+ij , where i denotes individual patients . since the latent variable y * does not completely equal the sum of the products of coefficients and variable values , the outcome y can not be determined with certainty from the data . assuming a logit distribution for the random error term makes it possible to model the probability that a patient is in a certain stage or worse ( y * > j ) . x i j ) , j = 1 , 2 ,  , m  1 . in contrast to the ppo model , however , the ordered logit model can only estimate one coefficient for each explanatory variable . this coefficient is assumed to be identical for all dichotomizations of the outcome variable ( in this case , stage iv versus ii / iii and stage iii / iv versus ii ) , and thus coefficient j = . this is the proportional odds or parallel regressions assumption . if this assumption is violated , which often happens in practice , estimates are invalid and important differences in the relationships at different thresholds may go unnoticed . the partial proportional odds model relaxes this assumption for variables where it does not hold.31 for these variables , a coefficient is estimated for each dichotomization . the final model was developed by stepwise backward elimination from the full model , which contained all variables . in each step the variable with the largest p - value was eliminated and the model was re - estimated . this process was repeated until all variables had at least a value of p  0.10 . we performed two sensitivity analyses in order to account for the fact that the proportion of patients with very severe copd ( gold iv ) in the dataset was smaller than the proportion of patients with moderate or severe copd . firstly , the final model was re - estimated in a sub - sample with a balanced distribution of patients across all three severity stages . this balance was achieved by using all gold stage iv patients and random draws from patients in stages ii and iii . secondly , the patients in gold stage iii and iv were combined into one group . the probability of being in  severe / very severe was then analyzed in a binary logit model with the same variables as in the final ppo model . all analyses were performed in stata 11.1 , ( statacorp lp , college station , tx ) using the gologit2 command for the ppo model.32 univariate tests were analysis of variance ( anova ) for continuous variables and chi - square for categorical variables . the regression results were used to predict the probabilities of being in each gold severity stage for each patient in the validation dataset . the agreement between the predicted and the observed stage was assessed using kappa statistics.33 for a kappa statistic , a value of 0 indicates that agreement has occurred by chance , whereas a value of 1 indicates perfect agreement . landis and koch suggested the following labels of agreement for ranges of values , which are often quoted : slight ( up to 0.20 ) , fair ( 0.210.40 ) , moderate ( 0.410.60 ) , substantial ( 0.610.80 ) , and almost perfect ( 0.811.0).34 we present kappa statistics unweighted and weighted . in the weighted kappa statistics , patients misclassified in the neighboring category ( ie , ii or iv instead of iii , or iii instead of ii or iv ) count as 0.5 agreement and patients misclassified in a non - neighboring category ( ii instead of iv or vice versa ) as no agreement ) . for the binary logit model , the c - statistic was computed as a measure of its discriminative performance . the c - statistic represents the area under the receiver operating characteristic ( roc ) curve , a plot of sensitivity against specificity . values can range from 0.5 ( no predictive ability ) to 1 ( perfect discrimination ) . lastly , in an analysis of correctly and incorrectly classified patients , the average values of fev1%-predicted were compared per severity stage . the uplift trial was a randomized , double - blind , placebo - controlled 4-year trial , investigating the effect of tiotropium on the yearly rate of decline in fev1 in patients with moderate to very severe copd according to gold ( stages ii to iv ) . 25 main inclusion criteria besides a diagnosis of copd included age of 40 years or over , a smoking history of at least 10 pack - years , a post - bronchodilator fev1-%predicted of 70% or less , and an fev1 of 70% or less of the forced vital capacity . key exclusion criteria were a history of asthma , a copd exacerbation or respiratory infection within 4 weeks before screening , use of supplemental oxygen for more than 12 hours per day . uplift collected information on patients demographics , comorbidities and co - medications as well as copd characteristics and smoking history , along with contacts with health care providers in various settings in the year preceding enrollment . in order to enhance the homogeneity of treatment patterns , the sample was limited to 3698 patients from western europe , the united states , australia and new zealand with copd severity stage ii , iii or iv . the sample was randomly split into two sets : one for developing the algorithm ( 2/3 of patients ) and the other for validation ( 1/3 of patients ) . potential variables for inclusion in the algorithm to approximate disease severity were those used in previous database research and those used in diagnostic and severity classifications.1417,2629 variables were then selected if available in uplift and likely to be available in routine databases . the final list of variables to be tested consisted of age ( continuously in years , squared , and in categories : 4049 , 5059 , 6069 , 7079 , 8090 , 90 + ) . bmi ( continuously and categorized : < 18.5 , 18.521.0 , 21.025.0 , 25.030.0 , > 30.0 ) , smoking status ( current / former smoker ) , pack - years of smoking ( continuously and categorized : < 20 , 2060 , > 60 ) , sex , type of respiratory maintenance medications ( short - acting and long - acting bronchodilators , inhaled corticosteroids ( ics ) , mucolytics , leukotriene modifiers , xanthine ) , number of medication types ( 0 to 5 ) , treatment of exacerbations ( number of courses of antibiotics or oral steroids , categorized : 0 , 1 , 2 , 3 in one year ) , resource use in one year ( number of emergency room visits with / without hospital admission , categorized ( 0 , 1 , 2 ) , number of scheduled and unscheduled general practitioner ( gp ) visits ( categorized : 0 , 1 , 2 , 34 , 5 and 0 , 1 , 2 , 3 , respectively ) , hospital admissions ( yes / no ) , use of oxygen at home ( yes / no ) , statin use ( yes / no ) , use of other cardiovascular medication ( yes / no ) , number of co - morbidities ( categorized : 0 , 1 , 2 , 35 , 69 , 10 ) , charlson comorbidity index30 ( categorized : 1 , 2 , 3 , 4 ) , presence of selected co - morbidities ( arrhythmia , coronary heart disease , vascular disease , hypertension , disorders of nervous system , stroke , diabetes , depression , anemia , and platelet disorders , osteoporosis and cataract ) , as well as interaction terms ( age  gender , age  bmi , bmi  gender ) . several partial proportional odds ( ppo ) ordered logit models were used to estimate the probabilities of being in gold stages ii ( fev1%-predicted > 50% ) , iii ( fev1%-predicted fev1%-predicted 30%50% ) , and iv ( fev1%-predicted < 30% ) , as copd severity categories are ordered . the ppo resembles the standard ordered logit model , which is the best - known ordered regression model . in our case , with three possible outcomes y , it estimates the probability that y > stage ii and the probability that y > stage iii . ordered regression models assume that the observed ordered outcome y is a function of a continuous and unobservable variable y*. two thresholds j are assumed to determine in which stage a patient is classified : y > yj if y * > j , where j denotes a specific outcome and tj is the upper limit of y * for this outcome . y * is related to the explanatory variables x : yi*=xij+ij , where i denotes individual patients . since the latent variable y * does not completely equal the sum of the products of coefficients and variable values , the outcome y can not be determined with certainty from the data . assuming a logit distribution for the random error term makes it possible to model the probability that a patient is in a certain stage or worse ( y * > j ) . j ) = exp( j + x i j)/1 + exp( j + x i j ) , j = 1 , 2 ,  , m  1 . in contrast to the ppo model , however , the ordered logit model can only estimate one coefficient for each explanatory variable . this coefficient is assumed to be identical for all dichotomizations of the outcome variable ( in this case , stage iv versus ii / iii and stage iii / iv versus ii ) , and thus coefficient j = . this is the proportional odds or parallel regressions assumption . if this assumption is violated , which often happens in practice , estimates are invalid and important differences in the relationships at different thresholds may go unnoticed . the partial proportional odds model relaxes this assumption for variables where it does not hold.31 for these variables , a coefficient is estimated for each dichotomization . the final model was developed by stepwise backward elimination from the full model , which contained all variables . in each step the variable with the largest p - value was eliminated and the model was re - estimated . this process was repeated until all variables had at least a value of p  0.10 . we performed two sensitivity analyses in order to account for the fact that the proportion of patients with very severe copd ( gold iv ) in the dataset was smaller than the proportion of patients with moderate or severe copd . firstly , the final model was re - estimated in a sub - sample with a balanced distribution of patients across all three severity stages . this balance was achieved by using all gold stage iv patients and random draws from patients in stages ii and iii . secondly , the patients in gold stage iii and iv were combined into one group . the probability of being in  severe / very severe was then analyzed in a binary logit model with the same variables as in the final ppo model . all analyses were performed in stata 11.1 , ( statacorp lp , college station , tx ) using the gologit2 command for the ppo model.32 univariate tests were analysis of variance ( anova ) for continuous variables and chi - square for categorical variables . the regression results were used to predict the probabilities of being in each gold severity stage for each patient in the validation dataset . the agreement between the predicted and the observed stage was assessed using kappa statistics.33 for a kappa statistic , a value of 0 indicates that agreement has occurred by chance , whereas a value of 1 indicates perfect agreement . landis and koch suggested the following labels of agreement for ranges of values , which are often quoted : slight ( up to 0.20 ) , fair ( 0.210.40 ) , moderate ( 0.410.60 ) , substantial ( 0.610.80 ) , and almost perfect ( 0.811.0).34 we present kappa statistics unweighted and weighted . in the weighted kappa statistics , patients misclassified in the neighboring category ( ie , ii or iv instead of iii , or iii instead of ii or iv ) count as 0.5 agreement and patients misclassified in a non - neighboring category ( ii instead of iv or vice versa ) as no agreement ) . in addition to the agreement per disease stage , the overall agreement was calculated . for the binary logit model , the c - statistic represents the area under the receiver operating characteristic ( roc ) curve , a plot of sensitivity against specificity . values can range from 0.5 ( no predictive ability ) to 1 ( perfect discrimination ) . lastly , in an analysis of correctly and incorrectly classified patients , the average values of fev1%-predicted were compared per severity stage . patient characteristics are summarized in table 1a , demonstrating that 47% were in gold stage ii , 44% in stage iii , and 9% in stage iv . patients in more severe disease states were more likely to have a lower bmi , to have quit smoking and to suffer from osteoporosis . the mean age of patients in stage iv was 23 years lower than in the less severe categories . the disease stage was not significantly associated with pack - years , the number of co - morbidities , charlson s comorbidity index , and various concomitant diseases . tables 1b and c show that disease stage was associated significantly with almost every type of medication and resource use . patients with more severe copd were more likely to use various types of pulmonary maintenance medication , oxygen at home , exacerbation medication , consult their gp more often ( on a scheduled and unscheduled basis ) , visit the emergency room with resulting admission more often , and were more likely to have been admitted for any reason to the hospital in the last year . the use of statins and other cardiovascular medication was more frequent in patients with less severe disease . there was no difference between the dataset used for developing the algorithm and the validation dataset . in the final model ( table 2 ) , the parallel regression assumption was violated for five variables ( age , gender , xanthine and oxygen use , and two categories of the bmi variable , overweight and obese ) . in these cases , iii / iv versus gold ii and gold iv versus gold ii / iii ) . for all other variables , one coefficient was estimated for both dichotomizations . in the multivariable analysis , a higher risk of more severe copd was significantly associated with being younger , being male , having a lower bmi , having quit smoking , suffering from osteoporosis , using oxygen , courses of oral steroids , having been hospitalized in the previous year , and certain types of respiratory maintenance medication ( long- and short - acting bronchodilators , xanthines , leukotriene modifiers , oral steroids ) . the thresholds were not statistically significantly different from 0 ( p - values : 0.466 and 0.338 , see table 2 ) . however , the first threshold was significantly lower than the second ( 0.25075 versus 0.40289 , p = 0.0422 ) . the following variables were not maintained in the final model : inhaled corticosteroids , mucolytics , cardiovascular medication , courses of antibiotics , most co - morbidities ( except osteoporosis ) , and healthcare resource use ( emergency room , gp consultations ) . sixty three percent of patients in stage ii , 53% of patients in stage iii and 8% of patients in stage iv were correctly identified with the final model ( see table 3a ) . the overall agreement was 53% ( unweighted ) to 77% ( weighted ) , leading to  slight to  fair kappa statistics of 0.151 ( unweighted ) and 0.230 ( weighted ) . the observed values of fev1%-predicted were slightly higher for copd stage ii patients who were correctly classified by the model than for patients who were incorrectly classified ( 59.4 versus 57.7% , respectively ) . in stage iii , the opposite was observed ( fev1%-predicted was 39.4% in the correctly classified versus 41.2% in the incorrectly classified ) , and in stage iv , fev1 was comparable among the correctly and incorrectly classified groups ( 24.9% and 25.6% ) . in the balanced sample ( regression results not presented ) , agreement was lower for stages ii and iii than in the unbalanced sample , but higher for stage iv : 53% , 33% and 57% respectively ( table 3b ) . fair with 0.215 and 0.290 , respectively . in the binary logit model ( table 3c ) , agreement for patients with moderate and severe / very severe copd was 54% and 68% . patient characteristics are summarized in table 1a , demonstrating that 47% were in gold stage ii , 44% in stage iii , and 9% in stage iv . patients in more severe disease states were more likely to have a lower bmi , to have quit smoking and to suffer from osteoporosis . the mean age of patients in stage iv was 23 years lower than in the less severe categories . the disease stage was not significantly associated with pack - years , the number of co - morbidities , charlson s comorbidity index , and various concomitant diseases . tables 1b and c show that disease stage was associated significantly with almost every type of medication and resource use . patients with more severe copd were more likely to use various types of pulmonary maintenance medication , oxygen at home , exacerbation medication , consult their gp more often ( on a scheduled and unscheduled basis ) , visit the emergency room with resulting admission more often , and were more likely to have been admitted for any reason to the hospital in the last year . the use of statins and other cardiovascular medication was more frequent in patients with less severe disease . there was no difference between the dataset used for developing the algorithm and the validation dataset . in the final model ( table 2 ) , the parallel regression assumption was violated for five variables ( age , gender , xanthine and oxygen use , and two categories of the bmi variable , overweight and obese ) . in these cases , a coefficient was estimated for each dichotomization ( gold iii / iv versus gold ii and gold iv versus gold ii / iii ) . for all other variables , one coefficient was estimated for both dichotomizations . in the multivariable analysis , a higher risk of more severe copd was significantly associated with being younger , being male , having a lower bmi , having quit smoking , suffering from osteoporosis , using oxygen , courses of oral steroids , having been hospitalized in the previous year , and certain types of respiratory maintenance medication ( long- and short - acting bronchodilators , xanthines , leukotriene modifiers , oral steroids ) . the thresholds were not statistically significantly different from 0 ( p - values : 0.466 and 0.338 , see table 2 ) . however , the first threshold was significantly lower than the second ( 0.25075 versus 0.40289 , p = 0.0422 ) . the following variables were not maintained in the final model : inhaled corticosteroids , mucolytics , cardiovascular medication , courses of antibiotics , most co - morbidities ( except osteoporosis ) , and healthcare resource use ( emergency room , gp consultations ) . sixty three percent of patients in stage ii , 53% of patients in stage iii and 8% of patients in stage iv were correctly identified with the final model ( see table 3a ) . the overall agreement was 53% ( unweighted ) to 77% ( weighted ) , leading to  slight to  fair kappa statistics of 0.151 ( unweighted ) and 0.230 ( weighted ) . the observed values of fev1%-predicted were slightly higher for copd stage ii patients who were correctly classified by the model than for patients who were incorrectly classified ( 59.4 versus 57.7% , respectively ) . in stage iii , the opposite was observed ( fev1%-predicted was 39.4% in the correctly classified versus 41.2% in the incorrectly classified ) , and in stage iv , fev1 was comparable among the correctly and incorrectly classified groups ( 24.9% and 25.6% ) . in the balanced sample ( regression results not presented ) , agreement was lower for stages ii and iii than in the unbalanced sample , but higher for stage iv : 53% , 33% and 57% respectively ( table 3b ) . fair with 0.215 and 0.290 , respectively . in the binary logit model ( table 3c ) , agreement for patients with moderate and severe / very severe copd was 54% and 68% . firstly , the variables that were independently related to more severe copd defined by lung function were a lower age , male gender , a lower bmi , being an ex - smoker , having osteoporosis , having been hospitalized in the previous year , using oxygen and certain types of respiratory maintenance medication ( long- and short - acting bronchodilators , xanthines , leukotrien modifiers , oral steroids ) . other variables expected to be associated with lung function impairment , such as other resource use variables , were not maintained in the final model . secondly , the performance of the final model was such that the confidence in using the selected variables to adjust for copd severity in the absence of lung function parameters was judged to be limited . of the variables in the final model , the impact of age may be partly due to a healthy survivor effect , whereas long - acting bronchodilators , multiple short - acting bronchodilators , oral steroids and xanthines are clearly indicated for treating more advanced stages of copd . a low bmi , often associated with loss of muscle mass is well known to be more frequent in severe copd and the higher prevalence of osteoporosis might be a side - effect of a long history of corticosteroid use . this made it possible to estimate different coefficients for the probability of being in stage iii / iv over stage ii than for the probability of being in stage iv , if this was required . the parallel regression assumption , which states that the coefficients are equal for both dichotomizations , was violated in five instances . in the final base case model , only 53% of patients were classified in the correct gold stage . especially for patients in stage iv , results from sensitivity analyses with a more balanced sample or a binary logit model were only slightly better . we chose to present both unweighted and weighted kappa statistics . in cases with more than two categories , it is customary to weight the kappa - statistics in order to penalize disagreements in terms of their seriousness  ie , a higher penalty for misclassifying a patient from stage ii as a stage iv patient than as a stage iii patient  whereas unweighted kappa treats all disagreements equally . however , in our study all disagreements were considered serious and the weighted overall agreement might give an overly favorable impression . based on all values of the kappa statistics , the agreement between predicted and observed gold stages can be characterized as slight to fair . firstly , the cut - off points between gold stages are inevitably somewhat arbitrary and artificial , especially as the decline in fev1 is a continuous process . a patient with a fev1%-predicted of 49% ( gold stage iii ) is probably less similar to a patient with a value of 31% ( also gold stage iii ) than to someone with a value of 51% ( gold stage ii ) . if this were an important explanation for our prediction results , misclassified patients in stage iv should have markedly better lung functions  closer to stage iii  than the correctly classified patients , while the opposite should be true for misclassified patients in stage ii . however , the actual differences in lung function between the correctly and incorrectly classified patients in our analysis were small ( 1.74%-point for stage ii and 0.66%-point for stage the second explanation concerns the source of our data , which was a clinical trial . it is conceivable that patients with very severe disease and many symptoms were less willing to participate in the trial . this problem would be expected to occur particularly when trying to predict gold stage iv . however , our models did not perform well at distinguishing moderate from severe patients either . in the balanced sample analysis , misclassifications occurred equally often in each gold stage . furthermore , the proportion of very severe patients in the sample does not appear to be low compared to the proportion in the general population of copd patients in finland,35 the netherlands,36 greece,37 the united kingdom,38 and a combination of european and north - american countries.39 moreover , a broad range of patients was allowed to participate . for examples , patients were permitted to use all respiratory medications concomitantly during the trial except inhaled anticholinergic drugs , thus closely resembling routine care . altogether , this protocol makes the uplift a suitable trial for this study . the third explanation would be the far from perfect association of resource use , which is often driven by symptoms and exacerbations , with gold stage . patients with a relatively good lung function do not necessarily experience fewer symptoms than patients with worse lung function . indeed , the uplift sample contains a non - negligible number of patients with very severe copd who do use little or no maintenance medication , as well as patients with moderate disease who use four or five different types of medications . overall , we observed an independent association of gold stage with several respiratory medications , but not with other types of resource use such as er visits and physician consults in our multivariable model , which was not expected beforehand . in univariate analyses , patients with more severe copd were more likely to have higher copd - related resource use ( almost all types of medication use , courses of oral steroids and antibiotics , scheduled and unscheduled gp visits and hospital admissions with and without visits to the emergency room ) , not to use statins and other cardiovascular medication , to be younger , male , underweight , ex - smoker , and suffer from osteoporosis . several database studies have attempted to adjust for possible differences in copd severity in the absence of lung function data . in an article comparing the assessment of copd patients in the uk general practice research database with the clinical opinion of the patient s gp , sin and tu20 assessed the effects of ics use on mortality and applied medication use , er visits and physician services as proxies of disease severity . they considered hospital admissions and medication use to be proxies for disease severity . in a study relating copd severity with cardiovascular disease , curkendall et al23 assumed that copd severity could be defined as the likelihood of being hospitalized , given the relationship with mortality . they concluded that this probability of hospitalization is associated with medication and oxygen use , previous hospitalizations , recent exacerbations , pneumonia and lung emphysema . griffin et al24 assessed the effects of tiotropium compared to combined ipratropium and salbutamol on exacerbations and hospitalizations . based on the current analyses , we conclude that the variables used in these database studies can not be relied upon to adequately adjust for copd severity in terms of lung function . data from a well - controlled trial setting indicated that copd severity defined by lung function thresholds can not be reliably predicted from patient characteristics and their previous healthcare use . this limitation should be considered when interpreting findings from database studies , and additional research should explore other methods allowing accounting for copd severity .
purposewhen comparing chronic obstructive lung disease ( copd ) interventions in database research , it is important to adjust for severity . global initiative for chronic obstructive lung disease ( gold ) guidelines grade severity according to lung function . 
 most databases lack data on lung function . 
 previous database research has approximated copd severity using demographics and healthcare utilization . 
 this study aims to derive an algorithm for copd severity using baseline data from a large respiratory trial ( uplift).methodspartial proportional odds logit models were developed for probabilities of being in gold stages ii , iii and iv . 
 concordance between predicted and observed stage was assessed using kappa - statistics . 
 models were estimated in a random selection of 2/3 of patients and validated in the remainder . 
 the analysis was repeated in a subsample with a balanced distribution across severity stages . 
 univariate associations of copd severity with the covariates were tested as well.resultsmore severe copd was associated with being male and younger , having quit smoking , lower bmi , osteoporosis , hospitalizations , using certain medications , and oxygen . after adjusting for these variables , co - morbidities , previous healthcare resource use ( eg , emergency room , hospitalizations ) and inhaled corticosteroids , xanthines , or mucolytics were no longer independently associated with copd severity , although they were in univariate tests . 
 the concordance was poor ( kappa = 0.151 ) and only slightly better in the balanced sample ( kappa = 0.215).conclusioncopd severity can not be reliably predicted from demographics and healthcare use . 
 this limitation should be considered when interpreting findings from database studies , and additional research should explore other methods to account for copd severity .