Constant coefficient tests for random coefficient regression

Random coefficient regression models have been applied in differentfields and they constitute a unifying setup for many statisticalproblems. The nonparametric study of this model started with Beranand Hall (1992) and it has become a fruitful framework. In thispaper we propose and study statistics for testing a basic hypothesisconcerning this model: the constancy of coefficients. The asymptoticbehavior of the statistics is investigated and bootstrapapproximations are used in order to determine the critical values ofthe test statistics. A simulation study illustrates the performanceof the proposals.

In vivo targeting of an anti-tumor antibody coupled to antigenic MHC class I complexes induces specific growth inhibition and regression of established syngeneic tumor grafts.

The concept of antibody-mediated targeting of antigenic MHC/peptide complexes on tumor cells in order to sensitize them to T-lymphocyte cytotoxicity represents an attractive new immunotherapy strategy. In vitro experiments have shown that an antibody chemically conjugated or fused to monomeric MHC/peptide can be oligomerized on the surface of tumor cells, rendering them susceptible to efficient lysis by MHC-peptide restricted specific T-cell clones. However, this strategy has not yet been tested entirely in vivo in immunocompetent animals. To this aim, we took advantage of OT-1 mice which have a transgenic T-cell receptor specific for the ovalbumin (ova) immunodominant peptide (257-264) expressed in the context of the MHC class I H-2K(b). We prepared and characterized conjugates between the Fab' fragment from a high-affinity monoclonal antibody to carcinoembryonic antigen (CEA) and the H-2K(b) /ova peptide complex. First, we showed in OT-1 mice that the grafting and growth of a syngeneic colon carcinoma line transfected with CEA could be specifically inhibited by systemic injections of the conjugate. Next, using CEA transgenic C57BL/6 mice adoptively transferred with OT-1 spleen cells and immunized with ovalbumin, we demonstrated that systemic injections of the anti-CEA-H-2K(b) /ova conjugate could induce specific growth inhibition and regression of well-established, palpable subcutaneous grafts from the syngeneic CEA-transfected colon carcinoma line. These results, obtained in a well-characterized syngeneic carcinoma model, demonstrate that the antibody-MHC/peptide strategy can function in vivo. Further preclinical experimental studies, using an anti-viral T-cell response, will be performed before this new form of immunotherapy can be considered for clinical use.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Solving higher-dimensional continuous time stochastic control problems by value function regression

The paper develops a method to solve higher-dimensional stochasticcontrol problems in continuous time. A finite difference typeapproximation scheme is used on a coarse grid of low discrepancypoints, while the value function at intermediate points is obtainedby regression. The stability properties of the method are discussed,and applications are given to test problems of up to 10 dimensions.Accurate solutions to these problems can be obtained on a personalcomputer.

Weighted Kernel regression

In the fixed design regression model, additional weights areconsidered for the Nadaraya--Watson and Gasser--M\"uller kernel estimators.We study their asymptotic behavior and the relationships between new andclassical estimators. For a simple family of weights, and considering theIMSE as global loss criterion, we show some possible theoretical advantages.An empirical study illustrates the performance of the weighted estimatorsin finite samples.

How wide is the gap? An investigation of gender wage differences using quantile regression

In this paper we examine the determinants of wages and decompose theobserved differences across genders into the "explained by differentcharacteristics" and "explained by different returns components"using a sample of Spanish workers. Apart from the conditionalexpectation of wages, we estimate the conditional quantile functionsfor men and women and find that both the absolute wage gap and thepart attributed to different returns at each of the quantiles, farfrom being well represented by their counterparts at the mean, aregreater as we move up in the wage range.

Fixed and random effects in Classical and Bayesian regression

This paper proposes a common and tractable framework for analyzingdifferent definitions of fixed and random effects in a contant-slopevariable-intercept model. It is shown that, regardless of whethereffects (i) are treated as parameters or as an error term, (ii) areestimated in different stages of a hierarchical model, or whether (iii)correlation between effects and regressors is allowed, when the sameinformation on effects is introduced into all estimation methods, theresulting slope estimator is also the same across methods. If differentmethods produce different results, it is ultimately because differentinformation is being used for each methods.

The importance of individual heterogeneity in the decomposition of measures of socioeconomic inequality in health: An approach based on quantile regression

This paper shows how recently developed regression-based methods for thedecomposition of health inequality can be extended to incorporateindividual heterogeneity in the responses of health to the explanatoryvariables. We illustrate our method with an application to the CanadianNPHS of 1994. Our strategy for the estimation of heterogeneous responsesis based on the quantile regression model. The results suggest that thereis an important degree of heterogeneity in the association of health toexplanatory variables which, in turn, accounts for a substantial percentageof inequality in observed health. A particularly interesting finding isthat the marginal response of health to income is zero for healthyindividuals but positive and significant for unhealthy individuals. Theheterogeneity in the income response reduces both overall health inequalityand income related health inequality.

Random measurement error and regression dilution bias.

Summary points: - The bias introduced by random measurement error will be different depending on whether the error is in an exposure variable (risk factor) or outcome variable (disease) - Random measurement error in an exposure variable will bias the estimates of regression slope coefficients towards the null - Random measurement error in an outcome variable will instead increase the standard error of the estimates and widen the corresponding confidence intervals, making results less likely to be statistically significant - Increasing sample size will help minimise the impact of measurement error in an outcome variable but will only make estimates more precisely wrong when the error is in an exposure variable

Estimation of percentage body fat in 6- to 13-year-old children by skinfold thickness, body mass index and waist circumference.

We evaluated the accuracy of skinfold thicknesses, BMI and waist circumference for the prediction of percentage body fat (PBF) in a representative sample of 372 Swiss children aged 6-13 years. PBF was measured using dual-energy X-ray absorptiometry. On the basis of a preliminary bootstrap selection of predictors, seven regression models were evaluated. All models included sex, age and pubertal stage plus one of the following predictors: (1) log-transformed triceps skinfold (logTSF); (2) logTSF and waist circumference; (3) log-transformed sum of triceps and subscapular skinfolds (logSF2); (4) log-transformed sum of triceps, biceps, subscapular and supra-iliac skinfolds (logSF4); (5) BMI; (6) waist circumference; (7) BMI and waist circumference. The adjusted determination coefficient (R² adj) and the root mean squared error (RMSE; kg) were calculated for each model. LogSF4 (R² adj 0.85; RMSE 2.35) and logSF2 (R² adj 0.82; RMSE 2.54) were similarly accurate at predicting PBF and superior to logTSF (R² adj 0.75; RMSE 3.02), logTSF combined with waist circumference (R² adj 0.78; RMSE 2.85), BMI (R² adj 0.62; RMSE 3.73), waist circumference (R² adj 0.58; RMSE 3.89), and BMI combined with waist circumference (R² adj 0.63; RMSE 3.66) (P < 0.001 for all values of R² adj). The finding that logSF4 was only modestly superior to logSF2 and that logTSF was better than BMI and waist circumference at predicting PBF has important implications for paediatric epidemiological studies aimed at disentangling the effect of body fat on health outcomes.

Effect of weight loss on resting energy expenditure in obese prepubertal children.

To assess the effect of weight loss on resting metabolic rate (RMR), the energy expenditure of eight obese prepubertal children (age 9 +/- 1 years; weight 48.7 +/- 9.1 kg; BMI 25.3 +/- 3.9) and of 14 age-matched children of normal body weight (age 9 +/- 1 years; weight 28.8 +/- 5.6 kg; BMI 16.5 +/- 1.7) was measured by indirect calorimetry. The obese children were reinvestigated after a mean weight loss of 5.4 +/- 1.2 kg induced by a six-months mixed hypocaloric diet. Before slimming, the obese group showed a higher daily energy intake than the control group (10.40 +/- 3.45 MJ/day vs 7.97 +/- 2.02 MJ/day respectively; P less than 0.05) but a similar value was observed per unit fat-free mass (FFM) (0.315 +/- 0.032 MJ/kgFFM/day vs 0.329 +/- 0.041 MJ/kgFFM/day respectively). The average RMR of the obese children was greater than that of the control group (5217 +/- 531 kJ/day vs 4477 +/- 506 kJ/day) but similar after adjusting for FFM (4728 +/- 3102 kJ/day vs 4899 +/- 3102 kJ/day). Weight loss resulted in a reduction in RMR (5217 +/- 531 kJ/day vs 4874 +/- 820 kJ/day), each kg of weight loss being accompanied by a decrease of RMR of 64 kJ (15.3 kcal) per day. The changes in RMR induced by weight loss paralleled the changes in FFM. No difference was found in average RQ in obese children vs controls (0.85 +/- 0.03 vs 0.87 +/- 0.03 respectively) and in the obese children before and after weight loss (0.87 +/- 0.02).(ABSTRACT TRUNCATED AT 250 WORDS)

Prediction of 18-month survival in patients with primary myelodysplastic syndrome. A regression model and scoring system based on the combination of chromosome findings and the Bournemouth score.

The predictive potential of six selected factors was assessed in 72 patients with primary myelodysplastic syndrome using univariate and multivariate logistic regression analysis of survival at 18 months. Factors were age (above median of 69 years), dysplastic features in the three myeloid bone marrow cell lineages, presence of chromosome defects, all metaphases abnormal, double or complex chromosome defects (C23), and a Bournemouth score of 2, 3, or 4 (B234). In the multivariate approach, B234 and C23 proved to be significantly associated with a reduction in the survival probability. The similarity of the regression coefficients associated with these two factors means that they have about the same weight. Consequently, the model was simplified by counting the number of factors (0, 1, or 2) present in each patient, thus generating a scoring system called the Lausanne-Bournemouth score (LB score). The LB score combines the well-recognized and easy-to-use Bournemouth score (B score) with the chromosome defect complexity, C23 constituting an additional indicator of patient outcome. The predicted risk of death within 18 months calculated from the model is as follows: 7.1% (confidence interval: 1.7-24.8) for patients with an LB score of 0, 60.1% (44.7-73.8) for an LB score of 1, and 96.8% (84.5-99.4) for an LB score of 2. The scoring system presented here has several interesting features. The LB score may improve the predictive value of the B score, as it is able to recognize two prognostic groups in the intermediate risk category of patients with B scores of 2 or 3. It has also the ability to identify two distinct prognostic subclasses among RAEB and possibly CMML patients. In addition to its above-described usefulness in the prognostic evaluation, the LB score may bring new insights into the understanding of evolution patterns in MDS. We used the combination of the B score and chromosome complexity to define four classes which may be considered four possible states of myelodysplasia and which describe two distinct evolutional pathways.

Método convencional e "bootstrap" para estimar o número de observações na determinação dos parâmetros da função K(Θ)

A variabilidade natural do solo dificulta a obtenção de valores que representem adequadamente as propriedades do solo em determinada área. O conhecimento do número mínimo de observações, que devem ser realizadas para representar, com um erro aceitável, uma propriedade ou característica do solo, é fundamental para que os resultados experimentais possam ser aplicados com segurança. No presente trabalho, comparou-se o método convencional (teste-t) com o método "bootstrap", com vistas em estimar o número de observações necessárias para calcular os parâmetros que caracterizam a relação entre a condutividade hidráulica e o teor de água do solo, determinada pelo método do perfil instantâneo. Realizou-se um experimento de drenagem num Latossolo Vermelho-Amarelo em Piracicaba (SP), numa parcela experimental com 45 pontos de observação distanciados de 1 m entre si. Observaram-se a umidade (com TDR) e o potencial mátrico (com tensiômetros) durante três semanas de redistribuição da água. Após processamento dos dados, o conjunto de valores mostrou uma distribuição não-normal, fazendo-se necessária a eliminação de "outliers" para a aplicação do método convencional, normalizando a distribuição. Assim, o uso do método tradicional só é recomendado após a confirmação da pertinência da eliminação dos "outliers". Ambos os métodos de análise requerem grande número de repetições, reafirmando que determinações da função condutividade hidráulica com poucas repetições não podem ser extrapoladas para áreas maiores.