On the consequences of misspecifing assumptions concerning residuals distribution in a repeated measures and nonlinear mixed modelling context

In this paper we describe the results of a simulation study performed to elucidate the robustness of the Lindstrom and Bates (1990) approximation method under non-normality of the residuals, under different situations. Concerning the fixed effects, the observed coverage probabilities and the true bias and mean square error values, show that some aspects of this inferential approach are not completely reliable. When the true distribution of the residuals is asymmetrical, the true coverage is markedly lower than the nominal one. The best results are obtained for the skew normal distribution, and not for the normal distribution. On the other hand, the results are partially reversed concerning the random effects. Soybean genotypes data are used to illustrate the methods and to motivate the simulation scenarios

Cation distribution and magnetization of BaFe12-2xCoxSnxO19 (x=0.9,1.28) single crystals

The distribution of Sn4+ cations within the five crystallographic sites of the magnetoplumbite (M) ‐like compound BaFe12−2xCoxSnxO19 has been analyzed using single‐crystal x‐ray‐diffraction data. The species Fe3+ and Co2+ cannot be distinguished using x rays because of their very similar atomic numbers; however, the calculation of the apparent valencies for the different sites allows an insight into the Co2+ cation segregation. The use of previous data from neutron powder diffraction allows a precise picture of the cation distribution, which indicates a pronounced site selectivity for both Sn4+ and Co2+ cations. The Sn4+ cations prefer the 4f2 sites and to a much lower extent the 12k sites, while they do not enter the octahedral 2a sites at all. Co2+ cations are distributed among tetrahedral and octahedral sites displaying a clear preference for the tetrahedral 4f1 sites. Magnetic measurements indicate that the compound still exhibits uniaxial anisotropy with the easy direction parallel to the c axis. Nevertheless, the magnetic structure shows a considerable degree of noncolinearity. A strong reduction of the magnetic anisotropy regarding that of the undoped compound is also detected.

Particle orientation distribution in Y-Fe2O3magnetic tapes by Mössbauer and hysteresis loop measurements

The particle orientation in several Y-Fe2O3magnetic tapes has been quantitatively evaluated by using the data of both Mössbauer and hysteresis loop measurements performed in the three orthogonal directions. A texture function has been obtained as a development of real harmonics. The profile of the texture function gives the quality of the different magnetic tapes. A different degree of particle orientation at the surface of the tape is evidenced by means of conversion electron Mössbauer spectra.

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

Testing the bivariate distribution of daily equity returns using copulas: an application to the Spanish stock market

In this paper we deal with the identification of dependencies between time series of equity returns. Marginal distribution functions are assumed to be known, and a bivariate chi-square test of fit is applied in a fully parametric copula approach. Several families of copulas are fitted and compared with Spanish stock market data. The results show that the t-copula generally outperforms other dependence structures, and highlight the difficulty in adjusting a significant number of bivariate data series

Spatial distribution and biomass of aquatic rooted macrophytes and their relevance in the metabolism of a Mediterranean coastal lagoon

This work aims to characterise the current autotrophic compartment of the Albufera des Grau coastal lagoon (Menorca, Balearic Islands) and to assess the relationship between the submerged macrophytes and the limnological parameters of the lagoon. During the study period the submerged vegetation was dominated by the macrophyte Ruppia cirrhosa, which formed dense extensive meadows covering 79% of the surface. Another macrophyte species, Potamogeton pectinatus, was also observed but only forming small stands near the rushing streams. Macroalgae were only occasionally observed. Macrophyte biomass showed a clear seasonal trend, with maximum values in July. The biomass of R. cirrhosa achieved 1760 g DW m-2, the highest biomass ever reported for this species in the literature. The seasonal production-decomposition cycle of the macrophyte meadows appears to drive the nutrient dynamics and carbon fluxes in the lagoon. Despite the significant biomass accumulation and the absence of a washout of nutrients and organic matter to the sea, the lagoon did not experience a dystrophic collapse. These results indicate that internal metabolism is more important than exchange processes in the lagoon.

Polydispersity index from linear viscoelastic data: Unimodal and bimodal linear polymer melts

This article describes a method for determining the polydispersity index Ip2=Mz/Mw of the molecular weight distribution (MWD) of linear polymeric materials from linear viscoelastic data. The method uses the Mellin transform of the relaxation modulus of a simple molecular rheological model. One of the main features of this technique is that it enables interesting MWD information to be obtained directly from dynamic shear experiments. It is not necessary to achieve the relaxation spectrum, so the ill-posed problem is avoided. Furthermore, a determinate shape of the continuous MWD does not have to be assumed in order to obtain the polydispersity index. The technique has been developed to deal with entangled linear polymers, whatever the form of the MWD is. The rheological information required to obtain the polydispersity index is the storage G′(ω) and loss G″(ω) moduli, extending from the terminal zone to the plateau region. The method provides a good agreement between the proposed theoretical approach and the experimental polydispersity indices of several linear polymers for a wide range of average molecular weights and polydispersity indices. It is also applicable to binary blends.

Finite mixture analysis of beauty-contest data using generalised beta distributions

This paper introduces a mixture model based on the beta distribution, without preestablishedmeans and variances, to analyze a large set of Beauty-Contest data obtainedfrom diverse groups of experiments (Bosch-Domenech et al. 2002). This model gives a bettert of the experimental data, and more precision to the hypothesis that a large proportionof individuals follow a common pattern of reasoning, described as iterated best reply (degenerate),than mixture models based on the normal distribution. The analysis shows thatthe means of the distributions across the groups of experiments are pretty stable, while theproportions of choices at dierent levels of reasoning vary across groups.

The normal distribution in some constrained sample spaces

Phenomena with a constrained sample space appear frequently in practice. This is the case e.g. with strictly positive data, or with compositional data, like percentages or proportions. If the natural measure of difference is not the absolute one, simple algebraic properties show that it is more convenient to work with a geometry different from the usual Euclidean geometry in real space, and with a measure different from the usual Lebesgue measure, leading to alternative models which better fit the phenomenon under study. The general approach is presented and illustrated using the normal distribution, both on the positive real line and on the D-part simplex. The original ideas of McAlister in his introduction to the lognormal distribution in 1879, are recovered and updated

Wave-height hazard analysis in Eastern Coast of Spain : Bayesian approach using generalized Pareto distribution

Standard practice of wave-height hazard analysis often pays little attention to the uncertainty of assessed return periods and occurrence probabilities. This fact favors the opinion that, when large events happen, the hazard assessment should change accordingly. However, uncertainty of the hazard estimates is normally able to hide the effect of those large events. This is illustrated using data from the Mediterranean coast of Spain, where the last years have been extremely disastrous. Thus, it is possible to compare the hazard assessment based on data previous to those years with the analysis including them. With our approach, no significant change is detected when the statistical uncertainty is taken into account. The hazard analysis is carried out with a standard model. Time-occurrence of events is assumed Poisson distributed. The wave-height of each event is modelled as a random variable which upper tail follows a Generalized Pareto Distribution (GPD). Moreover, wave-heights are assumed independent from event to event and also independent of their occurrence in time. A threshold for excesses is assessed empirically. The other three parameters (Poisson rate, shape and scale parameters of GPD) are jointly estimated using Bayes' theorem. Prior distribution accounts for physical features of ocean waves in the Mediterranean sea and experience with these phenomena. Posterior distribution of the parameters allows to obtain posterior distributions of other derived parameters like occurrence probabilities and return periods. Predictives are also available. Computations are carried out using the program BGPE v2.0

Analysis of sensitivity with respect to a compositional parameter : A comment on 'Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling' by L. Loosvelt and co-authors

A comment about the article “Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling” writen by L. Loosvelt and co-authors. The present comment is centered in three specific points. The first one is related to the fact that the authors avoid the use of ilr-coordinates. The second one refers to some generalization of sensitivity analysis when input parameters are compositional. The third tries to show that the role of the Dirichlet distribution in the sensitivity analysis is irrelevant

Predicción de la resistencia al flujo en ríos de montaña

Se han calibrado, validado y comparado tres modelos de resistencia al flujo de contorno granular: un modelo potencial y otros dos modelos desarrollados para condiciones de alta rugosidad relativa (uno basado en una modificación de la ley logarítmica de Prandtl-von Karman y otro fundamentado en un perfil de velocidad configurado en dos zonas: una uniforme en las proximidades de los elementos de rugosidad y otra superior que sigue una distribución logarítmica). Se ha empleado para ello un numeroso conjunto de 1.533 datos tomados en ríos y en canales de laboratorio, representativo de un amplio intervalo hidráulico y geomorfológico en el ámbito de ríos de grava y de montaña. Han resultado preferibles las ecuaciones ajustadas con los percentiles granulométricos mayores (d90 o d84) que las ajustadas con el diámetro mediano (d50), debido a la mayor capacidad explicativa alcanzada dado un modelo, la menor diferencia en la bondad de ajuste entre los diferentes modelos y la menor influencia del origen de los datos (río o canal de laboratorio). Las ecuaciones ajustadas de acuerdo con los modelos en donde se contemplan condiciones de alta rugosidad relativa presentan predicciones similares, exceptuando el intervalo macrorrugoso (y/d90 < 1), en el que es preferible la correspondiente al modelo fundamentado en el perfil de velocidad configurado en dos zonas. Se recomienda restringir la aplicación de la ecuación ajustada con arreglo a la ley potencial al intervalo de y/d90 comprendido entre uno y veinte, puesto que fuera de dicho intervalo tiende a infraestimar notablemente la resistencia al flujo.

Equivalent roughness of gravel-bed rivers

The relation between the equivalent roughness and different grain size percentiles of the sediment in gravel-bed rivers was determined under the hypothesis that the vertical distribution of the flow velocity follows a logarithmic law. A set of 954 data points was selected from rivers with gravel size sediment or larger, with a non-sinuous alignment and free of vegetation or obstacles. According to the results, the ks roughness is equivalent to approximately 2.4D90, 2.8D84, and 6.1D50. No correlation was detected between the sediment sorting and the sediment mobility index on one hand, and, on the other, the coefficient of proportionality of each grain size percentile.

Spectral analysis of the angular distribution function of back reflectors for thin film silicon solar cells

Nowadays, one of the most important challenges to enhance the efficiency of thin film silicon solar cells is to increase the short circuit intensity by means of optical confinement methods, such as textured back-reflector structures. In this work, two possible textured structures to be used as back reflectors for n-i-p solar cells have been optically analyzed and compared to a smooth one by using a system which is able to measure the angular distribution function (ADF) of the scattered light in a wide spectral range (350-1000 nm). The accurate analysis of the ADF data corresponding to the reflector structures and to the μc-Si:H films deposited onto them allows the optical losses due to the reflector absorption and its effectiveness in increasing light absorption in the μc-Si:H layer, mainly at long wavelengths, to be quantified.

25-years of biomonitoring in two mediterranean streams (Llobregat and Besos basins, NE Spain)

Biological water quality changes in two Mediterranean river basins from a network of 42 sampling sites assessed since 1979 are presented. In order to characterize the biological quality, the index FBILL, designed to characterize these rivers" quality using aquatic macroinvertebrates, is used. When comparing the data from recent years to older ones, only two headwater sites from the 42 had improved their water quality to good or very good conditions. In the middle or low river basin sites or even in headwater localities were river flow is reduced, the important investment to build up sewage water treatment systems and plants (more than 70 in 15 years) allowed for a small recovery from poor or very poor conditions to moderate water quality. Nevertheless still a significant number (25 %) of the localities remain in poor conditions. The evolution of the quality in several points of both basins shows how the main problems for the recovery of the biological quality is due to the water diverted for small hydraulic plants, the presence of saline pollution in the Llobregat River, and the insufficient water depuration. In the smaller rivers, and specially the Besòs the lack of dilution flows from the treatment plants is the main problem for water quality recovery.