A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

A Kalman filter based strategy for linear structural system identification based on multiple static and dynamic test data

The problem of identification of stiffness, mass and damping properties of linear structural systems, based on multiple sets of measurement data originating from static and dynamic tests is considered. A strategy, within the framework of Kalman filter based dynamic state estimation, is proposed to tackle this problem. The static tests consists of measurement of response of the structure to slowly moving loads, and to static loads whose magnitude are varied incrementally; the dynamic tests involve measurement of a few elements of the frequency response function (FRF) matrix. These measurements are taken to be contaminated by additive Gaussian noise. An artificial independent variable τ, that simultaneously parameterizes the point of application of the moving load, the magnitude of the incrementally varied static load and the driving frequency in the FRFs, is introduced. The state vector is taken to consist of system parameters to be identified. The fact that these parameters are independent of the variable τ is taken to constitute the set of ‘process’ equations. The measurement equations are derived based on the mechanics of the problem and, quantities, such as displacements and/or strains, are taken to be measured. A recursive algorithm that employs a linearization strategy based on Neumann’s expansion of structural static and dynamic stiffness matrices, and, which provides posterior estimates of the mean and covariance of the unknown system parameters, is developed. The satisfactory performance of the proposed approach is illustrated by considering the problem of the identification of the dynamic properties of an inhomogeneous beam and the axial rigidities of members of a truss structure.

Structural Relationships in 2,3-Bis-n-decyloxyanthracene and 12-Hydroxystearic Acid Molecular Gels and Aerogels Processed in Supercritical CO2

Supercritical carbon dioxide is used to prepare aerogels of two reference molecular organogelators, 2,3-bis-n-decyloxyanthracene (DDOA) (luminescent molecule) and 12-hydroxystearic acid (HSA). Electron microscopy reveals the fibrillar morphology of the aggregates generated by the protocol. SAXS and SANS measurements show that DDOA aerogels are crystalline materials exhibiting three morphs: (1) arrangements of the crystalline solid (2D p6m), (2) a second hexagonal morph slightly more compact, and (3) a packing specific of the fibers in the gel. Aggregates specific of the aerogel (volume fraction being typically phi approximate to 0.60) are developed over larger distances (similar to 1000 angstrom) and bear fewer defaults and residual strains than aggregates in the crystalline and gel phases. Porod, Scherrer and Debye-Bueche analyses of the scattering data have been performed. The first five diffraction peaks show small variations in position and intensity assigned to the variation of the number of fibers and their degree of vicinity within hexagonal bundles of the related SAFIN according to the Oster model. Conclusions are supported by the guidelines offered by the analysis of the situation in HSA aerogels for which the diffraction pattern can be described by two coexisting lamellar-like arrangements. The porosity of the aerogel, as measured by its specific surface extracted from the scattering invariant analysis, is only 1.8 times less than that of the swollen gel and is characteristic of a very porous material.

Morphology and localization of interstitial cells in the guinea pig bladder: structural relationships with smooth muscle and neurons.

PURPOSE: In the current study we examined the location of interstitial cell of Cajal (ICC)-like cells in the guinea pig bladder wall and studied their structural interactions with nerves and smooth muscle cells. MATERIALS AND METHODS: Whole mount samples and cryosections of bladder tissue were labeled with primary and fluorescent secondary antibodies, and imaged using confocal and multiphoton microscopy. RESULTS: Kit positive ICC-like cells were located below the urothelium, in the lamina propria region and throughout the detrusor. In the suburothelium they had a stellate morphology and appeared to network. They made connections with nerves, as shown by double labeling experiments with anti-kit and anti-protein gene product 9.5. A network of vimentin positive cells was also found, of which many but not all were kit positive. In the detrusor kit positive cells were most often seen at the edge of smooth muscle bundles. They were elongated with lateral branches, running in parallel with the bundles and closely associated with intramural nerves. Another population of kit positive cells was seen in the detrusor between muscle bundles. These cells had a more stellate-like morphology and made connections with each other. Kit positive cells were seen tracking nerve bundles and close to intramural ganglia. Vimentin positive cells were present in the detrusor, of which some were also kit positive. CONCLUSIONS: There are several populations of ICC-like cells throughout the guinea pig bladder wall. They differ in morphology and orientation but all make connections with intramural nerves and in the detrusor they are closely associated with smooth muscle cells.

Kit-positive interstitial cells in the rabbit urethra: structural relationships with nerves and smooth muscle.

OBJECTIVE: To identify interstitial cells (ICs) in the wall of the rabbit urethra using antibodies to the Kit receptor, and to examine their location, morphology and relationship with nerves and smooth muscle cells (SMCs), as studies of enzymatically isolated cells from the rabbit urethra have established that there are specialized cells that show spontaneous electrical activity and have morphological properties of ICs. MATERIALS AND METHODS: Urethral tissues from rabbits were fixed, labelled with antibodies and examined with confocal microscopy. Some specimens were embedded in paraffin wax and processed for histology. Histological sections from the most proximal third and mid-third region of rabbit urethra were stained with Masson's Trichrome to show their cellular arrangement. RESULTS: Sections from both regions had outer longitudinal and inner circular layers of SM, and a lamina propria containing connective tissue and blood vessels; the lumen was lined with urothelial cells. The mid-third region had a more developed circular SM layer than the most-proximal samples, and had extensive inner longitudinal SM bundles in the lamina propria. Labelling with anti-Kit revealed immunopositive cells within the wall of the rabbit urethra, in the circular and longitudinal layers of the muscularis. Double-labelling with an antibody to SM myosin showed Kit-positive cells on the boundary of the SM bundles, orientated parallel to the axis of the bundles. Others were in spaces between the bundles and often made contact with each other. Kit-positive cells were either elongated, with several lateral branches, or stellate with branches coming from a central soma. Similar cells could be labelled with vimentin antibodies. Their relationship with intramural nerves was examined by double immunostaining with an anti-neurofilament antibody. There were frequent points of contact between Kit-positive cells and nerves, with similar findings in specimens double-immunostained with anti-neuronal nitric oxide synthase (nNOS). CONCLUSION: Kit-positive ICs were found within the SM layers of the rabbit urethra, in association with nerves, on the edge of SM bundles and in the interbundle spaces. The contact with nNOS-containing neurones might imply participation in the nitrergic inhibitory neurotransmission of the urethra. PMID: 17212607 [PubMed - indexed for MEDLINE]

Projection-Based Statistical Inference in Linear Structural Models with Possibly Weak Instruments

It is well known that standard asymptotic theory is not valid or is extremely unreliable in models with identification problems or weak instruments [Dufour (1997, Econometrica), Staiger and Stock (1997, Econometrica), Wang and Zivot (1998, Econometrica), Stock and Wright (2000, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. One possible way out consists here in using a variant of the Anderson-Rubin (1949, Ann. Math. Stat.) procedure. The latter, however, allows one to build exact tests and confidence sets only for the full vector of the coefficients of the endogenous explanatory variables in a structural equation, which in general does not allow for individual coefficients. This problem may in principle be overcome by using projection techniques [Dufour (1997, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. AR-types are emphasized because they are robust to both weak instruments and instrument exclusion. However, these techniques can be implemented only by using costly numerical techniques. In this paper, we provide a complete analytic solution to the problem of building projection-based confidence sets from Anderson-Rubin-type confidence sets. The latter involves the geometric properties of “quadrics” and can be viewed as an extension of usual confidence intervals and ellipsoids. Only least squares techniques are required for building the confidence intervals. We also study by simulation how “conservative” projection-based confidence sets are. Finally, we illustrate the methods proposed by applying them to three different examples: the relationship between trade and growth in a cross-section of countries, returns to education, and a study of production functions in the U.S. economy.

Ternary erbium chromium sulfides : structural relationships and magnetic properties

Single crystals of four erbium-chromium sulfides have been grown by chemical vapor transport using iodine as the transporting agent. Single-crystal X-ray diffraction reveals that in Er(3)CrS(6) octahedral sites are occupied exclusively by Cr(3+) cations, leading to one-dimensional CrS(4)(5-) chains of edge-sharing octahedra, while in Er(2)CrS(4), Er(3+), and Cr(2+) cations occupy the available octahedral sites in an ordered manner. By contrast, in Er(6)Cr(2)S(11) and Er(4)CrS(7), Er(3+) and Cr(2+) ions are disordered over the octahedral sites. In Er(2)CrS(4), Er(6)Cr(2)S(11), and Er(4)CrS(7), the network of octahedra generates an anionic framework constructed from M(2)S(5) slabs of varying thickness, linked by one-dimensional octahedral chains. This suggests that these three phases belong to a series in which the anionic framework may be described by the general formula [M(2n+1)S(4n+3)](x-), with charge balancing provided by Er(3+) cations located in sites of high-coordination number within one-dimensional channels defined by the framework. Er(4)CrS(7), Er(6)Cr(2)S(11), and Er(2)CrS(4) may thus be considered as the n = 1, 2, and infinity members of this series. While Er(4)CrS(7) is paramagnetic, successive magnetic transitions associated with ordering of the chromium and erbium sub-lattices are observed on cooling Er(3)CrS(6) (T(C)(Cr) = 30 K; T(C)(Er) = 11 K) and Er(2)CrS(4) (T(N)(Cr) = 42 K, T(N)(Er) = 10 K) whereas Er(6)Cr(2)S(11) exhibits ordering of the chromium sub-lattice only (T(N) = 11.4 K).

Structural relationships in the solid state of the anti-chagas agent (E)-phenylethenylbenzofuroxan

The crystal structure and the vibrational spectrum of a potential drug for Chagas`s disease treatment, the (E)-isomer of phenylethenylbenzofuroxan 1 (5(6)(E)-[(2-phenylethenyl)]benzo[1,2-c]1,2,5-oxadiazole N-oxide), are reported. In order to provide insights into structural relationships, quantum mechanical calculations were employed starting from crystal structure. These results have given theoretical support to state interesting structural features, such as the effect of some intermolecular contacts on the molecule conformation and the electronic delocalization decreasing through atoms of the benzofuroxan moiety. Furthermore, the MOGUL comparative analysis in the Cambridge Structural Database provided additional evidences on these structural behaviors of compound 1. Intermolecular contacts interfere on the intramolecular geometry, as, for instance, on the phenyl group orientation, which is twisted by 12.32(6)A degrees from the ethenylbenzofuroxan plane. The experimental Raman spectrum of compound 1 presents unexpected frequency shift and also anomalous Raman activities. At last, the molecule skeleton deformation and the characteristic vibrational modes were correlated by matching the experimental Raman spectrum to the calculated one.

Mineral chemistry and structural relationships of inclusions in diamond samples

Diamant ist das härteste Mineral – und dazu ein Edelstein -, das unter höchstem Druck und hohen Temperaturen in tiefen kontinentalen Regionen der Erde kristallisiert. Die Mineraleinschlüsse in Diamanten werden durch die physikalische Stabilität und chemische Beständigkeit der umgebenden – eigentlich metastabilen -Diamant-Phase geschützt. Aufgrund der koexistierenden Phasenkombination ermöglichen sie, die Mineral-Entwicklung zu studieren, während deren der Einschlüssen und die Diamanten kristallisierten. rnDie Phasenkombinationen von Diamant und Chrom-Pyrop, Chrom-Diopsid, Chromit, Olivin, Graphit und Enstatit nebeneinander (teilweise in Berührungsexistenz) mit Chrom-Pyrop Einschlüssen wurden von neunundzwanzig Diamant-Proben von sechs Standorten in Südafrika (Premier, Koffiefontein, De Beers Pool, Finsch, Venetia und Koingnaas Minen) und Udachnaya (Sibirien/Russland) identifiziert und charakterisiert. Die Mineraleinschlüsse weisen z.T. kubo-oktaedrische Form auf, die unabhängig von ihren eigenen Kristallsystemen ausgebildet werden können. Das bedeutet, dass sie syngenetische Einschlüsse sind, die durch die sehr hohe Formenergie des umgebenden Diamanten morphologisch unter Zwang stehen. Aus zweidiemnsionalen Messungen der ersten Ordnung von charakteristischen Raman-Banden lassen sich relative Restdrucke in Diamanten zwischen Diamant und Einschlussmineral gewinnen; sie haben charakteristische Werte von ca. 0,4 bis 0,9 GPa um Chrom-Pyrop-Einschlüsse, 0,6 bis 2,0 GPa um Chrom-Diopsid-Einschlüsse, 0,3 bis 1,2 GPa um Olivin-Einschlüsse, 0,2 bis 1,0 GPa um Chromit-Einschlüsse, beziehungsweise 0,5 GPa um Graphit Einschlüsse.rnDie kristallstrukturellen Beziehung von Diamanten und ihren monomineralischen Einschlüssen wurden mit Hilfe der Quantifizierung der Winkelkorrelationen zwischen der [111] Richtung von Diamanten und spezifisch ausgewählten Richtungen ihrer mineralischen Einschlüsse untersucht. Die Winkelkorrelationen zwischen Diamant [111] und Chrom-Pyrop [111] oder Chromit [111] zeigen die kleinsten Verzerrungen von 2,2 bis zu 3,4. Die Chrom-Diopsid- und Olivin-Einschlüsse zeigen die Missorientierungswerte mit Diamant [111] bis zu 10,2 und 12,9 von Chrom-Diopsid [010] beziehungsweise Olivin [100].rnDie chemische Zusammensetzung von neun herausgearbeiteten (orientiertes Anschleifen) Einschlüssen (drei Chrom-Pyrop-Einschlüsse von Koffiefontein-, Finsch- und Venetia-Mine (zwei von drei koexistieren nebeneinander mit Enstatit), ein Chromit von Udachnaya (Sibirien/Russland), drei Chrom-Diopside von Koffiefontein, Koingnaas und Udachnaya (Sibirien/Russland) und zwei Olivin Einschlüsse von De Beers Pool und Koingnaas) wurden mit Hilfe EPMA und LA-ICP-MS analysiert. Auf der Grundlage der chemischen Zusammensetzung können die Mineraleinschlüsse in Diamanten in dieser Arbeit der peridotitischen Suite zugeordnet werden.rnDie Geothermobarometrie-Untersuchungen waren aufgrund der berührenden Koexistenz von Chrom-Pyrop- und Enstatit in einzelnen Diamanten möglich. Durchschnittliche Temperaturen und Drücke der Bildung sind mit ca. 1087 (± 15) C, 5,2 (± 0,1) GPa für Diamant DHK6.2 von der Koffiefontein Mine beziehungsweise ca. 1041 (± 5) C, 5,0 (± 0,1) GPa für Diamant DHF10.2 von der Finsch Mine zu interpretieren.rn

The proteolytic fragments generated by vertebrate proteasomes: structural relationships to major histocompatibility complex class I binding peptides.

Proteasomes are involved in the proteolytic generation of major histocompatibility complex (MHC) class I epitopes but their exact role has not been elucidated. We used highly purified murine 20S proteasomes for digestion of synthetic 22-mer and 41/44-mer ovalbumin partial sequences encompassing either an immunodominant or a marginally immunogenic epitope. At various times, digests were analyzed by pool sequencing and by semiquantitative electrospray ionization mass spectrometry. Most dual cleavage fragments derived from 22-mer peptides were 7-10 amino acids long, with octa- and nonamers predominating. Digestion of 41/44-mer peptides initially revealed major cleavage sites spaced by two size ranges, 8 or 9 amino acids and 14 or 15 amino acids, followed by further degradation of the latter as well as of larger single cleavage fragments. The final size distribution was slightly broader than that of fragments derived from 22-mer peptides. The majority of peptide bonds were cleaved, albeit with vastly different efficiencies. This resulted in multiple overlapping proteolytic fragments including a limited number of abundant peptides. The immunodominant epitope was generated abundantly whereas only small amounts of the marginally immunogenic epitope were detected. The frequency distributions of amino acids flanking proteasomal cleavage sites are correlated to that reported for corresponding positions of MHC class I binding peptides. The results suggest that proteasomal degradation products may include fragments with structural properties similar to MHC class I binding peptides. Proteasomes may thus be involved in the final stages of proteolytic epitope generation, often without the need for downstream proteolytic events.

Mimotope/anti-mimotope probing of structural relationships in platelet glycoprotein Ib alpha.

A bacteriophage library displaying random decapeptides was used to characterize the binding preference of C-34, a monoclonal antibody originally raised against platelet-type von Willebrand disease platelets heterozygous for the mutation 23OWKQ (G --> V)233V234 in the alpha chain of glycoprotein Ib (GPIb alpha). Three rounds of biopanning C-34 against the library resulted in striking convergence upon the sequence WNWRYREYV. Since no portion of this sequence corresponds to a recognizable peptide sequence within human platelet GPIb alpha, it may be considered a "mimotope" of the naturally occurring C-34 epitope, presumably bearing similarity to it in three-dimensional structure. Synthetic AWNWRYREYV peptide preincubated with C-34 fully neutralized the ability of C-34 to inhibit platelet aggregation, with an IC50 of approximately 6 microg/ml. When biotinylated AWNWRYREYV was subsequently bioparmed against the original decapeptide library, the sole clone demonstrating inhibitory activity above background level in a functional platelet assay displayed the sequence RHVAWWRQGV, and chemically synthesized peptide fully inhibited ristocetin-induced aggregation, with an IC50 of 200-400 microg/ml. Synthesized RHVAWWKQGV peptide exerted only slight inhibition, whereas RHVAWWKQVV peptide showed potent inhibitory activity. Moreover, whereas synthesized wild-type 228YVWKQGVDVK237 GPIb alpha peptide was virtually without inhibitory activity, the 228YVWKQ(G -->V) 233VDVK237 peptide fully inhibited ristocetin-induced aggregation, with an IC50 of approximately 400 microg/ml. These studies raise the possibility of an intramolecular association of peptide regions within GPIb alpha that may play a role in the regulation of von Willebrand factor-dependent platelet aggregation.

Inference for a skew extension of the Grubbs model

In this paper, we discuss inferential aspects for the Grubbs model when the unknown quantity x (latent response) follows a skew-normal distribution, extending early results given in Arellano-Valle et al. (J Multivar Anal 96:265-281, 2005b). Maximum likelihood parameter estimates are computed via the EM-algorithm. Wald and likelihood ratio type statistics are used for hypothesis testing and we explain the apparent failure of the Wald statistics in detecting skewness via the profile likelihood function. The results and methods developed in this paper are illustrated with a numerical example.

Multivariate measurement error models based on scale mixtures of the skew-normal distribution

Scale mixtures of the skew-normal (SMSN) distribution is a class of asymmetric thick-tailed distributions that includes the skew-normal (SN) distribution as a special case. The main advantage of these classes of distributions is that they are easy to simulate and have a nice hierarchical representation facilitating easy implementation of the expectation-maximization algorithm for the maximum-likelihood estimation. In this paper, we assume an SMSN distribution for the unobserved value of the covariates and a symmetric scale mixtures of the normal distribution for the error term of the model. This provides a robust alternative to parameter estimation in multivariate measurement error models. Specific distributions examined include univariate and multivariate versions of the SN, skew-t, skew-slash and skew-contaminated normal distributions. The results and methods are applied to a real data set.

On estimation and influence diagnostics for the Grubbs` model under heavy-tailed distributions

The Grubbs` measurement model is frequently used to compare several measuring devices. It is common to assume that the random terms have a normal distribution. However, such assumption makes the inference vulnerable to outlying observations, whereas scale mixtures of normal distributions have been an interesting alternative to produce robust estimates, keeping the elegancy and simplicity of the maximum likelihood theory. The aim of this paper is to develop an EM-type algorithm for the parameter estimation, and to use the local influence method to assess the robustness aspects of these parameter estimates under some usual perturbation schemes, In order to identify outliers and to criticize the model building we use the local influence procedure in a Study to compare the precision of several thermocouples. (C) 2008 Elsevier B.V. All rights reserved.