A note on the likelihood and moments of the skew-normal distribution

In this paper an alternative approach to the one in Henze (1986) is proposed for deriving the odd moments of the skew-normal distribution considered in Azzalini (1985). The approach is based on a Pascal type triangle, which seems to greatly simplify moments computation. Moreover, it is shown that the likelihood equation for estimating the asymmetry parameter in such model is generated as orthogonal functions to the sample vector. As a consequence, conditions for a unique solution of the likelihood equation are established, which seem to hold in more general setting.

LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

An extended BEM formulation for plates reinforced by rectangular beams

This article presents a BEM formulation developed particularly for analysis of plates reinforced by rectangular beams. This is an extended version of a Previous paper that only took into account bending effects. The problem is now re-formulated to consider bending and membrane force effects. The effects of the reinforcements are taken into account by using a simplified scheme that requires application of ail initial stress field to locally correct the bending and stretching stiffness of the reinforcement regions. The domain integrals due to the presence of the reinforcements are then transformed to the reinforcement/plate interface. To reduce the number of degrees of freedom related to the presence of the reinforcement, the proposed model was simplified to consider only bending and stretching rigidities in the direction of the beams. The complete model can be recovered by applying all six internal force correctors, corresponding to six degrees of freedom per node. Examples are presented to confirm the accuracy of the formulation and to illustrate the level of simplification introduced by this strong reduction in the number of degrees of freedom. (C) 2008 Elsevier Ltd. All rights reserved.

Extended Fault-Location Formulation for Power Distribution Systems

In this paper, an extended impedance-based fault-location formulation for generalized distribution systems is presented. The majority of distribution feeders are characterized by having several laterals, nonsymmetrical lines, highly unbalanced operation, and time-varying loads. These characteristics compromise traditional fault-location methods performance. The proposed method uses only local voltages and currents as input data. The current load profile is obtained through these measurements. The formulation considers load variation effects and different fault types. Results are obtained from numerical simulations by using a real distribution system from the Electrical Energy Distribution State Company of Rio Grande do Sul (CEEE-D), Southern Brazil. Comparative results show the technique robustness with respect to fault type and traditional fault-location problems, such as fault distance, resistance, inception angle, and load variation. The formulation was implemented as embedded software and is currently used at CEEE-D`s distribution operation center.

Modeling extended Petri nets compatible with GHENeSys IEC61131 for industrial automation

Petri net (PN) modeling is one of the most used formal methods in the automation applications field, together with programmable logic controllers (PLCs). Therefore, the creation of a modeling methodology for PNs compatible with the IEC61131 standard is a necessity of automation specialists. Different works dealing with this subject have been carried out; they are presented in the first part of this paper [Frey (2000a, 2000b); Peng and Zhou (IEEE Trans Syst Man Cybern, Part C Appl Rev 34(4):523-531, 2004); Uzam and Jones (Int J Adv Manuf Technol 14(10):716-728, 1998)], but they do not present a completely compatible methodology with this standard. At the same time, they do not maintain the simplicity required for such applications, nor the use of all-graphical and all-mathematical ordinary Petri net (OPN) tools to facilitate model verification and validation. The proposal presented here completes these requirements. Educational applications at the USP and UEA (Brazil) and the UO (Cuba), as well as industrial applications in Brazil and Cuba, have already been carried out with good results.

Lightning-Induced Voltages on Overhead Lines-Application of the Extended Rusck Model

Lightning-induced overvoltages have a considerable impact on the power quality of overhead distribution and telecommunications systems, and various models have been developed for the computation of the electromagnetic transients caused by indirect strokes. The most adequate has been shown to be the one proposed by Agrawal et al.; the Rusck model can be visualized as a particular case, as both models are equivalent when the lightning channel is perpendicular to the ground plane. In this paper, an extension of the Rusck model that enables the calculation of lightning-induced transients considering flashes to nearby elevated structures and realistic line configurations is tested against data obtained from both natural lightning and scale model experiments. The latter, performed under controlled conditions, can be used also to verify the validity of other coupling models and relevant codes. The so-called Extended Rusck Model, which is shown to be sufficiently accurate, is applied to the analysis of lightning-induced voltages on lines with a shield wire and/or surge arresters. The investigation conducted indicates that the ratio between the peak values of the voltages induced by typical first and subsequent strokes can be either greater or smaller than the unity, depending on the line configuration.

Time stability of soil water storage measured by neutron probe and the effects of calibration procedures in a small watershed

The knowledge of soil water storage (SWS) of soil profiles is crucial for the adoption of vegetation restoration practices. With the aim of identifying representative sites to obtain the mean SWS of a watershed, a time stability analysis of neutron probe evaluations of SWS was performed by the means of relative differences and Spearman rank correlation coefficients. At the same time, the effects of different neutron probe calibration procedures were explored on time stability analysis. mean SWS estimation. and preservation of the spatial variability of SWS. The selected watershed, with deep gullies and undulating slopes which cover an area of 20 ha, is characterized by an Ust-Sandiic Entisol and an Aeolian sandy soil. The dominant vegetation species are bunge needlegrass (Stipa bungeana Trim) and korshinsk peashrub (Carugano Korshinskii kom.). From June 11, 2007 to July 23,2008, SWS of the top1 m soil layer was evaluated for 20 dates, based on neutron probe data of 12 sampling sites. Three calibration procedures were employed: type 1, most complete, with each site having its own linear calibration equation (TrE); type II. with TrE equations extended over the whole field: and type III, with one single linear calibration curve for the whole field (UnE) and also correcting its intercept based on site specific relative difference analysis (RdE) and on linear fitting of data (RcE), both maintaining the same slope. A strong time stability of SWS estimated by TrE equations was identified. Soil particle size and soil organic matter content were recognized as the influencing factors for spatial variability of SWS. Land use influenced neither the spatial variability nor the time stability of SWS. Time stability analysis identified one site to represent the mean SWS of the whole watershed with mean absolute percentage errors of less than 10%, therefore. this site can be used as a predictor for the mean SWS of the watershed. Some equations of type II were found to be unsatisfactory to yield reliable mean SWS values or in preserving the associated soil spatial variability. Hence, it is recommended to be cautious in extending calibration equations to other sites since they might not consider the field variability. For the equations with corrected intercept (type III), which consider the spatial variability of calibration in a different way in relation to TrE, it was found that they can yield satisfactory means and standard deviation of SWS, except for the RdE equations, which largely leveled off the SWS values in the watershed. Correlation analysis showed that the neutron probe calibration was linked to soil bulk density and to organic matter content. Therefore, spatial variability of soil properties should be taken into account during the process of neutron probe calibration. This study provides useful information on the mean SWS observation with a time stable site and on distinct neutron probe calibration procedures, and it should be extended to soil water management studies with neutron probes, e.g., the process of vegetation restoration in wider area and soil types of the Loess Plateau in China. (C) 2009 Elsevier B.V. All rights reserved.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

Fate of Surface-Inoculated Escherichia coli O157:H7, Listeria monocytogenes, and Salmonella Typhimurium on Kippered Beef during Extended Storage at Refrigeration and Abusive Temperatures

The behavior of Escherichia coli O157:H7, Listeria monocytogenes, and Salmonella Typhimurium on kippered beef was evaluated. Individual pieces of the product were separately inoculated on the top and bottom surfaces with each three- to six-strain pathogen cocktail at ca. 6.0 log CFU per piece and stored at 4, 10, 21, or 30 degrees C for up to 28 days in each of two trials. When kippered beef was inoculated with E coli O157:H7, Salmonella Typhimurium, or L. monocytogenes and stored at 4, 10, 2 1, or 30 degrees C for up to 28 days, pathogen numbers decreased ca. 0.4 to 0.9, 1.0 to 1.8, 3.0 to >= 5.25, and >= 5.0 to 5.25 log CFU per piece, respectively. Average D-values for E. coli O157:H7, Salmonella Typhimurium, and L. monocytogenes stored at 4 to 30 degrees C for 28 days were ca. 41 to 4.6, 40.8 to 5.3, and 29.5 to 4.3 days, respectively. As expected, the higher the storage temperature, the greater the level and rate of inactivation for all three pathogens. These data establish that kippered beef does not provide an environment conducive to proliferation of these pathogens.

Predominance of CTX-M-type extended-spectrum beta-lactamase genes among enterobacterial isolates from outpatients in Brazil

Two hundred fifty-seven nalidixic acid-resistant enterobacterial isolates were collected in a Brazilian community from January 2000 to May 2005 to determine the prevalence of plasmid-encoded extended-spectrum beta-lactamases. The bla(CTX-M) genetic environment was determined by polymerase chain reaction and sequencing. Eleven isolates (4.2%) harbored a bla(CTX-M-2) gene, 3 isolates bla(CTX-M-9), 2 isolates bla(CTX-M-8), and 6 isolates bla(SHV-5). Two novel bla(CTX-M-2) variants, namely, bla(CTX-M-74) and bla(CTX-M-75), were identified. (C) 2009 Elsevier Inc. All rights reserved.

Investigation of the antinociceptive efficacy and relative potency of extended duration injectable 3-acylmorphine-6-sulfate prodrugs in rats

This investigation was designed to examine the antinociceptive activity in rats of 3-O-acyl prodrugs of M6S relative to the parent drug, after intravenous and intramuscular injection, using the tail flick latency test of antinociception. M6S, 3-acetylmorphine-6-sulfate (3AcM6S), 3-propionylmorphine-6-sulfate (3PrM6S), 3-butanoylmorphine-6-sulfate (3BuM6S) and 3-heptanoylmorphine-6-sulfate (3HpM6S) were administered by the IV route in a dose of 4.10 mu mol/kg. Relatively high levels of antinociception (>40% Maximum Possible Effect) were achieved following administration of M6S, 3AcM6S and 3PrM6S, whereas insignificant antinociception (<20%MPE) was achieved following administration of 3BuM6S or 3HpM6S. Although the mean duration of action for 3AcM6S (6 h) was longer than for M6S or 3PrM6S (4 h), the mean area (+/- S.E.M.) under the degree of antinociception versus time curve (AUG) for 3AcM6S (151.6 +/- 6.9%MPE h) was not significantly different (p <0.05) from that for M6S (120.8 +/- 32.7%MPE h) or for 3PrM6S (106.0 +/- 21.3%MPE h). The mean ED50 (range) doses for M6S, 3AcM6S and 3PrM6S were calculated to be 4.16 (3.61-4.48), 4.32 (3.55-5.09) and 4.54 (4.21-4.79) mu mol/kg, respectively. Preliminary studies were conducted on potential long-acting formulations containing 8 x ED50 doses of M6S and the 3-acetyl and 3-propionyl esters suspended in soybean oil. These showed that 3PrM6S gave a greater AUC (mean + S.E.M.) (1087.4 +/- 97.4%MPE h) and longer duration of action (20 h) than did M6S (613.1 +/- 155.9%MPE h; 10 h duration) or 3AcM6S (379.3 + 114.2%MPE h: 8 h duration). Further studies are needed to more fully investigate these findings. (C) 1998 Elsevier Science B.V. All rights reserved.

Extended delayed recall of AVLT word lists: Effects of age and sex on adult performance

An extension of a previous study of age and sex effects on verbal recall (Geffen, Moar, O'Hanlon, lark, & Geffen, 1990) examined forgetting of words over extended delays. The AVLT was administered to 201 normal adults (99 males, 102 females) ranging in age between 20 and 59 years. Recall was tested at intervals of 30 minutes, 24 hours, and 7 days after acquisition. Testing of the latter two intervals was conducted by telephone (Experiment 1, N = 177). After 30 minutes there was no significant loss of the 10 to 11 words retained from five acquisition trials. However, an overall mean of about one word was forgotten after 1 day and a further word after 7 days. The oldest age group (50-59 years) acquired fewer words and forgot more words than the younger groups. Females of all age groups performed slightly better than males at acquisition, at retention, and at recall after longer delays. A second experiment showed that telephone testing at the longer delay intervals was equivalent to testing face to face. These results extend the use of the AVLT by assessing memory decay beyond the immediate testing period. Telephone follow-up is a convenient and economical method of testing delayed recall.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.