Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Non-equilibrium energy and momentum accommodation coefficients of Ar atoms scattered from Ni(001) in the thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study energy and momentum transfer of low-energy Ar atoms scattered from the Ni(001) surface. The investigation concentrates on the dependence of these processes on incident energy, angles of incidence and surface temperature. Energy transfer exhibits a strong dependence on the surface temperature, at incident energies below 500 meV, and incident angles close to specular incidence. Above 500 meV, the surface temperature dependence vanishes, and a limiting value in the amount of energy transferred to the surface is attained. Momentum exchange is investigated in terms of tangential and normal components. Both components exhibit a weak surface temperature dependence, but they have opposite behaviours at all incidence angles. In each component, momentum can be lost or gained following the interaction with the surface. (C) 1997 Elsevier Science B.V.

Modelling HIV/AIDS length of stay in Portugal

Tese apresentada como requisito parcial para obtenção do grau de Doutor em Estatística e Gestão de Informação pelo Instituto Superior de Estatística e Gestão de Informação da Universidade Nova de Lisboa

Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.

Litoestratigrafia e deformação cenzóica na região de Icapuí, Ceará, e implicações para a estruturação de campos de petróleo na borda ocidental da bacia Potiguar(NE do Brasil)

This thesis deals with the sedimentological/stratigraphic and structural evolution of the sedimentary rocks that occur in the NW continental border of the Potiguar Basin. These rocks are well exposed along coastal cliffs between the localities of Lagoa do Mato and Icapuí, Ceará State (NE Brazil). The sedimentological/stratigraphic study involved, at the outcrop scale, detailed facies descriptions, profile mapping of the vertical succession of different beds, and columnar sections displaying inferred lateral relationships. The approach was complemented by granulometric and petrographic analyses, including the characterization of heavy mineral assemblages. The data set allowed to recognize two kinds of lithological units, a carbonate one of very restricted occurrence at the base of the cliffs, and three younger, distinct siliciclastic units, that predominate along the cliffs, in vertical and lateral extent. The carbonate rocks were correlated to the late Cretaceous Jandaíra Formation, which is covered by the siliciclastic Barreiras Formation. The Barreiras Formation occurs in two distinct structural settings, the usual one with nondeformed, subhorizontal strata, or as tilted beds, affected by strong deformation. Two lithofacies were recognized, vertically arranged or in fault contacts. The lower facies is characterized by silty-argillaceous sandstones with low-angle cross bedding; the upper facies comprises medium to coarse grained sandstones, with conglomeratic layers. The Tibau Formation (medium to coarse-grained sandstones with argillite intercalations) occurs at the NW side of the studied area, laterally interlayered with the Barreiras Formation. Eolic sediments correlated to the Potengi Formation overly the former units, either displaying an angular unconformity, or simply an erosional contact (stratigraphic unconformity). Outstanding structural features, identified in the Barreiras Formation, led to characterize a neocenozoic stress field, which generated faults and folds and/or reactivated older structures in the subjacent late cretaceous (to paleogene, in the offshore basin) section. The structures recognized in the Barreiras Formation comprise two distinct assemblages, namely a main extensional deformation between the localities of Ponta Grossa and Redonda, and a contractional style (succeeded by oblique extensional structures) at Vila Nova. In the first case, the structural assemblage is dominated by N-S (N±20°Az) steep to gently-dipping extensional faults, displaying a domino-style or listric geometry with associated roll-over structures. This deformation pattern is explained by an E-W/WNW extension, contemporaneous with deposition of the upper facies of the Barreiras Formation, during the time interval Miocene to Pleistocene. Strong rotation of blocks and faults generated low-angle distensional faults and, locally, subvertical bedding, allowing to estimate very high strain states, with extension estimates varying between 40% up to 200%. Numerous detachment zones, parallel to bedding, help to acommodate this intense deformation. The detachment surfaces and a large number of faults display mesoscopic features analoguous to the ones of ductile shear zones, with development of S-C fabrics, shear bands, sigmoidal clasts and others, pointing to a hydroplastic deformation regime in these cases. Local occurrences of the Jandaíra limestone are controled by extensional faults that exhume the pre-Barreiras section, including an earlier event with N-S extension. Finally, WNWtrending extensional shear zones and faults are compatible with the Holocene stress field along the present continental margin. In the Vila Nova region, close to Icapuí, gentle normal folds with fold hinges shallowly pluging to SSW affect the lower facies of the Barreiras Formation, displaying an incipient dissolution cleavage associated with an extension lineation at high rake (a S>L fabric). Deposition of the upper facies siliciclastics is controlled by pull-apart graben structures, bordered by N-NE-trending sinistral-normal shear zones and faults, characterizing an structural inversion. Microstructures are compatible with tectonic deformation of the sedimentary pile, burried at shallow depths. The observed features point to high pore fluid pressures during deformation of the sediments, producing hydroplastic structures through mechanisms of granular flow. Such structures are overprinted by microfractures and microfaults (an essentially brittle regime), tracking the change to microfracturing and frictional shear mechanisms accompanying progressive dewatering and sediment lithification. Correlation of the structures observed at the surface with those present at depth was tested through geophysical data (Ground Penetrating Radar, seismics and a magnetic map). EW and NE-trending lineaments are observed in the magnetic map. The seismic sections display several examples of positive flower structures which affect the base of the cretaceous sediments; at higher stratigraphic levels, normal components/slips are compatible with the negative structural inversion characterized at the surface. Such correlations assisted in proposing a structural model compatible with the regional tectonic framework. The strong neogenepleistocene deformation is necessarily propagated in the subsurface, affecting the late cretaceous section (Açu and Jandaíra formations), wich host the hydrocarbon reservoirs in this portion of the Potiguar Basin. The proposed structural model is related to the dextral transcurrent/transform deformation along the Equatorial Margin, associated with transpressive terminations of E-W fault zones, or at their intersections with NE-trending lineaments, such as the Ponta Grossa-Fazenda Belém one (the LPGFB, itself controlled by a Brasiliano-age strike-slip shear zone). In a first step (and possibly during the late Cretaceous to Paleogene), this lineament was activated under a sinistral transpressional regime (antithetic to the main dextral deformation in the E-W zones), giving way to the folds in the lower facies of the Barreiras Formation, as well as the positive flower structures mapped through the seismic sections, at depth. This stage was succeeded (or was penecontemporaneous) by the extensional structures related to a (also sinistral) transtensional movement stage, associated to volcanism (Macau, Messejana) and thermal doming processes during the Neogene-Pleistocene time interval. This structural model has direct implications to hydrocarbon exploration and exploitation activities at this sector of the Potiguar Basin and its offshore continuation. The structure of the reservoirs at depth (Açu Formation sandstones of the post-rift section) may be controlled (or at least, strongly influenced) by the deformation geometry and kinematics characterized at the surface. In addition, the deformation event recognized in the Barreiras Formation has an age close to the one postulated for the oil maturation and migration in the basin, between the Oligocene to the Miocene. In this way, the described structural cenario represents a valid model to understand the conditions of hydrocarbon transport and acummulation through space openings, trap formation and destruction. This model is potentially applicable to the NW region of the Potiguar Basin and other sectors with a similar structural setting, along the brazilian Equatorial Atlantic Margin

Tensor3D: A computer graphics program to simulate 3D real-time deformation and visualization of geometric bodies

Tensor3D is a geometric modeling program with the capacity to simulate and visualize in real-time the deformation, specified through a tensor matrix and applied to triangulated models representing geological bodies. 3D visualization allows the study of deformational processes that are traditionally conducted in 2D, such as simple and pure shears. Besides geometric objects that are immediately available in the program window, the program can read other models from disk, thus being able to import objects created with different open-source or proprietary programs. A strain ellipsoid and a bounding box are simultaneously shown and instantly deformed with the main object. The principal axes of strain are visualized as well to provide graphical information about the orientation of the tensor's normal components. The deformed models can also be saved, retrieved later and deformed again, in order to study different steps of progressive strain, or to make this data available to other programs. The shape of stress ellipsoids and the corresponding Mohr circles defined by any stress tensor can also be represented. The application was written using the Visualization ToolKit, a powerful scientific visualization library in the public domain. This development choice, allied to the use of the Tcl/Tk programming language, which is independent on the host computational platform, makes the program a useful tool for the study of geometric deformations directly in three dimensions in teaching as well as research activities. (C) 2007 Elsevier Ltd. All rights reserved.

Geochemistry and minerals at DSDP Leg 81 Holes

The Leg 81 basalts, drilled either from the margins ("dipping reflectors" sequence: Holes 552, 553A, and 554A) or from the "continental" side (Hole 555) of the Rockall Plateau microcontinent, are strongly light rare-earth element (LREE) depleted oceanic tholeiites. The basalts from the four holes are almost similar. Most of their primary characteristics have been preserved, although they have suffered alteration by seawater. From the petrological and mineralogical points of view, they resemble deep-ocean-floor basalts but show some peculiarities (occurrence of pigeonite and ilmenite as normal components of the groundmass differentiation sequences toward ferrobasalts). Their geochemical characteristics are dominated by their extreme depletion in the most hygromagmaphile elements (Th, Ta, La, and Nb), the concentrations of which are sometimes lower than the corresponding chondritic values. Leg 81 basalts are thus clearly different from continental tholeiites (flood basalts): Possible equivalents in the Thulean Tertiary Magmatic Province include the LREE-depleted tholeiites from the Upper Basaltic Series of the Faeroe Islands and the Preshal Mhor basalt type from the British Tertiary Province.

Mixture modelling for cluster analysis

Cluster analysis via a finite mixture model approach is considered. With this approach to clustering, the data can be partitioned into a specified number of clusters g by first fitting a mixture model with g components. An outright clustering of the data is then obtained by assigning an observation to the component to which it has the highest estimated posterior probability of belonging; that is, the ith cluster consists of those observations assigned to the ith component (i = 1,..., g). The focus is on the use of mixtures of normal components for the cluster analysis of data that can be regarded as being continuous. But attention is also given to the case of mixed data, where the observations consist of both continuous and discrete variables.

Short-Term Effects of Duodenocecostomy on Body Weight, Glucose Absorption, Serum Components, and Intestinal Histopathology in Four Normal Horses

The objective of this study was to evaluate duodenocecostomy in horses performed through a ventral midline laparotomy and report its influence oil body weight, glucose absorption, serum components, and characteristics of jejunum, cecum, and large colon histology. Four horses were submitted to the duodenocecostomy technique through a ventral midline laparotomy with animals in dorsal recumbency under inhalation anesthesia, followed by abdominal exploration. A side-to-side anastomosis was performed between the duodenojejunal flexure and the base of the cecum with two simple continuous suture lines of the serosal and muscular layers. The size of the opening created was approximately 2 cm in diameter. The mucosa layer was not Sutured. After 30 days, animals were submitted to a second laparotomy to check the patency of the duodenocaecal fistula. During both laparotomy procedures, excisional biopsies of different segments of the gastrointestinal tract were performed. Information on physical examination findings, results of hematologic and histopathologic evaluations, and oral glucose absorption test were recorded. The horses did not have significant weight loss from baseline, and absorption curve of glucose did not significantly vary from baseline. Only triglycerides had significant alterations. Histologic evaluation of jejunum, cecum, and large colon did not show alterations of intestinal structure and morphology. We concluded that the proposed technique, principally in relation to the fistula size and not suturing the mucosa layer, allowed partial or total Occlusion of the fistulae without the necessity of a second surgery and avoided the permanent bypass of ingesta and weight loss.

Identification of Trypanosoma brucei components involved in trypanolysis by normal human serum

Normal human serum (NHS) confers human resistance to infection by the parasite Trypanosoma brucei owing to the trypanolytic activity of apolipoprotein L1 (APOL1), present in two serum complexes termed Trypanolytic Factors (TLF-1 and -2). In order to identify parasite components involved in the intracellular trafficking and activity of TLFs, an inducible RNA interference (RNAi) genomic DNA library constructed in bloodstream form T. brucei was subjected to RNAi induction and selection for resistant parasites under NHS conditions favouring either TLF-1 or TLF-2 uptake. While TLF-1 conditions readily selected the haptoglobin-haemoglobin (HP-HB) surface receptor TbHpHbR as expected, given its known ability to bind TLF-1, under TLF-2 conditions no specific receptor for TLF-2 was identified. Instead, the screen allowed the identification of five distinct factors expected to be involved in the assembly of the vacuolar proton pump V-ATPase and consecutive endosomal acidification. These data confirm that lowering the pH during endocytosis is required for APOL1 toxic activity.

β-amyloid (Aβ) deposition in cognitively normal brain, dementia with Lewy bodies, and Alzheimer's disease:a study using principal components analysis

The densities of diffuse, primitive, and classic ß-amyloid (Aß) deposits were studied in the temporal lobe in cognitively normal brain, dementia with Lewy bodies (DLB), familial Alzheimer’s disease (FAD), and sporadic AD (SAD). Principal components analysis (PCA) was used to determine whether there were distinct differences between groups or whether Aß pathology was more continuously distributed from group to group. Three principal components (PC) were extracted from the data accounting for 56% of the total variance. Plots of cases in relation to the PC did not result in distinct groups but suggested overlap in Aß deposition between the groups. In addition, there were linear correlations between the densities of Aß deposits and the distribution of the cases along the PC in specific brain regions suggesting continuous variation from group to group. PC1 was associated with the degree of maturation of Aß deposits, PC2 with differences between FAD and SAD, and PC3 with the degree of spread of Aß pathology into the hippocampus. Apolipoprotein E (APOE) genotype was not associated with variation in Aß deposition between cases. PCA may be a useful method of studying the pathological interface between closely related neurodegenerative disorders.

Key Participants Of The Tumor Microenvironment Of The Prostate: An Approach Of The Structural Dynamic Of Cellular Elements And Extracellular Matrix Components During Epithelial-stromal Transition.

Cancer is a multistep process that begins with the transformation of normal epithelial cells and continues with tumor growth, stromal invasion and metastasis. The remodeling of the peritumoral environment is decisive for the onset of tumor invasiveness. This event is dependent on epithelial-stromal interactions, degradation of extracellular matrix components and reorganization of fibrillar components. Our research group has studied in a new proposed rodent model the participation of cellular and molecular components in the prostate microenvironment that contributes to cancer progression. Our group adopted the gerbil Meriones unguiculatus as an alternative experimental model for prostate cancer study. This model has presented significant responses to hormonal treatments and to development of spontaneous and induced neoplasias. The data obtained indicate reorganization of type I collagen fibers and reticular fibers, synthesis of new components such as tenascin and proteoglycans, degradation of basement membrane components and elastic fibers and increased expression of metalloproteinases. Fibroblasts that border the region, apparently participate in the stromal reaction. The roles of each of these events, as well as some signaling molecules, participants of neoplastic progression and factors that promote genetic reprogramming during epithelial-stromal transition are also discussed.

Gene expression relationship between prostate cancer cells of Gleason 3, 4 and normal epithelial cells as revealed by cell type-specific transcriptomes

Background: Prostate cancer cells in primary tumors have been typed CD10(-)/CD13(-)/CD24(hi)/CD26(+)/CD38(lo)/CD44(-)/CD104(-). This CD phenotype suggests a lineage relationship between cancer cells and luminal cells. The Gleason grade of tumors is a descriptive of tumor glandular differentiation. Higher Gleason scores are associated with treatment failure. Methods: CD26(+) cancer cells were isolated from Gleason 3+3 (G3) and Gleason 4+4 (G4) tumors by cell sorting, and their gene expression or transcriptome was determined by Affymetrix DNA array analysis. Dataset analysis was used to determine gene expression similarities and differences between G3 and G4 as well as to prostate cancer cell lines and histologically normal prostate luminal cells. Results: The G3 and G4 transcriptomes were compared to those of prostatic cell types of non-cancer, which included luminal, basal, stromal fibromuscular, and endothelial. A principal components analysis of the various transcriptome datasets indicated a closer relationship between luminal and G3 than luminal and G4. Dataset comparison also showed that the cancer transcriptomes differed substantially from those of prostate cancer cell lines. Conclusions: Genes differentially expressed in cancer are potential biomarkers for cancer detection, and those differentially expressed between G3 and G4 are potential biomarkers for disease stratification given that G4 cancer is associated with poor outcomes. Differentially expressed genes likely contribute to the prostate cancer phenotype and constitute the signatures of these particular cancer cell types.

Fluorescence and reflectance spectroscopy for identification of atherosclerosis in human carotid arteries using principal components analysis

Objectives: The aim of this work was to verify the differentiation between normal and pathological human carotid artery tissues by using fluorescence and reflectance spectroscopy in the 400- to 700-nm range and the spectral characterization by means of principal components analysis. Background Data: Atherosclerosis is the most common and serious pathology of the cardiovascular system. Principal components represent the main spectral characteristics that occur within the spectral data and could be used for tissue classification. Materials and Methods: Sixty postmortem carotid artery fragments (26 non-atherosclerotic and 34 atherosclerotic with non-calcified plaques) were studied. The excitation radiation consisted of a 488-nm argon laser. Two 600-mu m core optical fibers were used, one for excitation and one to collect the fluorescence radiation from the samples. The reflectance system was composed of a halogen lamp coupled to an excitation fiber positioned in one of the ports of an integrating sphere that delivered 5 mW to the sample. The photo-reflectance signal was coupled to a 1/4-m spectrograph via an optical fiber. Euclidean distance was then used to classify each principal component score into one of two classes, normal and atherosclerotic tissue, for both fluorescence and reflectance. Results: The principal components analysis allowed classification of the samples with 81% sensitivity and 88% specificity for fluorescence, and 81% sensitivity and 91% specificity for reflectance. Conclusions: Our results showed that principal components analysis could be applied to differentiate between normal and atherosclerotic tissue with high sensitivity and specificity.