Finite mixture regression model with random effects: application to neonatal hospital length of stay

A two-component mixture regression model that allows simultaneously for heterogeneity and dependency among observations is proposed. By specifying random effects explicitly in the linear predictor of the mixture probability and the mixture components, parameter estimation is achieved by maximising the corresponding best linear unbiased prediction type log-likelihood. Approximate residual maximum likelihood estimates are obtained via an EM algorithm in the manner of generalised linear mixed model (GLMM). The method can be extended to a g-component mixture regression model with the component density from the exponential family, leading to the development of the class of finite mixture GLMM. For illustration, the method is applied to analyse neonatal length of stay (LOS). It is shown that identification of pertinent factors that influence hospital LOS can provide important information for health care planning and resource allocation. (C) 2002 Elsevier Science B.V. All rights reserved.

Optimal debt, endogenous fertility, and human capital externalities in a model with altruistic bequests

and human capital externalities. Because of such externalities, education investment is too low and fertility is too high. While education subsidies are the conventional means to deal with these problems, we show that the optimal policy also comprises debt even when distortionary taxes are used. The reason is that debt tips the usual trade-off between children's quantity and quality in favor of the latter by increasing the bequest cost of children. The optimal debt-output ratio exceeds 10% for plausible parameterization. (C) 2002 Elsevier B.V. All rights reserved.

Optimal price regulation in a growth model with monopolistic suppliers of intermediate goods

In this paper we investigate the trade-off faced by regulators who must set a price for an intermediate good somewhere between the marginal cost and the monopoly price. We utilize a growth model with monopolistic suppliers of intermediate goods. Investment in innovation is required to produce a new intermediate good. Marginal cost pricing deters innovation, while monopoly pricing maximizes innovation and economic growth at the cost of some static inefficiency. We demonstrate the existence of a second-best price above the marginal cost but below the monopoly price, which maximizes consumer welfare. Simulation results suggest that substantial reductions in consumption, production, growth, and welfare occur where regulators focus on static efficiency issues by setting prices at or near marginal cost.

Exact solution of the XXZ Gaudin model with generic open boundaries

The XXZ Gaudin model with generic integrable boundaries specified by generic non-diagonal K-matrices is studied. The commuting families of Gaudin operators are diagonalized by the algebraic Bethe ansatz method. The eigenvalues and the corresponding Bethe ansatz equations are obtained. (C) 2004 Elsevier B.V. All rights reserved.

A carbon activation model with application to longan seed char gasification

In this paper a new structural model is presented to describe the evolution of porosity of char during the gasification process. The model assumes the char structure to be composed of bundles of parallel graphite layers, and the reactivities of each layer with the gasification agent are assumed to be different to represent the different degree of heterogeneity of each layer (i.e. each layer will react with the gasification agent at a different rate). It is this difference in the reactivity that allows micropores to be created during the course of gasification. This simple structural model enables the evolution of pore volume, pore geometrical surface area and the pore size distribution to be described with respect to the extent of char burn-off. The model is tested against the experimental data of gasification of longan seed-derived char with carbon dioxide and it is found that the agreement between the model and the data is reasonably satisfactory, especially the evolution of surface area and pore volume with burn-off.

A(n-1) Gaudin model with open boundaries

The A(n-1) Gaudin model with integrable boundaries specified by non-diagonal K-matrices is studied. The commuting families of Gaudin operators are diagonalized by the algebraic Bethe ansatz method. The eigenvalues and the corresponding Bethe ansatz equations are obtained. (c) 2005 Elsevier B.V. All rights reserved.

Stability and cycles in a cobweb model with heterogeneous expectations

We investigate the dynamics of a cobweb model with heterogeneous beliefs, generalizing the example of Brock and Hommes (1997). We examine situations where the agents form expectations by using either rational expectations, or a type of adaptive expectations with limited memory defined from the last two prices. We specify conditions that generate cycles. These conditions depend on a set of factors that includes the intensity of switching between beliefs and the adaption parameter. We show that both Flip bifurcation and Neimark-Sacker bifurcation can occur as primary bifurcation when the steady state is unstable.

A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

A scatterometer neural network sensor model with input noise

The ERS-1 satellite carries a scatterometer which measures the amount of radiation scattered back toward the satellite by the ocean's surface. These measurements can be used to infer wind vectors. The implementation of a neural network based forward model which maps wind vectors to radar backscatter is addressed. Input noise cannot be neglected. To account for this noise, a Bayesian framework is adopted. However, Markov Chain Monte Carlo sampling is too computationally expensive. Instead, gradient information is used with a non-linear optimisation algorithm to find the maximum em a posteriori probability values of the unknown variables. The resulting models are shown to compare well with the current operational model when visualised in the target space.

A scatterometer neural network sensor model with input noise

The ERS-1 satellite carries a scatterometer which measures the amount of radiation scattered back toward the satellite by the ocean's surface. These measurements can be used to infer wind vectors. The implementation of a neural network based forward model which maps wind vectors to radar backscatter is addressed. Input noise cannot be neglected. To account for this noise, a Bayesian framework is adopted. However, Markov Chain Monte Carlo sampling is too computationally expensive. Instead, gradient information is used with a non-linear optimisation algorithm to find the maximum em a posteriori probability values of the unknown variables. The resulting models are shown to compare well with the current operational model when visualised in the target space.

A flexible cost function model with risk

In examining bank cost efficiency in banking inclusion of risk-taking of banks is very important. In this paper we depart from the standard modeling approach and view risk intimately related to the technology. Thus, instead of controlling for risk by viewing them as covariates in the standard cost function we argue that the technology differs with risk, thereby meaning that the parameters of the parametric cost function changes with risk in a fully flexible manner. This is accomplished by viewing the parameters of the cost function as nonparametric functions of risk. We also control for country-specific effects in a fully flexible manner by using them as arguments of the nonparametric functions along with the risk variable. The resulting cost function then becomes semiparametric. The standard parametric model becomes a special case of our semiparametric model. We use the above modeling approach for banks in the EU countries. Actually, European financial integration is seen as a stepping stone for the development of a competitive single EU market that promotes efficiency and increases consumer welfare, changing the risk profile of the European banks. Particularly, financial integration allows more risk diversification and permits banks to use more advanced risk management instruments and systems, however it has at the same time increased the probability of systematic risks. Financial integration has increased the risk of contagion and changed its nature and scope. Consequently the bank’s risk seems to be an important issue to be investigated.

Asymptotically exact probability distribution for the Sinai model with finite drift

We obtain the exact asymptotic result for the disorder-averaged probability distribution function for a random walk in a biased Sinai model and show that it is characterized by a creeping behavior of the displacement moments with time, similar to v(mu n), where mu <1 is dimensionless mean drift. We employ a method originated in quantum diffusion which is based on the exact mapping of the problem to an imaginary-time Schrodinger equation. For nonzero drift such an equation has an isolated lowest eigenvalue separated by a gap from quasicontinuous excited states, and the eigenstate corresponding to the former governs the long-time asymptotic behavior.