Self-play: statistical significance

Heinz recently completed a comprehensive experiment in self-play using the FRITZ chess engine to establish the ‘decreasing returns’ hypothesis with specific levels of statistical confidence. This note revisits the results and recalculates the confidence levels of this and other hypotheses. These appear to be better than Heinz’ initial analysis suggests.

The infrared vibration-rotation spectrum of trans and cis nitrous acid

Infrared spectra of the trans and the cis isomers of nitrous acid, both HONO and DONO, have been observed in the gas phase using a Fourier transform interferometer with a resolution of about 0.05 cm−1 from 4000 to 500 cm−1. Rotational analyses are reported on eleven of the fundamentals and some overtones.

Statistical methods for examining genetic influences of resistance to anti-epileptic drugs

It is generally accepted that genetics may be an important factor in explaining the variation between patients’ responses to certain drugs. However, identification and confirmation of the responsible genetic variants is proving to be a challenge in many cases. A number of difficulties that maybe encountered in pursuit of these variants, such as non-replication of a true effect, population structure and selection bias, can be mitigated or at least reduced by appropriate statistical methodology. Another major statistical challenge facing pharmacogenetics studies is trying to detect possibly small polygenic effects using large volumes of genetic data, while controlling the number of false positive signals. Here we review statistical design and analysis options available for investigations of genetic resistance to anti-epileptic drugs.

Statistical downscaling of river flows

An extensive statistical ‘downscaling’ study is done to relate large-scale climate information from a general circulation model (GCM) to local-scale river flows in SW France for 51 gauging stations ranging from nival (snow-dominated) to pluvial (rainfall-dominated) river-systems. This study helps to select the appropriate statistical method at a given spatial and temporal scale to downscale hydrology for future climate change impact assessment of hydrological resources. The four proposed statistical downscaling models use large-scale predictors (derived from climate model outputs or reanalysis data) that characterize precipitation and evaporation processes in the hydrological cycle to estimate summary flow statistics. The four statistical models used are generalized linear (GLM) and additive (GAM) models, aggregated boosted trees (ABT) and multi-layer perceptron neural networks (ANN). These four models were each applied at two different spatial scales, namely at that of a single flow-gauging station (local downscaling) and that of a group of flow-gauging stations having the same hydrological behaviour (regional downscaling). For each statistical model and each spatial resolution, three temporal resolutions were considered, namely the daily mean flows, the summary statistics of fortnightly flows and a daily ‘integrated approach’. The results show that flow sensitivity to atmospheric factors is significantly different between nival and pluvial hydrological systems which are mainly influenced, respectively, by shortwave solar radiations and atmospheric temperature. The non-linear models (i.e. GAM, ABT and ANN) performed better than the linear GLM when simulating fortnightly flow percentiles. The aggregated boosted trees method showed higher and less variable R2 values to downscale the hydrological variability in both nival and pluvial regimes. Based on GCM cnrm-cm3 and scenarios A2 and A1B, future relative changes of fortnightly median flows were projected based on the regional downscaling approach. The results suggest a global decrease of flow in both pluvial and nival regimes, especially in spring, summer and autumn, whatever the considered scenario. The discussion considers the performance of each statistical method for downscaling flow at different spatial and temporal scales as well as the relationship between atmospheric processes and flow variability.

Incorporation of cis-9,trans-11 or trans-10,cis-12 conjugated linoleic acid into plasma and cellular lipids in healthy men

This study investigated the incorporation of cis-9,trans-11 conjugated linoleic acid (c9,t11 CLA) and trans-10,cis-12-CLA (t10,c12 CLA) into plasma and peripheral blood mononuclear cell (PBMC) lipids when consumed as supplements highly enriched in these isomers. Healthy men (n = 49, age 31 +/- 8 years) consumed one, two, and four capsules containing similar to600 mg of either c9,t11 CIA or t10,c12 CLA per capsule for sequential 8 week periods followed by a 6 week washout before consuming the alternative isomer. Both isomers were incorporated in a dosedependent manner into plasma phosphatidylcholine (PC) (c9,t11 CLA r = 0.779, t10,c12 CLA r = 0.738; P < 0.0001) and cholesteryl ester (CE) (c9,t11 CLA r = 0.706, t10,c12 CLA r = 0.788; P < 0.0001). Only t10,c12 CLA was enriched in plasma nonesterified fatty acids. Both c9,t11 CIA and t10,c12 CLA were incorporated linearly into PBMC total lipids (r = 0.285 and r = 0.273, respectively; P < 0.0005). The highest concentrations of c9,t11 CLA and t10,c12 CLA in PBMC lipids were 3- to 4-fold lower than those in plasma PC and CE. These data suggest that the level of intake is a major determinant of plasma and PBMC CLA content, although PBMCs appear to incorporate both CLA isomers less readily.

Incorporation of cis-9, trans-11 conjugated linoleic acid and vaccenic acid (trans-11 18 : 1) into plasma and leucocyte lipids in healthy men consuming dairy products naturally enriched in these fatty acids

The present study investigated whether consuming dairy products naturally enriched in cis-9, trans-11 (c9,t11) conjugated linoleic acid (CLA) by modification of cattle feed increases the concentration of this isomer in plasma and cellular lipids in healthy men. The study had a double-blind cross-over design. Subjects aged 34-60 years consumed dairy products available from food retailers for 1 week and then either control (0.17 g c9,t11 CLA/d; 0.31 g trans-vaccenic acid (tVA)/d) or CLA-enriched (1.43 g c9,t11 CLA/d; 4.71 g tVA/d) dairy products for 6 weeks. After 7 weeks washout, this was repeated with the alternate products. c9,t11 CLA concentration in plasma lipids was lower after consuming the control products, which may reflect the two-fold greater c9,t11 CLA content of the commercial products. Consuming the CLA-enriched dairy products increased the c9,t11 CLA concentration in plasma phosphatidylcholine (PC) (38 %; P=0.035), triacylglycerol (TAG) (22 %; P < 0.0001) and cholesteryl esters (205 %; P < 0.0001), and in peripheral blood mononuclear cells (PBMC) (238 %; P < 0.0001), while tVA concentration was greater in plasma PC (65 %; P=0.035), TAG (98 %; P=0.001) and PBMC (84 %; P=0.004). Overall, the present study shows that consumption of naturally enriched dairy products in amounts similar to habitual intakes of these foods increased the c9,t11 CLA content of plasma and cellular lipids.

Impact on CVD risk of modifying milk fat to decrease intake of SFA and increase intake of cis-MUFA

Despite the acknowledged benefits of reducing SFA intake few countries within the EU meet recognised targets. Milk and dairy products represent the single largest source of dietary SFA in most countries, yet epidemiological evidence indicates that milk has cardioprotective properties such that simply reducing consumption of dairy foods to meet SFA targets may not be a sound public health approach. The present paper explores the options for replacing some of the SFA in milk fat with cis-MUFA through alteration of the diet of the dairy cow, and the evidence that such changes can improve the indicators for CHD and CVD in general for the consumer. In addition, the outcome of such changes on risk factors for CHD and CVD at the population level is examined in the light of a modelling exercise involving data for eleven EU member states. Given the current and projected costs of health care, the results indicate that urgent consideration should be given to such a strategy.

Impact on CVD risk of modifying milk fat to decrease intake of SFA and increase intake of cis-MUFA

Despite the acknowledged benefits of reducing SFA intake few countries within the EU meet recognised targets. Milk and dairy products represent the single largest source of dietary SFA in most countries, yet epidemiological evidence indicates that milk has cardioprotective properties such that simply reducing consumption of dairy foods to meet SFA targets may not be a sound public health approach. The present paper explores the options for replacing some of the SFA in milk fat with cis-MUFA through alteration of the diet of the dairy cow, and the evidence that such changes can improve the indicators for CHD and CVD in general for the consumer. In addition, the outcome of such changes on risk factors for CHD and CVD at the population level is examined in the light of a modelling exercise involving data for eleven EU member states. Given the current and projected costs of health care, the results indicate that urgent consideration should be given to such a strategy.

Statistical evaluation of an acute dermal toxicity test using the dermal fixed dose procedure

The conventional method for the assessment of acute dermal toxicity (OECD Test Guideline 402, 1987) uses death of animals as an endpoint to identify the median lethal dose (LD50). A new OECD Testing Guideline called the dermal fixed dose procedure (dermal FDP) is being prepared to provide an alternative to Test Guideline 402. In contrast to Test Guideline 402, the dermal FDP does not provide a point estimate of the LD50, but aims to identify that dose of the substance under investigation that causes clear signs of nonlethal toxicity. This is then used to assign classification according to the new Globally Harmonised System of Classification and Labelling scheme (GHS). The dermal FDP has been validated using statistical modelling rather than by in vivo testing. The statistical modelling approach enables calculation of the probability of each GHS classification and the expected numbers of deaths and animals used in the test for imaginary substances with a range of LD50 values and dose-response curve slopes. This paper describes the dermal FDP and reports the results from the statistical evaluation. It is shown that the procedure will be completed with considerably less death and suffering than guideline 402, and will classify substances either in the same or a more stringent GHS class than that assigned on the basis of the LD50 value.

Statistical evaluation of the fixed concentration procedure for acute inhalation toxicity assessment

The conventional method for the assessment of acute inhalation toxicity (OECD Test Guideline 403, 1981) uses death of animals as an endpoint to identify the median lethal concentration (LC50). A new OECD Testing Guideline called the Fixed Concentration Procedure (FCP) is being prepared to provide an alternative to Test Guideline 403. Unlike Test Guideline 403, the FCP does not provide a point estimate of the LC50, but aims to identify an airborne exposure level that causes clear signs of nonlethal toxicity. This is then used to assign classification according to the new Globally Harmonized System of Classification and Labelling scheme (GHS). The FCP has been validated using statistical simulation rather than byin vivo testing. The statistical simulation approach predicts the GHS classification outcome and the numbers of deaths and animals used in the test for imaginary substances with a range of LC50 values and dose response curve slopes. This paper describes the FCP and reports the results from the statistical simulation study assessing its properties. It is shown that the procedure will be completed with considerably less death and suffering than Test Guideline 403, and will classify substances either in the same or a more stringent GHS class than that assigned on the basis of the LC50 value.

A statistical evaluation of the fixed dose procedure

The fixed-dose procedure (FDP) was introduced as OECD Test Guideline 420 in 1992, as an alternative to the conventional median lethal dose (LD50) test for the assessment of acute oral toxicity (OECD Test Guideline 401). The FDP uses fewer animals and causes less suffering than the conventional test, while providing information on the acute toxicity to allow substances to be ranked according to the EU hazard classification system. Recently the FDP has been revised, with the aim of providing further reductions and refinements, and classification according to the criteria of the Globally Harmonized Hazard Classification and Labelling scheme (GHS). This paper describes the revised FDP and analyses its properties, as determined by a statistical modelling approach. The analysis shows that the revised FDP classifies substances for acute oral toxicity generally in the same, or a more stringent, hazard class as that based on the LD50 value, according to either the GHS or the EU classification scheme. The likelihood of achieving the same classification is greatest for substances with a steep dose-response curve and median toxic dose (TD50) close to the LD50. The revised FDP usually requires five or six animals with two or fewer dying as a result of treatment in most cases.

Statistical design and analysis of pharmacogenetic trials

Pharmacogenetic trials investigate the effect of genotype on treatment response. When there are two or more treatment groups and two or more genetic groups, investigation of gene-treatment interactions is of key interest. However, calculation of the power to detect such interactions is complicated because this depends not only on the treatment effect size within each genetic group, but also on the number of genetic groups, the size of each genetic group, and the type of genetic effect that is both present and tested for. The scale chosen to measure the magnitude of an interaction can also be problematic, especially for the binary case. Elston et al. proposed a test for detecting the presence of gene-treatment interactions for binary responses, and gave appropriate power calculations. This paper shows how the same approach can also be used for normally distributed responses. We also propose a method for analysing and performing sample size calculations based on a generalized linear model (GLM) approach. The power of the Elston et al. and GLM approaches are compared for the binary and normal case using several illustrative examples. While more sensitive to errors in model specification than the Elston et al. approach, the GLM approach is much more flexible and in many cases more powerful. Copyright © 2005 John Wiley & Sons, Ltd.

Statistical methods for ordered categorical data based on a constrained odds model

The proportional odds model provides a powerful tool for analysing ordered categorical data and setting sample size, although for many clinical trials its validity is questionable. The purpose of this paper is to present a new class of constrained odds models which includes the proportional odds model. The efficient score and Fisher's information are derived from the profile likelihood for the constrained odds model. These results are new even for the special case of proportional odds where the resulting statistics define the Mann-Whitney test. A strategy is described involving selecting one of these models in advance, requiring assumptions as strong as those underlying proportional odds, but allowing a choice of such models. The accuracy of the new procedure and its power are evaluated.