Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t-distribution

Summary . A fairly general procedure is studied to perturb a multivariate density satisfying a weak form of multivariate symmetry, and to generate a whole set of non-symmetric densities. The approach is sufficiently general to encompass some recent proposals in the literature, variously related to the skew normal distribution. The special case of skew elliptical densities is examined in detail, establishing connections with existing similar work. The final part of the paper specializes further to a form of multivariate skew t -density. Likelihood inference for this distribution is examined, and it is illustrated with numerical examples.     

Estimation of noncentrality parameters

We propose some estimators of noncentrality parameters which improve upon usual unbiased estimators under quadratic loss. The distributions we consider are the noncentral chi-square and the noncentral F. However, we give more general results for the family of elliptically contoured distributions and propose a robust dominating estimator.     

Averaging Lorenz curves

A large number of functional forms has been suggested in the literature for estimating Lorenz curves that describe the relationship between income and population shares. The traditional way of overcoming functional-form uncertainty when estimating a Lorenz curve is to choose the function that best fits the data in some sense. In this paper we describe an alternative approach for accommodating functional-form uncertainty, namely, how to use Bayesian model averaging to average the alternative functional forms. In this averaging process, the different Lorenz curves are weighted by their posterior probabilities of being correct. Unlike a strategy of picking the best-fitting function, Bayesian model averaging gives posterior standard deviations that reflect the functional-form uncertainty. Building on our earlier work (Chotikapanich and Griffiths, 2002), we construct likelihood functions using the Dirichlet distribution and estimate a number of Lorenz functions for Australian income units. Prior information is formulated in terms of the Gini coefficient and the income shares of the poorest 10% and poorest 90% of the population. Posterior density functions for these quantities are derived for each Lorenz function and are averaged over all the Lorenz functions.     

Generalized additive models for location, scale and shape

Summary.  A general class of statistical models for a univariate response variable is presented which we call the generalized additive model for location, scale and shape (GAMLSS). The model assumes independent observations of the response variable y given the parameters, the explanatory variables and the values of the random effects. The distribution for the response variable in the GAMLSS can be selected from a very general family of distributions including highly skew or kurtotic continuous and discrete distributions. The systematic part of the model is expanded to allow modelling not only of the mean (or location) but also of the other parameters of the distribution of y , as parametric and/or additive nonparametric (smooth) functions of explanatory variables and/or random-effects terms. Maximum (penalized) likelihood estimation is used to fit the (non)parametric models. A Newton–Raphson or Fisher scoring algorithm is used to maximize the (penalized) likelihood. The additive terms in the model are fitted by using a backfitting algorithm. Censored data are easily incorporated into the framework. Five data sets from different fields of application are analysed to emphasize the generality of the GAMLSS class of models.     

基于单片机的粮仓温度测控系统设计

介绍了一种基于AT89C51单片机的温度测控系统。使用AD590作为温度传感器对粮仓内温度进行监测,实时显示温度。系统对温度的控制合理、有效,具有调试方便、精度高和可靠性好等优点。     

Bayesian analysis of a marked point process: Application in seismic hazard assessment

Point processes are the stochastic models most suitable for describing physical phenomena that appear at irregularly spaced times, such as the earthquakes. These processes are uniquely characterized by their conditional intensity, that is, by the probability that an event will occur in the infinitesimal interval (t, t+Δt), given the history of the process up tot. The seismic phenomenon displays different behaviours on different time and size scales; in particular, the occurrence of destructive shocks over some centuries in a seismogenic region may be explained by the elastic rebound theory. This theory has inspired the so-called stress release models: their conditional intensity translates the idea that an earthquake produces a sudden decrease in the amount of strain accumulated gradually over time along a fault, and the subsequent event occurs when the stress exceeds the strength of the medium. This study has a double objective: the formulation of these models in the Bayesian framework, and the assignment to each event of a mark, that is its magnitude, modelled through a distribution that depends at timet on the stress level accumulated up to that instant. The resulting parameter space is constrained and dependent on the data, complicating Bayesian computation and analysis. We have resorted to Monte Carlo methods to solve these problems.     

Real nonparametric regression using complex wavelets

Summary.  Wavelet shrinkage is an effective nonparametric regression technique, especially when the underlying curve has irregular features such as spikes or discontinuities. The basic idea is simple: take the discrete wavelet transform of data consisting of a signal corrupted by noise; shrink or remove the wavelet coefficients to remove the noise; then invert the discrete wavelet transform to form an estimate of the true underlying curve. Various researchers have proposed increasingly sophisticated methods of doing this by using real-valued wavelets. Complex-valued wavelets exist but are rarely used. We propose two new complex-valued wavelet shrinkage techniques: one based on multiwavelet style shrinkage and the other using Bayesian methods. Extensive simulations show that our methods almost always give significantly more accurate estimates than methods based on real-valued wavelets. Further, our multiwavelet style shrinkage method is both simpler and dramatically faster than its competitors. To understand the excellent performance of this method we present a new risk bound on its hard thresholded coefficients.     

国有企业的特殊性与我国国有企业的布局定位

国有企业特殊性的真实内涵表现为,它有着特殊的所有者和所有者代表、特殊的社会目标、特殊的社会经济功能和特殊的产权组织运作机制。因此,越是在以社会目标为主的产业领域,国有企业就越有存在的价值。尽管国有企业的布局是一种趋向性和开放式的动态布局,但它所传递的信息表明国有企业会在特殊产业领域中以合适的方式继续存在并扮演重要角色。     

Testing symmetry under a skew Laplace model

We develop tests of hypothesis about symmetry based on samples from possibly asymmetric Laplace distributions and present exact and limiting distribution of the test statistics. We postulate that the test statistic derived under the Laplace model is a rational choice as a measure of skewness and can be used in testing symmetry for other, quite general classes of skew distributions. Our results are applied to foreign exchange rates for 15 currencies.     

密度分布非均匀的反馈控制Logistic单种群模型平衡态的稳定性

本文对种群密度在非均匀分布情形下,考虑了具反馈控制的滞后 Logistic 生态模型平衡位置的稳定性;分别给出了在常时滞和弱连续时滞以及强连续时滞情况下的稳定性条件;其结果是对 Gopalsamy 在密度均匀分布情形下相应结果的推广.