in being responsible for about 350,000 new cases and 150,000 deaths yearly in france , cancer is a major health problem and its surveillance the utmost public health concern .
regarding surveillance , francim , the french network of cancer registries , is responsible for exhaustive collections of cancer cases in 10 to 14 french dpartements ( depending on the cancer type ) corresponding to 15% to 20% of the french population .
however , estimating epidemiological indicators at the scale of whole dpartements over the country is necessary not only to reveal etiological factors and geographical or social discrepancies but also to plan the needs in terms of medical resources ( prevention , treatment , and surveillance ) . over the previous ten years or so , francim and
the department of biostatistics of the hospices civils de lyon have been providing national estimations of cancer incidence [ 1 , 2 ] .
their usual approach to produce these estimations is to use registry incidence data together with cpidc mortality data ( centre d'epidmiologie sur les causes mdicales de dcs ) .
the principle is to calculate a mean ratio between incidence and mortality in an area covered by cancer registries then use that ratio with national mortality data to derive an estimation of national incidence .
while the mean ratio estimated in a registry area can be reasonably considered as representative of the ratio for the whole country , the same is not true at the level of a single dpartement because this ratio may be highly variable between dpartements and because identical incidence values do not necessarily lead to identical mortalities .
indeed , a great number of factors are able to affect patient survival and generate heterogeneity of the ratio between dpartements : differences in patient management ( diagnostic or therapeutic procedures ) , prevention or screening policies , or compliance of the population with these policies .
therefore , because incidence and mortality can not be used to provide dpartement - specific incidence estimations , a new approach should be sought for .
one interesting source of data with a nationwide coverage has been recently used together with registry data to estimate dpartement - specific cancer incidence : the hospital database of the programme de mdicalisation des systmes d'information mdicale ( pmsi ) [ 46 ] .
this medicoadministrative database is held to help manage health institutions and provide budget estimates .  
a previous work   has discussed the problem of using hospital stays from pmsi data to estimate dpartement - specific incidence of breast cancer .
however , the constant improvement of the quality of patient identification in pmsi data using a single - patient identifier makes it possible now to use patient - specific rather than stay - specific data   in the present paper , our objective was to estimate dpartement - specific incidence of colon - rectum , breast , kidney , and ovary cancers for 2007 using mean ratios of pmsi - extracted cases to registry - extracted ( incident ) cases .
the french agence technique de l'information sur l'hospitalisation ( atih ) made available its data on all short stays in all health institutions over france for 20022007 . in our analyses
, we kept the variables related to personal characteristics ( sex , age , and code of the residence area ) , hospital stay ( stay number and principal diagnosis according to the international classification of diseases ( icd-10 ) , and medical procedures according to the catalogue des actes mdicaux ( cdam ) until 2004   and to the new classification commune des actes mdicaux ( ccam )   from 2004 to 2007 , plus an anonymous alphanumerical patient identifier [ 11 , 12 ] to allow chaining of hospital stays of the same patient in successive institutions .
that identifier is systematically generated by procedure foin ( fonction d'occultation des informations nominatives ) in all french health institutions since 2001 .
the corresponding cim 10 codes were c18 to c21 for colon - rectum cancer , c50 for breast cancer , c56 and c57.0 to c57.4 for ovary cancer , and c64 to c66 plus c68 for kidney cancer .
the second algorithm extracted stays with cancer as principal diagnosis and with surgical procedures for cancer . using cdam and ccam codes , the latter extraction considered 95 procedures for colon - rectum cancer , 31 for breast cancer , 114 for ovary cancer , and 44 for kidney cancer .
after each type of extraction , hospital stays were ordered by their serial numbers to spot the first stay of each patient , then only these stays were kept for analysis .
these stays were then counted over each dpartement by age group : ten groups for colon - rectum cancer ( 1544 , 4549 , 5054 ,  , 8084 , and 85  yrs ) , eleven groups for kidney cancer in women ( 1539 , 4044 , 4549 ,  , 8084 , and 85  yrs ) , and thirteen groups for each of breast cancer , ovary cancer , and kidney cancer in men ( 1529 , 3034 , 3539 ,  , 8084 , 85  yrs ) .
francim network made available the data on incident cancer cases registered in 2004 ( the most recent and checked data set when the present study was initiated ) .
cancer sites were determined according to the international classification of diseases for oncology , third version ( icd - o-3 )   and corresponded to invasive tumors .
these codes were : c18 to c21 for colon - rectum cancers , c50 for breast cancer , c56 and c57.0 to c57.4 for ovary cancer ( excluding morphological codes 8442/3 , 8451/3 , 8461/3 , 8462/3 , 8472/3 , and 8473/3 ) , and c64 to c66 plus c68 for kidney cancer .
the cancer registries used were those of eleven dpartements : calvados , cte d'or , doubs , hrault , isre , loire - atlantique , bas - rhin , haut - rhin , sane - et - loire , somme , and tarn .
incidents cases of cancers were counted by the same age group as for hospital stays .
our approach was to model , in function of age , the pmsi cases / incident cases ratio ; that is , the ratio of the number of patients with hospital stays for cancer present in the pmsi database to the number of incident cases present in the registries .
this ratio was obtained from dpartements where the two sources of information exist ( i.e. , dpartements with a registry ) .
it was then applied to pmsi data of dpartements without registry , in order to estimate cancer incidence in these areas .
this ratio was obtained from registry and pmsi data of year 2004 ; it was then applied to pmsi data for 2007 to derive dpartement - specific incidence for 2007 .  
it is a calibration method where incidence , as obtained from cancer registries , is considered as reference or  true 
more precisely , the ratio is modelled as a function of age ( effects smoothed using cubic regression spline ) and dpartement ( considered as random effect ) .
further , to estimate the pmsi cases / incident cases ratio , a data quality criterion was required : that the chaining rate be greater or equal to 95% .
this chaining rate was defined as the ratio of the number of stays with personal identifier to the total number of stays .  
the analysis was carried out for the four cancer sites in women but only for colon - rectum cancer and kidney cancer in men . whenever the number of cancer cases was small , it was not possible to take into account between - dpartement variability in areas with registries .
we had then to sum up data from several dpartements and calculate the pmsi cases / incident cases ratio per age group over the entire zone .
the overall national incidence was calculated by summing all dpartement estimations . for validation , we compared these national incidence values to francim values previously obtained by modeling the incidence / mortality ratio [ 1 , 2 , 13 ] .
validation was carried out through three steps . in step 1 , the pmsi cases / incident cases ratio was calculated over all ages using each algorithm per cancer site - sex combination and dpartement .
indeed , dpartement - specific estimations stemming from this approach are invalid unless the ratio is homogeneous between dpartements ( always > 1 or always < 1 ) . in step 2 , the observed ratio for a given age group and dpartement was graphically compared to the modelled mean ratio over all dpartements with registries ( figure 1 ) .
however , to be applicable to all dpartements , the modelled mean ratio should not suffer a  dpartement effect  . in the presence of this effect , the observed ratios per age groups tend to be systematically higher or lower than the mean ratio whereas in its absence
, the observed ratios are distributed around the mean ratio . for cancer sites with high incidence ( breast , colon - rectum ) , step 3 was a cross - validation .
the incidence in a given dpartement with registry is estimated from pmsi data together with the pmsi cases / incident cases ratio obtained by a model from which the data of that dpartement were excluded . a comparison between the number of observed cases and the number of predicted cases yields a prediction error ( pe ) . under hypothesis h0 of a correct prediction ,
the pe obeys a  rule whose degree of freedom is equal to the number of age groups .
a 5%  -risk was adopted to set the critical value for rejecting h0 . for cancer with low incidence ( ovary , kidney ) , cross - validation
could not be used , because it was difficult to determine the statistical distribution of the pe and so only graphical validation was done .
in addition , a comparison between the total numbers of observed and predicted cases was carried out (  with one degree of freedom ) . a relative error ( re ) was calculated as the difference between the observed and predicted cases divided by the number of observed cases .
the results were mapped as standardized incidence ratios ( sir ) using software mapinfo ( version 7.0 ) .
in 2004 , the pmsi database was including 20,721,587 hospital stays of which 718,044 stays had cancer as principal diagnosis .
the chaining rate of those stays was 96.4% . in 2007 , there were 21,201,102 stays of which 721,823 had cancer as principal diagnosis and the chaining rate was 99% . to illustrate our validation procedure
the analyses relative to breast cancer were plainly detailed whereas those relative to the other cancer sites were less detailed .
six registries were selected to estimate the incidence of breast cancer ( table 1 ) .
four dpartements with registries were excluded because , in 2004 , their chaining rate was too low ; it ranged between 49.7% and 94.9% . 

table 1 shows the total number of cases by cancer site and case - extraction algorithm as well as the three steps of estimate validation : homogeneity of the pmsi cases / incident cases ratio through dpartements with ratio < 1 , the dpartement effect , and the cross - validation . regarding the homogeneity of the ratio ,
the results show that algorithm 2 was inadequate . indeed , whereas with algorithm 1 the number of pmsi cases was always higher than the incident cases in all dpartements , with algorithm 2 , the former number was sometimes higher ( 2 dpartements ) and sometimes lower ( 4 dpartements ) than the latter .   regarding the dpartement effect , figure 1 presents , for each dpartement with registries and for algorithm 1 , the observed ratio by age group as well as the modelled mean ratio over all dpartements with registries .
the absence of heterogeneity in that ratio between dpartements is illustrated by the fact that there was no dpartement in which the ratio per age group was systematically higher or lower than the modelled mean ratio ( though dpartement calvados tend to be systematically higher than the modelled mean ratio ) .
it can be therefore concluded that there is no dpartement effect for breast cancer with algorithm 1 ( the variance of the random effect was small : 0.023 ) .
regarding the cross - validation , table 2 presents the detailed results for breast cancer . with algorithm 1 ,
the differences between the observed and the predicted number of cases per age group ( pes ) were small whereas with algorithm 2 , two dpartements displayed large differences especially concerning the last age group .
furthermore , the difference between the total observed and the total predicted cases (  ) was small with algorithm 1 : thus , this algorithm 1 may be reliably used for breast cancer estimates . 

table 3 shows the national estimations obtained by adding the estimations obtained from all dpartements by algorithm 1 as well as the national estimations elaborated by francim   through the use of the incidence / mortality ratio [ 1 , 2 , 13 ] . comparing these two estimations
the world age - standardized incidence rate ( wasr ) was 98.6 cases per 100,000 and varied between dpartements from 71.4 to 127.1 ( table 4 ) .
to estimate dpartement - specific incidence of colon - rectum cancer , we used data from eight registries .
data from three dpartements with registries were excluded because the chaining rate in 2004 was too poor ( 72.6% to 94.9% ) .
irrespective of sex , the pmsi cases / incident cases ratios obtained with algorithm 1 were homogeneous ( the variance of the random effect was small : 0.024 ) , which was not the case with algorithm 2 for which there was a dpartement effect in both gender .
in addition , contrary to algorithm 1 , cross - validation invalidated the estimations made with algorithm 2 in two dpartements .
thus , we present only estimations made with algorithm 1 . the national incidence for 2007
the estimated national wasr was 38 cases per 100,000 in men and 24.8 in women . at the national level ,
our estimations were in high agreement with those of francim ( 40.8 cases per 100,000 in men and 24.8 cases per 100,000 in women ) . among dpartements ,
the wasr ranged from 23.1 to 48.9 in men and from 17.4 to 31.3 in women .
three dpartements with registries were excluded because the chaining rate was too poor ( 80% to 93.5% ) .
the pmsi cases / incident cases ratio carried out with algorithm 1 yielded homogeneous ratio between dpartements , and no dpartement effect was observed .
however , as already mentioned in section methods , cross - validation could not be carried out to validate estimate from algorithm 1 and thus , because of difficult formal validations of the dpartement - specific estimations , we only present national estimations based on this algorithm 1 .
the national incidence was 4,637 cases , which corresponds to a wasr of 8.6 cases per 100,000 .
thus , eight dpartement registries were considered for men ( one registry excluded because of a 92.5% chaining rate ) and five registries for women ( four registries excluded because of chaining rates ranging between 89.7% and 94.4% ) .   here
too , only a graphical validation of the estimations could be carried out and , because of difficult validations of dpartement - specific estimations , only national estimations are given .
as for ovary cancer and irrespective of sex , algorithm 1 performed better than algorithm 2 . in 2007 ,
the national incidence of kidney cancer was estimated at 6,775 in men and 3,273 in women ; the national wasr was 13.3 per 100,000 in men but much lower ( 5.2 per 100,000 ) in women .
dpartement - specific estimations for colon - rectum and breast cancer are shown in table 4 and sir maps of these cancers are shown in figures 2 , 3 ,   4 .
these maps were not constructed for ovary and kidney cancers because of difficult formal validations of dpartement - specific incidence .  
however , for colon - rectum cancer in men , the southwest was a low incidence area .
one well - marked low - incidence area for breast cancer was the southwest quadrant of france .
to estimate the incidence of the four cancers in each dpartement , case extraction from pmsi database used two algorithms .
algorithm 1 targeted all hospitalized cancer patients ; that is , those whose principal diagnosis is cancer because of positive laboratory tests , metastasis staging , cancer - related procedures , or sudden potentially fatal progression ( exacerbation or relapse ) .
this is common in incidence estimations based on hospital data [ 1517 ] and was confirmed here : the number found in pmsi data was higher than incident cases found in the registries ( truer with algorithm 1 than with algorithm 2 ) .
nevertheless , the pmsi cases / incident cases ratio obtained with algorithm 1 seemed somehow stable between dpartements , which allowed estimations of dpartement - specific incidence . in contrast , algorithm 2 that used initial surgical procedures was more selective ; it extracted a closer number of pmsi patients to the number of incident cases than algorithm 1 .
thus , the pmsi cases / incident cases ratio was more heterogeneous between dpartements ( at different degrees according to the cancer site ) .
besides , cross - validation revealed that incidence estimations in some dpartements with registries were not valid .
indeed , with algorithm 2 , the critical value of the prediction error was crossed for colon - rectum cancer in both sexes and for breast cancer in women . in sum
, the simpler and less selective algorithm 1 was more adequate than algorithm 2 to estimate dpartement - specific incidence . at the national level ,
our estimations were in agreement with francim projections , except for breast cancer ( our estimations were lower ) .
this was expected because estimations using incidence / mortality ratios do not take into account the recent trend towards a slow decline of breast cancer incidence in france   and in other countries [ 19 , 20 ] .
our estimation made in 2007 at 50,578 new yearly cases in france seemed thus more realistic than the 52,492 cases stemming from francim projections .
this comparison is interesting because it shows that in average , over all dpartements , our approach leads to reliable estimations .
finally , our graphical validation can not constitute a formal validation because , in borderline situations , one can not claim a  dpartement effect  .
the graphical method would be only as a tool to detect important departures from the model assumptions .
pmsi data concern only hospitalized patients ; thus , our method does not apply to cancers such as skin melanomas or basocellular cancers , which are usually treated early without hospitalization . besides , the method supposed identical treatment choices in all dpartements ; which motivated the choice of the four cancer sites . in two successive articles on colon and rectum cancers [ 22 , 23 ] , phelip et al .
have shown that surgical resection was performed in 90% of cases without significant geographical variation between dpartements .
a high variability in treatment choices leads to problems with case - extraction from the pmsi database .
for example , in aged men , prostate cancer can be treated by surgery , radiotherapy , or even hormone therapy alone ( without hospitalization ) and , in bladder cancer , the surgical treatment depends on the stage . if surgery is avoided , the lack of  surgery for cancer  in the pmsi database leads to missingness of pmsi cancer cases .
in addition , our method lays fundamentally on the hypothesis that within a given age group , the pmsi cases / incident cases ratio is constant between dpartements whereas several factors may affect that rate .
indeed , in a previous work , an insufficient chaining led to consider the number of stays rather than the number of patients .
consequently , the between - dpartement variability in the mean number of stays per patient  essentially due to very different hospitalization policies and coding practices  prevented correct estimations . in 2002 ,
the chaining rate was quite low ( 92% ) but improved up to 2004 ( nearly 96% ) because of the implementation of a 
the use of a single - patient identifier allows keeping a single record per patient whatever the number of hospital stays and insures a better homogeneity of pmsi data .
variability may also stem from various coding habits in different health institutions because of ignorance or misinterpretation of coding rules .
a wide interinstitution variability of the pmsi cases / incident cases ratio may lead to a wide between - dpartement variability , especially between dpartements with few health institutions . besides , the pmsi cases / incident cases ratio reflects the fact that for a given incidence level , the prevalence may vary between dpartements . indeed
, different rates of cancer - specific survivals between dpartements affect hospital prevalence ; thus , the number of pmsi cases .
in fact , survival in a given dpartement may vary with the presence / absence of systematic screening , the existence or not of a reference health institution , and the educational   and socioeconomic levels of the population .
if survival is high because of complete cures , hospital prevalence and pmsi cases will decrease , which will underestimate incidence , but if survival is high but associated with more procedures , hospital prevalence and pmsi cases will increase , which will overestimate dpartement - specific incidence .
the impact of different survival rates between dpartements on hospital prevalence is complex to seize and undoubtedly dependent on the cancer site under study .
another implicit hypothesis in our method is a constant pmsi cases / incident cases ratio from 2004 to 2007 .
this is plausible because data quality of both sources over that short period was deemed constant despite improvements in cancer therapies that would have changed the prevalence of cancer .
another effect of time is the change in the national standards of coding pmsi data .
for example , coding palliative care as  cancer  has been replaced by a specific icd-10 code ( z51.5 ) .
more recently , in 2009 , the rules for the choice of the principal diagnosis have changed ; the impact of that change should be evaluated .
nevertheless , that impact would be limited within the context of the four cancers under study here .
improvements of the present method are possible . a better followup of the same patients over several years would exclude a number of prevalent cases .
ideally , if all prevalent cases were excluded , the ratio would be interpreted as the proportion of hospitalized incident cases and would be no more affected by different prevalence in different dpartements .
another improvement would be to add data from the health insurance ( affections longue dure ( ald30 ) database of caisses d'assurance maladie ) ; a feasibility study of that possibility is underway .
using an adequate method , it seems now possible to estimate dpartement - specific incidence of some cancers for a given year .
nevertheless , this validation is only partial because dpartement - specific estimations will still suffer the basic assumption of similar coding practices in all hospitals .
 
 objective . use of cancer cases from registries and pmsi claims database to estimate dpartement - specific incidence of four major cancers . methods . case extraction used principal diagnosis then surgery codes . 
 pmsi cases / registry cases ratios for 2004 were modelled then dpartement - specific incidence for 2007 estimated using these ratios and 2007 pmsi cases . 
 results . for 2007 , 
 only colon - rectum and breast cancer estimations were satisfactorily validated for infranational incidence not ovary and kidney cancers . for breast , 
 the estimated national incidence was 50,578 cases and the incidence rate 98.6 cases per 100,000 person per year . for colon - rectum , incidence was 21,172 in men versus 18,327 in women and the incidence rate 38 per 100,000 versus 24.8 . for ovary , 
 the estimated incidence was 4,637 and the rate 8.6 per 100,000 . for kidney , incidence was 6,775 in men versus 3,273 in women and the rate 13.3 per 100.000 versus 5.2 . 
 conclusion . 
 incidence estimation using pmsi patient identifiers proved encouraging though still dependent on the assumption of uniform cancer treatments and coding . 
1. Introduction
2. Materials and Methods
3. Results
4. Discussion
5. Conclusion