Fourth order pseudo maximum likelihood methods

We extend PML theory to account for information on the conditional moments up to order four, but without assuming a parametric model, to avoid a risk of misspecification of the conditional distribution. The key statistical tool is the quartic exponential family, which allows us to generalize the PML2 and QGPML1 methods proposed in Gourieroux et al. (1984) to PML4 and QGPML2 methods, respectively. An asymptotic theory is developed. The key numerical tool that we use is the Gauss-Freud integration scheme that solves a computational problem that has previously been raised in several fields. Simulation exercises demonstrate the feasibility and robustness of the methods [Authors]

Notch1 and T-cell development: insights from conditional knockout mice.

Notch proteins influence cell-fate decisions in many developmental systems. Gain-of-function studies have suggested a crucial role for Notch1 signaling at several stages during lymphocyte development, including the B/T, alphabeta/gammadelta and CD4/CD8 lineage choices. Here, we critically re-evaluate these conclusions in the light of recent studies that describe inducible and tissue-specific targeting of the Notch1 gene.

The use of the likelihood ratio for evaluative and investigative purposes in comparative forensic handwriting examination

This paper extends previous research and discussion on the use of multivariate continuous data, which are about to become more prevalent in forensic science. As an illustrative example, attention is drawn here on the area of comparative handwriting examinations. Multivariate continuous data can be obtained in this field by analysing the contour shape of loop characters through Fourier analysis. This methodology, based on existing research in this area, allows one describe in detail the morphology of character contours throughout a set of variables. This paper uses data collected from female and male writers to conduct a comparative analysis of likelihood ratio based evidence assessment procedures in both, evaluative and investigative proceedings. While the use of likelihood ratios in the former situation is now rather well established (typically, in order to discriminate between propositions of authorship of a given individual versus another, unknown individual), focus on the investigative setting still remains rather beyond considerations in practice. This paper seeks to highlight that investigative settings, too, can represent an area of application for which the likelihood ratio can offer a logical support. As an example, the inference of gender of the writer of an incriminated handwritten text is forwarded, analysed and discussed in this paper. The more general viewpoint according to which likelihood ratio analyses can be helpful for investigative proceedings is supported here through various simulations. These offer a characterisation of the robustness of the proposed likelihood ratio methodology.

In vivo conditional Pax4 overexpression in mature islet β-cells prevents stress-induced hyperglycemia in mice.

OBJECTIVE To establish the role of the transcription factor Pax4 in pancreatic islet expansion and survival in response to physiological stress and its impact on glucose metabolism, we generated transgenic mice conditionally and selectively overexpressing Pax4 or a diabetes-linked mutant variant (Pax4R129 W) in β-cells. RESEARCH DESIGN AND METHODS Glucose homeostasis and β-cell death and proliferation were assessed in Pax4- or Pax4R129 W-overexpressing transgenic animals challenged with or without streptozotocin. Isolated transgenic islets were also exposed to cytokines, and apoptosis was evaluated by DNA fragmentation or cytochrome C release. The expression profiles of proliferation and apoptotic genes and β-cell markers were studied by immunohistochemistry and quantitative RT-PCR. RESULTS Pax4 but not Pax4R129 W protected animals against streptozotocin-induced hyperglycemia and isolated islets from cytokine-mediated β-cell apoptosis. Cytochrome C release was abrogated in Pax4 islets treated with cytokines. Interleukin-1β transcript levels were suppressed in Pax4 islets, whereas they were increased along with NOS2 in Pax4R129 W islets. Bcl-2, Cdk4, and c-myc expression levels were increased in Pax4 islets while MafA, insulin, and GLUT2 transcript levels were suppressed in both animal models. Long-term Pax4 expression promoted proliferation of a Pdx1-positive cell subpopulation while impeding insulin secretion. Suppression of Pax4 rescued this defect with a concomitant increase in pancreatic insulin content. CONCLUSIONS Pax4 protects adult islets from stress-induced apoptosis by suppressing selective nuclear factor-κB target genes while increasing Bcl-2 levels. Furthermore, it promotes dedifferentiation and proliferation of β-cells through MafA repression, with a concomitant increase in Cdk4 and c-myc expression.

The use of classification and regression trees to predict the likelihood of seasonal influenza.

Background Individual signs and symptoms are of limited value for the diagnosis of influenza. Objective To develop a decision tree for the diagnosis of influenza based on a classification and regression tree (CART) analysis. Methods Data from two previous similar cohort studies were assembled into a single dataset. The data were randomly divided into a development set (70%) and a validation set (30%). We used CART analysis to develop three models that maximize the number of patients who do not require diagnostic testing prior to treatment decisions. The validation set was used to evaluate overfitting of the model to the training set. Results Model 1 has seven terminal nodes based on temperature, the onset of symptoms and the presence of chills, cough and myalgia. Model 2 was a simpler tree with only two splits based on temperature and the presence of chills. Model 3 was developed with temperature as a dichotomous variable (≥38°C) and had only two splits based on the presence of fever and myalgia. The area under the receiver operating characteristic curves (AUROCC) for the development and validation sets, respectively, were 0.82 and 0.80 for Model 1, 0.75 and 0.76 for Model 2 and 0.76 and 0.77 for Model 3. Model 2 classified 67% of patients in the validation group into a high- or low-risk group compared with only 38% for Model 1 and 54% for Model 3. Conclusions A simple decision tree (Model 2) classified two-thirds of patients as low or high risk and had an AUROCC of 0.76. After further validation in an independent population, this CART model could support clinical decision making regarding influenza, with low-risk patients requiring no further evaluation for influenza and high-risk patients being candidates for empiric symptomatic or drug therapy.

The likelihood approach to compare populations : a study on DNA evidence and pitfalls of intuitions

The paper follows on from earlier work [Taroni F and Aitken CGG. Probabilistic reasoning in the law, Part 1: assessment of probabilities and explanation of the value of DNA evidence. Science & Justice 1998; 38: 165-177]. Different explanations of the value of DNA evidence were presented to students from two schools of forensic science and to members of fifteen laboratories all around the world. The responses were divided into two groups; those which came from a school or laboratory identified as Bayesian and those which came from a school or laboratory identified as non-Bayesian. The paper analyses these responses using a likelihood approach. This approach is more consistent with a Bayesian analysis than one based on a frequentist approach, as was reported by Taroni F and Aitken CGG. [Probabilistic reasoning in the law, Part 1: assessment of probabilities and explanation of the value of DNA evidence] in Science & Justice 1998.

Destruction of conditional insulinoma cell lines in NOD mice: a role for autoimmunity.

AIMS/HYPOTHESIS: betaTC-tet (H2(k)) is a conditional insulinoma cell line derived from transgenic mice expressing a tetracycline-regulated oncogene. Transgenic expression of several proteins implicated in the apoptotic pathways increase the resistance of betaTC-tet cells in vitro. We tested in vivo the sensitivity of the cells to rejection and the protective effect of genetic alterations in NOD mice. METHODS: betaTC-tet cells and genetically engineered lines expressing Bcl-2 (CDM3D), a dominant negative mutant of MyD88 or SOCS-1 were transplanted in diabetic female NOD mice or in male NOD mice with diabetes induced by high-dose streptozotocin. Survival of functional cell grafts in NOD-scid mice was also analyzed after transfer of splenocytes from diabetic NOD mice. Autoreactive T-cell hybridomas and splenocytes from diabetic NOD mice were stimulated by betaTC-tet cells. RESULTS: betaTC-tet cells and genetically engineered cell lines were all similarly rejected in diabetic NOD mice and in NOD-scid mice after splenocyte transfer. In 3- to 6-week-old male NOD mice treated with high-dose streptozotocin, the cells temporarily survived, in contrast with C57BL/6 mice treated with high-dose streptozotocin (indefinite survival) and untreated 3- to 6-week-old male NOD mice (rejection). The protective effect of high-dose streptozotocin was lost in older male NOD mice. betaTC-tet cells did not stimulate autoreactive T-cell hybridomas, but induced IL-2 secretion by splenocytes from diabetic NOD mice. CONCLUSION/INTERPRETATION: The autoimmune process seems to play an important role in the destruction of betaTC-tet cells in NOD mice. Genetic manipulations intended at increasing the resistance of beta cells were inefficient. Similar approaches should be tested in vivo as well as in vitro. High dose streptozotocin influences immune rejection and should be used with caution.

Dynamics of the recruitment of lateral amygdala neurons following successive fear learning to different conditional stimuli

In fear conditioning, an animal learns to associate an unconditioned stimulus (US), such as a shock, and a conditioned stimulus (CS), such as a tone, so that the presentation of the CS alone can trigger conditioned responses. Recent research on the lateral amygdala has shown that following cued fear conditioning, only a subset of higher-excitable neurons are recruited in the memory trace. Their selective deletion after fear conditioning results in a selective erasure of the fearful memory. I hypothesize that the recruitment of highly excitable neurons depends on responsiveness to stimuli, intrinsic excitability and local connectivity. In addition, I hypothesize that neurons recruited for an initial memory also participate in subsequent memories, and that changes in neuronal excitability affect secondary fear learning. To address these hypotheses, I will show that A) a rat can learn to associate two successive short-term fearful memories; B) neuronal populations in the LA are competitively recruited in the memory traces depending on individual neuronal advantages, as well as advantages granted by the local network. By performing two successive cued fear conditioning experiments, I found that rats were able to learn and extinguish the two successive short-term memories, when tested 1 hour after learning for each memory. These rats were equipped with a system of stable extracellular recordings that I developed, which allowed to monitor neuronal activity during fear learning. 233 individual putative pyramidal neurons could modulate their firing rate in response to the conditioned tone (conditioned neurons) and/or non- conditioned tones (generalizing neurons). Out of these recorded putative pyramidal neurons 86 (37%) neurons were conditioned to one or both tones. More precisely, one population of neurons encoded for a shared memory while another group of neurons likely encoded the memories' new features. Notably, in spite of a successful behavioral extinction, the firing rate of those conditioned neurons in response to the conditioned tone remained unchanged throughout memory testing. Furthermore, by analyzing the pre-conditioning characteristics of the conditioned neurons, I determined that it was possible to predict neuronal recruitment based on three factors: 1) initial sensitivity to auditory inputs, with tone-sensitive neurons being more easily recruited than tone- insensitive neurons; 2) baseline excitability levels, with more highly excitable neurons being more likely to become conditioned; and 3) the number of afferent connections received from local neurons, with neurons destined to become conditioned receiving more connections than non-conditioned neurons. - En conditionnement de la peur, un animal apprend à associer un stimulus inconditionnel (SI), tel un choc électrique, et un stimulus conditionné (SC), comme un son, de sorte que la présentation du SC seul suffit pour déclencher des réflexes conditionnés. Des recherches récentes sur l'amygdale latérale (AL) ont montré que, suite au conditionnement à la peur, seul un sous-ensemble de neurones plus excitables sont recrutés pour constituer la trace mnésique. Pour apprendre à associer deux sons au même SI, je fais l'hypothèse que les neurones entrent en compétition afin d'être sélectionnés lors du recrutement pour coder la trace mnésique. Ce recrutement dépendrait d'un part à une activation facilité des neurones ainsi qu'une activation facilité de réseaux de neurones locaux. En outre, je fais l'hypothèse que l'activation de ces réseaux de l'AL, en soi, est suffisante pour induire une mémoire effrayante. Pour répondre à ces hypothèses, je vais montrer que A) selon un processus de mémoire à court terme, un rat peut apprendre à associer deux mémoires effrayantes apprises successivement; B) des populations neuronales dans l'AL sont compétitivement recrutées dans les traces mnésiques en fonction des avantages neuronaux individuels, ainsi que les avantages consentis par le réseau local. En effectuant deux expériences successives de conditionnement à la peur, des rats étaient capables d'apprendre, ainsi que de subir un processus d'extinction, pour les deux souvenirs effrayants. La mesure de l'efficacité du conditionnement à la peur a été effectuée 1 heure après l'apprentissage pour chaque souvenir. Ces rats ont été équipés d'un système d'enregistrements extracellulaires stables que j'ai développé, ce qui a permis de suivre l'activité neuronale pendant l'apprentissage de la peur. 233 neurones pyramidaux individuels pouvaient moduler leur taux d'activité en réponse au son conditionné (neurones conditionnés) et/ou au son non conditionné (neurones généralisant). Sur les 233 neurones pyramidaux putatifs enregistrés 86 (37%) d'entre eux ont été conditionnés à un ou deux tons. Plus précisément, une population de neurones code conjointement pour un souvenir partagé, alors qu'un groupe de neurones différent code pour de nouvelles caractéristiques de nouveaux souvenirs. En particulier, en dépit d'une extinction du comportement réussie, le taux de décharge de ces neurones conditionné en réponse à la tonalité conditionnée est resté inchangée tout au long de la mesure d'apprentissage. En outre, en analysant les caractéristiques de pré-conditionnement des neurones conditionnés, j'ai déterminé qu'il était possible de prévoir le recrutement neuronal basé sur trois facteurs : 1) la sensibilité initiale aux entrées auditives, avec les neurones sensibles aux sons étant plus facilement recrutés que les neurones ne répondant pas aux stimuli auditifs; 2) les niveaux d'excitabilité des neurones, avec les neurones plus facilement excitables étant plus susceptibles d'être conditionnés au son ; et 3) le nombre de connexions reçues, puisque les neurones conditionné reçoivent plus de connexions que les neurones non-conditionnés. Enfin, nous avons constaté qu'il était possible de remplacer de façon satisfaisante le SI lors d'un conditionnement à la peur par des injections bilatérales de bicuculline, un antagoniste des récepteurs de l'acide y-Aminobutirique.

Uncertainty about the true source. A note on the likelihood ratio at the activity level.

This paper focuses on likelihood ratio based evaluations of fibre evidence in cases in which there is uncertainty about whether or not the reference item available for analysis - that is, an item typically taken from the suspect or seized at his home - is the item actually worn at the time of the offence. A likelihood ratio approach is proposed that, for situations in which certain categorical assumptions can be made about additionally introduced parameters, converges to formula described in existing literature. The properties of the proposed likelihood ratio approach are analysed through sensitivity analyses and discussed with respect to possible argumentative implications that arise in practice.

Conditional deletion of ferritin h in mice reduces B and T lymphocyte populations.

The immune system and iron availability are intimately linked as appropriate iron supply is needed for cell proliferation, while excess iron, as observed in hemochromatosis, may reduce subsets of lymphocytes. We have tested the effects of a ferritin H gene deletion on lymphocytes. Mx-Cre mediated conditional deletion of ferritin H in bone marrow reduced the number of mature B cells and peripheral T cells in all lymphoid organs. FACS analysis showed an increase in the labile iron pool, enhanced reactive oxygen species formation and mitochondrial depolarization. The findings were confirmed by a B-cell specific deletion using Fth(lox/lox) ; CD19-Cre mice. Mature B cells were strongly under-represented in bone marrow and spleen of the deleted mice, whereas pre-B and immature B cells were not affected. Bone marrow B cells showed increased proliferation as judged by the number of cells in S and G2/M phase as well as BrdU incorporation. Upon in vitro culture with B-cell activating factor of the tumor necrosis factor family (BAFF), ferritin H-deleted spleen B cells showed lower survival rates than wild type cells. This was partially reversed with iron-chelator deferiprone. The loss of T cells was also confirmed by a T cell-specific deletion in Fth(lox/lox) ;CD4-Cre mice. Our data show that ferritin H is required for B and T cell survival by actively reducing the labile iron pool. They further suggest that natural B and T cell maturation is influenced by intracellular iron levels and possibly deregulated in iron excess or deprivation.

Construction and analysis of a new conditional ferritin H knockout mouse strain/

Résumé Le fer joue un rôle important dans la plupart des fonctions biologiques mais sa présence excessive provoque la production de molécules réactives d'oxygène (ROS) qui peuvent contribuer à diverses maladies. La protéine de stockage du fer, la ferritine H, capte l'excès en fer et le stocke sous forme non-toxique, ce qui empêche des dommages potentiels. La délétion de la ferritine H dans des souris knock-out a été essayée antérieurement, mais ces souris mouraient au stade précoce du développement embryonnaire. Pour étudier l'importance du fer, et en particulier son stockage dans la ferritine, et pour pouvoir mieux comprendre les fonctions de la ferritine H, nous avons créé un modèle de souris knock-out conditionnelles de la ferritine H, selon le système classique de Cre-LoxP. Le premier exon et la région du promoteur du gène de la ferritine H ont été entourés de sites loxP. La mortalité embryonnaire provoquée par la délétion constitutive du gène de la ferritine H a été confirmée en croisant nos souris avec des souris exprimant nestin-Cre1. En croisant nos souris avec des souris transgéniques Mx-Cre, nous avons observé que l'induction de Cre par injection de polyI-polyC provoque la délétion presque complète de la ferritine H dans le foie (> 99%) et la rate (> 88%). Ces tissus ont également perdu une grande partie de leur réserve de fer. Cette observation apporte pour la première fois la preuve in vivo que la ferritine H est indispensable pour le stockage du fer, que les fonctions de la ferritine H et de la ferritine L ne sont pas équivalentes, et que la ferritine L ne peut pas assumer seule la fonction de stockage du fer. Dans le foie des souris knock-out, l'expression de l'ARN messager de l'hepcidine a été induite après 10 jours. En même temps, l'expression de l'ARN messager des gènes codant pour des protéines de l'absorption de fer (DMT1, ferroportin, Dcytb1 et hephaestin) a été réprimée mais dans le duodénum seulement. L'expression d'hepcidine est inversément corrélée avec celle des gènes liés à l'absorption de fer. Cette observation corrobore des études antérieures. Mais, en plus, elle montre également que cette répression se produit seulement dans l'intestin. Nous pouvons ainsi tirer la conclusion suivante : ou bien l'hepcidine a un récepteur spécifique dans le duodénum ou bien les gènes liés à l'absorption de fer dans le duodénum ont un facteur spécifique de transcription sensible à l'hepcidine. Aucune répression de DMT1 et de ferroportin n'a été observée dans les macrophages de la rate après l'induction d'hepcidine. La délétion de ferritine H a entraîné une augmentation du taux de mortalité des cellules hépatiques, ainsi que des altérations dans l'architecture normale du tissu de la rate. Vu par l'immunohistologie, le nombre de lymphocytes B et T était réduit dans la rate, tendant à démontrer que la ferritine H et l'homéostase du fer jouent un rôle dans l'immunité. En conclusion, le modèle de souris knock-out conditionnelles de la ferritine H nous fournit un outil précieux pour l'étude in vivo du rôle joué par la ferritine dans l'homéostase du fer, dans les dommages créés par les ROS, ainsi que dans l'apoptose et l'immunité. Summary Iron plays an important role in most biological functions. However, excess of iron results in production of reactive oxygen species (ROS) which could substantially contribute to pathology of various diseases. Ferritin H scavenges excess of iron and stores it in non-toxic form and potentially prevents the damage. Fenitin H targeting in mice has been attempted before, however, straight knockout was lethal in early embryonic stage. To study the role of iron and its storage protein ferritin and to further elucidate ferritin H functions, we aimed at creating a conditional ferritin H knockout mouse model by classical Cre-LoxP system. First exon along with promoter region of the ferritin H gene was foxed. Embryonic lethality of the constitutive ferritin H deletion was confirmed by crossing the foxed mice with mice expressing nestin Cre-1 as transgene. Almost complete deletion was observed in liver (> 99%) and spleen (>88%) upon induction of Cre by injecting polyI-polyC in Fth Lox/Lox; MxCre mice. These tissues also lost substantial fraction of their iron stores. This provides first in vivo evidence that ferritin H is required for iron storage, ferritin H and L functions are not redundant and that ferritin L cannot perform iron storage function alone. Hepcidin mRNA expression was induced after 10 days in the livers of deleted mice and, simultaneously, mRNA expression of iron absorption related genes (DMT 1, ferroportin, Dcytb1 and hephaestin) was repressed in duodenum only. Hepcidin expression is inversely correlated with that of duodenal iron absorption related genes. This is in agreement with previous studies. However, we also show that this repression happens only in intestine. This leads to the conclusion that either hepcidin has a specific receptor in duodenum or the iron absorption related genes have duodenum specific transcription factor that is responsive to hepcidin. No repression of DMT1 and ferroportin was observed in spleen macrophages upon hepcidin induction. Ferritin H deletion showed increased cell death in liver and disruption of normal architecture of spleen. B lymphocytes were reduced in spleen on immunohistology which point towards a role of ferritin H and iron homeostasis in immunity. In conclusion, ferritin H conditional knockout mouse model provides us with an invaluable tool to study the in vivo role of ferritin H in iron homeostasis, ROS mediated damage, apoptosis and immunity.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.