Bayesian estimation of regression parameters in elliptical measurement error models

The main object of this paper is to discuss the Bayes estimation of the regression coefficients in the elliptically distributed simple regression model with measurement errors. The posterior distribution for the line parameters is obtained in a closed form, considering the following: the ratio of the error variances is known, informative prior distribution for the error variance, and non-informative prior distributions for the regression coefficients and for the incidental parameters. We proved that the posterior distribution of the regression coefficients has at most two real modes. Situations with a single mode are more likely than those with two modes, especially in large samples. The precision of the modal estimators is studied by deriving the Hessian matrix, which although complicated can be computed numerically. The posterior mean is estimated by using the Gibbs sampling algorithm and approximations by normal distributions. The results are applied to a real data set and connections with results in the literature are reported. (C) 2011 Elsevier B.V. All rights reserved.

Minimum size for the occurrence of vortex matter in a square mesoscopic superconductor

The study of superconducting samples in mesoscopic scale presented a remarkable improvement during the last years. Certainly, such interest is based on the fact that when the size of the samples is close to the order of the temperature dependent coherence length xi(T), and/or the size of the penetration depth lambda(T), there are some significant modifications on the physical properties of the superconducting state. This contribution tests the square cross-section size limit for the occurrence (or not) of vortices in mesoscopic samples of area L-2, where L varies discretely from 1 xi(0) to 8 xi(0).The time dependent Ginzburg-Landau (TDGL) equations approach is used upon taking the order parameter and the local magnetic field invariant along the z-direction. The vortex configurations at the equilibrium can be obtained from the TDGL equations for superconductivity as the system relaxes to the stationary state.The obtained results show that the limit of vortex penetration is for the square sample of size 3 xi(0) x 3 xi(0) in which only a single vortex are allowed into the sample. For smaller specimens, no vortex can be formed and the field entrance into the sample is continuous and the total flux penetration occurs at higher values of H/H-c2(0), where H-c2(T) is the upper critical field. Otherwise, for larger samples different vortices patterns can be observed depending on the sample size. (c) 2007 Elsevier B.V. All rights reserved.

Distribution of stream macroalgae in three tropical drainage basins of southeastern Brazil

Three drainage basins belonging to different drainage systems, but in close proximity, were evaluated to compare the distribution of macroalgal communities in the northwest region of São Paulo State, southeastern Brazil. Monthly samplings were carried out from September 1992 through September 1993 in three sites along the main river of each basin. 10 m length cross segments were evaluated for species per cent cover and richness, on both the population and community levels. Selected stream variables were measured: specific conductance, temperature, turbidity, mean cross-sectional area of the reach, current velocity, pH, and oxygen saturation (%). Principal Coordinate Analysis (PCO), chi-square good-of-fitness, and Pearson Product-Moment correlation coefficient were applied to evaluate the distribution of the macroalgal species. The survey resulted in 36 taxa of macroalgae, of which Cyanophyta was the dominant group (17 taxa or 41.7% of the total), followed by Chlorophyta (15 taxa or 41.7%), Rhodophyta (3 taxa or 8.3%) and Chrysophyta (1 taxon or 2.8%). Stigeoclonium helveticum, 'Chantransia' stage of Batrachospermum spp., and B. delicatulum were the most widespread and frequent macroalgae throughout the basins. The analyses showed that conductance and current velocity were the factors most closely related to the distribution of the macroalgal species. Positive correlation between richness and percent cover was determined, which reinforces the patchiness of stream macroalgal distribution.

Spatial and temporal distribution of the zoobenthos community during the filling up period of Porto Primavera Reservoir (Paraná River, Brazil)

Este estudo é parte do monitoramento limnológico empreendido pela Companhia Energética do Estado de São Paulo (CESP) durante o processo de enchimento do reservatório de Porto Primavera (Usina Hidrelétrica Engenheiro Sérgio Motta). Este reservatório, localizado no alto rio Paraná, entre os Estados de São Paulo e Mato Grosso, é o quarto maior do país. A primeira etapa de enchimento do lago começou em dezembro de 1998 e a segunda em março de 2001. Amostras para a análise da comunidade bentônica e das características sedimentológicas foram coletadas trimestralmente entre agosto de 1999 e novembro de 2001 e também em agosto de 2002 (11 campanhas). As coletas foram feitas em 13 estações de amostragem distribuídas no reservatório e em uma localizada a jusante da barragem. 128 táxons de invertebrados foram encontrados, sendo Mollusca, Annelida, Insecta e Nematoda os grupos dominantes durante praticamente todos os meses analisados. A classe Insecta foi a melhor representada, com 9 diferentes ordens, dentro das quais os Diptera contribuíram com a ocorrência de 63 táxons. A espécie exótica de bivalve Corbicula fluminea foi registrada em todas as estações de amostragem mostrando sua grande capacidade para colonizar novos habitats em regiões neotropicais. Variações consideráveis na densidade da fauna foram observadas para os diferentes períodos e locais analisados. A densidade máxima (média de 7812 ind.m-2) foi registrada no centro do reservatório enquanto que as densidades mínimas foram registradas na zona lacustre próxima à barragem (média de 9 ind.m-2). A maior riqueza de espécies por local/período (24 táxons) foi encontrada no trecho superior do reservatório (trecho fluvial). A diversidade máxima foi observada nas zonas superior e central do reservatório, com valores de 3.82 e 3.86 (bits.ind-1) no início (agosto/1999) e final (agosto/2002) do processo de enchimento, respectivamente. Não foi encontrado um padrão de distribuição dos grupos faunísticos que pudesse estar associado com a textura granulométrica dos diferentes locais amostrados. Por outro lado, constatou-se a diminuição, ou mesmo a não ocorrência de organismos, nas estações com elevada concentração de matéria orgânica (>40%) em baixo estado de degradação (grandes detritos vegetais). Tal fato pode estar relacionado com a falta de depósitos de sedimentos, dificultando a fixação de organismos da fauna bentônica, bem como com condições químicas mais redutoras em função da intensidade dos processos de decomposição da fitomassa inundada.

Minimum blood lactate and muscle protein of rats during swimming exercise

Few studies dealing with effort intensity during swimming exercise in rats have been reported in the literature. Recently, with the use of the lactate minimum test (LMT), our group estimated the minimum blood lactate (MBL) of rats during swimming exercises. This information allowed accurate evaluation of the effort intensity developed by rats during swimming exercise. The present study was designed to evaluate the effects of swimming exercise sessions in below, equivalent and above intensities to MBL, on protein metabolism of rats. Adult (90 days) sedentary male Wistar rats were used in the present study. Mean values of MBL, in the present study, were obtained at blood concentration of 6.7 +/- 0.4 mmol/L with a load of 5% bw. The animals were sacrificed at rest (R) or immediately after a single swimming session (30 min) supporting loads below (3.5% bw), equivalent (5.0% bw) and high load (6.5% bw) to AT. Blood samples were collected each 5 min of exercise for lactate determination. Soleus muscle protein synthesis (amount of L-[C-14] fenil alanyn incorporation to protein) and breakdown (tyrosin release) rates were evaluated. Blood lactate concentrations (mmol/L) stabilized with the below (5.4 +/- 0.01) and equivalent (6.4 +/- 0.006) to MBL but increased, progressively, with the high load. There were no differences in protein synthesis (pmol/mg.h) among rest values (65.2 +/- 3.4) and after-exercise supporting the loads below (61.5 +/- 1.3) and the equivalent (60.7+/-1.7) to MBL but there was a decrease with the high load (36.6+/-2.0). Protein breakdown rates (pmol/g.h) increase after exercise supporting the loads below (227.0 +/- 6.1), equivalent (227.9 +/- 6.0) and high (363.6 +/- 7.1) to MBL in relation to the rest (214.3 +/- 6.0). The results indicate the viability of the application of LMT in studies with rats since it detected alterations imposed by exercise.

Control of thermal damage in grinding by digital signal processing of raw acoustic emission

Grinding process is usually the last finishing process of a precision component in the manufacturing industries. This process is utilized for manufacturing parts of different materials, so it demands results such as low roughness, dimensional and shape error control, optimum tool-life, with minimum cost and time. Damages on the parts are very expensive since the previous processes and the grinding itself are useless when the part is damaged in this stage. This work aims to investigate the efficiency of digital signal processing tools of acoustic emission signals in order to detect thermal damages in grinding process. To accomplish such a goal, an experimental work was carried out for 15 runs in a surface grinding machine operating with an aluminum oxide grinding wheel and ABNT 1045 e VC131 steels. The acoustic emission signals were acquired from a fixed sensor placed on the workpiece holder. A high sampling rate acquisition system at 2.5 MHz was used to collect the raw acoustic emission instead of root mean square value usually employed. In each test AE data was analyzed off-line, with results compared to inspection of each workpiece for burn and other metallurgical anomaly. A number of statistical signal processing tools have been evaluated.

A Galactic ring of minimum stellar density near the solar orbit radius

We analyse the secular effects of a long-lived Galactic spiral structure on the stellar orbits with mean radii close to the corotation resonance. By test-particle simulations and different spiral potential models with parameters constrained on observations, we verified the formation of a minimum with amplitude ∼30–40 per cent of the background disc stellar density at corotation. Such a minimum is formed by the secular angular momentum transfer between stars and the spiral density wave on both sides of corotation. We demonstrate that the secular loss (gain) of angular momentum and decrease (increase) of mean orbital radius of stars just inside (outside) corotation can counterbalance the opposite trend of exchange of angular momentum shown by stars orbiting the librational points L4/5 at the corotation circle. Such secular processes actually allow steady spiral waves to promote radial migration across corotation. We propose some pieces of observational evidence for the minimum stellar density in the Galactic disc, such as its direct relation to the minimum in the observed rotation curve of the Galaxy at the radius r ∼ 9 kpc (for R0 = 7.5 kpc), as well as its association with a minimum in the distribution of Galactic radii of a sample of open clusters older than 1Gyr. The closeness of the solar orbit adius to the corotation resonance implies that the solar orbit lies inside a ring of minimum surface density (stellar + gas). This also implies a correction to larger values for the estimated total mass of the Galactic disc, and consequently, a greater contribution of the disc componente to the inner rotation curve of the Galaxy.

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Individualization of 5-fluorocytosine therapy

Using a convenient and fast HPLC procedure we determined serum concentrations of the fungistatic agent 5-fluorocytosine (5-FC) in 375 samples from 60 patients treated with this drug. The mean trough concentration (n = 127) was 64.3 mg/l (range: 11.8-208.0 mg/l), the mean peak concentration (n = 122) was 99.9 mg/l (range: 25.6-263.8 mg/l), the mean nonpeak/nontrough concentration (n = 126) was 80.1 mg/l (range: 10.5-268.0 mg/l). Totally 134 (35.7%) samples were outside the therapeutic range (25-100 mg/l), 108 (28.8%) being too high, 26 (6.9%) being too low. Forty-four (73%) patients showed 5-FC serum concentrations outside the therapeutic range at least once during the treatment course. In a prospective study we performed 65 dosage predictions on 30 patients by use of a 3-point method previously developed for aminoglycoside dosage adaptation. The mean absolute prediction error of the dosage adaptation was +0.7 mg/l (range: -26.0 to +28.0 mg/l). The root mean square prediction error was 10.7 mg/l. The mean predicted concentration (65.3 mg/l) agreed very well with the mean measured concentration (64.6 mg/l). The frequency distribution of 5-FC serum concentrations indicates that 5-FC monitoring is important. The applied pharmacokinetic method allows individual adaptations of 5-FC dosage with a clinically acceptable prediction error.

Determination of the minimum alveolar concentration of isoflurane in Shetland ponies using constant current or constant voltage electrical stimulation

OBJECTIVE: To determine the minimum alveolar concentration (MAC) of isoflurane in Shetland ponies using a sequence of three different supramaximal noxious stimulations at each tested concentration of isoflurane rather than a single stimulation. STUDY DESIGN: Prospective, experimental trial. ANIMALS: Seven 4-year-old, gelding Shetland ponies. METHODS: The MAC of isoflurane was determined for each pony. Three different modes of electrical stimulation were applied consecutively (2 minute intervals): two using constant voltage (90 V) on the gingiva via needle- (CVneedle) or surface-electrodes (CVsurface) and one using constant current (CC; 40 mA) via surface electrodes applied to the skin over the digital nerve. The ability to clearly interpret the responses as positive, the latency of the evoked responses and the inter-electrode resistance were recorded for each stimulus. RESULTS: Individual isoflurane MAC (%) values ranged from 0.60 to 1.17 with a mean (+/-SD) of 0.97 (+/-0.17). The responses were more clearly interpreted with CC, but did not reach statistical significance. The CVsurface mode produced responses with a longer delay. The CVneedle mode was accompanied by variable inter-electrode resistances resulting in uncontrolled stimulus intensity. At 0.9 MAC, the third stimulation induced more positive responses than the first stimulation, independent of the mode of stimulation used. CONCLUSIONS: The MAC of isoflurane in the Shetland ponies was lower than expected with considerable variability among individuals. Constant current surface electrode stimulations were the most repeatable. A summation over the sequence of three supramaximal stimulations was observed around 0.9 MAC. CLINICAL RELEVANCE: The possibility that Shetland ponies require less isoflurane than horses needs further investigation. Constant current surface-electrode stimulations were the most repeatable. Repetitive supramaximal stimuli may have evoked movements at isoflurane concentrations that provide immobility when single supramaximal stimulation was applied.

Ultraviolet radiation intensity predicts the relative distribution of dermatomyositis and anti-Mi-2 autoantibodies in women.

OBJECTIVE: Because studies suggest that ultraviolet (UV) radiation modulates the myositis phenotype and Mi-2 autoantigen expression, we conducted a retrospective investigation to determine whether UV radiation may influence the relative prevalence of dermatomyositis and anti-Mi-2 autoantibodies in the US. METHODS: We assessed the relationship between surface UV radiation intensity in the state of residence at the time of onset with the relative prevalence of dermatomyositis and myositis autoantibodies in 380 patients with myositis from referral centers in the US. Myositis autoantibodies were detected by validated immunoprecipitation assays. Surface UV radiation intensity was estimated from UV Index data collected by the US National Weather Service. RESULTS: UV radiation intensity was associated with the relative proportion of patients with dermatomyositis (odds ratio [OR] 2.3, 95% confidence interval [95% CI] 0.9-5.8) and with the proportion of patients expressing anti-Mi-2 autoantibodies (OR 6.0, 95% CI 1.1-34.1). Modeling of these data showed that these associations were confined to women (OR 3.8, 95% CI 1.3-11.0 and OR 17.3, 95% CI 1.8-162.4, respectively) and suggests that sex influences the effects of UV radiation on autoimmune disorders. Significant associations were not observed in men, nor were UV radiation levels related to the presence of antisynthetase or anti-signal recognition particle autoantibodies. CONCLUSION: This first study of the distribution of myositis phenotypes and UV radiation exposure in the US showed that UV radiation may modulate the clinical and immunologic expression of autoimmune disease in women. Further investigation of the mechanisms by which these effects are produced may provide insights into pathogenesis and suggest therapeutic or preventative strategies.

Estimating the mean and standard deviation of concentration distributions using censored samples

Environmental data sets of pollutant concentrations in air, water, and soil frequently include unquantified sample values reported only as being below the analytical method detection limit. These values, referred to as censored values, should be considered in the estimation of distribution parameters as each represents some value of pollutant concentration between zero and the detection limit. Most of the currently accepted methods for estimating the population parameters of environmental data sets containing censored values rely upon the assumption of an underlying normal (or transformed normal) distribution. This assumption can result in unacceptable levels of error in parameter estimation due to the unbounded left tail of the normal distribution. With the beta distribution, which is bounded by the same range of a distribution of concentrations, $\rm\lbrack0\le x\le1\rbrack,$ parameter estimation errors resulting from improper distribution bounds are avoided. This work developed a method that uses the beta distribution to estimate population parameters from censored environmental data sets and evaluated its performance in comparison to currently accepted methods that rely upon an underlying normal (or transformed normal) distribution. Data sets were generated assuming typical values encountered in environmental pollutant evaluation for mean, standard deviation, and number of variates. For each set of model values, data sets were generated assuming that the data was distributed either normally, lognormally, or according to a beta distribution. For varying levels of censoring, two established methods of parameter estimation, regression on normal ordered statistics, and regression on lognormal ordered statistics, were used to estimate the known mean and standard deviation of each data set. The method developed for this study, employing a beta distribution assumption, was also used to estimate parameters and the relative accuracy of all three methods were compared. For data sets of all three distribution types, and for censoring levels up to 50%, the performance of the new method equaled, if not exceeded, the performance of the two established methods. Because of its robustness in parameter estimation regardless of distribution type or censoring level, the method employing the beta distribution should be considered for full development in estimating parameters for censored environmental data sets. ^

DISTRIBUTION FUNCTIONS, MEAN SQUARED ERRORS, AND CONFIDENCE LIMITS FOR RIDGE REGRESSION ESTIMATORS

A large number of ridge regression estimators have been proposed and used with little knowledge of their true distributions. Because of this lack of knowledge, these estimators cannot be used to test hypotheses or to form confidence intervals.^ This paper presents a basic technique for deriving the exact distribution functions for a class of generalized ridge estimators. The technique is applied to five prominent generalized ridge estimators. Graphs of the resulting distribution functions are presented. The actual behavior of these estimators is found to be considerably different than the behavior which is generally assumed for ridge estimators.^ This paper also uses the derived distributions to examine the mean squared error properties of the estimators. A technique for developing confidence intervals based on the generalized ridge estimators is also presented. ^

Profiles of diatoms, pollen, XRD, XRF and organic composition of two sediment cores (PG2022 and PG2023) from Lake Kyutyunda in Yakutia, NE Siberia, Russia

Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between ~13,500 and ~8,900 cal. years BP and possibly during the 8,200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8°C during the Holocene Thermal Maximum (HTM) between ~8,900 and ~4,500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3,000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice-sheet vanished 7,000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions.

(Table 2) Distribution of alkenone unsaturation index UK37 of the eastern North Atlantic

This paper reports the concentrations and within-class distributions of long-chain alkenones and alkyl alkenoates in the surface waters (0-50 m) of the eastern North Atlantic, and correlates their abundance and distribution with those of source organisms and with water temperature and other environmental variables. We collected these samples of >0.8 µm particulate material from the euphotic zone along the JGOFS 20°W longitude transect, from 61°N to 24°N, during seven cruises of the UK-JGOFS Biogeochemical Ocean Flux Study (BOFS) in 1989-1991; the biogeographical range of our 53 samples extends from the cold (<10°C), nutrient-rich and highly productive subarctic waters of the Iceland Basin to the warm (>25°C) oligotrophic subtropical waters off Africa. Surface water concentrations of total alkenone and alkenoates ranged from <50 ng/l in oligotrophic waters below 40°N to 2000-4500 ng/l in high latitude E. huxleyi blooms, and were well correlated with E. huxleyi cell densities, supporting the assumption that E. huxleyi is the predominant source of these compounds in the present day North Atlantic. The within-class distribution of the C37 and C38 alkenones and C36 alkenoates varied strongly as a function of temperature, and was largely unaffected by nutrient concentration, bloom status and other surface water properties. The biosynthetic response of the source organisms to growth temperature differed between the cold (<16°C) waters above 47°N and the warmer waters to the south. In cold (<16°C) waters above 47°N, the relative amounts of alkenoates and C38 alkenones synthesized was a strong function of growth temperature, while the unsaturation ratio of the alkenones (C37 and C38) was uncorrelated with temperature. Conversely, in warm (>16°C) waters below 47°N, the relative proportions of alkenoates and alkenones synthesized remained constant with increasing temperature while the unsaturation ratios of the C37 and C38 methyl alkenones (Uk37 and Uk38Me, respectively) increased linearly. The fitted regressions of Uk37 and Uk38Me versus temperature for waters >16°C were both highly significant (r**2 > 0.96) and had identical slopes (0.057) that were 50% higher than the slope (0.034) of the temperature calibration of Uk37 reported by Prahl and Wakeham (1987; doi:10.1038/330367a0) over the same temperature range. These observations suggest either a physiological adjustment in biochemical response to growth temperature above a 16-17°C threshold and/or variation between different E. huxleyi strains and/or related species inhabiting the cold and warm water regions of the eastern North Atlantic. Using our North Atlantic data set, we have produced multivariate temperature calibrations incorporating all major features of the alkenone and alkenoate data set. Predicted temperatures using multivariate calibrations are largely unbiased, with a standard error of approximately ±1°C over the entire data range. In contrast, simpler calibration models cannot adequately incorporate regional diversity and nonlinear trends with temperature. Our results indicate that calibrations based upon single variables, such as Uk37, can be strongly biased by unknown systematic errors arising from natural variability in the biosynthetic response of the source organisms to growth temperature. Multivariate temperature calibration can be expected to give more precise estimates of Integrated Production Temperatures (IPT) in the sedimentary record over a wider range of paleoenvironmental conditions, when derived using a calibration data set incorporating a similar range of natural variability in biosynthetic response.