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

Impact of incorrect assumptions on the covariance structure of random effects and/or residuals in nonlinear mixed models for repeated measures data

In this paper we analyse, using Monte Carlo simulation, the possible consequences of incorrect assumptions on the true structure of the random effects covariance matrix and the true correlation pattern of residuals, over the performance of an estimation method for nonlinear mixed models. The procedure under study is the well known linearization method due to Lindstrom and Bates (1990), implemented in the nlme library of S-Plus and R. Its performance is studied in terms of bias, mean square error (MSE), and true coverage of the associated asymptotic confidence intervals. Ignoring other criteria like the convenience of avoiding over parameterised models, it seems worst to erroneously assume some structure than do not assume any structure when this would be adequate.

Effects of new nonlinear couplings in relativistic effective field theory

We extend the relativistic mean field theory model of Sugahara and Toki by adding new couplings suggested by modern effective field theories. An improved set of parameters is developed with the goal to test the ability of the models based on effective field theory to describe the properties of finite nuclei and, at the same time, to be consistent with the trends of Dirac-Brueckner-Hartree-Fock calculations at densities away from the saturation region. We compare our calculations with other relativistic nuclear force parameters for various nuclear phenomena.

Transient and preparation colored-noise effects: The nonlinear relaxation-time approach

We develop a singular perturbation approach to the problem of the calculation of a characteristic time (the nonlinear relaxation time) for non-Markovian processes driven by Gaussian colored noise with small correlation time. Transient and initial preparation effects are discussed and explicit results for prototype situations are obtained. New effects on the relaxation of unstable states are predicted. The approach is compared with previous techniques.

Decay of unstable states in the presence of colored noise and random initial conditions. I. Theory of nonlinear relaxation times

The general theory of nonlinear relaxation times is developed for the case of Gaussian colored noise. General expressions are obtained and applied to the study of the characteristic decay time of unstable states in different situations, including white and colored noise, with emphasis on the distributed initial conditions. Universal effects of the coupling between colored noise and random initial conditions are predicted.

Nonlinear Saffman-Taylor instability

We show, both theoretically and experimentally, that the interface between two viscous fluids in a Hele-Shaw cell can be nonlinearly unstable before the Saffman-Taylor linear instability point is reached. We identify the family of exact elastica solutions [Nye et al., Eur. J. Phys. 5, 73 (1984)] as the unstable branch of the corresponding subcritical bifurcation which ends up at a topological singularity defined by interface pinchoff. We devise an experimental procedure to prepare arbitrary initial conditions in a Hele-Shaw cell. This is used to test the proposed bifurcation scenario and quantitatively asses its practical relevance.

Mapping, organic matter mass and water volume of a peatland in Serra do Espinhaço Meridional

Peatlands form in areas where net primary of organic matter production exceeds losses due to the decomposition, leaching or disturbance. Due to their chemical and physical characteristics, bogs can influence water dynamics because they can store large volumes of water in the rainy season and gradually release this water during the other months of the year. In Diamantina, Minas Gerais, Brazil, a peatland in the environmental protection area of Pau-de-Fruta ensures the water supply of 40,000 inhabitants. The hypothesis of this study is that the peat bogs in Pau-de-Fruta act as an environment for carbon storage and a regulator of water flow in the Córrego das Pedras basin. The objective of this study was to estimate the water volume and organic matter mass in this peatland and to study the influence of this environment on the water flow in the Córrego das Pedras basin. The peatland was mapped using 57 transects, at intervals of 100 m. Along all transects, the depth of the peat bog, the Universal Transverse Mercator (UTM) coordinates and altitude were recorded every 20 m and used to calculate the area and volume of the peatland. The water volume was estimated, using a method developed in this study, and the mass of organic matter based on samples from 106 profiles. The peatland covered 81.7 hectares (ha), and stored 497,767 m³ of water, representing 83.7 % of the total volume of the peat bog. The total amount of organic matter (OM) was 45,148 t, corresponding to 552 t ha-1 of OM. The peat bog occupies 11.9 % of the area covered by the Córrego das Pedras basin and stores 77.6 % of the annual water surplus, thus controlling the water flow in the basin and consequently regulating the water course.

Systematic weakly nonlinear analysis of interfacial instabilities in Hele-Shaw flows

We develop a systematic method to derive all orders of mode couplings in a weakly nonlinear approach to the dynamics of the interface between two immiscible viscous fluids in a Hele-Shaw cell. The method is completely general: it applies to arbitrary geometry and driving. Here we apply it to the channel geometry driven by gravity and pressure. The finite radius of convergence of the mode-coupling expansion is found. Calculation up to third-order couplings is done, which is necessary to account for the time-dependent Saffman-Taylor finger solution and the case of zero viscosity contrast. The explicit results provide relevant analytical information about the role that the viscosity contrast and the surface tension play in the dynamics of the system. We finally check the quantitative validity of different orders of approximation and a resummation scheme against a physically relevant, exact time-dependent solution. The agreement between the low-order approximations and the exact solution is excellent within the radius of convergence, and is even reasonably good beyond this radius.

Systematic weakly nonlinear analysis of radial viscous fingering

We present a weakly nonlinear analysis of the interface dynamics in a radial Hele-Shaw cell driven by both injection and rotation. We extend the systematic expansion introduced in [E. Alvarez-Lacalle et al., Phys. Rev. E 64, 016302 (2001)] to the radial geometry, and compute explicitly the first nonlinear contributions. We also find the necessary and sufficient condition for the uniform convergence of the nonlinear expansion. Within this region of convergence, the analytical predictions at low orders are compared satisfactorily to exact solutions and numerical integration of the problem. This is particularly remarkable in configurations (with no counterpart in the channel geometry) for which the interplay between injection and rotation allows that condition to be satisfied at all times. In the case of the purely centrifugal forcing we demonstrate that nonlinear couplings make the interface more unstable for lower viscosity contrast between the fluids.

The nonlinear relaxation time and quasideterministic-theory approaches to characterize the decay of unestable states

We present the relationship between nonlinear-relaxation-time (NLRT) and quasideterministic approaches to characterize the decay of an unstable state. The universal character of the NLRT is established. The theoretical results are applied to study the dynamical relaxation of the Landau model in one and n variables and also a laser model.

Digital soil mapping: strategy for data pre-processing

The region of greatest variability on soil maps is along the edge of their polygons, causing disagreement among pedologists about the appropriate description of soil classes at these locations. The objective of this work was to propose a strategy for data pre-processing applied to digital soil mapping (DSM). Soil polygons on a training map were shrunk by 100 and 160 m. This strategy prevented the use of covariates located near the edge of the soil classes for the Decision Tree (DT) models. Three DT models derived from eight predictive covariates, related to relief and organism factors sampled on the original polygons of a soil map and on polygons shrunk by 100 and 160 m were used to predict soil classes. The DT model derived from observations 160 m away from the edge of the polygons on the original map is less complex and has a better predictive performance.

Confinement of discrete breathers in inhomogeneously profiled nonlinear chains

We investigate numerically the scattering of a moving discrete breather on a pair of junctions in a Fermi-Pasta-Ulam chain. These junctions delimit an extended region with different masses of the particles. We consider (i) a rectangular trap, (ii) a wedge shaped trap, and (iii) a smoothly varying convex or concave mass profile. All three cases lead to DB confinement, with the ease of trapping depending on the profile of the trap. We also study the collision and trapping of two DBs within the profile as a function of trap width, shape, and approach time at the two junctions. The latter controls whether one or both DBs are trapped.

From lattice-gas models to nonlinear diffusion models: A derivation of macroscopic equations and fluctuations

We derive nonlinear diffusion equations and equations containing corrections due to fluctuations for a coarse-grained concentration field. To deal with diffusion coefficients with an explicit dependence on the concentration values, we generalize the Van Kampen method of expansion of the master equation to field variables. We apply these results to the derivation of equations of phase-separation dynamics and interfacial growth instabilities.

Nonlinear relaxation time and the detection of weak signals

Laser systems can be used to detect very weak optical signals. The physical mechanism is the dynamical process of the relaxation of a laser from an unstable state to a steady stable state. We present an analysis of this process based on the study of the nonlinear relaxation time. Our analytical results are compared with numerical integration of the stochastic differential equations that model this process.