 Research article
 Open Access
 Open Peer Review
 Published:
BonferroniHolm and permutation tests to compare health data: methodological and applicative issues
BMC Medical Research Methodology volume 18, Article number: 81 (2018)
Abstract
Background
Statistical methodology is a powerful tool in the health research; however, there is wide accord that statistical methodologies are not usually used properly. In particular when multiple comparisons are needed, it is necessary to check the rate of false positive results and the potential inflation of type I errors. In this case, permutation testing methods are useful to check the simultaneous significance level and identify the most significant factors.
Methods
In this paper an application of permutation tests, in the medical context of Inflammatory Bowel Diseases, is performed. The main goal is to assess the existence of significant differences between Crohn’s Disease (CD) and Ulcerative Colitis (UC). The Sequentially Rejective Multiple Test (BonferroniHolm procedure) is used to find which of the partial tests are effectively significant and solve the problem of the multiplicity control.
Results
Applying NonParametric Combination (NPC) Test for partial and combined tests we conclude that Crohn’s Disease patients and Ulcerative Colitis patients differ between them for most examined variables. UC patients compared with the CD patients, have a higher diagnosis age, not show smoking status, proportion of patients treated with immunosuppressants or with biological drugs is lower than the CD patients, even if the duration of such therapies is longer. CD patients have a higher rate of rehospitalization. Diabetes is more present in the subpopulation of UC patients. Analyzing the Charlson score we can highlight that UC patients have a more severe clinical situation than CD patients. Finally, CD patients are more frequently subject to surgery compared to UC. Appling of the Bonferroni Holm procedure, which provided adjusted pvalues, we note that only nine of the examined variables are statistically significant: Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Rehospitalization, Gender and Duration of Immunosoppressive Therapy. Therefore, we can conclude that these are the specific variables that can discriminate effectively the Crohn’s Disease and Ulcerative Colitis groups.
Conclusions
We identified significant variables that discriminate the two groups, satisfying the multiplicity problem, in fact we can affirm that Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Hospitalization, Gender and Duration of Immunosoppressive Therapy are the effectively significant variables.
Background
Statistical methodology is an useful and powerful tool in the medical scientific research; therefore, an important increase in the use of statistical methods having been documented in most medical journals [23, 32, 46].
In recent years permutation tests increased in applications to solve complex multivariate problems. Permutation tests are essentially of an exact nonparametric nature in a conditional context, where conditioning is on the pooled observed data set which is often a set of sufficient statistics in the null hypothesis. Whereas, the reference null distribution of most parametric tests is only known asymptotically [39].
There are, however, many complex multivariate problems (quite common in biostatistics, clinical trials, engineering, the environment, epidemiology, experimental data, industrial statistics, pharmacology, psychology, social sciences, etc.) that are difficult to solve outside the conditional framework and in particular outside the method of Non Parametric Combination (NPC) of dependent permutation tests [38].
Permutation tests and bootstrap methods have very wideranging applications, both share a common potential drawback: as dataintensive resampling methods, both can be runtime prohibitive when applied to large or even mediumsized datasets. The data explosion over the past few decades has made this a common occurrence and it highlights the increasing need for faster and more efficient permutation tests and bootstrap algorithms [31]. The permutation test essentially works by combining two important principles: exchangeability and conditioning.
The main goal of this paper is, applying the NPC test methodology, to study a specific medical problem with a large amount of patients (about 1700) in order to assess the existence of significant differences between subjects affected by two Inflammatory Bowel Diseases (IBD); in particular, Crohn’s Disease (CD) and Ulcerative Colitis (UC), with reference to a great number of variables. In this case we are in presence of an authentic real complex problem to be solved; for its solution, the permutation methods are better than the ordinary parametric methods because do not require strong assumptions that are extremely difficult to justify. Since several variables are considered, we also propose an application of the BonferroniHolm procedure for the multiplicity control. In the paper theoretical, methodological and applied aspects [44] have been fruitfully integrated with specific competences from medicine field [33].
The medical context: IBD
The inflammatory bowel diseases (IBD) are chronic inflammatory diseases of the intestinal mucosa; they include only Crohn’s Disease (CD) and Ulcerative Colitis (UC).

Crohn’s Disease can affect the entire gastrointestinal tract, from mouth to anus. In about 90% of cases, the disease mostly affects the last part of the small intestine (ileum) and the colon. It is characterized by intestinal ulcers, often alternating with stretches of healthy gut, which, if not properly treated, can lead to complications (such as stenosis or fistula) that may require surgery. Immunosuppressive therapy and regular monitoring are used to control the disease and its progression in most cases.

Ulcerative Colitis primarily affects the rectum and may involve part or all of the colon. The main clinical symptoms are diarrhoea, often with blood and mucus, and abdominal pain. The course of the disease is characterized by the alternation of acute episodes followed by periods of clinical remission. The medical therapy of this disease is based on administration of antiinflammatory drugs and immunosuppressants. If not properly treated, chronic inflammation can lead over time to irreversible alterations of intestinal cells with the possible development of cancerous lesions. In rare cases (refractory to medical therapy) it is necessary to make a total colectomy surgery.
The causes of IBD are not yet clear. However, most experts agree that several factors may play a causal role in the disease: genetics and, therefore, familiarity is clearly implicated in the disease; in fact, in 20% of cases, individuals with IBD have a first degree relative (up to first cousins) who suffers from ulcerative colitis or Crohn’s disease; other causes are abnormal reactions of the immune system and, as last, environmental factors. Although the exact cause of IBD is not clear, there are certain triggering factors that can cause a worsening of symptoms.
These include

stress (in some subjects the emotional stress can lead to an exacerbation of symptoms);

recently exposure to some types of antiinflammatory drugs (FANS) or antibiotics;

intake of some foods;

smoke.
It is estimated that in Italy about 200,000 people are now suffering from these diseases. The diagnosis of new cases in the last 10 years and the number of patients increased by about 20 times. IBD hit with the same frequency the two sexes, with a clinical onset that is placed between 15 and 45 years. It is important to emphasize that neither the UC or CD are contagious. The two diseases are different, even if they affect the same apparatus. Therefore, a statistical comparison between patients affected by CD and UC is very interesting, from a medical and scientific point of view, in order to assess the differences between them.
Methods
Permutation tests: The reasons
Parametric tests usually imply an approach to the hypothesis test problem that require a series of stringent hypotheses, which are often in practice difficult to justify, particularly in medical research [49, 50]. These assumptions are sometimes arbitrarily established. Generally, without any justification, biomedical studies assume:

a)
multivariate normality;

b)
random sampling;

c)
homoschedasticity;

d)
allocation to treatment is independent.
In other words, the concept that “all models are wrong but some are useful” is often adopted without an adequate critical spirit so that one can be confident that the resulting approximation can be considered acceptable for the specific problem. Conversely, nonparametric statistical tests try to keep assumptions at a lower level, possibly avoiding those that are hard to justify. By doing so, they rely on less stringent and more realistic foundations and are intrinsically robust.
Permutation tests: The methodology
In this section we introduce the theoretical aspects of Non Parametric Combination (NPC) test, based on permutation solution [36]. Permutation tests [12] represent an effective solution for problems concerning the testing of multiple hypotheses, that are difficult or even impossible to face in a parametric context. This multivariate procedure allows to reach effective solutions concerning problems of multidimensional hypotheses verifying by nonparametric permutation inference [34]; it is used in different application fields that concern verifying of multidimensional hypotheses with a complexity that cannot be managed in parametric context [43].
In comparison to the classical approach, NPC test is characterized by several advantages:

it does not require normality and homoschedasticity assumptions ([28]; Janssen A. [27];

it draws any type of variable [35];

it assumes a good behaviour also in presence of missing data; “without relevant loss of information we may remove from the permutation sample space, associated with the whole data set, all data permutation in which the actual sample size of really observed data are not sufficient for approximation. We must establish a kind of restriction on the permutation space, provided that this restriction does not imply biased effects on inferential conclusion”. The missing data can be missing at random (MAR) or not missing at random (NMAR). “The missing data are missing at random (MAR), if the conditional probability of the observed pattern of missing data given the missing data and the value of the observed data is the same for all possible values of the missing data. If the missing data are missing not at random (MNAR), then in order to make valid parametric inferences, the missing data process must be properly specified. The specification of a model which correctly represents the missing data process seems the only way to eliminate the inferencial bias caused by nonresponses in a parametric framework. In the literature, various models have been proposed, most of which concern cases in which nonresponses are confined to a single variable.” ([36], pp. 232–243). We can state that the permutation analysis can be run when there is missingness and is valid when we have missing completely at random (MCAR) data. So, NPC test allows to ignore missingness by removing all unobserved units from the data set and to obtain exact permutation solutions;

it is powerful in presence of low sampling size [9];

it resolves multivariate problems without the necessity to specify the dependence structure among variables [5, 6, 20];

it allows stratified analyses;

it allows to test multivariate restricted alternative hypothesis (to verify the directionality for a specific alternative hypothesis);

it solves problems in which the number of observed subjects is smaller than that of variables [17].
The NPC method is optimal when you want to identify any different patterns between the layers. It allows to realize the control of possible confounding factors using data poststratification techniques. For the control of these factors, which is performed by randomization in clinical trials, an observational context is used in the socalled poststratification. Furthermore, this methodology can also be used with heterogeneous response variables. The NPC method has proven to be robust in the presence of heterogeneity [3].
All these properties make NPC test very flexible [2, 24, 54] and widely applied in several fields; in particular we cite recent applications in medical context [4, 7, 8, 10, 16, 25, 48, 53] and in genetics [14].
By means of the mentioned procedure, it is preliminarily possible to define a set of K onedimensional permutation test, denominated partial test, through which the marginal contribution of every responsevariable can be examined while comparing groups.
The partial tests are nonparametrically combined through CMC (Conditional Monte Carlo) procedure in combined tests, using an opportune combination function (generally Fisher, Tippett or Liptak); these tests globally verify the existence of differences among the multivariate distributions in the groups.
Let us suppose that K variables are observed on two groups (c = 1,2)of n_{c} subjects each. So, the observed data are X = (X_{icu}, i = 1,...,K; c = 1,2; u = 1,...,n_{c}).
According to Roy’s UnionIntersection notation [45], the null hypothesis states the distributional equality in of two Kdimensional variables, that is
where a breakdown into K subnull hypotheses is emphasized. Indeed, global H_{0} is true if all K subnull are jointly true. The alternative is.
which is true when at least one subalternative is true.
The distributional equality stated by H_{0} implies that the observed data vectors are exchangeable between two groups. Without loss of generality, we suppose that for each subhypothesis H_{0i} against H_{1i} there is a suitable partial permutation test T_{i} assumed to be significant for large values.
The system of hypotheses is set in such a way that the related partial tests are jointly processed, so that they can be combined nonparametrically by taking into account their underlying dependence structure within the nonparametric combination method (NPC). We notice that, especially when the number of variables is large, the underlying dependence structure can be more complex than pairwise linear, as it is common described by multivariate Gaussian distribution. So, it is impossible to deal with it by proper estimators of all related regression coefficients, the number and type of which are typically unknown. Thus, it must be worked out nonparametrically. This implies turning to the permutation testing principle and specifically to the NPC.
It is worth noting to observe that permutation tests enjoy several important properties. Among these we underline:

a)
the similarity, that is the rate of rejection of H_{0}, when it is true, is α uniformly for all possible sample data and independently whichever the underling distribution;

b)
under the alternative, the rejection rate of H_{0} is not smaller than α uniformly for all sample data and all underlying distributions, which imply a form of uniform unbiasedness.
The analysis was performed using Methodologica Srl (2001) NPC Test: Statistical Software for Multivariate Permutation Tests (Methodologica Copyright). In the calculation of raw pvalue 10,000 permutations were implemented.
The BonferroniHolm procedure
The Bonferroni  Holm procedure [26] allows to solve the problem of multiple comparisons [1]; it provides control are the family wise error rate (the probability of witnessing one or more Type I errors), by adjusting the rejection criteria for each hypothesis, and offers a simple method, uniformly more powerful than the classical Bonferroni correction. It works as follows:

1.
all pvalues are sorted from smallest to largest. Let’s indicate with K the number of the pvalues;

2.
if the first pvalue is greater than or equal to α/K, the procedure is stopped and no pvalues are significant. Otherwise, we go on.

3.
the first pvalue is declared significant and afterwards the second pvalue is compared to α/(K1). If the second pvalue is greater than or equal to α/(K1), the procedure is stopped and no further pvalues are significant. Otherwise, we go on until the ith ordered pvalue is such that:
p _{ (i) ≥} α /(Ki + 1).
BonferroniHolm procedure is the most widely recommended way to reduce the apparent significance of effects. The great advantage with the sequentially rejective Bonferroni test (as well as with the classical Bonferroni test) is its flexibility [47]. There are no restrictions on the type of tests, the only requirement being that it should be possible to calculate the obtained level for each separate test.
Results
In a multicenter retrospective observational study, we investigated the disease occurrence and course in the first three years in 1722 patients followed the Gastroenterology Unit of several Hospitals, located in the Italian territory: Bari, Cagliari, Catania, Desio (Monza and Brianza), Florence, Messina, Milan, Naples, Padua, Palermo, Rome, San Giovanni Rotondo (Foggia).
The data, deriving from the various hospital centers, were organized in a single dataset by Prof. Walter Fries, Director of Gastroenterology Unit of the University Hospital “G. Martino” in Messina (see [21, 22]). The distribution of outcomes does not vary by center because we verified the condition of equality among the means of the covariates in the different centers, applying the Analysis of Variance (ANOVA) test; it provided nonsignificant results for all variables, denoting the existence of similarities among the means. Before applying the NPC test methodology we also assessed possible heterogeneity or homogeneity in the data, deriving from the different centers through the application of Levene’s test; it was used to assess if 12 samples, deriving from the twelve hospital centers had equal variances. Since the test was not significant for all the examined variables, the condition of “homogeneity of variances” in the data coming from the different centers was established.The analysis was performed in order to assess the existence of significant differences between patients affected by CD and UC, in the context of the IBD. Specifically, we examined data concerning 631 CD patients (36.6%) and 1091 UC patients (63.4%). Disease patterns, medical and surgical therapies, and risk factors for disease outcomes were analyzed. In particular, for each patient (in the respect of privacy) we acquired information about twentytwo variables: diagnosis age, gender, smoking habit (yes or no), use of immunosuppressive therapy (yes or no) and its duration, treatment with biological drugs (yes or no) and its duration, rehospitalization (yes or no), adverse events (yes or no), infections (yes or no), cancers (yes or no), diabetes (yes or no), hypertension (yes or no), heart failure (yes or no), kidney failure (yes or no), pulmonary failure (yes or no), neuropathy (yes or no), liver disease (yes or no), Charlson Index (the most widely used index to predict the tenyear mortality for a patient who may have comorbidity conditions; its score are 1, 2, 3 or 6, depending on the risk of death), surgery (yes or no), final exitus (survivor or died) and followup time. The hypotheses system is the following:
where 1 and 2 are the two examined inflammatory bowel diseases.
We used “Likelihood Ratio test” for categorical variables and “differences for two means test” for numerical ones. The used statistical package was NPC test, version 2.0, Statistical Software for Multivariate Nonparametric Permutation Test, Copyright 2001, Methodologica s.r.l.
In Table 1 we report, for both groups of patients, mean ± standard deviations (for numerical variables) and percentages (for categorical variables). The last column of the Table 1 shows the partial pvalues obtained by the application of NPC Test for analyzing the differences between the two examined groups; the last row shows the combined pvalue, referred to all twentytwo variables.
Examining the results achieved by applying NPC tests for partial and combined tests, we have to notice the high significance of the combined test, that provides guarantee affirming that patients with Crohn’s Disease and Ulcerative Colitis significantly differ between them, in relation to the examined variables. Focusing our attention on raw pvalues of partial tests, we can see that some variables significantly discriminate the two different subpopulations; in particular the UC patients, in compared to the CD patients, have a higher diagnosis age, do not show a marked smoking status, the proportion of patients treated with immunosuppressants or with biological drugs is lower than the CD patients, even if the duration of such therapies is longer. CD patients have a higher rate of rehospitalization; probably this is related to the significant greater occurrence of adverse events (rather than UC). Diabetes is more present in the subpopulation of UC patients. Analyzing the Charlson score we can highlight that UC patients have a more severe clinical situation than CD patients. Finally, the CD patients are more frequently subject to surgery compared to UC.
Since we are in presence of a high number of variables, we applied the Sequentially Rejective Multiple Test to determine which of the partial tests are effectively significant into discrimination between CD and UC patients.
In Table 2 we report, for each variable, the raw pvalues, the iindex (number expressing the ascending sort of raw pvalues) and the adjusted pvalues. Examining the raw pvalue (obtained from the NPC test), we note that twelve variables are apparently significant. After application of the Bonferroni Holm procedure, which provided adjusted pvalues, we can note that only nine of these variables were statistically significant; in accordance to iindex, they are: Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Rehospitalization, Gender and Duration of Immunosoppressive Therapy. So, with our data we can conclude that they are the only variables that significantly discriminate the Crohn’s Disease and Ulcerative Colitis groups.
Discussion
In general, IBDs affect 2.2 million people in Europe [15] and in Italy the estimated incidence of ulcerative colitis is 5.2 cases per 100,000 inhabitants per year, with a prevalence of approximately 70,150 cases / 100,000, and for Crohn’s disease 2.3 cases per 100,000 inhabitants per year, with a prevalence of 20–40 cases / 100,000 [42, 51].
In particular, we know that Crohn’s disease is spread all over the world and reaches the highest prevalence in Western nations. The ratio of affected females and males is around 1.35: 1 and many studies show that smokers are twice as likely to develop Crohn’s disease compared to nonsmokers [11, 13]. Our study, in line with previous literature, shows how CD patients, when compared with UC patients, do not exhibit a marked “smoker status” in the sense that smoking is more a cause of Crohn’s disease than of Ulcerative Colitis. Avoiding smoking, in a way, helps reduce the likelihood of contracting the disease.
IBDs can lead to various complications within the intestine, including obstruction, fistula and abscess development, as well as increase the risk of cancer in the inflammation area. For example, individuals with Crohn’s disease involving the small intestine are at greater risk for intestinal cancer. There is no certainty care yet [52].
Unfortunately, the IBD cannot be promptly prevented [29], even if complications and evolution can be prevented. Our analysis has made it possible to show more accurately the variables that most cause this disease. For this reason it is recommended to focus on the latter for prevention purposes.
As recalled several times, the use of nonparametric tests makes it possible to narrow the range of significant variables to focus on those of the most critical for preventive purposes.Ultimately, our analysis has made it possible to outline the variables that most discriminate these diseases.
From the statistical point of view, in this paper, one of the purposes was to examine and critically discuss the theoretical and practical relevance of permutation tests, demonstrating their effectiveness and ease of use in medical research. In literature NPC permutation tests have been successfully applied in many biomedical and epidemiologic fields, including gastroenterology [18, 19].
For statistical properties, the permutation tests have interesting property; in particular they are exact for any, even very small, sample size. This means that their null distributions, which are used to compute the pvalue, are known for each data set and for each sample size and this implies to controlling I and II error types. On the contrary, nonparametric tests are asymptotically guaranteed only for large sample sizes.
Besides, considering simultaneously different hypotheses, the problem of multiplicity or multiple testing problem arises. An incorrect approach is to test each hypothesis separately, using some level of significance α; in this case the real α level is bigger than nominal fixed level. Besides, the multiple testing approach consist to test simultaneous the set of hypoteses null and to use some appropriate correction to reached the desidered α level.
Specifically, the HolmBonferroni method is an approach that controls the probability that one or more type I errors will be adjusted, using adequate criteria for rejecting each of the individual hypotheses or comparisons. The comparison between groups is complex for the presence of multiple variables. This problem with parametric methods cannot be solved because the assumptions are too stringent.
NPC tests outweigh some of the limitations that traditional multivariate hypothesis verification procedures have, such as the ability to include a large number of variables. At the same time NPC tests offer a large number of advantages:

a)
this is an exact inferential procedure for any finite size of the sample;

b)
the solution is robust compared to the actual random distribution below the data (or error);

c)
the NPC procedure implicitly takes into account the underlying dependency structure of the response variables;

d)
it is not affected by the problem of loss of degrees of freedom when the number of variables increases.
Indeed, in contrast to traditional methods, increasing the number of information outputs also the power of the NPC test increases, i.e. the probability of detecting a true effect also increases monotonously [37]. In this sense, the NPC methodology can provide an effective and robust tool for statistical analysis of both experimental and observational medical studies.
In particular, in this paper we tried to show as the permutation tests are helpful for largesized data analysis in many applications contexts. In large data sets consisting of 1000 or more observations, performance of the permutation test appears equivalent to that of the asymptotic test; on the other hand, the NPC test, based on permutation solution, can be appropriately applied when the assumption for asymptotic tests are fulfilled [30]. In addition, unlike the classical nonparametric tests, the NPC method entails testing a global null hypothesis consisting of the intersection of K > 1 partial subhypotheses. In essence, the global null states that all of its constituent subhypotheses are true. \par The global alternative hypothesis is the union of K subalternatives. In this way NPC provides in multivariate context the combined pvalue, by means of an adequate combining function.
From the application point of view, we have great interest in evaluating this combined pvalue because it provides a result that takes into account the contribution of all examined variables; on the other hand, no other nonparametric test provides the advantage of a combined pvalue [41].
This particular feature justifies our choice of the NPC test as methodologically appropriate solution. In particular we applied permutation tests to perform comparison between a large number of patients affected by Crohn’s Disease and Ulcerative Colitis. Both of these illness are inflammatory bowel diseases, involving more than 100,000 people in Italy; they often arise in young people, go on for a lifetime and manifest alterations of the intestinal canal, causing relationship and working problems.
The results achieved applying NPC tests underline the high significance of the combined test, that shows that patients with Crohn’s Disease significantly differ from Ulcerative Colitis patients. Looking at the partial tests, we can notice that the differences between groups are referable to most of the examined variables; in particular the UC patients have a higher diagnosis age than CD patients, not showing a marked smoking status, the proportion of patients treated with immunosuppressants or with biological drugs is lower than the CD patients, even if the duration of these therapies is longer. On the other hand, CD patients have a higher rate of rehospitalization; probably it is related to the significant greater occurrence of adverse events (rather than UC). Diabetes is more recurrent in UC group. Moreover, UC patients have a more severe clinical profile, such as defined by Charlson score. Finally, the CD patients are more frequently subjected to surgery.
The findings of the study have a limit, which is represented by the sampling plan. Since the patients followed in the different hospital centers were examined and enrolled in the analyzed sample, we must admit that a sampling of convenience was chosen; it provides for the selection of the sample on the basis of criteria of convenience or practicality; it does not offer to all units of the population the same possibility of becoming part of the sample.
Conclusions
From a methodological point of view, thanks to BonferroniHolm procedure we were able to identify the really significant variables that discriminate the groups in exam, satisfying the multiplicity problem [40]. On the bases of the results we can affirm that Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Rehospitalization, Gender and Duration of Immunosoppressive Therapy are the variables effectively significant.
It is notable that the BonferroniHolm procedure leaves unchanged the original data information and allows a better interpretation of the results.
Until a few years ago the use of largesized data did not receive particular attention from researchers. Today the conspicuous availability of large amounts of data and the need of their analysis requires an adjustment of data processing methodologies, with careful attention to all the sources of variation in data. In this context, the nonparametric procedures, such as permutation tests, are widely applicable because of the numerous optimal properties of which they are characterized.
In the end, we can argue that the causes of IBD are not yet clear. In this paper we have identified the really significant variables that discriminate the groups under exam, satisfying the multiplicity problem. In fact we can affirm that Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Rehospitalization, Gender and Duration of Immunosuppressive Therapy are the effectively significant variables which can explain the occurrence of these diseases.
This work does not intend to provide a contribution in the clinical field of the IBD literature, but wants to allow a reflection on the possibility of using the NPC methodology to compare two chronic diseases (CD and UC) that affect the intestine (just the IBD) but which differ in some specific aspects.
It seems that the incidence of CD affects males and females with the same frequency (even if several studies allow to affirm that the female sex, especially if under the age of 45, presents a 20–30% greater risk of Crohn’s disease compared to males).
Among the environmental factors, the most important is the smoke that, curiously, seems to predispose to the CD rather than to the UC. IBDs are diseases that require medical therapy, close clinical surveillance and an appropriate therapeutic regimen. Medical therapy is based on the use of drugs such as immunosuppressant and biological drugs, but patients with CD are more frequently being subjected to such forms of therapy.
Medical therapy aims to induce clinical remission of the disease and keep patients free from relapses of the disease. In fact, the statistical comparison reveals that patients with CD (rather than UC) more frequently report cases of hospitalization due to IBD; probably it is related to the significant greater occurrence of adverse events.
Diabetes is more recurrent in UC group; this is already known because type 1 diabetes is the third most common comorbidity in patients with UC (after psoriasis and rheumatoid arthritis). However, diabetes can also complicate postoperative recovery in patients suffering from ulcerative colitis.
Moreover, UC patients have a more severe clinical profile, such as defined by Charlson score, the comorbidity of ulcerative colitis with other disorders of an EXTRAintestinal nature is very frequent.
Finally, the CD patients are more frequently subjected to surgery: surgery is an almost obligatory stage in the natural history of Crohn’s disease (about 70% of cases). The surgical intervention of intestinal resection, however, is almost invariably followed by recurrence of lesions (endoscopic relapse) and symptoms (clinical relapse).
Abbreviations
 ANOVA:

Analysis of variance
 CD:

Crohn’s disease
 IBD:

Inflammatory bowel diseases
 MAR:

Missing at random
 MCAR:

Missing completely at random
 NMAR:

Not missing at random
 NPC:

NonParametric combination
 UC:

Ulcerative colitis
References
 1.
Aickin M, Gensler H. Adjusting for multiple testing when reporting research results: the Bonferroni vs. Holm methods. Am J Public Health 1996;86(5): 726–728.
 2.
Alibrandi A, Giacalone M, Zirilli A, Moleti M. NPC to assess effects of maternal iodine nutrition and thyroid status on children cognitive development. In Proceedings of Compstat 2016, 22nd International Conference on Computatìonal Statistics 2016. ISBN/EAN: 978–90–73592360.
 3.
Antolini L, Bolzan M, Salmaso L. Metodi non parametrici per la verifica di ipotesi in indagini multicentriche. Statistica. 2007;62(3):523–33.
 4.
Arboretti Giancristofaro R, Marozzi M, Salmaso L. Repeated measures designs: a permutation approach for testing for active effects. Far East J Theoret Stat. 2005;16(2):303–25.
 5.
Arboretti Giancristofaro R, Brombin C. Overview of NonParametric combinationbased permutation tests for multivariate multisample problems. Statistica. 2014;74(3):233–46.
 6.
Basso D, Chiarandini M, Salmaso L. Synchronized permutation tests in I×J designs. J Stat Plan Inference. 2007;137(8):2564–78.
 7.
Bonnini S, Pesarin F, Salmaso L. Statistical Analysis in biomedical studies: an application of NPC Test to a clinical trial on a respiratory drug. In Congresso Nazionale della Società Italiana di Biometria. Società Italiana di Biometria; 2003. p. 10710.
 8.
Bonnini S, Corain L, Munaò F, Salmaso L. Neurocognitive effects in welders exposed to Aluminium: an application of the NPC test and NPC ranking methods. JISS. 2006;15(2):191–208.
 9.
Brombin C, Salmaso L. Multiaspect permutation tests in shape analysis with small sample size. Comput Stat Data Anal. 2009;53(12):3921–31. https://doi.org/10.1016/j.csda.2009.05.010.
 10.
Callegaro A, Pesarin R, Salmaso L. Test di permutazione per il confronto di curve di sopravvivenza. Statistica Applicata. 2003;15(2):241–61.
 11.
Cosnes J. Tobacco and IBD: relevance in the understanding of disease mechanisms and clinical practice. Best Pract Res Clin Gastroenterol. 2004;18(3):481–96. https://doi.org/10.1016/j.bpg.2003.12.003. PMID 15157822
 12.
Corain L, Salmaso L. Multivariate and multistrata nonparametric tests: the NPC method. J Modern Appl Stat Methods. 2004;3(2):443–61.
 13.
Corrao G, Tragnone A, Caprilli R, Trallori G, Papi C, Andreoli A, Di Paolo M, Riegler G, 2Rigo GP, Ferraù O, Mansi C, Ingrosso M, Valpiani D. Risk of inflammatory bowel disease attributable to smoking, oral contraception and breastfeeding in Italy: a nationwide casecontrol study. Int J Epidemiol. 1998;27(3):397–404.
 14.
Di Castelnuovo A, Mazzaro D, Pesarin R, Salmaso L. Test di permutazione multidimensionali in problemi d'inferenza isotonica: un'applicazione alla genetica. Statistica. 2000;60(4):691–700.
 15.
Edward V, Loftus JR. Clinical epidemiology of inflammatory bowel disease: incidence, prevalence and environmental influences. Gastroenterology. 2004;126:1504–17.
 16.
Finos L., Pesarin R, Salmaso L., Solari A. Nonparametric iterated combined tests for genetic differentiation. In Atti XLIH Riunione Scientica SIS 2004; CLEUP, Padova.
 17.
Finos L, Salmaso L. Weighted methods controlling the multiplicity when the number of variables is much higher than the number of observations. J Nonparametr Stat. 2006;18(2):245–61.
 18.
Floreani A, Caroli D, Variola A. A 35year followup of a large cohort of patients with primary biliary cirrhosis seen at a single Centre. Liver Int. 2011;31:361–8.
 19.
Floreani A, Cazzagon N, Franceschet I, Canesso F, Salmaso L, Baldo V. Metabolic syndrome associated with primary biliary cirrhosis. J Clin Gastroenterol. 2015;49:57–60.
 20.
Friedrich S, Brunner E, Pauly M. Permuting longitudinal data in spite of the dependencies. J Multivar Anal. 2017;153:255–65.
 21.
Fries W, Viola A, Manetti N, Frankovic I, Pugliese D, Monterubbianesi R, Samperi L. Disease patterns in lateonset ulcerative colitis: results from the IGIBD “AGED study”. Dig Liver Dis. 2017a;49(1):17–23.
 22.
Fries W, Viola A, Manetti N, Frankovic I, Pugliese D, Monterubbianesi R, Scalisi G, Aratari A, Cantoro L, Cappello M, Samperi L, Saibeni S, Casella G, Mocci G, Rea M, Furfaro F, Contaldo A, Magarotto A, Calella F, Manguso F, Inserra G, Privitera AC, Principi M, Castiglione F, Caprioli F, Ardizzone S, Danese S, Papi C, Bossa F, Kohn A, Armuzzi A, D’Incà R, Annese V, Alibrandi A, Bonovas S, Fiorino G, Italian Group for the study of Inflammatory Bowel Disease (IGIBD). Disease patterns in lateonset ulcerative colitis: results from the IGIBD "AGED study". Digestive Liver Diseas. 2017b;49(1):17–23.
 23.
Galimberti S, Valsecchi MG. Multivariate permutation test to compare survival curves for matched data. BMC Med Res Methodol. 2013;13(1):16.
 24.
Giacalone M., Zirilli A., Alibrandi A.. The use of permutation tests on largesized datasets. In Proceedings of the 48th Scientific Meeting of the Italian Statistical Society; Università degli Studi di Salerno, Monica Pratesi and Cira Perna Editors;2016. ISBN: 9788861970618.
 25.
Giacalone M, Zirilli A, Moleti M, Alibrandi A. Does the iodized salt therapy of pregnant mothers increase the children IQ? Empirical evidence of a statistical study based on permutation tests. Qual Quant. 2018;52:1423–35.
 26.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6(2):65–70.
 27.
Janssen A. Studentized permutation tests for noniid hypotheses and the generalized Behrensfisher problem. Stat Probabil Lett. 1997;36(1):9–21.
 28.
Klingenberg B, Solari A, Salmaso L, Pesarin F. Testing marginal homogeneity against stochastic order in multivariate ordinai data. Biometrics. 2008;65(2):452–62. https://doi.org/10.1111/j.15410420.2008.01067.X.
 29.
Kobashi G, Hata A, Uchida K, Ishige T, Abukawa D, Tajiri H, Uchiyama K, Hirota Y, Nagai M, T. J. P. I. B. D. Research Group. A casecontrol study to detect genetic and acquired risk factors for pediatric inflammatory bowel disease. Int J Epidemiol. 2015;44(1):232.
 30.
Ludbrook J, Dudley H. Why permutation tests are superior to t and F tests in biomedicai research. Am Stat. 1998;52(2):127–32.
 31.
Opdyke JD. Bootstraps, permutation tests and sampling orders of magnitude faster using SAS, Computational StatisticsWIREs 2013;5(5):390–405.
 32.
Pajouheshnia R, Pestman WR, Teerenstra S, Groenwold RHH. A computational approach to compare regression modelling strategies in prediction research. BMC Med Res Methodol. 2016;16:1.
 33.
Peek N, Holmes JH, Sun J. Technical challenges for big data in biomedicine and health: data sources, infrastructure, and analytics. Yearbook of Medicai Informatica. 2014;9(1):42–7. https://doi.org/10.15265/IY20140018.
 34.
Pesarin F. Multivariate Permutation Test. Chichester: Wiley and Sons; 2001.
 35.
Pesarin F, Salmaso L. Permutation Tests For Univariate And Multivariate Ordered Categorical Data. Aust J Statistica. 2006;35(2):315–4.
 36.
Pesarin F, Salmaso L. Permutation Tests for Complex Data. Theory, Applications and Software (a). Chichester: Wiley and Sons; 2010.
 37.
Pesarin F, Salmaso L. Finitesample consistency of combinationbased permutation tests with application to repeated measures designs. J Nonparametr Stat. 2010b;22:669–84.
 38.
Pesarin F, Salmaso L. The permutation testing approach: a review. Statistica. 2010;70(4):481–509.
 39.
Pesarin F, Salmaso L. Stat Comput. 2012;22:639. https://doi.org/10.1007/s1122201192610.
 40.
Pesarin F. Permutation tests: multivariate. Wiley StatsRef: Statistics Reference Online; 2014. p. 1–15.
 41.
Racioppi M, Salmaso L, Brombin C, Arboretti R, D'Agostino D, Colombo R, Serretta V, Brausi M, Casetta G, Gontero P, Hurle R, Tenaglia R, Altieri V, Bartoletti R, Maffezzini M, Siracusano S, Morgia G, Bassi PF. The clinical use of statistical permutation test methodology: a tool for identifying predective variables of outcome. Urol Int. 2015;94(3):262–9.
 42.
Ranzi T, Bodini P, Zambelli A, Politi P, Lupinacci G, Campanini MC, Dal Lago AL, Lisciandrano D, Bianchi PA. Epidemiological aspects of inflammatory bowel disease in a north Italian population: a 4 year prospective study. Eur J Gastroenterol Hepatol. 1996;8:657–61.
 43.
Reiss PT, Lei H, Maarten M. Fast functiononscalar regression with penalized basis expansions. Int J Biostat. 2010;6(1)
 44.
Rezzani A. Big Data. Architettura, tecnologie e metodi per l'utilizzo di grandi basi di dati; Apogeo education, Maggioli Editore, Milano; 2013.
 45.
Roy SN. On heuristic method of test construction and its use in multivariate analysis. Ann Math Stat. 1953;24:220–8.
 46.
Royston P, Altman DG. External validation of a cox prognostic model: principles and methods. BMC Med Res Methodol. 2013;13(1):33.
 47.
Rubin DB. Evaluations of the optimal discovery procedure for multiple testing. Int J Biostat. 2016;12(1):21–9.
 48.
Salmaso L. Permutation tests in screening twolevel factorial experiments. Adv App Stat. 2005;5(1):91–110.
 49.
Seibold H, Zeileis A, Hothron T. Modelbased recursive partitioning for subgroup analyses. Int J Biostat. 2016;12(1):45–3.
 50.
Sturino J, Zorych I, Mallick B, Pokusaeva K, Chang YY, Carroll RJ, Bliznuyk N. Statistical methods for comparative phenomics using highthroughput phenotype microarrays. Int J Biostat. 2010;6(1).
 51.
Tragnone A, Corrao G, Miglio F, Caprilli R, Lanfranchi GA. Incidence of inflammatory bowel disease in Italy: a nationwide populationbased study. Int J Epidemiol. 1996;25:1044–52.
 52.
Ueno F, Nakayama Y, Hagiwara E, Kurimoto S, Hibi T. Impact of inflammatory bowel disease on Japanese patients’ quality of life: results of a patient questionnaire survey. J Gastroenterol. 2017;52(5):555–67.
 53.
Zirilli A, Alibrandi A. A permutation solution to compare two hepatocellular carcinoma markers. JP J Biostat. 2011;5:2,109–19.
 54.
Zirilli A, Alibrandi A. The alteration of t,tmuconic acid and sphenilmercapturic acid levels due to benzene exposure: an application of NPC test. JP J of Biostat. 2012;7(2):91–104.
Acknowledgments
We have the pleasure of thanking Professor Walter Fries (Director of Gastroenterology Unit of the “G. Martino” University Hospital in Messina) for providing data and for his medical support into the realization of this paper.
Funding
No funding was available for this research.
Availability of data and materials
Professor Walter Fries (Director of Gastroenterology Unit of the “G. Martino” University Hospital in Messina) provided the data for this research.
Author information
Affiliations
Contributions
All the authors have developed the idea and designed the methodological aspects of the paper, showing a novel and innovative empirical study. They have given the same effort in the data analysis and statistical methods. GM gave his contribution in the data analysis, developing the comparison of the nonparametric techniques, ensuring the adequacy of the methodologies to data; ZA made contributions to the conception and paper design, ensuring that questions related to the accuracy and integrity of any part of the work were appropriately investigated and resolved; CPC has been involved in drafting of the manuscript, revising it critically for important intellectual contents and statistical methodology evaluations; AA gave her contribution in the acquisition of data working with medical teams in the results interpretation, providing also a general supervision of the research group. All the authors approved the final manuscript.
Corresponding author
Correspondence to Massimiliano Giacalone.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
Even if the Gastroenterology Unit of the “G. Martino” University Hospital in Messina provided the data, it had no input into the data selection, analysis design and interpretation of results.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Permutation tests
 BonferroniHolm procedure
 Multiplicity control
 Inflammatory bowel diseases
 Comparative analysis