Principal differential analysis with a continuous covariate: low-dimensional approximations for functional data

Given a collection of n curves that are independent realizations of a functional variable, we are interested in finding patterns in the curve data by exploring low-dimensional approximations to the curves. It is assumed that the data curves are noisy samples from the vector space span <texlscub>f ₁, …, f _m </texlscub>, where f ₁, …, f _m are unknown functions on the real interval (0, T) with square-integrable derivatives of all orders m or less, and m<n. Ramsay [Principal differential analysis: Data reduction by differential operators, J. R. Statist. Soc. Ser. B 58 (1996), pp. 495–508] first proposed the method of regularized principal differential analysis (PDA) as an alternative to principal component analysis for finding low-dimensional approximations to curves. PDA is based on the following theorem: there exists an annihilating linear differential operator (LDO) ? of order m such that ?f _i =0, i=1, …, m [E.A. Coddington and N. Levinson, Theory of Ordinary Differential Equations, McGraw-Hill, New York, 1955, Theorem 6.2]. PDA specifies m, then uses the data to estimate an annihilating LDO. Smooth estimates of the coefficients of the LDO are obtained by minimizing a penalized sum of the squared norm of the residuals. In this context, the residual is that part of the data curve that is not annihilated by the LDO. PDA obtains the smooth low dimensional approximation to the data curves by projecting onto the null space of the estimated annihilating LDO; PDA is thus useful for obtaining low-dimensional approximations to the data curves whether or not the interpretation of the annihilating LDO is intuitive or obvious from the context of the data. This paper extends PDA to allow for the coefficients in the LDO to smoothly depend upon a single continuous covariate. The estimating equations for the coefficients allowing for a continuous covariate are derived; the penalty of Eilers and Marx [Flexible smoothing with B-splines and penalties, Statist. Sci. 11(2) (1996), pp. 89–121] is used to impose smoothness. The results of a small computer simulation study investigating the bias and variance properties of the estimator are reported.     

Foraging in the tropics: relationships among species’ abundances,niche asymmetries and body condition in an urban avian assemblage

Two main theories attempt to explain species coexistence: the neutral theory considers all the species as equivalents so biodiversity is mainly regarded as a function of total available resources (i.e. niche expansion), while the niche theory stresses the relevance of differences in niche use between species (i.e. niche packing). The relative importance of these forces is under discussion and has been largely tested in natural ecosystems. However, few studies have addressed this issue in tropical-urban environments. In this work we studied niche overlap asymmetries among the most common urban resident birds, and the effect of habitat type on this pattern, in a subtropical location of South China (Nanning, Guangxi). We found differences in abundances and niche use among species and urban habitats (parks, streets and orchards). We also recorded strong asymmetries in niche use between species, which we divided into three categories: species with positive asymmetries, which showed highly specific niches and were able to exploit other species’ niches; species with neutral asymmetries, which showed high levels of niche overlap with the rest of species; and species with negative asymmetries, which showed low specificity in their niches and were unable to exploit other species’ niches. These differences in niche use correlated with differences in fat scores. Species from the first group showed higher fat scores than other two groups. Ultimately, niche asymmetries correlated with species’ body condition and mediated their differences in abundance, which supports the view that in this urban context niche theory is more appropriate than neutral theory.     

Methods for estimating relative risk in studies of common binary outcomes

Studying the effect of exposure or intervention on a dichotomous outcome is very common in medical research. Logistic regression (LR) is often used to determine such association which provides odds ratio (OR). OR often overestimates the effect size for prevalent outcome data. In such situations, use of relative risk (RR) has been suggested. We propose modifications in Zhang and Yu and Diaz-Quijano methods. These methods were compared with stratified Mantel Haenszel method, LR, log binomial regression (LBR), Zhang and Yu method, Poisson/Cox regression, modified Poisson/Cox regression, marginal probability method, COPY method, inverse probability of treatment weighted LBR, and Diaz-Quijano method. Our proposed modified Diaz-Quijano (MDQ) method provides RR and its confidence interval similar to those estimated by modified Poisson/Cox and LBRs. The proposed modifications in Zhang and Yu method provides better estimate of RR and its standard error as compared to Zhang and Yu method in a variety of situations with prevalent outcome. The MDQ method can be used easily to estimate the RR and its confidence interval in the studies which require reporting of RRs. Regression models which directly provide the estimate of RR without convergence problems such as the MDQ method and modified Poisson/Cox regression should be preferred.