Using discriminant analysis to determine the breaking criterion for an ISW propagating over a ridge

This study aims to develop methods that are capable of deciding the breaking criterion for an internal solitary wave (ISW) propagating over a submarine ridge. Laboratory experiments were conducted in a wave tank to measure ISWs propagating over a submarine ridge. The results suggest that the ISW-ridge interaction can be grouped according to three degrees of magnitude based on the blockage parameter ζ and the degree of blocking B. For classification reasons, we first present an alternative decision model for evaluating the interaction of ISWs with an underwater ridge in a two-layer system. This approach is based on a multivariate statistical method and discriminant analysis. Information obtained from the eigenvalues is used to combine different ratio measures which are defined according to every single input and output. The discriminant model effectively classifies units into distinct predefined groups. An experimental simulation is conducted to demonstrate the practical implementation of the ISW-ridge interaction. The results of the method applied in this example are statistically significant, demonstrating the effectiveness of the ISW-ridge interaction classification method.     

Geographic assignment of seabirds to their origin: combining morphologic, genetic, and biogeochemical analyses.

Longline fisheries, oil spills, and offshore wind farms are some of the major threats increasing seabird mortality at sea, but the impact of these threats on specific populations has been difficult to determine so far. We tested the use of molecular markers, morphometric measures, and stable isotope (delta15N and delta13C) and trace element concentrations in the first primary feather (grown at the end of the breeding period) to assign the geographic origin of Calonectris shearwaters. Overall, we sampled birds from three taxa: 13 Mediterranean Cory's Shearwater (Calonectris diomedea diomedea) breeding sites, 10 Atlantic Cory's Shearwater (Calonectris diomedea borealis) breeding sites, and one Cape Verde Shearwater (C. edwardsii) breeding site. Assignment rates were investigated at three spatial scales: breeding colony, breeding archipelago, and taxa levels. Genetic analyses based on the mitochondrial control region (198 birds from 21 breeding colonies) correctly assigned 100% of birds to the three main taxa but failed in detecting geographic structuring at lower scales. Discriminant analyses based on trace elements composition achieved the best rate of correct assignment to colony (77.5%). Body measurements or stable isotopes mainly succeeded in assigning individuals among taxa (87.9% and 89.9%, respectively) but failed at the colony level (27.1% and 38.0%, respectively). Combining all three approaches (morphometrics, isotopes, and trace elements on 186 birds from 15 breeding colonies) substantially improved correct classifications (86.0%, 90.7%, and 100% among colonies, archipelagos, and taxa, respectively). Validations using two independent data sets and jackknife cross-validation confirmed the robustness of the combined approach in the colony assignment (62.5%, 58.8%, and 69.8% for each validation test, respectively). A preliminary application of the discriminant model based on stable isotope delta15N and delta13C values and trace elements (219 birds from 17 breeding sites) showed that 41 Cory's Shearwaters caught by western Mediterranean long-liners came mainly from breeding colonies in Menorca (48.8%), Ibiza (14.6%), and Crete (31.7%). Our findings show that combining analyses of trace elements and stable isotopes on feathers can achieve high rates of correct geographic assignment of birds in the marine environment, opening new prospects for the study of seabird mortality at sea.     

Sensitivity analysis on the ecological bias for Seoul tuberculosis data

In ecological studies, researchers often try to convey the analysis results to individual level based on aggregate data. In order to do this correctly, the possibility of ecological bias should be studied and addressed. One of the key ideas used to address the ecological bias issue is to derive the ecological model from the individual model and to check whether the parameter of interest in the individual model is identifiable in the ecological model. However, the procedure depends on unverifiable assumptions, and we recommend checking how sensitive the results are to these unverifiable assumptions. We analyzed the tuberculosis data that was collected in Seoul in 2005 using a spatial ecological regression model for the aggregate count data with spatial correlation, and found that the deprivation index is likely to have a small positive effect on the occurrence risk of tuberculosis in individual level in Seoul. We considered this finding in various aspects by performing in depth sensitivity analyses. In particular, our findings are shown to be robust to the distribution assumptions for the individual exposure and missing binary covariate across various scenarios.     

Environmental data management: The identification of outliers

Increased emphasis on pollution control and abatement has often necessitated development of large scale data bases. While sophisticated techniques have been developed and employed for data storage and manipulation, parallel developments in analyzing the accuracy and reliability of the data have been absent (see, for instance, the broad spectrum of requirements outlined in the Federal Water Pollution Control Act Amendments of 1972).This paper centers on the latter and sets out a procedure for data editing and “outlier” identification based on an application of discriminant analysis. A hypothetical example is included along with some suggested applications.     

Detecting an orientation component in animal paths when the preferred direction is individual-dependent

An orientation component leads to directionally biased paths, with major consequences in animal population redistribution. Classical orientation analyses, which focus on the overall direction of motion, are useless for detecting such a component when the preferred direction is not common to the whole population, but differs from one path to another. In-depth path analyses are required in this case. They consist of determining whether paths are more suitably represented as biased or unbiased random walks. The answer is not easy because most animals' paths show some forward persistence propensity that acts as a purely local directional bias and, hence, blurs the possible occurrence of an additional, consistent bias in a preferred direction. I highlight the key differences between biased and unbiased random walks and the different ways orientation mechanisms can generate a consistent directional bias. I then examine the strength and weakness of the available methods likely to detect it. Finally, I introduce a new procedure based on the backward evolution of the beeline distance, from the end of the path, which might correspond to a goal toward which the animal orients itself, to each of the animal's preceding locations. This new procedure proves to be very efficient, as it requires only a small sample of short paths for detecting a possible orientation component.     

Cross-validation of species distribution models: removing spatial sorting bias and calibration with a null model

Species distribution models are usually evaluated with cross-validation. In this procedure evaluation statistics are computed from model predictions for sites of presence and absence that were not used to train (fit) the model. Using data for 226 species, from six regions, and two species distribution modeling algorithms (Bioclim and MaxEnt), I show that this procedure is highly sensitive to "spatial sorting bias": the difference between the geographic distance from testing-presence to training-presence sites and the geographic distance from testing-absence (or testing-background) to training-presence sites. I propose the use of pairwise distance sampling to remove this bias, and the use of a null model that only considers the geographic distance to training sites to calibrate cross-validation results for remaining bias. Model evaluation results (AUC) were strongly inflated: the null model performed better than MaxEnt for 45% and better than Bioclim for 67% of the species. Spatial sorting bias and area under the receiver-operator curve (AUC) values increased when using partitioned presence data and random-absence data instead of independently obtained presence-absence testing data from systematic surveys. Pairwise distance sampling removed spatial sorting bias, yielding null models with an AUC close to 0.5, such that AUC was the same as null model calibrated AUC (cAUC). This adjustment strongly decreased AUC values and changed the ranking among species. Cross-validation results for different species are only comparable after removal of spatial sorting bias and/or calibration with an appropriate null model.     

Assessing the Geographic Representativeness of Genebank Collections: the Case of Bolivian Wild Potatoes

Abstract: Genebank collection databases can be used for ecogeographical studies under the assumption that the accessions are a geographically unbiased sample. We evaluated the representativeness of a collection of wild potatoes from Bolivia and defined and assessed four types of bias: species, species-area, hotspot, and infrastructure. Species bias is the sampling of some species more often than others. Species-area bias is a sampling that is disproportionate to the total area in which a species is found. Hotspot bias is the disproportionate sampling of areas with high levels of diversity. Infrastructure bias is the disproportionate sampling of areas near roads and towns. Each of these biases is present in the Bolivian wild potato collection. The infrastructure bias was strong: 60% of all wild potato accessions were collected within 2 km of a road, as opposed to 22%, if collections had been made randomly. This analysis can serve as a guide for future collecting trips. It can also provide baseline information for the application of genebank data in studies based on geographic information systems.     

Assessing desertification risk using system stability condition analysis

This paper describes a procedure for evaluating the desertification risk in threatened areas. The procedure is based on an eight-equation dynamic model of a generic human–resource system that can be applied to different desertification syndromes. For each application, interest focuses on finding all the possible long-term final states of the system and on defining the conditions that mark out sustainability and long-term desertification by means of unambiguous specific parameter relations. The procedure is applied to three typified cases in Spain: (A) rainfed crops in areas with high soil erosion risk; (B) irrigated intensive agricultural systems; and (C) commercial rangelands. Results show that, in case A, high profit scenarios are responsible for the final extension of desertification but do not determine the specific threshold between sustainability and desertification. They do, however, in cases B and C.     

干旱区绿洲土壤盐渍化程度遥感定量评价

土壤盐渍化是制约渭干河-库车河三角洲绿洲内部植被生长最主要的生态环境地质问题,也是影响区域农业生产的第一障碍性问题.利用Landsat TM遥感数据,构建了多维向量空间下的11个遥感定量指标.首先,采用因子分析方法分析不同程度的盐渍化与这11个指标之间的相关性,并建立因子得分模型.其次,采用Fisher逐步判别分析法从这些指标中筛选出8个与盐渍化程度密切相关的指标:氧化铁、归一化盐分指数、地表反照率、植被指数、TM1、TM5、地表温度、湿度指标,并对这些指标建立判别方程,然后对渭干河—库车河三角洲绿洲典型样区盐渍化程度进行评价,将所有的实地样本分别代入判别方程,通过分析,发现判别方程的精度较高,与实际情况较吻合.轻度盐渍化的判别精度为80.3%,中度盐渍化的判别精度为68.75%,重度盐渍化的判别精度为83.65%.中度盐渍化土地的判别精度与轻度盐渍化、重度盐渍化土地的判别精度相比较低,这与中度盐渍化的光谱特征与其它地物易混淆有关.总之,所选取的指标简单,易于获取,有利于盐渍化的定量分析与评价.     

Sensitivity of species-distribution models to error,bias, and model design: An application to resource selection functions for woodland caribou

Models that predict distribution are now widely used to understand the patterns and processes of plant and animal occurrence as well as to guide conservation and management of rare or threatened species. Application of these methods has led to corresponding studies evaluating the sensitivity of model performance to requisite data and other factors that may lead to imprecise or false inferences. We expand upon these works by providing a relative measure of the sensitivity of model parameters and prediction to common sources of error, bias, and variability. We used a one-at-a-time sample design and GPS location data for woodland caribou (Rangifer tarandus caribou) to assess one common species-distribution model: a resource selection function. Our measures of sensitivity included change in coefficient values, prediction success, and the area of mapped habitats following the systematic introduction of geographic error and bias in occurrence data, thematic misclassification of resource maps, and variation in model design. Results suggested that error, bias and model variation have a large impact on the direct interpretation of coefficients. Prediction success and definition of important habitats were less responsive to the perturbations we introduced to the baseline model. Model coefficients, prediction success, and area of ranked habitats were most sensitive to positional error in species locations followed by sampling bias, misclassification of resources, and variation in model design. We recommend that researchers report, and practitioners consider, levels of error and bias introduced to predictive species-distribution models. Formal sensitivity and uncertainty analyses are the most effective means for evaluating and focusing improvements on input data and considering the range of values possible from imperfect models.     

A performance indicator for ranked set sampling using ranking error probability matrix

Ranked set sampling (RSS) is a sampling procedure that has been shown to provide more efficient procedures than simple random sampling, in particular the Mann-Whitney-Wilcoxon (MWW) statistic and the empirical distribution function (EDF). We briefly review the work of Bohn (1992) and Stokes and Sager (1988) on the effect of imperfect ranking on the RSS-based MWW test and on the RSS-based EDF, respectively. We propose a model for a ranking error probability matrix which we hope will become a useful tool for evaluating RSS-based statistical procedures     

Variances and variance estimators of the improved ratio estimators under adaptive cluster sampling

Although not design-unbiased, the ratio estimator is recognized as more efficient when a certain degree of correlation exists between the variable of primary interest and the auxiliary variable. Meanwhile, the Rao–Blackwell method is another commonly used procedure to improve estimation efficiency. Various improved ratio estimators under adaptive cluster sampling (ACS) that make use of the auxiliary information together with the Rao–Blackwellized univariate estimators have been proposed in past research studies. In this article, the variances and the associated variance estimators of these improved ratio estimators are proposed for a thorough framework of statistical inference under ACS. Performance of the proposed variance estimators is evaluated in terms of the absolute relative percentage bias and the empirical mean-squared error. As expected, results show that both the absolute relative percentage bias and the empirical mean-squared error decrease as the initial sample size increases for all the variance estimators. To evaluate the confidence intervals based on these variance estimators and the finite-population Central Limit Theorem, the coverage rate and the interval width are used. These confidence intervals suffer a disadvantage similar to that of the conventional ratio estimator. Hence, alternative confidence intervals based on a certain type of adjusted variance estimators are constructed and assessed in this article.     

Site-selection bias and apparent population declines in long-term studies

Detecting population declines is a critical task for conservation biology. Logistical difficulties and the spatiotemporal variability of populations make estimation of population declines difficult. For statistical reasons, estimates of population decline may be biased when study sites are chosen based on abundance of the focal species. In this situation, apparent population declines are likely to be detected even if there is no decline. This site-selection bias is mentioned in the literature but is not well known. We used simulations and real population data to examine the effects of site-selection biases on inferences about population trends. We used a left-censoring method to detect population-size patterns consistent with site-selection bias. The site-selection bias is an important consideration for conservation biologists, and we offer suggestions for minimizing or mitigating it in study design and analysis. Article impact statement: Estimates of population declines are biased if studies begin in large populations, and time-series data show a signature of such an effect.     

Using unobtrusive sensors to measure and minimize Hawthorne effects: Evidence from cookstoves

People act differently when they know they are being observed. This phenomenon, the Hawthorne effect, can bias estimates of program impacts. Unobtrusive sensors substituting for human observation can alleviate this bias. To demonstrate this potential, we used temperature loggers to measure fuel-efficient cookstoves as a replacement for three-stone fires. We find a large Hawthorne effect: when in-person measurement begins, participants increase fuel-efficient stove use approximately three hours/day (53%) and reduce three-stone fire use by approximately two hours/day (29%). When in-person measurement ends, participants reverse those changes. We then examine how this Hawthorne effect biases estimates of fuel use and pollution concentrations. Our results reinforce concerns about Hawthorne effects, especially in policy-relevant impact evaluations. By measuring the Hawthorne effect we permit researchers to correct for the bias it introduces.     

Spatial optimization of protected area placement incorporating ecological, social and economical criteria

We describe two approaches for spatial optimization of protected area placement, both based on maximizing an objective function that incorporates ecological, social, and economical criteria. Of these, a seed cell selection procedure works by evaluating potential cells for protection one by one, picking the one that maximizes the objective function, adding seed cells. This continues to full protection of the project area. The other is a Monte Carlo approach, which uses a likelihood sampling procedure based on weighted importance layers of conservation interest to evaluate alternative protected area sizing and placement. This is similar to the objective function of Marxan, a priority-selection decision-support tool based on optimization algorithms using geographic information system data. The two approaches are alternative options in a common spatial optimization module, which uses the time- and spatial-dynamic Ecospace model for the evaluations. The optimizations are implemented as components of the Ecopath with Ecosim approach and software. In a case study, we find that there can be protected area zoning that will accommodate economical and social factors, without causing ecological deterioration. We also find a tradeoff between including cells of special conservation interest, and the economic and social interests. While this does not need to be a general feature, it emphasizes the need to use modeling techniques to evaluate the tradeoff.     

Value of time: Speeding behavior and gasoline prices

‘Value of Time’ (VOT) is a key parameter in economics and policy. This paper presents an alternative method to estimate VOT by analyzing an hourly dataset on drivers speeding behavior as a function of the gasoline price. Our identification strategy is novel as it is based on the intensive margin. In comparison, previous studies reveal VOT on the extensive margin, but choice alternatives have multiple attributes thereby potentially confounding estimates. Consistent with the range of the prior literature, we find a VOT of about 50% of the wage rate and analyze sources of bias from accidents and traffic tickets. These bias functions suggest that previous stated preference VOT estimates are likely downward whereas previous revealed preference estimates are likely upward biased.     

Sources of bias in ecological studiesof non-rare events

Ecological studies investigate relationships at the level of the group, rather than at the level of the individual. Although such studies are a common design in epidemiology, it is well-known that estimates may be subject to ecological bias. Most discussion of ecological bias has focused on rare disease events, where the tractability of the loglinear model allows some characterization of the nature of different biases. This paper concentrates on non-rare events, where the Poisson approximation to the binomial distribution is not appropriate. We limit the discussion to bias that arises from within-area variability in exposures and confounders. Our aims are to investigate the likely sizes and directions of bias and, where possible, to suggest methods for controlling the bias or for addressing the sensitivity of inference to assumptions on the nature of the bias. We illustrate that for non-rare events it is much more difficult to characterize the direction of bias than in the rare case. A series of simple numerical examples based on a chronic study of respiratory health illustrate the ideas of the paper.     

Use of opportunistic sightings and expert knowledge to predict and compare Whooping Crane stopover habitat

Predicting a species’ distribution can be helpful for evaluating management actions such as critical habitat designations under the U.S. Endangered Species Act or habitat acquisition and rehabilitation. Whooping Cranes (Grus americana) are one of the rarest birds in the world, and conservation and management of habitat is required to ensure their survival. We developed a species distribution model (SDM) that could be used to inform habitat management actions for Whooping Cranes within the state of Nebraska (U.S.A.). We collated 407 opportunistic Whooping Crane group records reported from 1988 to 2012. Most records of Whooping Cranes were contributed by the public; therefore, developing an SDM that accounted for sampling bias was essential because observations at some migration stopover locations may be under represented. An auxiliary data set, required to explore the influence of sampling bias, was derived with expert elicitation. Using our SDM, we compared an intensively managed area in the Central Platte River Valley with the Niobrara National Scenic River in northern Nebraska. Our results suggest, during the peak of migration, Whooping Crane abundance was 262.2 (90% CI 40.2?3144.2) times higher per unit area in the Central Platte River Valley relative to the Niobrara National Scenic River. Although we compared only 2 areas, our model could be used to evaluate any region within the state of Nebraska. Furthermore, our expert‐informed modeling approach could be applied to opportunistic presence‐only data when sampling bias is a concern and expert knowledge is available.