Tests for Comparing Mark-Specific Hazards and Cumulative Incidence Functions

It is of interest in some applications to determine whether there is a relationship between a hazard rate function (or a cumulative incidence function) and a mark variable which is only observed at uncensored failure times. We develop nonparametric tests for this problem when the mark variable is continuous. Tests are developed for the null hypothesis that the mark-specific hazard rate is independent of the mark versus ordered and two-sided alternatives expressed in terms of mark-specific hazard functions and mark-specific cumulative incidence functions. The test statistics are based on functionals of a bivariate test process equal to a weighted average of differences between a Nelson--Aalen-type estimator of the mark-specific cumulative hazard function and a nonparametric estimator of this function under the null hypothesis. The weight function in the test process can be chosen so that the test statistics are asymptotically distribution-free.Asymptotically correct critical values are obtained through a simple simulation procedure. The testing procedures are shown to perform well in numerical studies, and are illustrated with an AIDS clinical trial example. Specifically, the tests are used to assess if the instantaneous or absolute risk of treatment failure depends on the amount of accumulation of drug resistance mutations in a subject's HIV virus. This assessment helps guide development of anti-HIV therapies that surmount the problem of drug resistance.

Validation of the LittlEARS((R)) Auditory Questionnaire in children with normal hearing

OBJECTIVES: With more children receiving cochlear implants during infancy, there is a need for validated assessments of pre-verbal and early verbal auditory skills. The LittlEARS Auditory Questionnaire is presented here as the first module of the LittlEARS test battery. The LittlEARS Auditory Questionnaire was developed and piloted to assess the auditory behaviour of normal hearing children and hearing impaired children who receive a cochlear implant or hearing aid prior to 24 months of age. This paper presents results from two studies: one validating the LittlEARS Auditory Questionnaire on children with normal hearing who are German speaking and a second validating the norm curves found after adaptation and administration of the questionnaire to children with normal hearing in 15 different languages. METHODS: Scores from a group of 218 German and Austrian children with normal hearing between 5 days and 24 months of age were used to create a norm curve. The questionnaire was adapted from the German original into English and then 15 other languages to date. Regression curves were found based on parental responses from 3309 normal hearing infants and toddlers. Curves for each language were compared to the original German validation curve. RESULTS: The results of the first study were a norm curve which reflects the age-dependence of auditory behaviour, reliability and homogeneity as a measure of auditory behaviour, and calculations of expected and critical values as a function of age. Results of the second study show that the regression curves found for all the adapted languages are essentially equal to the German norm curve, as no statistically significant differences were found. CONCLUSIONS: The LittlEARS Auditory Questionnaire is a valid, language-independent tool for assessing the early auditory behaviour of infants and toddlers with normal hearing. The results of this study suggest that the LittlEARS Auditory Questionnaire could also be very useful for documenting children's progress with their current amplification, providing evidence of the need for implantation, or highlighting the need for follow-up in other developmental areas.

Assessment of interrater agreement using a binomial-based statistic: K(bin)

This project develops K(bin), a relatively simple, binomial based statistic for assessing interrater agreement in which expected agreement is calculated a priori from the number of raters involved in the study and number of categories on the rating tool. The statistic is logical in interpretation, easily calculated, stable for small sample sizes, and has application over a wide range of possible combinations from the simplest case of two raters using a binomial scale to multiple raters using a multiple level scale.^ Tables of expected agreement values and tables of critical values for K(bin) which include power to detect three levels of the population parameter K for n from 2 to 30 and observed agreement $\ge$.70 calculated at alpha =.05,.025, and.01 are included.^ An example is also included which describes the use of the tables for planning and evaluating an interrater reliability study using the statistic, K(bin). ^

Burrowing in Marine Muds by Crack Propagation: Kinematics and Forces

The polychaete Nereis virens burrows through muddy sediments by exerting dorsoventral forces against the walls of its tongue-depressor- shaped burrow to extend an oblate hemispheroidal crack. Stress is concentrated at the crack tip, which extends when the stress intensity factor (K-I) exceeds the critical stress intensity factor (K-Ic). Relevant forces were measured in gelatin, an analog for elastic muds, by photoelastic stress analysis, and were 0.015 +/- 0.001 N (mean +/- s.d.;N= 5). Measured elastic moduli (E) for gelatin and sediment were used in finite element models to convert the forces in gelatin to those required in muds to maintain the same body shapes observed in gelatin. The force increases directly with increasing sediment stiffness, and is 0.16 N for measured sediment stiffness of E=2.7x10(4) Pa. This measurement of forces exerted by burrowers is the first that explicitly considers the mechanical behavior of the sediment. Calculated stress intensity factors fall within the range of critical values for gelatin and exceed those for sediment, showing that crack propagation is a mechanically feasible mechanism of burrowing. The pharynx extends anteriorly as it everts, extending the crack tip only as far as the anterior of the worm, consistent with wedge-driven fracture and drawing obvious parallels between soft-bodied burrowers and more rigid, wedge-shaped burrowers(i.e. clams). Our results raise questions about the reputed high energetic cost of burrowing and emphasize the need for better understanding of sediment mechanics to quantify external energy expenditure during burrowing.

Anion exchange membrane soil nitrate predicts turfgrass color and yield.

Desirable nitrogen (N) management practices for turfgrass supply sufficient N for high quality turf while limiting excess soil N. Previous studies suggested the potential of anion exchange membranes (AEMs) for predicting turfgrass color, quality, or yield. However, these studies suggested a wide range of critical soil nitrate-nitrogen (NO3-N) values across sample dates. A field experiment, in randomized complete block design with treatments consisting of nine N application rates, was conducted on a mixed species cool-season turfgrass lawn across two growing seasons. Every 2 wk from May to October, turfgrass color was assessed with three different reflectance meters, and soil NO3-N was measured with in situ AEMs. Cate-Nelson models were developed comparing relative reflectance value and yield to AEM desorbed soil NO3-N pooled across all sample dates. These models predicted critical AEM soil NO3-N values from 0. 45 to 1.4 micro g cm-2 d-1. Turf had a low probability of further positive response to AEM soil NO3-N greater than these critical values. These results suggest that soil NO3-N critical values from AEMs may be applicable across sample dates and years and may serve to guide N fertilization to limit excess soil NO3-N.

Rare earth and other element contents in surface layer bottom sediments and ferromanganese nodules along the Transatlantic profile

Behavior of rare earth elements (REE) and Th is studied along the Transatlantic transect at 22°N. It is shown that both REE and Th contents relative to Al (the most lithogenic element) increase toward the pelagic region. The increasing trend becomes more complicated due to variations in content of biogenic calcium carbonate that acts as a diluting component in sediments. REE composition varies symmetrically relative to the Mid-Atlantic Ridge (MAR) emphasizing weak hydrothermal influence on sediments of the ridge axis, although the well-known criteria for hydrothermal contribution, such as Al/(Al+Mn+Fe) and (Fe+Mn)/Ti, do not reach critical values. Variations in REE content and composition allowed to distinguish the following five sediment zones in the transect: (I) terrigenous sediments of the Nares abyssal plain; (II) pelagic sediments of the North American Basin; (III) carbonate ooze of the MAR axis; (IV) pelagic sediments of the Canary Basin; and (V) terrigenous clay and calcareous mud of the African continental slope and slope base. Ferromanganese nodules of the hydrogenous type with extremely high Ce (up to 1801 ppm) and Th (up to 138 ppm) contents occur in pelagic sediments. It is ascertained that P, REE, and Th contents depend on Fe content in Atlantic sediments. Therefore, one can suggest that only minor amount of phosphorus is bound with bone debris. Low concentration of bone debris phosphorus is a result of relatively high sedimentation rates in the Atlantic Ocean, as compared with those in pelagic regions of the Pacific Ocean.

Surface time series of nitrate, silicate and chlorophyll a (1999-2011) and CTD profiles (1999) at the PhytoCly station, Bay of Calvi, Corsica, Mediterranean Sea

This work is based on a long time series of data collected in the well-preserved Bay of Calvi (Corsica island, Ligurian Sea, NW Mediterranean) between 1979 and 2011, which include physical characteristics (31 years), chlorophyll a (chl a, 15 years), and inorganic nutrients (13 years). Because samples were collected at relatively high frequencies, which ranged from daily to biweekly during the winter-spring period, it was possible to (1) evidence the key role of two interacting physical variables, i.e. water temperature and wind intensity, on nutrient replenishment and phytoplankton dynamics during the winter-spring period, (2) determine critical values of physical factors that explained interannual variability in the replenishment of surface nutrients and the winter-spring phytoplankton bloom, and (3) identify previously unrecognized characteristics of the planktonic ecosystem. Over the >30 year observation period, the main driver of nutrient replenishment and phytoplankton (chl a) development was the number of wind events (mean daily wind speed >5 m s-1) during the cold-water period (subsurface water <13.5°C). According to winter intensity, there were strong differences in both the duration and intensity of nutrient fertilization and phytoplankton blooms (chl a). The trophic character of the Bay of Calvi changed according to years, and ranged from very oligotrophic (i.e. subtropical regime, characterized by low seasonal variability) to mesotrophic (i.e. temperate regime, with a well-marked increase in nutrient concentrations and chl a during the winter-spring period) during mild and moderate winters, respectively. A third regime occurred during severe winters characterized by specific wind conditions (i.e. high frequency of northeasterly winds), when Mediterranean "high nutrient - low chlorophyll" conditions occurred as a result of enhanced crossshore exchanges and associated offshore export of the nutrient-rich water. There was no long-term trend (e.g. climatic) in either nutrient replenishment or the winter-spring phytoplankton bloom between 1979 and 2011, but both nutrients and chl a reflected interannual and decadal changes in winter intensity.

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.

Edge-to-edge matching and its applications - Part I. Application to the simple HCP/BCC system

The edge-to-edge matching crystallographic model has been used to predict all the orientation relationships (OR) between crystals that have simple hexagonal close packed (HCP) and body-centered cubic (BCC) structures. Using the critical values for the interatomic spacing misfit along the matching directions and the cl-value mismatch between matching planes, the model predicted all the four common ORs, namely the Burgers OR, the Potter OR, the Pitsch-Schrader OR and the Rong Dunlop OR, together with the corresponding habit planes. Taking the c(H)/a(H) and a(H)/a(B) ratios as variables, where H and B denote the HCP and BCC structures respectively, the model also predicted the relationship between these variables and the four ORs. These predictions are perfectly consistent with the published experimental results. As was the case in the FCC/BCC system, the edge-to-edge matching model has been shown to be a powerful tool for predicting the crystallographic features of diffusion-controlled phase transformations. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Finite deformation model of simple shear of fault with microrotations: apparent strain localisation and en-echelon fracture pattern

Strain localisation is a widespread phenomenon often observed in shear and compressive loading of geomaterials, for example, the fault gouge. It is believed that the main mechanisms of strain localisation are strain softening and mismatch between dilatancy and pressure sensitivity. Observations show that gouge deformation is accompanied by considerable rotations of grains. In our previous work as a model for gouge material, we proposed a continuum description for an assembly of particles of equal radius in which the particle rotation is treated as an independent degree of freedom. We showed that there exist critical values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-surface layers of the fault, even in the absence of inelasticity. Here, we generalise the model to the case of finite deformations characteristic for the gouge deformation. We derive objective constitutive relationships relating the Jaumann rates of stress and moment stress to the relative strain and curvature rates, respectively. The model suggests that the pattern of localisation remains the same as in the linear case. However, the presence of the Jaumann terms leads to the emergence of non-zero normal stresses acting along and perpendicular to the shear layer (with zero hydrostatic pressure), and localised along the mid-line of the gouge; these stress components are absent in the linear model of simple shear. These additional normal stresses, albeit small, cause a change in the direction in which the maximal normal stresses act and in which en-echelon fracturing is formed.

Are European equity markets efficient? New evidence from fractal analysis

We report an empirical analysis of long-range dependence in the returns of eight stock market indices, using the Rescaled Range Analysis (RRA) to estimate the Hurst exponent. Monte Carlo and bootstrap simulations are used to construct critical values for the null hypothesis of no long-range dependence. The issue of disentangling short-range and long-range dependence is examined. Pre-filtering by fitting a (short-range) autoregressive model eliminates part of the long-range dependence when the latter is present, while failure to pre-filter leaves open the possibility of conflating short-range and long-range dependence. There is a strong evidence of long-range dependence for the small central European Czech stock market index PX-glob, and a weaker evidence for two smaller western European stock market indices, MSE (Spain) and SWX (Switzerland). There is little or no evidence of long-range dependence for the other five indices, including those with the largest capitalizations among those considered, DJIA (US) and FTSE350 (UK). These results are generally consistent with prior expectations concerning the relative efficiency of the stock markets examined. © 2011 Elsevier Inc.

Unifractality and multifractality in the Italian stock market

Tests for random walk behaviour in the Italian stock market are presented, based on an investigation of the fractal properties of the log return series for the Mibtel index. The random walk hypothesis is evaluated against alternatives accommodating either unifractality or multifractality. Critical values for the test statistics are generated using Monte Carlo simulations of random Gaussian innovations. Evidence is reported of multifractality, and the departure from random walk behaviour is statistically significant on standard criteria. The observed pattern is attributed primarily to fat tails in the return probability distribution, associated with volatility clustering in returns measured over various time scales. © 2009 Elsevier Inc. All rights reserved.

Performance measurement with multiple interrelated variables and threshold target levels:evidence from retail firms in the US

In this study, we developed a DEA-based performance measurement methodology that is consistent with performance assessment frameworks such as the Balanced Scorecard. The methodology developed in this paper takes into account the direct or inverse relationships that may exist among the dimensions of performance to construct appropriate production frontiers. The production frontiers we obtained are deemed appropriate as they consist solely of firms with desirable levels for all dimensions of performance. These levels should be at least equal to the critical values set by decision makers. The properties and advantages of our methodology against competing methodologies are presented through an application to a real-world case study from retail firms operating in the US. A comparative analysis between the new methodology and existing methodologies explains the failure of the existing approaches to define appropriate production frontiers when directly or inversely related dimensions of performance are present and to express the interrelationships between the dimensions of performance.

Avaliação e correções de testes de hipóteses em modelos de sobrevivência com fração de cura

Survival models deals with the modelling of time to event data. In certain situations, a share of the population can no longer be subjected to the event occurrence. In this context, the cure fraction models emerged. Among the models that incorporate a fraction of cured one of the most known is the promotion time model. In the present study we discuss hypothesis testing in the promotion time model with Weibull distribution for the failure times of susceptible individuals. Hypothesis testing in this model may be performed based on likelihood ratio, gradient, score or Wald statistics. The critical values are obtained from asymptotic approximations, which may result in size distortions in nite sample sizes. This study proposes bootstrap corrections to the aforementioned tests and Bartlett bootstrap to the likelihood ratio statistic in Weibull promotion time model. Using Monte Carlo simulations we compared the nite sample performances of the proposed corrections in contrast with the usual tests. The numerical evidence favors the proposed corrected tests. At the end of the work an empirical application is presented.

An Alternative Goodness-of-fit Test for Normality with Unknown Parameters

Goodness-of-fit tests have been studied by many researchers. Among them, an alternative statistical test for uniformity was proposed by Chen and Ye (2009). The test was used by Xiong (2010) to test normality for the case that both location parameter and scale parameter of the normal distribution are known. The purpose of the present thesis is to extend the result to the case that the parameters are unknown. A table for the critical values of the test statistic is obtained using Monte Carlo simulation. The performance of the proposed test is compared with the Shapiro-Wilk test and the Kolmogorov-Smirnov test. Monte-Carlo simulation results show that proposed test performs better than the Kolmogorov-Smirnov test in many cases. The Shapiro Wilk test is still the most powerful test although in some cases the test proposed in the present research performs better.