What Belongs Where? Variable Selection for Zero-Inflated Count Models with an Application to the Demand for Health Care

This paper develops stochastic search variable selection (SSVS) for zero-inflated count models which are commonly used in health economics. This allows for either model averaging or model selection in situations with many potential regressors. The proposed techniques are applied to a data set from Germany considering the demand for health care. A package for the free statistical software environment R is provided.

On Identification of Bayesian DSGE Models

In recent years there has been increasing concern about the identification of parameters in dynamic stochastic general equilibrium (DSGE) models. Given the structure of DSGE models it may be difficult to determine whether a parameter is identified. For the researcher using Bayesian methods, a lack of identification may not be evident since the posterior of a parameter of interest may differ from its prior even if the parameter is unidentified. We show that this can even be the case even if the priors assumed on the structural parameters are independent. We suggest two Bayesian identification indicators that do not suffer from this difficulty and are relatively easy to compute. The first applies to DSGE models where the parameters can be partitioned into those that are known to be identified and the rest where it is not known whether they are identified. In such cases the marginal posterior of an unidentified parameter will equal the posterior expectation of the prior for that parameter conditional on the identified parameters. The second indicator is more generally applicable and considers the rate at which the posterior precision gets updated as the sample size (T) is increased. For identified parameters the posterior precision rises with T, whilst for an unidentified parameter its posterior precision may be updated but its rate of update will be slower than T. This result assumes that the identified parameters are pT-consistent, but similar differential rates of updates for identified and unidentified parameters can be established in the case of super consistent estimators. These results are illustrated by means of simple DSGE models.

Model Uncertainty in Panel Vector Autoregressive Models.

We develop methods for Bayesian model averaging (BMA) or selection (BMS) in Panel Vector Autoregressions (PVARs). Our approach allows us to select between or average over all possible combinations of restricted PVARs where the restrictions involve interdependencies between and heterogeneities across cross-sectional units. The resulting BMA framework can find a parsimonious PVAR specification, thus dealing with overparameterization concerns. We use these methods in an application involving the euro area sovereign debt crisis and show that our methods perform better than alternatives. Our findings contradict a simple view of the sovereign debt crisis which divides the euro zone into groups of core and peripheral countries and worries about financial contagion within the latter group.

Model Uncertainty in Panel Vector Autoregressive Models.

We develop methods for Bayesian model averaging (BMA) or selection (BMS) in Panel Vector Autoregressions (PVARs). Our approach allows us to select between or average over all possible combinations of restricted PVARs where the restrictions involve interdependencies between and heterogeneities across cross-sectional units. The resulting BMA framework can find a parsimonious PVAR specification, thus dealing with overparameterization concerns. We use these methods in an application involving the euro area sovereign debt crisis and show that our methods perform better than alternatives. Our findings contradict a simple view of the sovereign debt crisis which divides the euro zone into groups of core and peripheral countries and worries about financial contagion within the latter group.

Computing Markov-Perfect Optimal Policies in Business-Cycle Models

Time-inconsistency is an essential feature of many policy problems (Kydland and Prescott, 1977). This paper presents and compares three methods for computing Markov-perfect optimal policies in stochastic nonlinear business cycle models. The methods considered include value function iteration, generalized Euler-equations, and parameterized shadow prices. In the context of a business cycle model in which a scal authority chooses government spending and income taxation optimally, while lacking the ability to commit, we show that the solutions obtained using value function iteration and generalized Euler equations are somewhat more accurate than that obtained using parameterized shadow prices. Among these three methods, we show that value function iteration can be applied easily, even to environments that include a risk-sensitive scal authority and/or inequality constraints on government spending. We show that the risk-sensitive scal authority lowers government spending and income-taxation, reducing the disincentive households face to accumulate wealth.

Proteomic analysis of cytokine induced proteins in human intestinal epithelial cells: implications for inflammatory bowel diseases.

A role for cytokine regulated proteins in epithelial cells has been suggested in the pathogenesis of inflammatory bowel diseases (IBD). The aim of this study was to identify such cytokine regulated targets using a proteomic functional approach. Protein patterns from (35)S-radiolabeled homogenates of cultured colon epithelial cells were compared before and after exposure to interferon-gamma, interleukin-1beta and interleukin-6. Proteins were separated by two-dimensional polyacrylamide gel electrophoresis. Both autoradiographies and silver stained gels were analyzed. Proteins showing differential expression were identified by tryptic in-gel digestion and mass spectrometry. Metabolism related proteins were also investigated by Western blot analysis. Tryptophanyl-tRNA synthetase, indoleamine-2,3-dioxygenase, heterogeneous nuclear ribonucleoprotein JKTBP, interferon-induced 35kDa protein, proteasome subunit LMP2 and arginosuccinate synthetase were identified as cytokine modulated proteins in vitro. Using purified epithelial cells from patients, overexpression of indoleamine-2,3-dioxygenase, an enzyme involved in tryptophan metabolism, was confirmed in Crohn's disease as well as in ulcerative colitis, as compared to normal mucosa. No such difference was found in diverticulitis. Potentially, this observation opens new avenues in the treatment of IBD.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit.

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Stochastic demography and the neutral substitution rate in class-structured populations.

The neutral rate of allelic substitution is analyzed for a class-structured population subject to a stationary stochastic demographic process. The substitution rate is shown to be generally equal to the effective mutation rate, and under overlapping generations it can be expressed as the effective mutation rate in newborns when measured in units of average generation time. With uniform mutation rate across classes the substitution rate reduces to the mutation rate.

Estimation of dynamic latent variable models using simulated nonparametric moments

Abstract. Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Because conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. It is shown that as the number of simulations diverges, the estimator is consistent and a higher-order expansion reveals the stochastic difference between the infeasible GMM estimator based on the same moment conditions and the simulated version. In particular, we show how to adjust standard errors to account for the simulations. Monte Carlo results show how the estimator may be applied to a range of dynamic latent variable (DLV) models, and that it performs well in comparison to several other estimators that have been proposed for DLV models.

Three Essays on the Formation and Impact of Economic Policy

Introduction In my thesis I argue that economic policy is all about economics and politics. Consequently, analysing and understanding economic policy ideally has at least two parts. The economics part, which is centered around the expected impact of a specific policy on the real economy both in terms of efficiency and equity. The insights of this part point into which direction the fine-tuning of economic policies should go. However, fine-tuning of economic policies will be most likely subject to political constraints. That is why, in the politics part, a much better understanding can be gained by taking into account how the incentives of politicians and special interest groups as well as the role played by different institutional features affect the formation of economic policies. The first part and chapter of my thesis concentrates on the efficiency-related impact of economic policies: how does corporate income taxation in general, and corporate income tax progressivity in specific, affect the creation of new firms? Reduced progressivity and flat-rate taxes are in vogue. By 2009, 22 countries are operating flat-rate income tax systems, as do 7 US states and 14 Swiss cantons (for corporate income only). Tax reform proposals in the spirit of the "flat tax" model typically aim to reduce three parameters: the average tax burden, the progressivity of the tax schedule, and the complexity of the tax code. In joint work, Marius Brülhart and I explore the implications of changes in these three parameters on entrepreneurial activity, measured by counts of firm births in a panel of Swiss municipalities. Our results show that lower average tax rates and reduced complexity of the tax code promote firm births. Controlling for these effects, reduced progressivity inhibits firm births. Our reading of these results is that tax progressivity has an insurance effect that facilitates entrepreneurial risk taking. The positive effects of lower tax levels and reduced complexity are estimated to be significantly stronger than the negative effect of reduced progressivity. To the extent that firm births reflect desirable entrepreneurial dynamism, it is not the flattening of tax schedules that is key to successful tax reforms, but the lowering of average tax burdens and the simplification of tax codes. Flatness per se is of secondary importance and even appears to be detrimental to firm births. The second part of my thesis, which corresponds to the second and third chapter, concentrates on how economic policies are formed. By the nature of the analysis, these two chapters draw on a broader literature than the first chapter. Both economists and political scientists have done extensive research on how economic policies are formed. Thereby, researchers in both disciplines have recognised the importance of special interest groups trying to influence policy-making through various channels. In general, economists base their analysis on a formal and microeconomically founded approach, while abstracting from institutional details. In contrast, political scientists' frameworks are generally richer in terms of institutional features but lack the theoretical rigour of economists' approaches. I start from the economist's point of view. However, I try to borrow as much as possible from the findings of political science to gain a better understanding of how economic policies are formed in reality. In the second chapter, I take a theoretical approach and focus on the institutional policy framework to explore how interactions between different political institutions affect the outcome of trade policy in presence of special interest groups' lobbying. Standard political economy theory treats the government as a single institutional actor which sets tariffs by trading off social welfare against contributions from special interest groups seeking industry-specific protection from imports. However, these models lack important (institutional) features of reality. That is why, in my model, I split up the government into a legislative and executive branch which can both be lobbied by special interest groups. Furthermore, the legislative has the option to delegate its trade policy authority to the executive. I allow the executive to compensate the legislative in exchange for delegation. Despite ample anecdotal evidence, bargaining over delegation of trade policy authority has not yet been formally modelled in the literature. I show that delegation has an impact on policy formation in that it leads to lower equilibrium tariffs compared to a standard model without delegation. I also show that delegation will only take place if the lobby is not strong enough to prevent it. Furthermore, the option to delegate increases the bargaining power of the legislative at the expense of the lobbies. Therefore, the findings of this model can shed a light on why the U.S. Congress often practices delegation to the executive. In the final chapter of my thesis, my coauthor, Antonio Fidalgo, and I take a narrower approach and focus on the individual politician level of policy-making to explore how connections to private firms and networks within parliament affect individual politicians' decision-making. Theories in the spirit of the model of the second chapter show how campaign contributions from lobbies to politicians can influence economic policies. There exists an abundant empirical literature that analyses ties between firms and politicians based on campaign contributions. However, the evidence on the impact of campaign contributions is mixed, at best. In our paper, we analyse an alternative channel of influence in the shape of personal connections between politicians and firms through board membership. We identify a direct effect of board membership on individual politicians' voting behaviour and an indirect leverage effect when politicians with board connections influence non-connected peers. We assess the importance of these two effects using a vote in the Swiss parliament on a government bailout of the national airline, Swissair, in 2001, which serves as a natural experiment. We find that both the direct effect of connections to firms and the indirect leverage effect had a strong and positive impact on the probability that a politician supported the government bailout.

Agglutinating Secretory IgA Preserves Intestinal Epithelial Cell Integrity during Apical Infection by Shigella flexneri.

Shigella flexneri, by invading intestinal epithelial cells (IECs) and inducing inflammatory responses of the colonic mucosa, causes bacillary dysentery. Although M cells overlying Peyer's patches are commonly considered the primary site of entry of S. flexneri, indirect evidence suggests that bacteria can also use IECs as a portal of entry to the lamina propria. Passive delivery of secretory IgA (SIgA), the major immunoglobulin secreted at mucosal surfaces, has been shown to protect rabbits from experimental shigellosis, but no information exists as to its molecular role in maintaining luminal epithelial integrity. We have established that the interaction of virulent S. flexneri with the apical pole of a model intestinal epithelium consisting of polarized Caco-2 cell monolayers resulted in the progressive disruption of the tight junction network and actin depolymerization, eventually resulting in cell death. The lipopolysaccharide (LPS)-specific agglutinating SIgAC5 monoclonal antibody (MAb), but not monomeric IgAC5 or IgGC20 MAbs of the same specificity, achieved protective functions through combined mechanisms, including limitation of the interaction between S. flexneri and epithelial cells, maintenance of the tight junction seal, preservation of the cell morphology, reduction of NF-κB nuclear translocation, and inhibition of proinflammatory mediator secretion. Our results add to the understanding of the function of SIgA-mediated immune exclusion by identifying a mode of action whereby the formation of immune complexes translates into maintenance of the integrity of epithelial cells lining the mucosa. This novel mechanism of protection mediated by SIgA is important to extend the arsenal of effective strategies to fight against S. flexneri mucosal invasion.

Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

Oxidative stress and inflammation response after nanoparticle exposure : differences between human lung cell monocultures and an advanced three-dimensional model of the human epithelial airways

Combustion-derived and manufactured nanoparticles (NPs) are known to provoke oxidative stress and inflammatory responses in human lung cells; therefore, they play an important role during the development of adverse health effects. As the lungs are composed of more than 40 different cell types, it is of particular interest to perform toxicological studies with co-cultures systems, rather than with monocultures of only one cell type, to gain a better understanding of complex cellular reactions upon exposure to toxic substances. Monocultures of A549 human epithelial lung cells, human monocyte-derived macrophages and monocyte-derived dendritic cells (MDDCs) as well as triple cell co-cultures consisting of all three cell types were exposed to combustion-derived NPs (diesel exhaust particles) and to manufactured NPs (titanium dioxide and single-walled carbon nanotubes). The penetration of particles into cells was analysed by transmission electron microscopy. The amount of intracellular reactive oxygen species (ROS), the total antioxidant capacity (TAC) and the production of tumour necrosis factor (TNF)-a and interleukin (IL)-8 were quantified. The results of the monocultures were summed with an adjustment for the number of each single cell type in the triple cell co-culture. All three particle types were found in all cell and culture types. The production of ROS was induced by all particle types in all cell cultures except in monocultures of MDDCs. The TAC and the (pro-)inflammatory reactions were not statistically significantly increased by particle exposure in any of the cell cultures. Interestingly, in the triple cell co-cultures, the TAC and IL-8 concentrations were lower and the TNF-a concentrations were higher than the expected values calculated from the monocultures. The interplay of different lung cell types seems to substantially modulate the oxidative stress and the inflammatory responses after NP exposure. [Authors]

What do we gain from simplicity versus complexity in species distribution models?

Species distribution models (SDMs) are widely used to explain and predict species ranges and environmental niches. They are most commonly constructed by inferring species' occurrence-environment relationships using statistical and machine-learning methods. The variety of methods that can be used to construct SDMs (e.g. generalized linear/additive models, tree-based models, maximum entropy, etc.), and the variety of ways that such models can be implemented, permits substantial flexibility in SDM complexity. Building models with an appropriate amount of complexity for the study objectives is critical for robust inference. We characterize complexity as the shape of the inferred occurrence-environment relationships and the number of parameters used to describe them, and search for insights into whether additional complexity is informative or superfluous. By building 'under fit' models, having insufficient flexibility to describe observed occurrence-environment relationships, we risk misunderstanding the factors shaping species distributions. By building 'over fit' models, with excessive flexibility, we risk inadvertently ascribing pattern to noise or building opaque models. However, model selection can be challenging, especially when comparing models constructed under different modeling approaches. Here we argue for a more pragmatic approach: researchers should constrain the complexity of their models based on study objective, attributes of the data, and an understanding of how these interact with the underlying biological processes. We discuss guidelines for balancing under fitting with over fitting and consequently how complexity affects decisions made during model building. Although some generalities are possible, our discussion reflects differences in opinions that favor simpler versus more complex models. We conclude that combining insights from both simple and complex SDM building approaches best advances our knowledge of current and future species ranges.

Wnt secretion and gradient formation.

Concentration gradients formed by the lipid-modified morphogens of the Wnt family are known for their pivotal roles during embryogenesis and adult tissue homeostasis. Wnt morphogens are also implicated in a variety of human diseases, especially cancer. Therefore, the signaling cascades triggered by Wnts have received considerable attention during recent decades. However, how Wnts are secreted and how concentration gradients are formed remains poorly understood. The use of model organisms such as Drosophila melanogaster has provided important advances in this area. For instance, we have previously shown that the lipid raft-associated reggie/flotillin proteins influence Wnt secretion and spreading in Drosophila. Our work supports the notion that producing cells secrete Wnt molecules in at least two pools: a poorly diffusible one and a reggie/flotillin-dependent highly diffusible pool which allows morphogen spreading over long distances away from its source of production. Here we revise the current views of Wnt secretion and spreading, and propose two models for the role of the reggie/flotillin proteins in these processes: (i) reggies/flotillins regulate the basolateral endocytosis of the poorly diffusible, membrane-bound Wnt pool, which is then sorted and secreted to apical compartments for long-range diffusion, and (ii) lipid rafts organized by reggies/flotillins serve as "dating points" where extracellular Wnt transiently interacts with lipoprotein receptors to allow its capture and further spreading via lipoprotein particles. We further discuss these processes in the context of human breast cancer. A better understanding of these phenomena may be relevant for identification of novel drug targets and therapeutic strategies.