Robust formulations for clustering-based large-scale classification

Chebyshev-inequality-based convex relaxations of Chance-Constrained Programs (CCPs) are shown to be useful for learning classifiers on massive datasets. In particular, an algorithm that integrates efficient clustering procedures and CCP approaches for computing classifiers on large datasets is proposed. The key idea is to identify high density regions or clusters from individual class conditional densities and then use a CCP formulation to learn a classifier on the clusters. The CCP formulation ensures that most of the data points in a cluster are correctly classified by employing a Chebyshev-inequality-based convex relaxation. This relaxation is heavily dependent on the second-order statistics. However, this formulation and in general such relaxations that depend on the second-order moments are susceptible to moment estimation errors. One of the contributions of the paper is to propose several formulations that are robust to such errors. In particular a generic way of making such formulations robust to moment estimation errors is illustrated using two novel confidence sets. An important contribution is to show that when either of the confidence sets is employed, for the special case of a spherical normal distribution of clusters, the robust variant of the formulation can be posed as a second-order cone program. Empirical results show that the robust formulations achieve accuracies comparable to that with true moments, even when moment estimates are erroneous. Results also illustrate the benefits of employing the proposed methodology for robust classification of large-scale datasets.

Entgiftungsstoffwechsel der Nordsee-Kliesche (Limanda limanda)

In the framework of monitoring programmes organized under several sea protection conventions (HELSINKI Conv., OSPAR Conv.) the contracting parties are requested to develop appropriate techniques for Biological-Effect- Monitoring. In following these recommendations the Institut for Fisheries Ecology studies the 7-ethoxyresorufin- O-deethylase (EROD) activity in the liver of dab. EROD represents one enzyme of the cytochrome P-450 species, also called mixed function oxygenases (MFO), which is induced by certain organic contaminants, e.g. PCBs. On the other hand, an influence of natural factors like season, temperature or spawning on the EROD activity may be possible. The present study represents an insight into the status of the EROD activity in North Sea dab. Ultimately, we intend to decide if EROD activity is an appropriate tool to detect effects of contaminants. The EROD activity in the liver of 687 dabs, caught in the North Sea at different seasons in 1995 and 1996 with the fishery research vessel “Walther Herwig III”, has been determined and the data obtained have been statistically evaluated. The logarithmically transformed values of the EROD activity are following approximately a normal distribution. Due to the wide variation of the enzyme activities and due to the small number of samples minor differences between samples are not detectable. Nevertheless, comparing the enzyme activities at different sites of the North Sea, some significant differences have been identified. A model for the discription of seasonal variations of EROD activity, developed at the Biologische Anstalt Helgoland, could be helpful for interpretation.

An evaluation of the survival of Atlantic salmon (Salmo salar L.)fry stocked in eight streams in the North West of England.

Eight streams from the North West of England were stocked with Atlantic salmon (Salmo salar L.) fed fry at densities ranging from 1 to 4/m2 over a period of up to three years to evaluate survival to the end of the first an d second growing periods and hence assess the value of stocking as a management practice. Survival to the end of the first growin g period (mean duration of 108 days) was found to vary between 7.8 and 41.3% with a mean of 22% and CV of 0.44. Survival from the end of the first growing period to the end of the second growing period (mean duration of 384 days) ranged from 19.9 to 34.1% with a mean of 26.3% and CV of 0.21. Survival was found to be positively related to 0+ trout density (P < 0.05) and negatively related to altitude (P < 0.05). A comparison of the raw survival data (non standardised with respect to duration of experiments) with that from other studies in relation to stocking densities revealed a negative relationship between fry survival and stocking density (P < 0.05). Densities in excess of 5/m2 tended to result in lower levels of survival. Post stocking fry dispersal patterns were examined for the 1991 data. On average 86.7% of the number of fry surviving remained within the stocked zone by the end of the first growing period. With the exception of one stream there was little in the way of dispersal beyond the stocked zone. The dispersal pattern approximated to the normal distribution (P < 0.05). It was estimated that stocking can result in a net gain of fish to a river system compared with natural productivity, however the numerical significance of this gain and its cost effectiveness need to be determined on a river specific basis.

Length-weight relationship of mudskippers (Gobiidae: Oxudercinae) in the coastal areas of Selangor, Malaysia

Parameters a and b of the length-weight relationship (LWR) were estimated for eleven species of mudskippers caught in the coastal areas of Selangor, Malaysia. The values of b ranged from 2.56 to 3.50 with the mean b equal to 2.95 (n=11; sd=0.302). A normal distribution of the calculated LWR exponent (b) was obtained.

Length-weight relationship of marine fish species off Sao Sebastiao system, Sao Paulo, Southeastern Brazil

The parameters a and b of the length-weight relationship of the form W = a . L super(b) were estimated for 57 fish species sampled in Sao Sebastiao Channel and shelf system in 1997, Sao Paulo, Brazil. The b values ranged from 2.746 to 3.617. The Student's t-test revealed that mot (44) species had b values significantly different from 3. A normal distribution of the calculated LWR exponents (b) was obtained.

Reverberações da Lei n 8.213/91 à luz das concepções de deficiência: estudo de caso sobre o programa de sensibilização gerencial em uma organização privada de ensino profissionalizante

Com base nos rebatimentos da Lei n 8.213/1991, que prevê a obrigatoriedade legal de empresas privadas brasileiras contratarem de 2% a 5% de beneficiários da Previdência Social reabilitados ou de pessoas com deficiência (PcDs) habilitadas em seus quadros funcionais, o objetivo da presente pesquisa foi avaliar as repercussões da Lei de Cotas, tendo como referência as concepções de deficiência preponderantes dentro de uma organização privada de ensino profissionalizante por meio do seu Projeto de Sensibilização Gerencial. O projeto de sensibilização gerencial objetivou consolidar as etapas iniciais, na implementação de um programa institucional de valorização da diversidade. A hipótese principal formulada para a pesquisa, dentre outras relativas às concepções de deficiência, por parte de gestores, é que as reverberações positivas produzidas pela Lei de Cotas atingem, inclusive, as organizações empresariais supostamente distanciadas do movimento inclusionista. Para fins de avaliação das concepções de deficiência foram utilizados delineamentos de pesquisa estatística, compreendendo o teste de Shapiro-Wilk para determinar se as respostas às perguntas configuravam ou não uma distribuição normal, além do Coeficiente de Correlação de Pearson para medir a intensidade da relação linear entre as variáveis estudadas. Serviram como instrumentos da pesquisa um questionário sociodemográfico e um inventário de concepções de deficiência (ICD), sendo este direcionado para o objeto atitudinal, considerando-se as percepções sociais favoráveis e desfavoráveis no grupo pesquisado. Esse inventário de concepções de deficiência é composto de sete blocos de asserções e de uma escala do tipo Likert de seis pontos, que foi aplicada em um universo que contempla estrategicamente 60 participantes selecionados em três (3) grupos (Grupo piloto 1 envolvendo 30 participantes das áreas do Comitê Gestor do projeto na empresa; Grupo piloto 2 envolvendo 12 trainees; e Grupo Gerencial envolvendo 18 participantes, incluindo Gerentes de Área e Gerentes de Equipe da Superintendência de Produtos Educacionais). De posse dos resultados da avaliação das concepções de deficiência pelo ICD, foram realizados workshops de sensibilização com os participantes com o objetivo de sensibilizá-los e disseminar o conhecimento sobre inclusão social e laboral de PcDs, as ações de políticas públicas na atualidade, a natureza das deficiências, considerando-se os aspectos sociais da profissionalização, empregabilidade de pessoas com deficiência na empresa. Para tanto, foram empregadas técnicas e procedimentos lúdicos, além de debates para fins de reflexão crítica por parte dos participantes. A avaliação de reação foi conduzida ao término desses workshops. O conjunto dos dados levantados até então possibilitou proceder-se a um diagnóstico das concepções que prevalecem sobre PcDs na organização alvo da pesquisa. Os resultados evidenciaram a coexistência de concepções distintas da deficiência, indicando que, embora concepções negativas se perpetuem, as reverberações da Lei de Cotas têm apresentado também repercussões visivelmente positivas valorizando, assim, as ações corporativas apontadas para a diversidade humana, no contexto do trabalho. Em termos conclusivos, considera-se, no entanto, que o processo de inclusão laboral deva ser percebido por parte dos gestores como contínuo e em direção à mudança do comportamento humano nas organizações frente à profissionalização de PcDs. Posteriormente, tendo-se como suporte os resultados da presente pesquisa, um plano de ação institucional será implementado, como proposta de um programa balizado em 10 projetos sintéticos que servirão de modelo para empresas brasileiras interessadas em incluir a diversidade e reter talentos com deficiência em seus postos de trabalho, de modo a garantir-lhes o direito de exercício pleno da cidadania.

Qualidade do debridamento obtido pela técnica de instrumentação da lima única de NiTi em canais ovais e circulares

O presente avalia a qualidade da remoção de tecido pulpar após o preparo químico-cirúrgico realizado com a técnica da lima única, descrita por Ghassan Yared. Ainda não há publicado pesquisa sobre os resultados desta técnica. Este estudo compara o percentual de tecido pulpar remanescente em canais radiculares ovais e circulares de incisivos inferiores recém-extraídos que possuíssem polpa viva e armazenados em formol a 10%. Foram comparadas duas técnicas: ProTaper Universal e a técnica da lima única F2. Após uma rigorosa seleção, quarenta e oito dentes com polpa viva que possuíam indicação de extração, foram preparados, classificados em canais ovais e circulares, separados aleatoriamente em 4 grupos e instrumentados com as duas técnicas. O grupo controle, com 12 espécimes, não recebeu nenhum tipo de intervenção. G1 (n=12), canais ovais, instrumentados com a técnica ProTaper Universal; G2 (n=12), canais ovais instrumentados com a técnica da lima única F2; G3 (n=12), canais circulares instrumentados com a técnica ProTaper Universal; G4 (n=12), canais circulares instrumentados com a técnica da lima única F2. Então, seções transversais foram preparadas para avaliação histológica. A análise da quantidade de tecido pulpar remanescente foi avaliada digitalmente. A análise preliminar dos dados brutos em conjunto de todos os grupos experimentais revelou um padrão de distribuição normal por meio do uso do teste Kolmogorov-Smirnov. A análise foi realizada, e os dados brutos foram avaliados através de métodos não-paramétricos: Teste H Kruskal-Wallis. O valor percentual mínimo de tecido remanescente foi de 0% e o máximo de 37,78% entre todos os grupos. Os valores relativos a quantidade de tecido pulpar remanescente variaram entre 0 a 43.47% m2. Os resultados do Teste H Kruskal-Wallis não revelaram diferenças entre as seções mais apicais (p > 0.05). Entretanto, foi encontrada diferença significante entre as seções mais apicais e a seção do terço médio (p < 0.05). Também foram encontradas diferenças significantes quando canais circulares foram comparados com canais ovais independente da técnica de instrumentação utilizada (p < 0.05). Porém, entre as duas técnicas de instrumentação estudadas, tanto nos canais ovais quanto para os os canais circulares, não houve diferença estatística significante (p > 0.05). A proposta deste estudo é a de fazer uma reflexão sobre a real necessidade de um grande número de instrumentos para o total preparo de canais radiculares, uma vez que nenhuma das técnicas foi capaz de debridar por completo o espaço do canal radicular.

A stochastic model of temporal variations in monthly temperature, precipitation, snowfall, and resulting snowpack

We report a Monte Carlo representation of the long-term inter-annual variability of monthly snowfall on a detailed (1 km) grid of points throughout the southwest. An extension of the local climate model of the southwestern United States (Stamm and Craig 1992) provides spatially based estimates of mean and variance of monthly temperature and precipitation. The mean is the expected value from a canonical regression using independent variables that represent controls on climate in this area, including orography. Variance is computed as the standard error of the prediction and provides site-specific measures of (1) natural sources of variation and (2) errors due to limitations of the data and poor distribution of climate stations. Simulation of monthly temperature and precipitation over a sequence of years is achieved by drawing from a bivariate normal distribution. The conditional expectation of precipitation. given temperature in each month, is the basis of a numerical integration of the normal probability distribution of log precipitation below a threshold temperature (3°C) to determine snowfall as a percent of total precipitation. Snowfall predictions are tested at stations for which long-term records are available. At Donner Memorial State Park (elevation 1811 meters) a 34-year simulation - matching the length of instrumental record - is within 15 percent of observed for mean annual snowfall. We also compute resulting snowpack using a variation of the model of Martinec et al. (1983). This allows additional tests by examining spatial patterns of predicted snowfall and snowpack and their hydrologic implications.

A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning.

The statistical properties of recreational catch rate data for some fish stocks off the northeast U.S. coast

Recreational fisheries in the waters off the northeast U.S. target a variety of pelagic and demersal fish species, and catch and effort data sampled from recreational fisheries are a critical component of the information used in resource evaluation and management. Standardized indices of stock abundance developed from recreational fishery catch rates are routinely used in stock assessments. The statistical properties of both simulated and empirical recreational fishery catch-rate data such as those collected by the National Marine Fisheries Service (NMFS) Marine Recreational Fishery Statistics Survey (MRFSS) are examined, and the potential effects of different assumptions about the error structure of the catch-rate frequency distributions in computing indices of stock abundance are evaluated. Recreational fishery catch distributions sampled by the MRFSS are highly contagious and overdispersed in relation to the normal distribution and are generally best characterized by the Poisson or negative binomial distributions. The modeling of both the simulated and empirical MRFSS catch rates indicates that one may draw erroneous conclusions about stock trends by assuming the wrong error distribution in procedures used to developed standardized indices of stock abundance. The results demonstrate the importance of considering not only the overall model fit and significance of classification effects, but also the possible effects of model misspecification, when determining the most appropriate model construction.

Seasonal growth of King George whiting (Sillaginodes punctata) estimated from length-at-age samples of the legal-size harvest

Samples of 11,000 King George whiting (Sillaginodes punctata) from the South Australian commercial and recreational catch, supplemented by research samples, were aged from otoliths. Samples were analyzed from three coastal regions and by sex. Most sampling was undertaken at fish processing plants, from which only fish longer than the legal minimum length were obtained. A left-truncated normal distribution of lengths at monthly age was therefore employed as model likelihood. Mean length-at-monthly-age was described by a generalized von Bertalanffy formula with sinusoidal seasonality. Likelihood standard deviation was modeled to vary allometrically with mean length. A range of related formulas (with 6 to 8 parameters) for seasonal mean length at age were compared. In addition to likelihood ratio tests of relative fit, model selection criteria were a minimum occurrence of high uncertainties (>20% SE), of high correlations (>0.9, >0.95, and >0.99) and of parameter estimates at their biological limits, and we sought a model with a minimum number of parameters. A generalized von Bertalanffy formula with t0 fixed at 0 was chosen. The truncated likelihood alleviated the overestimation bias of mean length at age that would otherwise accrue from catch samples being restricted to legal sizes.

东灵山地区辽东栎群中植物种群的时空格局研究

在北京东灵山地区以辽东栎（Quercus liaotungensis Koidz.）林群落为主要对象，研究了辽东栎种群和主要伴生种的空间格局、种间联结以及辽东栎种子雨的时间格局，井对空间格局的表达方法进行尝试性的探讨。 对辽东栎林群落中主要乔木树种辽东栎、大叶白蜡（Fraxinus thynchophylla）和五角枫（Acer mono）种群空间格局与动态的考察结果显示：各种群大多表现出聚集分布的趋势;辽东栎种群在各个龄级层均显示出较强优势，但种群内竞争激烈，龄级结构相对稳定;人叶白蜡的幼苗对环境适应能力较强，在成熟林中呈现出较五角枫种群更明显的增长型趋势：幼林中3个种群均呈衰退型趋势。 种间联结性的测定结果较好地反映了植物对环境的适应性，用关联分析法作出的树状图利用临界种将样方划分为若干组，体现出小生境之间的差异，反映了植物对其生境亦具有一定指示作用。 辽东栎种子雨的时间格局研究结果表明：种群种子雨密度随时间的变化呈正态分布趋势，种群水平利个体水平的种子雨累积随时间变化均符台Logistic曲线，具有很高的相关性。 空间格局表达方法的探讨显示：可以用集合论的方法对数学模型进行形象表达。

Predicting surface vibration from underground railways through inhomogeneous soil

Noise and vibration from underground railways is a major source of disturbance to inhabitants near subways. To help designers meet noise and vibration limits, numerical models are used to understand vibration propagation from these underground railways. However, the models commonly assume the ground is homogeneous and neglect to include local variability in the soil properties. Such simplifying assumptions add a level of uncertainty to the predictions which is not well understood. The goal of the current paper is to quantify the effect of soil inhomogeneity on surface vibration. The thin-layer method (TLM) is suggested as an efficient and accurate means of simulating vibration from underground railways in arbitrarily layered half-spaces. Stochastic variability of the soils elastic modulus is introduced using a KL expansion; the modulus is assumed to have a log-normal distribution and a modified exponential covariance kernel. The effect of horizontal soil variability is investigated by comparing the stochastic results for soils varied only in the vertical direction to soils with 2D variability. Results suggest that local soil inhomogeneity can significantly affect surface velocity predictions; 90 percent confidence intervals showing 8 dB averages and peak values up to 12 dB are computed. This is a significant source of uncertainty and should be considered when using predictions from models assuming homogeneous soil properties. Furthermore, the effect of horizontal variability of the elastic modulus on the confidence interval appears to be negligible. This suggests that only vertical variation needs to be taken into account when modelling ground vibration from underground railways. © 2012 Elsevier Ltd. All rights reserved.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.