The systematic profiling of false identity documents: method validation and performance evaluation using seizures known to originate from common and different sources

False identity documents constitute a potential powerful source of forensic intelligence because they are essential elements of transnational crime and provide cover for organized crime. In previous work, a systematic profiling method using false documents' visual features has been built within a forensic intelligence model. In the current study, the comparison process and metrics lying at the heart of this profiling method are described and evaluated. This evaluation takes advantage of 347 false identity documents of four different types seized in two countries whose sources were known to be common or different (following police investigations and dismantling of counterfeit factories). Intra-source and inter-sources variations were evaluated through the computation of more than 7500 similarity scores. The profiling method could thus be validated and its performance assessed using two complementary approaches to measuring type I and type II error rates: a binary classification and the computation of likelihood ratios. Very low error rates were measured across the four document types, demonstrating the validity and robustness of the method to link documents to a common source or to differentiate them. These results pave the way for an operational implementation of a systematic profiling process integrated in a developed forensic intelligence model.

Development of Performance Properties of Ternary Mixtures and Concrete Pavement Mixture Design and Analysis (MDA): Effect of Paste Quality on Fresh and Hardened Properties of Ternary Mixtures, TPF-5(117), 2012

The purpose of this study was to investigate the effect of cement paste quality on the concrete performance, particularly fresh properties, by changing the water-to-cementitious materials ratio (w/cm), type and dosage of supplementary cementitious materials (SCM), and airvoid system in binary and ternary mixtures. In this experimental program, a total matrix of 54 mixtures with w/cm of 0.40 and 0.45; target air content of 2%, 4%, and 8%; a fixed cementitious content of 600 pounds per cubic yard (pcy), and the incorporation of three types of SCMs at different dosages was prepared. The fine aggregate-to- total aggregate ratio was fixed at 0.42. Workability, rheology, air-void system, setting time, strength, Wenner Probe surface resistivity, and shrinkage were determined. The effects of paste variables on workability are more marked at the higher w/cm. The compressive strength is strongly influenced by the paste quality, dominated by w/cm and air content. Surface resistivity is improved by inclusion of Class F fly ash and slag cement, especially at later ages. Ternary mixtures performed in accordance with their ingredients. The data collected will be used to develop models that will be part of an innovative mix proportioning procedure.

Development of Performance Properties of Ternary Mixtures: Field Demonstrations and Project Summary, TPF-5(117), 2012

Supplementary cementitious materials (SCM) have become common parts of modern concrete practice. The blending of two or three cementitious materials to optimize durability, strength, or economics provides owners, engineers, materials suppliers, and contractors with substantial advantages over mixtures containing only portland cement. However, these advances in concrete technology and engineering have not always been adequately captured in specifications for concrete. Users need specific guidance to assist them in defining the performance requirements for a concrete application and the selection of optimal proportions of the cementitious materials needed to produce the required durable concrete. The fact that blended cements are currently available in many regions increases options for mixtures and thus can complicate the selection process. Both Portland and blended cements have already been optimized by the manufacturer to provide specific properties (such as setting time, shrinkage, and strength gain). The addition of SCMs (as binary, ternary, or even more complex mixtures) can alter these properties, and therefore has the potential to impact the overall performance and applications of concrete. This report is the final of a series of publications describing a project aimed at addressing effective use of ternary systems. The work was conducted in several stages and individual reports have been published at the end of each stage.

Adherence to oral anticancer treatment: focus on quality of execution during continuous schema of treatment

Background and objective: Oral anti-cancer treatments have expanded rapidly over the last years. While taking oral tablets at home ensures a better quality of life, it also exposes patients to the risk of sub-optimal adherence. The objective of this study is to assess how well ambulatory cancer patients execute their prescribed dosing regimen while they are engaged with continuous anti-cancer treatments. Design: This is an on-going longitudinal study. Consecutive patients starting an oral treatment are proposed to enter the study by the oncologist. Then they are referred to the pharmacy, where their oral anticancer treatment is dispensed in a Medication Event Monitoring System (MEMSTM), which records date and time of each opening of the drug container. Electronically compiled dosing history data from the MEMS are summarized and used as feedback during semistructured interviews with the pharmacist, which are dedicated to prevention and management of side effects. Interviews are scheduled before each medical visit. Report of the interview is available to the oncologist via an on-line secured portal. Setting: Seamless care approach between a Multidisciplinary Oncology Center and the Pharmacy of an Ambulatory Care and Community Medicine Department. Main outcome measures: For each patient, the comparison between the electronically compiled dosing history and the prescribed regimen was summarized using a daily binary indicator indicating whether yes or no the patient has taken the medication as prescribed. Results: Study started in March 2008. Among 22 eligible patients, 19 were included (11 men, median age 63 years old) and 3 (14%) refused to participate. 15 patients were prescribed a QD regimen, 3 patients a BID and 1 patient switched from QD to BID during follow-up. Median follow up was 182 days (IQR 72-252). Early discontinuation happened in four patients: side effects (n = 1), psychiatric reasons (n = 1), cancer progression (n = 1) and death (n = 1). On average, the daily number of medications was taken as prescribed in 99% of the follow-up days. Conclusions: Execution of the prescribed dosing regimens was almost perfect during the first 6 months. Maintaining this high degree of regimen execution and persistence over time might however be challenging in this population and need therefore to be confirmed in larger and longer follow-up cohort studies.

Direct Effect of Birth Weight on Childhood Blood Pressure: a Causal Mediation Analysis

Background: Numerous studies have shown a negative association between birth weight (BW) and blood pressure (BP) later in life. To estimate the direct effect of BW on BP, it is conventional to condition on current weight (CW). However, such conditioning can induce collider stratification bias in the estimate of the direct effect. Objective: To bound the potential bias due to U, an unmeasured common cause of CW and BP, on the estimate of the (controlled) direct effect of BW on BP. Methods: Data from a school based study in Switzerland were used (N = 4,005; 2,010 B/1,995 G; mean age: 12.3 yr [range: 10.1-14.9]). Measured common causes of BW-BP (SES, smoking, body weight, and hypertension status of the mother) and CW-BP (breastfeeding and child's physical activity and diet) were identified with DAGs. Linear regression models were fitted to estimate the association between BW and BP. Sensitivity analyses were conducted to assess the potential effect of U on the association between BW and BP. U was assumed 1) to be a binary variable that affected BP by the same magnitude in low BWand in normal BW children and 2) to have a different prevalence in low BW children and in normal BW children for a given CW. Results: A small negative association was observed between BW and BP [beta: -0.3 mmHg/kg (95% CI: -0.9 to 0.3)]. The association was strengthened upon conditioning for CW [beta: -1.5 mmHg/kg (95% CI: -2.1 to -0.9)]. Upon further conditioning on common causes of BW-BP and CW-BP, the association did not change substantially [beta: -1.4 mmHg/kg (95% CI: -2.0 to -0.8)]. The negative association could be explained by U only if U was strongly associated with BP and if there was a large difference in the prevalence of U between low BWand normal BW children. Conclusion: The observed negative association between BW and BP upon adjustment for CW was not easily explained by an unmeasured common cause of CWand BP.

Polydispersity index from linear viscoelastic data: Unimodal and bimodal linear polymer melts

This article describes a method for determining the polydispersity index Ip2=Mz/Mw of the molecular weight distribution (MWD) of linear polymeric materials from linear viscoelastic data. The method uses the Mellin transform of the relaxation modulus of a simple molecular rheological model. One of the main features of this technique is that it enables interesting MWD information to be obtained directly from dynamic shear experiments. It is not necessary to achieve the relaxation spectrum, so the ill-posed problem is avoided. Furthermore, a determinate shape of the continuous MWD does not have to be assumed in order to obtain the polydispersity index. The technique has been developed to deal with entangled linear polymers, whatever the form of the MWD is. The rheological information required to obtain the polydispersity index is the storage G′(ω) and loss G″(ω) moduli, extending from the terminal zone to the plateau region. The method provides a good agreement between the proposed theoretical approach and the experimental polydispersity indices of several linear polymers for a wide range of average molecular weights and polydispersity indices. It is also applicable to binary blends.

Determinants of methicillin-susceptible Staphylococcus aureus native bone and joint infection treatment failure: a retrospective cohort study.

BACKGROUND: Although methicillin-susceptible Staphylococcus aureus (MSSA) native bone and joint infection (BJI) constitutes the more frequent clinical entity of BJI, prognostic studies mostly focused on methicillin-resistant S. aureus prosthetic joint infection. We aimed to assess the determinants of native MSSA BJI outcomes. METHODS: Retrospective cohort study (2001-2011) of patients admitted in a reference hospital centre for native MSSA BJI. Treatment failure determinants were assessed using Kaplan-Meier curves and binary logistic regression. RESULTS: Sixty-six patients (42 males [63.6%]; median age 61.2 years; interquartile range [IQR] 45.9-71.9) presented an acute (n = 38; 57.6%) or chronic (n = 28; 42.4%) native MSSA arthritis (n = 15; 22.7%), osteomyelitis (n = 19; 28.8%) or spondylodiscitis (n = 32; 48.5%), considered as "difficult-to-treat" in 61 cases (92.4%). All received a prolonged (27.1 weeks; IQR, 16.9-36.1) combined antimicrobial therapy, after surgical management in 37 cases (56.1%). Sixteen treatment failures (24.2%) were observed during a median follow-up period of 63.3 weeks (IQR, 44.7-103.1), including 13 persisting infections, 1 relapse after treatment disruption, and 2 super-infections. Independent determinants of treatment failure were the existence of a sinus tract (odds ratio [OR], 5.300; 95% confidence interval [CI], 1.166-24.103) and a prolonged delay to infectious disease specialist referral (OR, 1.134; 95% CI 1.013-1.271). CONCLUSIONS: The important treatment failure rate pinpointed the difficulty of cure encountered in complicated native MSSA BJI. An early infectious disease specialist referral is essential, especially in debilitated patients or in presence of sinus tract.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

Polymorphism of the grasshopper Rhammatocerus schistocercoides populations revealed by RAPD

The objective of this work was to study the genetic variability of the grasshopper Rhammatocerus schistocercoides (Orthoptera: Acrididae) using RAPD analysis among individuals from three populations, one from Colombia and two from Brazil (Goiás and Mato Grosso States). Ninety scorable binary markers were obtained by fingerprinting with 11 oligonucleotide primers. Most of the polymorphism was attributed to 42 markers with variable frequency among the different populations. Although the existence of significant difference among populations (P<0.0001), most of the genetic variability was found within populations (87.7% of total variation). Pairwise distances between Colombian and Brazilian populations were 0.12 (P<0.0001) and 0.18 (P<0.0001) for Goiás and Mato Grosso, respectively. The pairwise distance between Goiás and Mato Grosso populations was 0.06 (P<0.0001). These data indicated that the phenotypic differences among populations are associated mainly with the geographical distances between the Brazilian and Colombian populations.

SCRIB and PUF60 Are Primary Drivers of the Multisystemic Phenotypes of the 8q24.3 Copy-Number Variant.

Copy-number variants (CNVs) represent a significant interpretative challenge, given that each CNV typically affects the dosage of multiple genes. Here we report on five individuals with coloboma, microcephaly, developmental delay, short stature, and craniofacial, cardiac, and renal defects who harbor overlapping microdeletions on 8q24.3. Fine mapping localized a commonly deleted 78 kb region that contains three genes: SCRIB, NRBP2, and PUF60. In vivo dissection of the CNV showed discrete contributions of the planar cell polarity effector SCRIB and the splicing factor PUF60 to the syndromic phenotype, and the combinatorial suppression of both genes exacerbated some, but not all, phenotypic components. Consistent with these findings, we identified an individual with microcephaly, short stature, intellectual disability, and heart defects with a de novo c.505C>T variant leading to a p.His169Tyr change in PUF60. Functional testing of this allele in vivo and in vitro showed that the mutation perturbs the relative dosage of two PUF60 isoforms and, subsequently, the splicing efficiency of downstream PUF60 targets. These data inform the functions of two genes not associated previously with human genetic disease and demonstrate how CNVs can exhibit complex genetic architecture, with the phenotype being the amalgam of both discrete dosage dysfunction of single transcripts and also of binary genetic interactions.

Decision curve analysis revisited: overall net benefit, relationships to ROC curve analysis, and application to case-control studies.

ABSTRACT: BACKGROUND: Decision curve analysis has been introduced as a method to evaluate prediction models in terms of their clinical consequences if used for a binary classification of subjects into a group who should and into a group who should not be treated. The key concept for this type of evaluation is the "net benefit", a concept borrowed from utility theory. METHODS: We recall the foundations of decision curve analysis and discuss some new aspects. First, we stress the formal distinction between the net benefit for the treated and for the untreated and define the concept of the "overall net benefit". Next, we revisit the important distinction between the concept of accuracy, as typically assessed using the Youden index and a receiver operating characteristic (ROC) analysis, and the concept of utility of a prediction model, as assessed using decision curve analysis. Finally, we provide an explicit implementation of decision curve analysis to be applied in the context of case-control studies. RESULTS: We show that the overall net benefit, which combines the net benefit for the treated and the untreated, is a natural alternative to the benefit achieved by a model, being invariant with respect to the coding of the outcome, and conveying a more comprehensive picture of the situation. Further, within the framework of decision curve analysis, we illustrate the important difference between the accuracy and the utility of a model, demonstrating how poor an accurate model may be in terms of its net benefit. Eventually, we expose that the application of decision curve analysis to case-control studies, where an accurate estimate of the true prevalence of a disease cannot be obtained from the data, is achieved with a few modifications to the original calculation procedure. CONCLUSIONS: We present several interrelated extensions to decision curve analysis that will both facilitate its interpretation and broaden its potential area of application.

Transcranial Doppler Pulsatility Index: What it is and What it Isn't.

BACKGROUND: Transcranial Doppler (TCD) pulsatility index (PI) has traditionally been interpreted as a descriptor of distal cerebrovascular resistance (CVR). We sought to evaluate the relationship between PI and CVR in situations, where CVR increases (mild hypocapnia) and decreases (plateau waves of intracranial pressure-ICP). METHODS: Recordings from patients with head-injury undergoing monitoring of arterial blood pressure (ABP), ICP, cerebral perfusion pressure (CPP), and TCD assessed cerebral blood flow velocities (FV) were analyzed. The Gosling pulsatility index (PI) was compared between baseline and ICP plateau waves (n = 20 patients) or short term (30-60 min) hypocapnia (n = 31). In addition, a modeling study was conducted with the "spectral" PI (calculated using fundamental harmonic of FV) resulting in a theoretical formula expressing the dependence of PI on balance of cerebrovascular impedances. RESULTS: PI increased significantly (p < 0.001) while CVR decreased (p < 0.001) during plateau waves. During hypocapnia PI and CVR increased (p < 0.001). The modeling formula explained more than 65% of the variability of Gosling PI and 90% of the variability of the "spectral" PI (R = 0.81 and R = 0.95, respectively). CONCLUSION: TCD pulsatility index can be easily and quickly assessed but is usually misinterpreted as a descriptor of CVR. The mathematical model presents a complex relationship between PI and multiple haemodynamic variables.

Computational problems in perfect phylogeny haplotyping: Typing without calling the allele

A haplotype is an m-long binary vector. The XOR-genotype of two haplotypes is the m-vector of their coordinate-wise XOR. We study the following problem: Given a set of XOR-genotypes, reconstruct their haplotypes so that the set of resulting haplotypes can be mapped onto a perfect phylogeny (PP) tree. The question is motivated by studying population evolution in human genetics and is a variant of the PP haplotyping problem that has received intensive attention recently. Unlike the latter problem, in which the input is '' full '' genotypes, here, we assume less informative input and so may be more economical to obtain experimentally. Building on ideas of Gusfield, we show how to solve the problem in polynomial time by a reduction to the graph realization problem. The actual haplotypes are not uniquely determined by the tree they map onto and the tree itself may or may not be unique. We show that tree uniqueness implies uniquely determined haplotypes, up to inherent degrees of freedom, and give a sufficient condition for the uniqueness. To actually determine the haplotypes given the tree, additional information is necessary. We show that two or three full genotypes suffice to reconstruct all the haplotypes and present a linear algorithm for identifying those genotypes.

Multiscale assessment of plant diversity in wooded pastures of the Jura mountains

Résumé La diminution de la biodiversité, à toutes les échelles spatiales et sur l'ensemble de la planète, compte parmi les problèmes les plus préoccupants de notre époque. En terme de conservation, il est aujourd'hui primordial de mieux comprendre les mécanismes qui créent et maintiennent la biodiversité dans les écosystèmes naturels ou anthropiques. La présente étude a pour principal objectif d'améliorer notre compréhension des patrons de biodiversité végétale et des mécanismes sous jacents, dans un écosystème complexe, riche en espèces et à forte valeur patrimoniale, les pâturages boisés jurassiens. Structure et échelle spatiales sont progressivement reconnues comme des dimensions incontournables dans l'étude des patrons de biodiversité. De plus, ces deux éléments jouent un rôle central dans plusieurs théories écologiques. Toutefois, peu d'hypothèses issues de simulations ou d'études théoriques concernant le lien entre structure spatiale du paysage et biodiversité ont été testées de façon empirique. De même, l'influence des différentes composantes de l'échelle spatiale sur les patrons de biodiversité est méconnue. Cette étude vise donc à tester quelques-unes de ces hypothèses et à explorer les patrons spatiaux de biodiversité dans un contexte multi-échelle, pour différentes mesures de biodiversité (richesse et composition en espèces) à l'aide de données de terrain. Ces données ont été collectées selon un plan d'échantillonnage hiérarchique. Dans un premier temps, nous avons testé l'hypothèse élémentaire selon laquelle la richesse spécifique (le nombre d'espèces sur une surface donnée) est liée à l'hétérogénéité environnementale quelque soit l'échelle. Nous avons décomposé l'hétérogénéité environnementale en deux parties, la variabilité des conditions environnementales et sa configuration spatiale. Nous avons montré que, en général, la richesse spécifique augmentait avec l'hétérogénéité de l'environnement : elle augmentait avec le nombre de types d'habitats et diminuait avec l'agrégation spatiale de ces habitats. Ces effets ont été observés à toutes les échelles mais leur nature variait en fonction de l'échelle, suggérant une modification des mécanismes. Dans un deuxième temps, la structure spatiale de la composition en espèces a été décomposée en relation avec 20 variables environnementales et 11 traits d'espèces. Nous avons utilisé la technique de partition de la variation et un descripteur spatial, récemment développé, donnant accès à une large gamme d'échelles spatiales. Nos résultats ont montré que la structure spatiale de la composition en espèces végétales était principalement liée à la topographie, aux échelles les plus grossières, et à la disponibilité en lumière, aux échelles les plus fines. La fraction non-environnementale de la variation spatiale de la composition spécifique avait une relation complexe avec plusieurs traits d'espèces suggérant un lien avec des processus biologiques tels que la dispersion, dépendant de l'échelle spatiale. Dans un dernier temps, nous avons testé, à plusieurs échelles spatiales, les relations entre trois composantes de la biodiversité : la richesse spécifique totale d'un échantillon (diversité gamma), la richesse spécifique moyenne (diversité alpha), mesurée sur des sous-échantillons, et les différences de composition spécifique entre les sous-échantillons (diversité beta). Les relations deux à deux entre les diversités alpha, beta et gamma ne suivaient pas les relations attendues, tout du moins à certaines échelles spatiales. Plusieurs de ces relations étaient fortement dépendantes de l'échelle. Nos résultats ont mis en évidence l'importance du rapport d'échelle (rapport entre la taille de l'échantillon et du sous-échantillon) lors de l'étude des patrons spatiaux de biodiversité. Ainsi, cette étude offre un nouvel aperçu des patrons spatiaux de biodiversité végétale et des mécanismes potentiels permettant la coexistence des espèces. Nos résultats suggèrent que les patrons de biodiversité ne peuvent être expliqués par une seule théorie, mais plutôt par une combinaison de théories. Ils ont également mis en évidence le rôle essentiel joué par la structure spatiale dans la détermination de la biodiversité, quelque soit le composant de la biodiversité considéré. Enfin, cette étude souligne l'importance de prendre en compte plusieurs échelles spatiales et différents constituants de l'échelle spatiale pour toute étude relative à la diversité spécifique. Abstract The world-wide loss of biodiversity at all scales has become a matter of urgent concern, and improving our understanding of local drivers of biodiversity in natural and anthropogenic ecosystems is now crucial for conservation. The main objective of this study was to further our comprehension of the driving forces controlling biodiversity patterns in a complex and diverse ecosystem of high conservation value, wooded pastures. Spatial pattern and scale are central to several ecological theories, and it is increasingly recognized that they must be taken -into consideration when studying biodiversity patterns. However, few hypotheses developed from simulations or theoretical studies have been tested using field data, and the evolution of biodiversity patterns with different scale components remains largely unknown. We test several such hypotheses and explore spatial patterns of biodiversity in a multi-scale context and using different measures of biodiversity (species richness and composition), with field data. Data were collected using a hierarchical sampling design. We first tested the simple hypothesis that species richness, the number of species in a given area, is related to environmental heterogeneity at all scales. We decomposed environmental heterogeneity into two parts: the variability of environmental conditions and its spatial configuration. We showed that species richness generally increased with environmental heterogeneity: species richness increased with increasing number of habitat types and with decreasing spatial aggregation of those habitats. Effects occurred at all scales but the nature of the effect changed with scale, suggesting a change in underlying mechanisms. We then decomposed the spatial structure of species composition in relation to environmental variables and species traits using variation partitioning and a recently developed spatial descriptor, allowing us to capture a wide range of spatial scales. We showed that the spatial structure of plant species composition was related to topography at the coarsest scales and insolation at finer scales. The non-environmental fraction of the spatial variation in species composition had a complex relationship with several species traits, suggesting a scale-dependent link to biological processes, particularly dispersal. Finally, we tested, at different spatial scales, the relationships between different components of biodiversity: total sample species richness (gamma diversity), mean species .richness (alpha diversity), measured in nested subsamples, and differences in species composition between subsamples (beta diversity). The pairwise relationships between alpha, beta and gamma diversity did not follow the expected patterns, at least at certain scales. Our result indicated a strong scale-dependency of several relationships, and highlighted the importance of the scale ratio when studying biodiversity patterns. Thus, our results bring new insights on the spatial patterns of biodiversity and the possible mechanisms allowing species coexistence. They suggest that biodiversity patterns cannot be explained by any single theory proposed in the literature, but a combination of theories is sufficient. Spatial structure plays a crucial role for all components of biodiversity. Results emphasize the importance of considering multiple spatial scales and multiple scale components when studying species diversity.

FAM-MDR: a flexible family-based multifactor dimensionality reduction technique to detect epistasis using related individuals

We propose a novel multifactor dimensionality reduction method for epistasis detection in small or extended pedigrees, FAM-MDR. It combines features of the Genome-wide Rapid Association using Mixed Model And Regression approach (GRAMMAR) with Model-Based MDR (MB-MDR). We focus on continuous traits, although the method is general and can be used for outcomes of any type, including binary and censored traits. When comparing FAM-MDR with Pedigree-based Generalized MDR (PGMDR), which is a generalization of Multifactor Dimensionality Reduction (MDR) to continuous traits and related individuals, FAM-MDR was found to outperform PGMDR in terms of power, in most of the considered simulated scenarios. Additional simulations revealed that PGMDR does not appropriately deal with multiple testing and consequently gives rise to overly optimistic results. FAM-MDR adequately deals with multiple testing in epistasis screens and is in contrast rather conservative, by construction. Furthermore, simulations show that correcting for lower order (main) effects is of utmost importance when claiming epistasis. As Type 2 Diabetes Mellitus (T2DM) is a complex phenotype likely influenced by gene-gene interactions, we applied FAM-MDR to examine data on glucose area-under-the-curve (GAUC), an endophenotype of T2DM for which multiple independent genetic associations have been observed, in the Amish Family Diabetes Study (AFDS). This application reveals that FAM-MDR makes more efficient use of the available data than PGMDR and can deal with multi-generational pedigrees more easily. In conclusion, we have validated FAM-MDR and compared it to PGMDR, the current state-of-the-art MDR method for family data, using both simulations and a practical dataset. FAM-MDR is found to outperform PGMDR in that it handles the multiple testing issue more correctly, has increased power, and efficiently uses all available information.