Estimation of Gaussian process regression model using probability distance measures

A new class of parameter estimation algorithms is introduced for Gaussian process regression (GPR) models. It is shown that the integration of the GPR model with probability distance measures of (i) the integrated square error and (ii) Kullback–Leibler (K–L) divergence are analytically tractable. An efficient coordinate descent algorithm is proposed to iteratively estimate the kernel width using golden section search which includes a fast gradient descent algorithm as an inner loop to estimate the noise variance. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Combining satellite observations and reanalysis energy transports to estimate global net surface energy fluxes 1985-2012

Two methods are developed to estimate net surface energy fluxes based upon satellite-based reconstructions of radiative fluxes at the top of atmosphere and the atmospheric energy tendencies and transports from the ERA-Interim reanalysis. Method 1 applies the mass adjusted energy divergence from ERA-Interim while method 2 estimates energy divergence based upon the net energy difference at the top of atmosphere and the surface from ERA-Interim. To optimise the surface flux and its variability over ocean, the divergences over land are constrained to match the monthly area mean surface net energy flux variability derived from a simple relationship between the surface net energy flux and the surface temperature change. The energy divergences over the oceans are then adjusted to remove an unphysical residual global mean atmospheric energy divergence. The estimated net surface energy fluxes are compared with other data sets from reanalysis and atmospheric model simulations. The spatial correlation coefficients of multi-annual means between the estimations made here and other data sets are all around 0.9. There are good agreements in area mean anomaly variability over the global ocean, but discrepancies in the trend over the eastern Pacific are apparent.

Uniform factorial decay estimate for the remainder of rough Taylor expansion

We establish an uniform factorial decay estimate for the Taylor approximation of solutions to controlled differential equations. Its proof requires a factorial decay estimate for controlled paths which is interesting in its own right.

Using Africa’s protected area network to estimate the global population of a threatened and declining species: a case study of the critically endangered White-headed Vulture Trigonoceps occipitalis

The White-headed Vulture Trigonoceps occipitalis (WhV) is uncommon and largely restricted to protected areas across its range in sub-Saharan Africa. We used the World Database on Protected Areas to identify protected areas (PAs) likely to contain White-headed Vultures. Vulture occurrence on road transects in Southern, East, and West Africa was adjusted to nests per km2 using data from areas with known numbers of nests and corresponding road transect data. Nest density was used to calculate the number of WhV nests within identified PAs and from there extrapolated to estimate the global population. Across a fragmented range, 400 PAs are estimated to contain 1893 WhV nests. Eastern Africa is estimated to contain 721 nests, Central Africa 548 nests, Southern Africa 468 nests, and West Africa 156 nests. Including immature and nonbreeding birds, and accounting for data deficient PAs, the estimated global population is 5475 - 5493 birds. The identified distribution highlights are alarming: over 78% (n = 313) of identified PAs contain fewer than five nests. A further 17% (n = 68) of PAs contain 5 - 20 nests and 4% (n = 14) of identified PAs are estimated to contain >20 nests. Just 1% (n = 5) of PAs are estimated to contain >40 nests; none is located in West Africa. Whilst ranging behavior of WhVs is currently unknown, 35% of PAs large enough to hold >20 nests are isolated by more than 100 km from other PAs. Spatially discrete and unpredictable mortality events such as poisoning pose major threats to small localized vulture populations and will accelerate ongoing local extinctions. Apart from reducing the threat of poisoning events, conservation actions promoting linkages between protected areas should be pursued. Identifying potential areas for assisted re-establishment via translocation offers the potential to expand the range of this species and alleviate risk.

An improved estimate of leaf area index based on the histogram analysis of hemispherical photographs

Leaf area index (LAI) is a key parameter that affects the surface fluxes of energy, mass, and momentum over vegetated lands, but observational measurements are scarce, especially in remote areas with complex canopy structure. In this paper we present an indirect method to calculate the LAI based on the analyses of histograms of hemispherical photographs. The optimal threshold value (OTV), the gray-level required to separate the background (sky) and the foreground (leaves), was analytically calculated using the entropy crossover method (Sahoo, P.K., Slaaf, D.W., Albert, T.A., 1997. Threshold selection using a minimal histogram entropy difference. Optical Engineering 36(7) 1976-1981). The OTV was used to calculate the LAI using the well-known gap fraction method. This methodology was tested in two different ecosystems, including Amazon forest and pasturelands in Brazil. In general, the error between observed and calculated LAI was similar to 6%. The methodology presented is suitable for the calculation of LAI since it is responsive to sky conditions, automatic, easy to implement, faster than commercially available software, and requires less data storage. (C) 2008 Elsevier B.V. All rights reserved.

Virtual Pole from Magnetic Anomaly (VPMA): A procedure to estimate the age of a rock from its magnetic anomaly only

Virtual Pole from Magnetic Anomaly (VPMA) is a new multi-disciplinary methodology that estimates the age of a source rock from its magnetic anomaly, taken directly from available aeromagnetic data. The idea is to use those anomalies in which a strong remanent magnetic component is likely to occur. Once the total magnetization of the anomaly is computed through any of the currently available methods, the line that connects all virtual paleogeographic poles is related with the position, on a paleogeographic projection, of the appropriate age fragment of the APWT curve. We applied this procedure to five (5) well-known magnetic anomalies of the South American plate in SE Brazil, all of them associated to alkaline complexes of Mesozoic age. The apparent ages obtained from VPMA on three of the anomalies where the radiometric age of the source rock is known - Tapira, Araxa and Juquia were inside the error interval of the published ages. The VPMA apparent ages of the other two, where the age of the source rock is not known (Registro and Pariqueracu magnetic anomalies) were geologically coherent. We expect that the application of the VPMA methodology as a reconnaissance geochronological tool may contribute to geological knowledge over continental areas, especially when the source rocks of the magnetic anomalies am unknown or buried below superficial sediments. (C) 2009 Elsevier B.V. All rights reserved.

The Diffusion Kernel Filter

A particle filter method is presented for the discrete-time filtering problem with nonlinear ItA ` stochastic ordinary differential equations (SODE) with additive noise supposed to be analytically integrable as a function of the underlying vector-Wiener process and time. The Diffusion Kernel Filter is arrived at by a parametrization of small noise-driven state fluctuations within branches of prediction and a local use of this parametrization in the Bootstrap Filter. The method applies for small noise and short prediction steps. With explicit numerical integrators, the operations count in the Diffusion Kernel Filter is shown to be smaller than in the Bootstrap Filter whenever the initial state for the prediction step has sufficiently few moments. The established parametrization is a dual-formula for the analysis of sensitivity to gaussian-initial perturbations and the analysis of sensitivity to noise-perturbations, in deterministic models, showing in particular how the stability of a deterministic dynamics is modeled by noise on short times and how the diffusion matrix of an SODE should be modeled (i.e. defined) for a gaussian-initial deterministic problem to be cast into an SODE problem. From it, a novel definition of prediction may be proposed that coincides with the deterministic path within the branch of prediction whose information entropy at the end of the prediction step is closest to the average information entropy over all branches. Tests are made with the Lorenz-63 equations, showing good results both for the filter and the definition of prediction.

Consistency in approximate entropy given by a volumetric estimate

Non-linear methods for estimating variability in time-series are currently of widespread use. Among such methods are approximate entropy (ApEn) and sample approximate entropy (SampEn). The applicability of ApEn and SampEn in analyzing data is evident and their use is increasing. However, consistency is a point of concern in these tools, i.e., the classification of the temporal organization of a data set might indicate a relative less ordered series in relation to another when the opposite is true. As highlighted by their proponents themselves, ApEn and SampEn might present incorrect results due to this lack of consistency. In this study, we present a method which gains consistency by using ApEn repeatedly in a wide range of combinations of window lengths and matching error tolerance. The tool is called volumetric approximate entropy, vApEn. We analyze nine artificially generated prototypical time-series with different degrees of temporal order (combinations of sine waves, logistic maps with different control parameter values, random noises). While ApEn/SampEn clearly fail to consistently identify the temporal order of the sequences, vApEn correctly do. In order to validate the tool we performed shuffled and surrogate data analysis. Statistical analysis confirmed the consistency of the method. (C) 2008 Elsevier Ltd. All rights reserved.

First molecular estimate of sex-ratio of southern right whale calves, Eubalaena australis, for Brazilian waters

The southern right whale (Eubalaena australis) was one of the most intensively hunted whales between the 17th and 20th centuries in the southern hemisphere. Recent estimates indicate that today there are around 7000 whales, representing 5 to 10% Of its original population. On the other hand, recent studies estimated that the population that migrates to the Brazilian coast grew by 14% from 1987 to 2003. However, there is no information about sex-ratio for adults or for calves in this region, which is an important parameter for understanding the biology of the species. We present here the first estimate Of calves` sex-ratio of southern right whales found along the southern Brazilian coast, one of the most important wintering grounds for the species. Sex was molecularly indentified for 21 biopsies collected from calves between 1998 and 2002, along the coast of Rio Grande do Sul and Santa Catarina States, in southern Brazil. The sex-ratio was two females for one male, however, it was not statistically different (chi(2) test, alpha = 0.05; df = 1) from the expected ratio of 1:1. This result is in accordance with the sex-ratio estimated for the species of all ages using external morphology (and behaviour in formation), (is well as for most species of baleen whales.

Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

The univalent polynomial of Suffridge as a summability kernel

A positive summability trigonometric kernel {K(n)(theta)}(infinity)(n=1) is generated through a sequence of univalent polynomials constructed by Suffridge. We prove that the convolution {K(n) * f} approximates every continuous 2 pi-periodic function f with the rate omega(f, 1/n), where omega(f, delta) denotes the modulus of continuity, and this provides a new proof of the classical Jackson`s theorem. Despite that it turns out that K(n)(theta) coincide with positive cosine polynomials generated by Fejer, our proof differs from others known in the literature.

In vitro basal cytotoxicity assay applied to estimate acute oral systemic toxicity of grandisin and its major metabolite

Preclinical investigations can start with preliminary in vitro studies before using animal models. Following this approach, the number of animals used in preclinical acute toxicity testing can be reduced. In this study, we employed an in-house validated in vitro cytotoxicity test based on the Spielmann approach for toxicity evaluation of the lignan grandisin, a candidate anticancer agent, and its major metabolite. the 4-O-demethylgrandisin, by neutral red uptake (NRU) assay, on mouse fibroblasts Balb/c 3T3 cell line. Using different concentrations of grandisin and its major metabolite (2.31; 1.16; 0.58; 0.29; 0.14; 0.07; 0.04; 0.002 mu M) in Balb/c 3T3-A31 NRU cytotoxicity assay, after incubation for 48 h, we obtained IC(50) values for grandisin and its metabolite of 0.078 and 0.043 mu M, respectively. The computed LD(50) of grandisin and 4-O-demethylgrandisin were 617.72 and 429.95 mg/kg, respectively. Both were classified under the Globally Harmonized System as category 4. Since pharmacological and toxicological data are crucial in the developmental stages of drug discovery, using an in vitro assay we demonstrated that grandisin and its metabolite exhibit distinct toxicity profiles. Furthermore, results presented in this work can contribute to reduce the number of animals required in subsequent pharmacological/toxicological studies. (C) 2010 Elsevier GmbH. All rights reserved.

A Consumer Surplus Estimate of Peace & Love festival in Borlänge : A Travel Cost Approach

This thesis uses zonal travel cost method (ZTCM) to estimate consumer surplus of Peace & Love festival in Borlänge, Sweden. The study defines counties as zones of origin of the visitors. Visiting rates from each zone are estimated based on survey data. The study is novel due to the fact that mostly TCM has been applied in the environmental and recreational sector, not for short term events, like P&L festival. The analysis shows that travel cost has a significantly negative effect on visiting rate as expected. Even though income has previously shown to be significant in similar studies, it turns out to be insignificant in this study. A point estimate for the total consumer surplus of P&L festival is 35.6 million Swedish kronor. However, this point estimate is associated with high uncertainty since a 95 % confidence interval for it is (17.9, 53.2). It is also important to note that the estimated value only represents one part of the total economic value, the other values of the festival's totaleconomic value have not been estimated in this thesis.

Likelihood estimate of treatment effects under selection bias

We consider methods for estimating causal effects of treatment in the situation where the individuals in the treatment and the control group are self selected, i.e., the selection mechanism is not randomized. In this case, simple comparison of treated and control outcomes will not generally yield valid estimates of casual effects. The propensity score method is frequently used for the evaluation of treatment effect. However, this method is based onsome strong assumptions, which are not directly testable. In this paper, we present an alternative modeling approachto draw causal inference by using share random-effect model and the computational algorithm to draw likelihood based inference with such a model. With small numerical studies and a real data analysis, we show that our approach gives not only more efficient estimates but it is also less sensitive to model misspecifications, which we consider, than the existing methods.