Reconstruction of polygonal images

Let T be a two-dimensional region, and let X be a surface dejined on T. The values of X on T, constitute an image, or pattern. The true value of X at any point on T cannot be directly observed, but data can be recorded which provide information about X. The aim is to reconstruct X using the prior knowledge that X will vary smoothly over most of T, but may exhibit jump discontinuities over line segments. This information can be incorporated via Bayes' theorem, using a polygonal Markov random field on T as prior distribution. Under this continuum model, X may in principle be estimated according to standard criteria. In practice, the techniques rely on simulation of the posterior distribution. A natural family of conjugate priors is identified, and a class of spatial-temporal Markov processes is constructed on the uncountable state space; simulation then proceeds by a method of analogous to the Gibbs sampler.     

Bayesian tests for matched proportions

A sample of n subjects is observed in each of two states, S₁-and S₂. In each state, a subject is in one of two conditions, X or Y. Thus, a subject may be recorded as showing a change if its condition in the two states is ‘Y,X’ or ‘X,Y’ and, otherwise, the condition is unchanged. We consider a Bayesian test of the null hypothesis that the probability of an ‘X,Y’ change exceeds that of a ‘Y,X’ change by amount k_O. That is, we develop the posterior distribution of k_O, the difference between the two probabilities and reject the null hypothesis if k lies outside the appropriate posterior probability interval. The performance of the method is assessed by Monte Carlo and other numerical studies and brief tables of exact critical values are presented     

Nonparametric estimation of the hazard function by using a model selection method: estimation of cancer deaths in Hiroshima atomic bomb survivors

Summary.  Controversy has intensified regarding the death-rate from cancer that is induced by a dose of radiation. In the models that are usually considered the hazard function is an increasing function of the dose of radiation. Such models can mask local variations. We consider the models of excess relative risk and of absolute risk and propose a nonparametric estimation of the effect of the dose by using a model selection procedure. This estimation deals with stratified data. We approximate the function of the dose by a collection of splines and select the best one according to the Akaike information criterion. In the same way between the models of excess relative risk or excess absolute risk, we choose the model that best fits the data. We propose a bootstrap method for calculating a pointwise confidence interval of the dose function. We apply our method for estimating the solid cancer and leukaemia death hazard functions to Hiroshima.     

An alternate version of the conceptual predictive statistic based on a symmetrized discrepancy measure

The conceptual predictive statistic, C_p, is a widely used criterion for model selection in linear regression. C_p serves as an estimator of a discrepancy, a measure that reflects the disparity between the generating model and a fitted candidate model. This discrepancy, based on scaled squared error loss, is asymmetric: an alternate measure is obtained by reversing the roles of the two models in the definition of the measure. We propose a variant of the C_p statistic based on estimating a symmetrized version of the discrepancy targeted by C_p. We claim that the resulting criterion provides better protection against overfitting than C_p, since the symmetric discrepancy is more sensitive towards detecting overspecification than its asymmetric counterpart. We illustrate our claim by presenting simulation results. Finally, we demonstrate the practical utility of the new criterion by discussing a modeling application based on data collected in a cardiac rehabilitation program at University of Iowa Hospitals and Clinics.     

Estimation with truncated inverse binomial sampling

Consider truncated samples taken from an infinite population with a fixed uumber n of observations recorded. A randon number X of items must be sampled in order to observe after trunca-tion. adified maximun liicel.ihooa estimators of X (assumed un-'mown) and of the population parameters are develo~cd, and a computing scheme is given for the exponential distribution. On the basis of asymptotu ,operties, some estimators are singled out and compared with the usual maximum likelihood estimators.     

Nonparametric estimation of survival functions for doubly censored observations when the lifetime hazard rate is proportionally related to the censoring hazard rates

Let T, X and Y be non-negative random variables, where T is the time of occurrence of an event of interest, X and Y being the lefl and right censoring variables respectively.

In this paper we propose a nonparametric estimator of the survival function, ST, when T, X and Y are supposed to be independent and their corresponding hazard rates are proportionally related. In this way, our results extend Ebrahimi's work (1985) to the doubly censored data case.     

Adaptive estimation of the conditional cumulative distribution function from current status data

Consider a positive random variable of interest Y depending on a covariate X, and a random observation time T independent of Y given X. Assume that the only knowledge available about Y is its current status at time T  : δ=_I{Y≤T}$\delta ={\mathbb{I}}_{\{Y\le T\}}$ with I$\mathbb{I}$ the indicator function. This paper presents a procedure to estimate the conditional cumulative distribution function F of Y given X   from an independent identically distributed sample of (X,T,δ)$(X,T,\delta )$.     

Pseudolikelihood estimation in a class of problems with response-related missing covariates

Many practical situations involve a response variable Y and covariates X , where data on (Y, X ) are incomplete for some portion of a sample of individuals. We consider two general types of pseudolikelihood estimation for problems in which missingness may be response-related. These are typically simpler to implement than ordinary maximum likelihood, which in this context is semiparametric. Asymptotics for the pseudolikelihood methods are presented, and simulations conducted to investigate the methods for an important class of problems involving lifetime data. Our results indicate that for these problems the two methods are effective and comparable with respect to efficiency.     

Covariate effect on the life expectancy and percentile residual life functions under the proportional hazards and the accelerated life models

The present paper is concerned with statistical models for the dependence of survival time or time to occurrence of an event, such as time to tumor, on a vector X of covariates or prognostic variables such as age, sex, blood pressure, length of exposure to a toxic material, etc., measured on a group of individuals in biomedical investigations. It is assumed that the covariates influence the distribution of time to tumor only through a linear predictor μ =β^′X.

The object of our paper is to investigate the effect due to the covariates on the Life Expectancy and the Percentile Residual Life (PRL) function of a family of organisms under the proportional hazards and the accelerated life models. The key result is that the families of survival distributions under these models have the 'setting the clock back to zero' property if the family of baseline survival distributions does. This property is a generalization of the lack of memory property of the exponential distribution. Simple examples of the members of this family are the linear hazard exponential, Pareto and Gompertz life distributions.

As a simple application of the main results obtained in the present paper, we have considered a stochastic survival model recently proposed by Chiang and Conforti (1989) for the time-to-tumor distribution in the context of a large-scale serial sacrifice experiment by the National Center of Toxicological Research (NCTR). This involves some mice that were fed 2-AAF from infancy and those that developed bladder and/or liver neoplasms, see Farmer et al (1980). It is shown that their stochastic model for tumor incidence intensity at time t leads to a family of survival models that has the setting the clock back to zero property. The survival functions and the effect of the vector X of covariates on the PRL and the tumor-free life expectancies are evaluated for the proportional hazards and accelerated life models.     

Bivariate Contours of Copula

Quantile functions associated with bivariate copulas are considered. Some of their structural properties are studied. Quantile functions allow one to express the cdf of the random variable C(X, Y), where (X, Y) is distributed as C(x, y) and where C is a copula. Quantile functions provide also a simple algorithm for simulating random observations.     

Effect on probabilities and quantiles of adding a quantity with small variance

This paper gives simple approximations for the distribution function and quantiles of the sum X + Y when X is a continuous variable and Y is an independent variable with variance small compared to that of X . The approximations are based around the distribution function or quantiles of X and require only the first two or three moments of Y to be known. Example evaluations with X having a normal, Student's t or chi-squared distribution suggest that the approximations are good in unbounded tail regions when the ratio of variances is less than 0.2.     

Large sample tests for testing symmetry and independence in some bivariate exponential models

Bivariate Exponential Distribution (BVED) were introduced by Freund (1961), Marshall and Olkin (1967) and Block and Basu (1974) as models for the distributions of (X,Y) the failure times of dependent components (C₁,C₂). We study the structure of these models and observe that Freund model leads to a regular exponential family with a four dimensional orthogonal parameter. Marshall-Olkin model involving three parameters leads to a conditional or piece wise exponential family and Block-Basu model which also depends on three parameters is a sub-model of the Freund model and is a curved exponential family. We obtain a large sample tests for symmetry as well as independence of (X,Y) in each of these models by using the Generalized Likelihood Ratio Tests (GLRT) or tests basesd on MLE of the parameters and root n consistent estimators of their variance-covariance matrices.     

On Censored Bivariate Random Variables: Copula,Characterization, and Estimation

Let (X, Y) be a bivariate random vector with joint distribution function F_{X, Y}(x, y) = C(F(x), G(y)), where C is a copula and F and G are marginal distributions of X and Y, respectively. Suppose that (X_i, Y_i), i = 1, 2, …, n is a random sample from (X, Y) but we are able to observe only the data consisting of those pairs (X_i, Y_i) for which X_i ? Y_i. We denote such pairs as (X*_i, Y_i*), i = 1, 2, …, ν, where ν is a random variable. The main problem of interest is to express the distribution function F_{X, Y}(x, y) and marginal distributions F and G with the distribution function of observed random variables X* and Y*. It is shown that if X and Y are exchangeable with marginal distribution function F, then F can be uniquely determined by the distributions of X* and Y*. It is also shown that if X and Y are independent and absolutely continuous, then F and G can be expressed through the distribution functions of X* and Y* and the stress–strength reliability P{X ? Y}. This allows also to estimate P{X ? Y} with the truncated observations (X*_i, Y_i*). The copula of bivariate random vector (X*, Y*) is also derived.     

Characterization of generalized markov-polya and generalized polya-eggenberger distributions

A discrete model is considered where the original observation is subjected to partial destruction according to the Generalized Markov-Polya (GMP) damage model. A characterization of the Generalized Polya-Eggenberger distribution (GPED) is given in the context of the Rao-Rubin condition. More specifically, if the probability that an observation n of a non-negative integer valued r.v.X is reduced to an integer k during a damage, process is given by the GMPD, and if the resulting r.v.Y is such thatrit satisfies the RR-conditlon, then X has a GPED. Secondly, if N = A + B, where B is the missing part and A is the recorded part such that the conditional distribution P(A= x|N=n) is the GMPD, then the r.v.'s A and B are independent if, and only if, N has a GPED. Several other characterizations are also given for these two distributions. The results of Rao-Rubin ‘1964’, Patil-Ratnaparkhi (1977) and Consul (1975) follow as special cases.     

Marginal Regression Model with Time-Varying Coefficients for Panel Data

In this article, we formulate a class of semiparametric marginal means models with a mixture of time-varying and time-independent parameters for analyzing panel data. For inference about the regression parameters, an estimation procedure is developed and asymptotic properties of the proposed estimators are established. In addition, some tests are presented for investigating whether or not covariate effects vary with time. The finite-sample behavior of the proposed methods is examined in simulation studies, and the data from an AIDS clinical trial study are used to illustrate the methodology.     

Probability density function of quotient of order statistics from the pareto,power and weibull distributions

This paper deals with the probability density functions of quotient of order statistics. We use the Mellin transform technique, to find the distribution of the quotient Z= X/Xwhere X.,X(i < j) are the ith and jth order statistics from the Pareto, Power and Weibull distributions     

城镇化发展与财政政策相关关系的实证分析

近些年来的实践经验表明,我国城镇化水平的逐步上升与财政政策的支持有着密切的关系,因此建立合理的城镇化财政政策支持体系是我国面临的重大现实课题。本文突破了规范性研究的范畴,将财政政策对城镇化发展的支持进行量化,通过建立在向量自回归模型的基础上,运用脉冲响应函数和预测方差分解来研究财政政策和城镇化发展之间的动态相关性以及三者之间的交互响应情况和响应路径。     

Decision-theoretic estimation of generalized variance and generalized precision

Consider a random data matrix X=(X₁,...,X_k):pXk with independent columns [sathik] and an independent p X p Wishart matrix [sathik]. Estimators dominating the best affine equivariant estimators of [sathik] are obtained under four types of loss functions. Improved estimators (Testimators) of generalized variance and generalized precision are also considered under convex entropy loss (CEL).     

MINIMAX ESTIMATION OF A MULTIVARIATE NORMAL MEAN UNDER A CONVEX LOSS FUNCTION

Let X = (X₁, - X_p)prime; ˜ N_p (μ, Σ) where μ= (μ₁, -, μ_p)' and Σ= diag (Σ²₁, -, Σ²_p) are both unknown and p3. Let (n_i - 2) w_i/Σ²_i! X²_ni, independent. of w_i (I ≠ j = 1, -, p). Assume that (w₁, -, w_p) and X are independent. Define W = diag (w₁, -, w_p) and ¶ X ¶²_w= X'W^-1Q^-1W^-1X where Q = diag (q₁, -,n q_p), q_i > 0, i = 1, -, p. In this paper, the minimax estimator of Berger & Bock (1976), given by δ (X, W) = [I_p - r(X, W) ¶ X ¶^-2_w Q^-1W^-1] X, is shown to be minimax relative to the convex loss (δ - μ)'[αQ + (1 - α) Σ^-1] δ - μ)/C, where C =α tr (Σ) + (1 - α)p and 0 α 1, under certain conditions on r(X, W). This generalizes the above mentioned result of Berger & Bock.