Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completely absent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involved parts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method is introduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that the theoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approach has reasonable properties from a compositional point of view. In particular, it is “natural” in the sense that it recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in the same paper a substitution method for missing values on compositional data sets is introduced

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades <-1 to >9) were separated, and each fraction was analysed for its chemical composition. The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø <8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar). To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend. Key words: sediment, geochemistry, grain size, regression, step function

Speculative parallelization of multipath radiosity algorithm

In computer graphics, global illumination algorithms take into account not only the light that comes directly from the sources, but also the light interreflections. This kind of algorithms produce very realistic images, but at a high computational cost, especially when dealing with complex environments. Parallel computation has been successfully applied to such algorithms in order to make it possible to compute highly-realistic images in a reasonable time. We introduce here a speculation-based parallel solution for a global illumination algorithm in the context of radiosity, in which we have taken advantage of the hierarchical nature of such an algorithm

An Interval Intelligent-based Approach for Fault Detection and Modelling

Not considered in the analytical model of the plant, uncertainties always dramatically decrease the performance of the fault detection task in the practice. To cope better with this prevalent problem, in this paper we develop a methodology using Modal Interval Analysis which takes into account those uncertainties in the plant model. A fault detection method is developed based on this model which is quite robust to uncertainty and results in no false alarm. As soon as a fault is detected, an ANFIS model is trained in online to capture the major behavior of the occurred fault which can be used for fault accommodation. The simulation results understandably demonstrate the capability of the proposed method for accomplishing both tasks appropriately

Prediction of the bulking phenomenon in wastewater treatment plants

The control and prediction of wastewater treatment plants poses an important goal: to avoid breaking the environmental balance by always keeping the system in stable operating conditions. It is known that qualitative information — coming from microscopic examinations and subjective remarks — has a deep influence on the activated sludge process. In particular, on the total amount of effluent suspended solids, one of the measures of overall plant performance. The search for an input–output model of this variable and the prediction of sudden increases (bulking episodes) is thus a central concern to ensure the fulfillment of current discharge limitations. Unfortunately, the strong interrelation between variables, their heterogeneity and the very high amount of missing information makes the use of traditional techniques difficult, or even impossible. Through the combined use of several methods — rough set theory and artificial neural networks, mainly — reasonable prediction models are found, which also serve to show the different importance of variables and provide insight into the process dynamics

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Comparing meta-analysis and ecological-longitudinal analysis in time-series studies: a case study of the effects of air pollution on mortality in three Spanish cities

The objective of this paper is to introduce a diVerent approach, called the ecological-longitudinal, to carrying out pooled analysis in time series ecological studies. Because it gives a larger number of data points and, hence, increases the statistical power of the analysis, this approach, unlike conventional ones, allows the complementation of aspects such as accommodation of random effect models, of lags, of interaction between pollutants and between pollutants and meteorological variables, that are hardly implemented in conventional approaches. Design—The approach is illustrated by providing quantitative estimates of the short-termeVects of air pollution on mortality in three Spanish cities, Barcelona,Valencia and Vigo, for the period 1992–1994. Because the dependent variable was a count, a Poisson generalised linear model was first specified. Several modelling issues are worth mentioning. Firstly, because the relations between mortality and explanatory variables were nonlinear, cubic splines were used for covariate control, leading to a generalised additive model, GAM. Secondly, the effects of the predictors on the response were allowed to occur with some lag. Thirdly, the residual autocorrelation, because of imperfect control, was controlled for by means of an autoregressive Poisson GAM. Finally, the longitudinal design demanded the consideration of the existence of individual heterogeneity, requiring the consideration of mixed models. Main results—The estimates of the relative risks obtained from the individual analyses varied across cities, particularly those associated with sulphur dioxide. The highest relative risks corresponded to black smoke in Valencia. These estimates were higher than those obtained from the ecological-longitudinal analysis. Relative risks estimated from this latter analysis were practically identical across cities, 1.00638 (95% confidence intervals 1.0002, 1.0011) for a black smoke increase of 10 μg/m3 and 1.00415 (95% CI 1.0001, 1.0007) for a increase of 10 μg/m3 of sulphur dioxide. Because the statistical power is higher than in the individual analysis more interactions were statistically significant,especially those among air pollutants and meteorological variables. Conclusions—Air pollutant levels were related to mortality in the three cities of the study, Barcelona, Valencia and Vigo. These results were consistent with similar studies in other cities, with other multicentric studies and coherent with both, previous individual, for each city, and multicentric studies for all three cities

Desenvolupament d'un sistema de previsió de risc d'infecció per Stemphylium vesicarium. Avaluació, validació i implementació en parcel·les experimentals en camps comercials de perera

L'estemfiliosi de la perera és una malaltia fúngica ocasionada pel fong Deuteromicet Stemphylium vesicarium. Els símptomes característics d'aquesta malaltia són lesions necròtiques en fulles i fruits, provocant el depreciament econòmic d'aquests últims. Si les condicions climàtiques són favorables aquesta malaltia pot arribar a provocar la pèrdua total de la producció. La malaltia es localitza bàsicament en la conca mediterrània, sent les zones fructícoles de Girona a Catalunya i Emilia-Romagna a Itàlia els llocs on les pèrdues de producció són més grans. Actualment l'únic mètode de control és la cobertura permanent dels òrgans de la planta amb fungicides. Per obtenir aquesta protecció és necessari realitzar aplicacions preventives amb una cadència fixa, setmanal o quinzenal, de diferents fungicides en els arbres. Aquesta cadència d'aplicació comporta que al llarg del període vegetatiu de la perera es realitzin entre 20 i 25 aplicacions de fungicida, per cadències setmanals, o entre 10 i 13 aplicacions si la cadència és quinzenal. Aquest elevat nombre de tractaments comporta uns efectes toxicològics a nivell mediambiental i per a la salut humana, un augment del cost econòmic de la producció i en alguns casos afavorir l'aparició de resistències del patogen a determinats fungicides. Molts d' aquests tractaments es realitzen sense que les condicions ambientals siguin favorables a l'inici d'infeccions del fong, i són per tant innecessàries. Amb l'objectiu d'obtenir una eina per determinar el risc d'infeccions en funció de les condicions climàtiques i realitzar d'aquesta manera els tractaments amb fungicides només quan són necessaris, es va elaborar el Sistema pel Tractament Racionalitzat de l'Estemfiliosi de la Perera (STREP). El model consisteix en relacionar quantitativament amb la severitat de la malaltia els paràmetres climàtics de durada del període d'humectació i temperatura mitjana durant aquests períodes. En aquest treball s'ha determinat l'eficàcia en el control de la malaltia per diferents fungicides aplicats de manera preventiva amb una cadència fixa en diferents finques comercials i durant varis anys. S'ha constatat que no hi ha cap producte totalment eficaç en el control de la malaltia, sent l'eficàcia obtinguda mitjana. Aquesta eficàcia de control esta en funció del nivell de pressió de la malaltia, obtenint eficàcies inferiors quan el nivell en els testimonis no tractats és elevat. En assajos en condicions controlades s'ha determinat que els fungicides analitzats aplicats de manera curativa o postinfecciosa, quan ja s'han iniciat les infeccions, presenten una eficàcia més baixa que quan s'apliquen preventivament. El tiram és el producte amb una eficàcia més elevada en la major part de les condicions experimentals. S' ha analitzat durant 5 anys, 1992-1996, en diferents plantacions de perera de Girona, la dinàmica dels paràmetres ambientals de durada deis períodes amb humectació, la temperatura mitjana de l'aire, la temperatura mitjana durant els períodes d'humectació, la durada dels períodes amb humitat relativa superior al 90% i la pluviometria. S 'ha constatat que la humectació deguda a la rosada és la principal causa d'humectació i que és un procés lent i continu. Les durades més freqüents d'humectació han estat períodes inferiors a 12 hores i les temperatures més freqüents durant aquests períodes han estat entre 11 i 21°C. S'ha determinat també que existeixen períodes amb humitat relativa elevada però sense formació de pel·lícula d'aigua i s'ha observat que ocasionalment existeixen períodes breus d'interrupció de la humectació en que poden existir humanitats relatives elevades o baixes. La capacitat predictiva del model STREP ha estat avaluada mitjançant 42 assajos en condicions naturals i controlades, determinant que en més del 95 % dels casos avaluats les prediccions de malaltia s'han correspost amb les observades, i que només en dos casos el model ha infraestimat el nivell de risc d'infecció. Demostrant així que el model STREP és fiable en les seves prediccions. S'ha estudiat l'efecte de la humitat relativa durant els períodes interromputs d'humectació i les durades d'aquests períodes, en les infeccions produïdes pel fong. S'ha constatat que el model STREP hauria d'incorporar com a períodes separats els períodes d'humectació interromputs amb relatives elevades (>98%) s'haurien de considerar com a períodes continus. Per una altra banda s'ha comprovat que els períodes amb humitats relatives elevades sense formació d'humectació en condicions naturals no són suficients per provocar l'inici d'infeccions de S. Vesicarium. S'ha utilitzat el model STREP com a eina per guiar els tractaments fungicides en diferents parcel·les de finques comercials durant dos anys. D'aquests assajos es pot concloure que utilitzant valors de risc acumulat SA:0.4 com a llindar per a realitzar els tractaments es controla la malaltia amb una eficàcia semblant als tractaments setmanals, però amb una reducció del nombre d'aplicacions necessàries entre un 25 i un 50%. Utilitzar un llindar SA:0.6 en finques amb varietats sensibles i una pressió d'inòcul del patogen elevada presenta bons resultats però inferiors als obtinguts amb tractaments setmanals, no obstant l'estalvi de tractaments ha estat molt important (50-70%). De tots aquests resultats es conclou que el model STREP és adequat per ésser utilitzat en programes pilot a escala comercial en el control racional de l'estemfiliosi de la perera. En funció d'aquests resultats es considerarà la seva implementació en Estacions d'Avisos Fitosanitaris.

Innovació i desestacionalització de la demanda hotelera en un districte turístic: el paper de les xarxes relacionals

La tesi es planteja en quina mesura l'estructura de la xarxa social existent en un districte turístic entre els diferents actors que la conformen, i la posició que aquests ocupen en la mateixa, afecta a la estacionalitat de la demanda hotelera del districte. S'assumeix que els fluxos d'informació i coneixement entre els actors del districte, i també entre aquests i d'altres actors externs determinaran la seva capacitat d'innovació, i per tant, la seva capacitat de creació de nous productes i ofertes turístiques capaces d'atraure turistes fora de la temporada alta. Així, es pretén analitzar si els hotels millor posicionats dins de l'estructura de relacions d'un districte turístic d'àmbit comarcal es corresponen amb els que presenten una millor estacionalitat en termes d'allargament de la temporada, de mesos d'obertura i de millora de l'ocupació en els dos darrers períodes. Addicionalment, s'analitza en quina mesura els segments de mercat turístic als que els hotels dirigeixen els seus esforços expliquen les diferencies d'estacionalitat entre ells.

La presión laboral tendenciosa (Mobbing)

El trabajo de investigación surge en el año 2001, ante la necesidad de hacer frente a una nueva realidad jurídica, el mobbing. Para ello fue decisivo el estudio de lo publicado (básicamente de ramas ajenas al Derecho) pero sobre todo las entrevistas con las víctimas de mobbing y sus asociaciones; este extremo, unido a la ausencia de un tratamiento internacional, obligó a un camino autodidacta para definir mobbing jurídicamente. La Tesis define mobbing como la presión laboral tendente a la autoeliminación de un trabajador mediante su denigración (presión laboral tendenciosa), y con ello por primera vez se tiene una definición de mobbing en línea y media, con plena validez jurídica, que es susceptible de ser memorizada y por lo tanto divulgada, para corregir el problema. El denominado "concepto uniformado de mobbing" recalca la denigración como mecanismo frente a los tratos degradantes y recalca la autoeliminación como finalidad de un comportamiento doloso. El trabajo aporta fórmulas para deslindar casos de mobbing de otras figuras próximas, y en este sentido debe citarse "la regla del 9" para saber si hay mobbing; en sede de estadísticas se critican metodológicamente muchas de ellas presentadas hasta el momento y se aporta alguna en sede de Tribunales; pero sobre todo se advierte de los riesgos jurídicos de una previsible regulación específica antimobbing, mediante el examen de las distintas definiciones que se han esgrimido hasta el momento. La segunda parte de la Tesis profundiza sobre el grado de sensibilización de nuestro ordenamiento jurídico y Tribunales, a cuyo fin se ha trabajado con más de un centernar y medio de sentencias dictadas sobre la materia, y por supuesto la totalidad de las recogidas en las bases de datos de las editoriales. El análisis sirve para apreciar la bondad de la sistemática aquí defendida, poniendo en evidencia errores, y contradicciones. La Tesis advierte que la presión laboral tendenciosa más allá de vulnerar el derecho constitucional al trabajo, o los derechos fundamentales a la integridad moral y el honor, es una transgresión a todo un "espíritu constitucional", y en este sentido se analiza con detalle tanto la posibilidad de recurrir en amparo, como el derecho a la indemnidad para quien se enfrenta a esta situación. Advirtiendo de las ventajas de efectuar esta reacción mediante la modalidad procesal de tutela de los derechos fundamentales, se analiza la recurrida acción del art.50 ET, donde se realizan aportaciones sugerentes como el plazo prescripción o la "doctrina de los antecedentes", y se otorgan respuestas a las preguntas sobre obligación de seguir trabajando y ejecución provisional. En sede de acciones de Seguridad Social, la Tesis distingue entre la incapacidad temporal y permanente (depresiones) y la muerte y supervivencia, aportándose sobre la primera la técnica denominada "interpretación en tres niveles" y descartando la posibilidad de considerar accidente de trabajo el suicidio tras un mobbing por imperativo legal, pero aportando un sucedáneo bastante razonable como es el accidente no laboral. Junto a ello se razona por la viabilidad del recargo del art.123 LGSS. Civilmente, la Tesis se posiciona de "lege ferenda" por reconducir este tipo de acciones resarcitorias del daño psíquico y moral al orden civil, por una mayor explicación sobre el origen del quantum, pero sobre todo considera inadmisible la STS 11-3-04, y ello por una pluralidad de argumentos, pero sobre todo por cuanto viene a autorizar "de facto" este tipo de conductas. La posibilidad de accionar administrativamente frente a este riesgo psicosocial se analiza en un doble terreno, la empresa y la Administración. Si bien el cauce sobre el primero tiene algunos meandros que se desbelan, la situación es radicalmente frustrante en la Administración -donde se encuentra el mayor caldo de cultivo del mobbing- , y ello por el RD 707/2002, pero todavía en mayor medida por el Criterio Técnico 34/2003 mediante el cual la interpretación del Director General de la Inspección de Trabajo y Seguridad Social ha venido tácitamente a derogar parcialmente la Ley de Prevención de Riesgos Laborales para la Administración. En materia penal, la Tesis se decanta "a priori" por dos tipos penales, los delitos contra los derechos de los trabajadores, y el delito de trato degradante; sin embargo, en la práctica sólo este segundo es el camino que puede alcanzar buen puerto. Finalmente se realiza un estudio detallado de la Ley 62/2003, ley que se divulgó como reguladora del acoso moral, y que después se defiende como un avance frente al mobbing. La Tesis advierte que no es cierto ni lo uno, ni lo otro, habiendo creado un "espejismo legal" que puede perjudicar a las víctimas de mobbing, además de no servir su estructura para una futura regulación explícita antimobbing.

Quantified real constraint solving using modal intervals with applications to control

Les restriccions reals quantificades (QRC) formen un formalisme matemàtic utilitzat per modelar un gran nombre de problemes físics dins els quals intervenen sistemes d'equacions no-lineals sobre variables reals, algunes de les quals podent ésser quantificades. Els QRCs apareixen en nombrosos contextos, com l'Enginyeria de Control o la Biologia. La resolució de QRCs és un domini de recerca molt actiu dins el qual es proposen dos enfocaments diferents: l'eliminació simbòlica de quantificadors i els mètodes aproximatius. Tot i això, la resolució de problemes de grans dimensions i del cas general, resten encara problemes oberts. Aquesta tesi proposa una nova metodologia aproximativa basada en l'Anàlisi Intervalar Modal, una teoria matemàtica que permet resoldre problemes en els quals intervenen quantificadors lògics sobre variables reals. Finalment, dues aplicacions a l'Enginyeria de Control són presentades. La primera fa referència al problema de detecció de fallades i la segona consisteix en un controlador per a un vaixell a vela.