Subsampling the mean of heavy-tailed dependent observations

We establish the validity of subsampling confidence intervals for themean of a dependent series with heavy-tailed marginal distributions.Using point process theory, we study both linear and nonlinear GARCH-liketime series models. We propose a data-dependent method for the optimalblock size selection and investigate its performance by means of asimulation study.

Asymptotic robust inferences in the analysis of mean and covariance structures

Structural equation models are widely used in economic, socialand behavioral studies to analyze linear interrelationships amongvariables, some of which may be unobservable or subject to measurementerror. Alternative estimation methods that exploit different distributionalassumptions are now available. The present paper deals with issues ofasymptotic statistical inferences, such as the evaluation of standarderrors of estimates and chi--square goodness--of--fit statistics,in the general context of mean and covariance structures. The emphasisis on drawing correct statistical inferences regardless of thedistribution of the data and the method of estimation employed. A(distribution--free) consistent estimate of $\Gamma$, the matrix ofasymptotic variances of the vector of sample second--order moments,will be used to compute robust standard errors and a robust chi--squaregoodness--of--fit squares. Simple modifications of the usual estimateof $\Gamma$ will also permit correct inferences in the case of multi--stage complex samples. We will also discuss the conditions under which,regardless of the distribution of the data, one can rely on the usual(non--robust) inferential statistics. Finally, a multivariate regressionmodel with errors--in--variables will be used to illustrate, by meansof simulated data, various theoretical aspects of the paper.

Option pricing under stochastic volatility and stochastic interest rate in the Spanish case

Among the underlying assumptions of the Black-Scholes option pricingmodel, those of a fixed volatility of the underlying asset and of aconstantshort-term riskless interest rate, cause the largest empirical biases. Onlyrecently has attention been paid to the simultaneous effects of thestochasticnature of both variables on the pricing of options. This paper has tried toestimate the effects of a stochastic volatility and a stochastic interestrate inthe Spanish option market. A discrete approach was used. Symmetricand asymmetricGARCH models were tried. The presence of in-the-mean and seasonalityeffectswas allowed. The stochastic processes of the MIBOR90, a Spanishshort-terminterest rate, from March 19, 1990 to May 31, 1994 and of the volatilityofthe returns of the most important Spanish stock index (IBEX-35) fromOctober1, 1987 to January 20, 1994, were estimated. These estimators wereused onpricing Call options on the stock index, from November 30, 1993 to May30, 1994.Hull-White and Amin-Ng pricing formulas were used. These prices werecomparedwith actual prices and with those derived from the Black-Scholesformula,trying to detect the biases reported previously in the literature. Whereasthe conditional variance of the MIBOR90 interest rate seemed to be freeofARCH effects, an asymmetric GARCH with in-the-mean and seasonalityeffectsand some evidence of persistence in variance (IEGARCH(1,2)-M-S) wasfoundto be the model that best represent the behavior of the stochasticvolatilityof the IBEX-35 stock returns. All the biases reported previously in theliterature were found. All the formulas overpriced the options inNear-the-Moneycase and underpriced the options otherwise. Furthermore, in most optiontrading, Black-Scholes overpriced the options and, because of thetime-to-maturityeffect, implied volatility computed from the Black-Scholes formula,underestimatedthe actual volatility.

Shrews in suburbia: an application of Goodman's extinction model

An equation is applied for calculating the expected persistence time of an unstructured population of the white-toothed shrew Crocidura russula from Preverenges, a suburban area in western Switzerland. Population abundance data from March and November between 1977 and 1988 were fit to the logistic density dependence model to estimate mean population growth rate as a function of population density. The variance in mean growth rate was approximated with two different models. The largest estimated persistence time was less than a few decades, the smallest less than 10 years. The results are sensitive to the magnitude of variance in population growth rate. Deviations from the logistic density dependence model in November are quite well explained by weather variables but those in March are uncorrelated with weather variables. Variability in population growth rates measured in winter months may be better explained by behavioural mechanisms. Environmental variability, dispersal of juveniles and refugia within the range of the population may contribute to its long-term survival.

25-hydroxyvitamin D concentration is inversely associated with serum MMP-9 in a cross-sectional study of African American ESRD patients.

BACKGROUND: Circulating 25-hydroxyvitamin D [25(OH)D] concentration is inversely associated with peripheral arterial disease and hypertension. Vascular remodeling may play a role in this association, however, data relating vitamin D level to specific remodeling biomarkers among ESRD patients is sparse. We tested whether 25(OH)D concentration is associated with markers of vascular remodeling and inflammation in African American ESRD patients.METHODS: We conducted a cross-sectional study among ESRD patients receiving maintenance hemodialysis within Emory University-affiliated outpatient hemodialysis units. Demographic, clinical and dialysis treatment data were collected via direct patient interview and review of patients records at the time of enrollment, and each patient gave blood samples. Associations between 25(OH)D and biomarker concentrations were estimated in univariate analyses using Pearson's correlation coefficients and in multivariate analyses using linear regression models. 25(OH) D concentration was entered in multivariate linear regression models as a continuous variable and binary variable (<15 ng/ml and =15 ng/ml). Adjusted estimate concentrations of biomarkers were compared between 25(OH) D groups using analysis of variance (ANOVA). Finally, results were stratified by vascular access type.RESULTS: Among 91 patients, mean (standard deviation) 25(OH)D concentration was 18.8 (9.6) ng/ml, and was low (<15 ng/ml) in 43% of patients. In univariate analyses, low 25(OH) D was associated with lower serum calcium, higher serum phosphorus, and higher LDL concentrations. 25(OH) D concentration was inversely correlated with MMP-9 concentration (r = -0.29, p = 0.004). In multivariate analyses, MMP-9 concentration remained negatively associated with 25(OH) D concentration (P = 0.03) and anti-inflammatory IL-10 concentration positively correlated with 25(OH) D concentration (P = 0.04).CONCLUSIONS: Plasma MMP-9 and circulating 25(OH) D concentrations are significantly and inversely associated among ESRD patients. This finding may suggest a potential mechanism by which low circulating 25(OH) D functions as a cardiovascular risk factor.

The disturbing 'rise' of global income inequality

We use aggregate GDP data and within-country income shares for theperiod 1970-1998 to assign a level of income to each person in theworld. We then estimate the gaussian kernel density function for theworldwide distribution of income. We compute world poverty rates byintegrating the density function below the poverty lines. The $1/daypoverty rate has fallen from 20% to 5% over the last twenty five years.The $2/day rate has fallen from 44% to 18%. There are between 300 and500 million less poor people in 1998 than there were in the 70s.We estimate global income inequality using seven different popularindexes: the Gini coefficient, the variance of log-income, two ofAtkinson s indexes, the Mean Logarithmic Deviation, the Theil indexand the coefficient of variation. All indexes show a reduction in globalincome inequality between 1980 and 1998. We also find that most globaldisparities can be accounted for by across-country, not within-country,inequalities. Within-country disparities have increased slightly duringthe sample period, but not nearly enough to offset the substantialreduction in across-country disparities. The across-country reductionsin inequality are driven mainly, but not fully, by the large growth rateof the incomes of the 1.2 billion Chinese citizens. Unless Africa startsgrowing in the near future, we project that income inequalities willstart rising again. If Africa does not start growing, then China, India,the OECD and the rest of middle-income and rich countries diverge awayfrom it, and global inequality will rise. Thus, the aggregate GDP growthof the African continent should be the priority of anyone concerned withincreasing global income inequality.

Principal curves and principal oriented points

Principal curves have been defined Hastie and Stuetzle (JASA, 1989) assmooth curves passing through the middle of a multidimensional dataset. They are nonlinear generalizations of the first principalcomponent, a characterization of which is the basis for the principalcurves definition.In this paper we propose an alternative approach based on a differentproperty of principal components. Consider a point in the space wherea multivariate normal is defined and, for each hyperplane containingthat point, compute the total variance of the normal distributionconditioned to belong to that hyperplane. Choose now the hyperplaneminimizing this conditional total variance and look for thecorresponding conditional mean. The first principal component of theoriginal distribution passes by this conditional mean and it isorthogonal to that hyperplane. This property is easily generalized todata sets with nonlinear structure. Repeating the search from differentstarting points, many points analogous to conditional means are found.We call them principal oriented points. When a one-dimensional curveruns the set of these special points it is called principal curve oforiented points. Successive principal curves are recursively definedfrom a generalization of the total variance.

Resting metabolic rate in six- to ten-year-old obese and nonobese children.

The resting metabolic rate (RMR) and body composition of 130 obese and nonobese prepubertal children, aged 6 to 10 years, were assessed by indirect calorimetry and skin-fold thickness, respectively. The mean (+/- SD) RMR was 4619 +/- 449 kJ.day-1 (164 +/- 31 kJ.kg body weight-1 x day-1) in the 62 boys and 4449 +/- 520 kJ.day-1 (147 +/- 32 kJ.kg body weight-1 x day-1) in the 68 girls. Fat-free mass was the best single predictor of RMR (R2 = 0.64; p < 0.001). Step-down multiple regression analysis, with independent variables such as age, gender, weight, and height, allowed several RMR predictive equations to be developed. An equation for boys is as follows: RMR (kJ.day-1) = 1287 + 28.6 x Weight(kg) + 23.6 x Height(cm) - 69.1 x Age(yr) (R2 = 0.58; p < 0.001). An equation for girls is as follows: RMR (kJ.day-1 = 1552 + 35.8 x Weight (kg) + 15.6 x Height (cm) - 36.3 x Age (yr) (R2 = 0.69; p < 0.001). Comparison between the measured RMR and that predicted by currently used formulas showed that most of these equations tended to overestimate the RMR of both genders, especially in overweight children.

Breeding system and genetic variance in the monogamous, semi-social shrew, Crocidura russula

The population-genetic consequences of monogamy and male philopatry (a rare breeding system in mammals) were investigated using microsatellite markers in the semisocial and anthropophilic shrew Crocidura russula. A hierarchical sampling design over a 16-km geographical transect revealed a large genetic diversity (h = 0.813) with significant differentiation among subpopulations (F-ST = 5-6%), which suggests an exchange of 4.4 migrants per generation. Demic effective-size estimates were very high, due both to this limited gene inflow and to the inner structure of subpopulations. These were made of 13-20 smaller units (breeding groups), comprising an estimate of four breeding pairs each. Members of the same breeding groups displayed significant coancestries (F-LS = 9-10%), which was essentially due to strong male kinship: syntopic males were on average related at the half-sib level. Female dispersal among breeding groups was not complete (similar to 39%), and insufficient to prevent inbreeding. From our results, the breeding strategy of C. russula seems less efficient than classical mammalian systems (polygyny and male dispersal) in disentangling coancestry from inbreeding, but more so in retaining genetic variance.

Variance reduction methods for simulation of densities on Wiener space

We develop a general error analysis framework for the Monte Carlo simulationof densities for functionals in Wiener space. We also study variancereduction methods with the help of Malliavin derivatives. For this, wegive some general heuristic principles which are applied to diffusionprocesses. A comparison with kernel density estimates is made.

A two-mean reverting-factor model of the term structure of interest rates

This paper presents a two--factor model of the term structure ofinterest rates. We assume that default free discount bond prices aredetermined by the time to maturity and two factors, the long--term interestrate and the spread (difference between the long--term rate and theshort--term (instantaneous) riskless rate). Assuming that both factorsfollow a joint Ornstein--Uhlenbeck process, a general bond pricing equationis derived. We obtain a closed--form expression for bond prices andexamine its implications for the term structure of interest rates. We alsoderive a closed--form solution for interest rate derivatives prices. Thisexpression is applied to price European options on discount bonds andmore complex types of options. Finally, empirical evidence of the model'sperformance is presented.

An empirical evaluation of small area estimators

This paper investigates the comparative performance of five small areaestimators. We use Monte Carlo simulation in the context of boththeoretical and empirical populations. In addition to the direct andindirect estimators, we consider the optimal composite estimator withpopulation weights, and two composite estimators with estimatedweights: one that assumes homogeneity of within area variance andsquare bias, and another one that uses area specific estimates ofvariance and square bias. It is found that among the feasibleestimators, the best choice is the one that uses area specificestimates of variance and square bias.

Heidelberg Retina Tomograph in Glaucoma Case-finding

Purpose: To compare the disease classification using optic nerve head imaging against clinical examination combined with perimetry in the case-finding environment.Methods: During a glaucoma screening event at the University Hospital in Lausanne 148 attendees opted to undergo an extended screening exam consisting of perimetry with Octopus (Haag-Streit, Koeniz, Switzerland), tonometry, gonioscopy, slit-lamp optic nerve head examination and Heidelberg Retinal Tomorgraph (HRT, Heidelberg engineering). Classification criteria of an abnormal clinical exam consisted of a cup to disc ratio > 0.6, other disc abnormality (e.g. notching), closed angles, Van Herrick < ¼ and intraocular pressure >21mmHg. In perimetry a square root loss variance (sLV) >3.0 dB was considered abnormal. Those cases with both an abnormal clinical exam and abnormal perimetric results were considered to be glaucoma suspects. An abnormal HRT exam was a global Moorfields Regression Analysis (MRA) result of "outside normal limits". The results from the worse eye for each instrument are reported.Results: The mean age of the patients was 59.9 years (SD ±14.8years). 46 subjects failed the clinical exam. 55 subjects had sLV>3.0dB. 86 subjects failed one or more of the testing components, 21 subjects failed both the clinical exam and perimetry. Two cases of advanced glaucoma were diagnosed on the day of the event. 20 subjects were referred due to other ocular pathology (10 cataract; 2 AMD; 8 other). 16 subjects were outside normal limits with the HRT, of which 11 failed the clinical exam, 9 failed the perimetric exam and 7 failed all three. HRT had a sensitivity of 33% CI [15%,57%] with a specificity of 93% CI [87%,97%].Conclusions: HRT shows good specificity, however the low sensitivity makes it of limited use in the proposed case finding scenario.

Ambulatory blood pressure monitoring: a mean to stratify cardiovascular risk.

OBJECTIVE: Current hypertension guidelines stress the importance to assess total cardiovascular risk but do not describe precisely how to use ambulatory blood pressures in the cardiovascular risk stratification. METHOD: We calculated here global cardiovascular risk according to 2003 European Society of Hypertension/European Society of Cardiology guidelines in 127 patients in whom daytime ambulatory blood pressures were recorded and carotid/femoral ultrasonography performed. RESULTS: The presence of ambulatory blood pressures >or =135/85 mmHg shifted cardiovascular risk to higher categories, as did the presence of hypercholesterolemia and, even more so, the presence of atherosclerotic plaques. CONCLUSION: Further studies are, however, needed to define the position of ambulatory blood pressures in the assessment of cardiovascular risk.