共查询到20条相似文献,搜索用时 15 毫秒
1.
P. Economou 《Journal of applied statistics》2009,36(1):53-66
Biased sampling from an underlying distribution with p.d.f. f(t), t>0, implies that observations follow the weighted distribution with p.d.f. f w (t)=w(t)f(t)/E[w(T)] for a known weight function w. In particular, the function w(t)=t α has important applications, including length-biased sampling (α=1) and area-biased sampling (α=2). We first consider here the maximum likelihood estimation of the parameters of a distribution f(t) under biased sampling from a censored population in a proportional hazards frailty model where a baseline distribution (e.g. Weibull) is mixed with a continuous frailty distribution (e.g. Gamma). A right-censored observation contributes a term proportional to w(t)S(t) to the likelihood; this is not the same as S w (t), so the problem of fitting the model does not simply reduce to fitting the weighted distribution. We present results on the distribution of frailty in the weighted distribution and develop an EM algorithm for estimating the parameters of the model in the important Weibull–Gamma case. We also give results for the case where f(t) is a finite mixture distribution. Results are presented for uncensored data and for Type I right censoring. Simulation results are presented, and the methods are illustrated on a set of lifetime data. 相似文献
2.
Consider an ergodic Markov chain X(t) in continuous time with an infinitesimal matrix Q = (qij) defined on a finite state space {0, 1,…, N}. In this note, we prove that if X(t) is skip-free positive (negative, respectively), i.e., qij, = 0 for j > i+ 1 (i > j+ 1), then the transition probability pij(t) = Pr[X(t)=j | X(0) =i] can be represented as a linear combination of p0N(t) (p(m)(N0)(t)), 0 ≤ m ≤N, where f(m)(t) denotes the mth derivative of a function f(t) with f(0)(t) =f(t). If X(t) is a birth-death process, then pij(t) is represented as a linear combination of p0N(m)(t), 0 ≤m≤N - |i-j|. 相似文献
3.
This paper presents a methodology for model fitting and inference in the context of Bayesian models of the type f(Y
|
X,θ)f(X|θ)f(θ), where Y is the (set of) observed data, θ is a set of model parameters and X is an unobserved (latent) stationary stochastic process induced by the first order transition model f(X
(t+1)|X
(t),θ), where X
(t) denotes the state of the process at time (or generation) t. The crucial feature of the above type of model is that, given θ, the transition model f(X
(t+1)|X
(t),θ) is known but the distribution of the stochastic process in equilibrium, that is f(X|θ), is, except in very special cases, intractable, hence unknown. A further point to note is that the data Y has been assumed to be observed when the underlying process is in equilibrium. In other words, the data is not collected
dynamically over time.
We refer to such specification as a latent equilibrium process (LEP) model. It is motivated by problems in population genetics
(though other applications are discussed), where it is of interest to learn about parameters such as mutation and migration
rates and population sizes, given a sample of allele frequencies at one or more loci. In such problems it is natural to assume
that the distribution of the observed allele frequencies depends on the true (unobserved) population allele frequencies, whereas
the distribution of the true allele frequencies is only indirectly specified through a transition model.
As a hierarchical specification, it is natural to fit the LEP within a Bayesian framework. Fitting such models is usually
done via Markov chain Monte Carlo (MCMC). However, we demonstrate that, in the case of LEP models, implementation of MCMC
is far from straightforward. The main contribution of this paper is to provide a methodology to implement MCMC for LEP models.
We demonstrate our approach in population genetics problems with both simulated and real data sets. The resultant model fitting
is computationally intensive and thus, we also discuss parallel implementation of the procedure in special cases. 相似文献
4.
Gilles Stupfler 《Journal of nonparametric statistics》2014,26(4):721-735
Since its introduction, the pointwise asymptotic properties of the kernel estimator f?n of a probability density function f on ?d, as well as the asymptotic behaviour of its integrated errors, have been studied in great detail. Its weak convergence in functional spaces, however, is a more difficult problem. In this paper, we show that if fn(x)=(f?n(x)) and (rn) is any nonrandom sequence of positive real numbers such that rn/√n→0 then if rn(f?n?fn) converges to a Borel measurable weak limit in a weighted Lp space on ?d, with 1≤p<∞, the limit must be 0. We also provide simple conditions for proving or disproving the existence of this Borel measurable weak limit. 相似文献
5.
Let f(x) and g(x) denote two probability density functions and g(x)≠0. There are two ways to estimate the density ratio f(x)/g(x). One is to estimate f(x) and g(x) first and then the ratio, the other is to estimate f(x)/g(x) directly. In this paper, we derive asymptotic mean square errors and central limit theorems for both estimators. 相似文献
6.
Antonio Cuevas Wenceslao González-Manteiga Alberto Rodríguez-Casal 《Australian & New Zealand Journal of Statistics》2006,48(1):7-19
Given an unknown function (e.g. a probability density, a regression function, …) f and a constant c, the problem of estimating the level set L(c) ={f≥c} is considered. This problem is tackled in a very general framework, which allows f to be defined on a metric space different from . Such a degree of generality is motivated by practical considerations and, in fact, an example with astronomical data is analyzed where the domain of f is the unit sphere. A plug‐in approach is followed; that is, L(c) is estimated by Ln(c) ={fn≥c} , where fn is an estimator of f. Two results are obtained concerning consistency and convergence rates, with respect to the Hausdorff metric, of the boundaries ?Ln(c) towards ?L(c) . Also, the consistency of Ln(c) to L(c) is shown, under mild conditions, with respect to the L1 distance. Special attention is paid to the particular case of spherical data. 相似文献
7.
M. C. Jones 《统计学通讯:理论与方法》2013,42(10):1835-1843
Knowledge concerning the family of univariate continuous distributions with density function f and distribution function F defined through the relation f(x) = F α(x)(1 ? F(x))β, α, β ? , is reviewed and modestly extended. Symmetry, modality, tail behavior, order statistics, shape properties based on the mode, L-moments, and—for the first time—transformations between members of the family are the general properties considered. Fully tractable special cases include all the complementary beta distributions (including uniform, power law and cosine distributions), the logistic, exponential and Pareto distributions, the Student t distribution on 2 degrees of freedom and, newly, the distribution corresponding to α = β = 5/2. The logistic distribution is central to some of the developments of the article. 相似文献
8.
AbstractThis paper considers a partially non linear model E(Y|X, z, t) = f(X, β) + zTg(t) and gives its T-type estimate, which is a weighted quasi-likelihood estimate using sieve method and can be obtained by EM algorithm. The influence functions and asymptotic properties of T-type estimate (consistency and asymptotic normality) are discussed, and convergence rate of both parametric and non parametric components are obtained. Simulation results show the shape of influence functions and prove that the T-type estimate performs quite well. The proposed estimate is also applied to a data set and compared with the least square estimate and least absolute deviation estimate. 相似文献
9.
《随机性模型》2013,29(2):147-156
We consider a population of n individuals. Each of these individuals generates a discrete time branching stochastic process. We study the number of ancestors S(n,t) whose offspring at time t exceeds level θ(t), where θ(t) is some positive valued function. It is proved that S(n,t) may be approximated as t → ∞ and n → ∞ by some stochastic processes with independent increments. 相似文献
10.
Odile Pons 《Statistics》2013,47(5):377-388
Nonparametric estimators of the survival function S(t) = P(T ≥ t) for a partially observed time variable T have been defined by several methods, in particular, by integral self-consistency equations. The author establishes explicit expressions of the estimators in an additive form and extend this approach to several cases: a left-truncated and right-censored variable and the left-censored or left-truncated sojourn times of a right-censored semi-Markov process. These estimators are always identical to the product-limit estimators if hazard functions may be defined. 相似文献
11.
《Journal of Statistical Computation and Simulation》2012,82(1-4):341-356
This paper investigates two “non-exact” t-type tests, t( k2) and t(k2), of the individual coefficients of a linear regression model, based on two ordinary ridge estimators. The reported results are built on a simulation study covering 84 different models. For models with large standard errors, the ridge-based t-tests have correct levels with considerable gain in powers over those of the least squares t-test, t(0). For models with small standard errors, t(k1) is found to be liberal and is not safe to use while, t(k2) is found to slightly exceed the nominal level in few cases. When tie two ridge tests art: not winners, the results indicate that they don't loose much against t(0). 相似文献
12.
M. Afshari 《统计学通讯:理论与方法》2013,42(22):4156-4169
In this article, we consider a sample point (t j , s j ) including a value s j = f(t j ) at height s j and abscissa (time or location) t j . We apply wavelet decomposition by using shifts and dilations of the basic Häar transform and obtain an algorithm to analyze a signal or function f. We use this algorithm in practical to approximating function by numerical example. Some relationships between wavelets coefficients and asymptotic distribution of wavelet coefficients are investigated. At the end, we illustrate the results on simulated data by using MATLAB and R software. 相似文献
13.
We consider a stochastic fluid inventory model based on a (s, k, S) policy. The content level W = {W(t): t ≥ 0} increases or decreases according to a fluid-flow rate modulated by an n-state continuous time Markov chain (CTMC). W starts at W(0) = S; whenever W(t) drops to level s, an order is placed to take the inventory back to level S, which the supplier will carry out after an exponential leadtime. However, if during the leadtime the content level reaches k, the order is suppressed. We obtain explicit formulas for the expected discounted costs. The derivations are based on the optional sampling theorem (OST) to the multidimensional martingale and on fluid flow techniques. 相似文献
14.
A. Boyarsky 《Revue canadienne de statistique》1976,4(1):143-149
Let Wt be a one-dimensional Brownian motion on the probability space (Ω,F,P), and let dxt = a(xt)dt + b(xt)dwt, b2(x) > 0, be a one-dimensional Ito stochastic differential equation. For a(x) = a0 + a1x + … + anxn on a bounded interval we obtain a lower bound for p(t,x,y), the transition density function of the homogeneous Markov process xt, depending directly on the coefficients a0,a1, …, an, and b(x). 相似文献
15.
Jing-Jing Zhang 《统计学通讯:理论与方法》2013,42(13):3993-4017
ABSTRACTConsider the heteroscedastic partially linear errors-in-variables (EV) model yi = xiβ + g(ti) + εi, ξi = xi + μi (1 ? i ? n), where εi = σiei are random errors with mean zero, σ2i = f(ui), (xi, ti, ui) are non random design points, xi are observed with measurement errors μi. When f( · ) is known, we derive the Berry–Esseen type bounds for estimators of β and g( · ) under {ei,?1 ? i ? n} is a sequence of stationary α-mixing random variables, when f( · ) is unknown, the Berry–Esseen type bounds for estimators of β, g( · ), and f( · ) are discussed under independent errors. 相似文献
16.
P.C. Consul 《统计学通讯:理论与方法》2013,42(20):2039-2046
The class of Lagrangian probability distributions ‘LPD’, given by the expansion of a probability generating function f‘t’ under the transformation u = t/g‘t’ where g ‘t’ is also a p.g.f., has been substantially widened by removing the restriction that the defining functions g ‘t’ and f‘t’ be probability generating functions. The class of modified power series distributions defined by Gupta ‘1974’ has been shown to be a sub-class of the wider class of LPDs 相似文献
17.
18.
《随机性模型》2013,29(1):25-37
For a shot-noise process X(t) with Poisson arrival times and exponentially diminishing shocks of i.i.d. sizes, we consider the first time T b at which a given level b > 0 is exceeded. An integral equation for the joint density of T b and X(T b ) is derived and, for the case of exponential jumps, solved explicitly in terms of Laplace transforms (LTs). In the general case we determine the ordinary LT of the function ? P(T b > t) in terms of certain LTs derived from the distribution function H(x; t) = P(X(t) ≤ x), considered as a function of both variables x and t. Moreover, for G(t, u) = P(T b > t, X(t) < u), that is the joint distribution function of sup0 ≤ s ≤ t X(s) and X(t), an integro-differential equation is presented, whose unique solution is G(t, u). 相似文献
19.
Let {ξi} be an absolutely regular sequence of identically distributed random variables having common density function f(x). Let Hk(x,y) (k=1, 2,…) be a sequence of Borel-measurable functions and fn(x)=n?1(Hn(x,ξ1)+…+Hn(x,ξn)) the empirical density function. In this paper, the asymptotic property of the probability P(supx|fn(x)?f(x)|>ε) (n→∞) is studied. 相似文献
20.
Consider a planner choosing treatments for observationally identical persons who vary in their response to treatment. There are two treatments with binary outcomes. One is a status quo with known population success rate. The other is an innovation for which the data are the outcomes of an experiment. Karlin and Rubin [1956. The theory of decision procedures for distributions with monotone likelihood ratio. Ann. Math. Statist. 27, 272–299] assumed that the objective is to maximize the population success rate and showed that the admissible rules are the KR-monotone rules. These assign everyone to the status quo if the number of experimental successes is below a specified threshold and everyone to the innovation if experimental success exceeds the threshold. We assume that the objective is to maximize a concave-monotone function f(·) of the success rate and show that admissibility depends on the curvature of f(·). Let a fractional monotone rule be one where the fraction of persons assigned to the innovation weakly increases with the number of experimental successes. We show that the class of fractional monotone rules is complete if f(·) is concave and strictly monotone. Define an M-step monotone rule to be a fractional monotone rule with an interior fractional treatment assignment for no more than M consecutive values of the number of experimental successes. The M-step monotone rules form a complete class if f(·) is differentiable and has sufficiently weak curvature. Bayes rules and the minimax-regret rule depend on the curvature of the welfare function. 相似文献