Causality in temporal systems: Characterization and a survey

A time series {Y_t} ‘causes’ another time series {Y_t}, in the sense defined by C.W.J. Granger, if present Y can be predicted better by using past values of X than by not doing so, other relevant information (including the past of Y) being used in either case. In this paper we (1) classify the possible causality relationships between two series X and Y, using an analogy to events in a sample space; (2) review existing work and present some new results on alternative characterizations of the more important causality events; and (3) compare several recent procedures for the empirical detection of causality.     

Estimating the conditional mode of a stationary stochastic process from noisy observations

Let {(X _i, Y _i,)}, i≥1, be a strictly stationary process from noisy observations. We examine the effect of the noise in the response Y and the covariates X on the nonparametric estimation of the conditional mode function. To estimate this function we are using deconvoluting kernel estimators. The asymptotic behavior of these estimators depends on the smoothness of the noise distribution, which is classified as either ordinary smooth or super smooth. Uniform convergence with almost sure convergence rates is established for strongly mixing stochastic processes, when the noise distribution is ordinary smooth. Received: April 1998     

Characterization of the multivariate normal distribution by conditional normal distributions

Summary LetX andY be two random vectors with values in ℝ ^k and ℝ∝, respectively. IfZ=(X ^T,Y ^T) ^T is multivariate normal thenX givenY=y andY givenX=x are (multivariate) normal; the converse is wrong. In this paper simple additional conditions are stated such that the converse is true, too. Furthermore, the case is treated that the random vectorZ=(X ₁ ^T , …,X _t ^T ) ^T is splitted intot≥3 partsX ₁, …,X _t.     

Strong uniform convergence for the estimator of the regression function under ϕ-mixing conditions

Suppose the observations (X _i,Y _i), i=1,…, n, are ϕ-mixing. The strong uniform convergence and convergence rate for the estimator of the regression function was studied by serveral authors, e.g. G. Collomb (1984), L. Gy?rfi et al. (1989). But the optimal convergence rates are not reached unless the Y _i are bounded or the E exp (a|Y _i|) are bounded for some a>0. Compared with the i.i.d. case the convergence of the Nadaraya-Watson estimator under ϕ-mixing variables needs strong moment conditions. In this paper we study the strong uniform convergence and convergence rate for the improved kernel estimator of the regression function which has been suggested by Cheng P. (1983). Compared with Theorem A in Y. P. Mack and B. Silverman (1982) or Theorem 3.3.1 in L. Gy?rfi et al. (1989), we prove the convergence for this kind of estimators under weaker moment conditions. The optimal convergence rate for the improved kernel estimator is attained under almost the same conditions of Theorem 3.3.2 in L. Gy?rfi et al. (1989). Received: September 1999     

Strong consistency of the MLE under random censoring

LetX ₁, ...,X _n be an i.i.d. sample from some parametric family {_θ :θ (Θ} of densities. In the random censorship model one observesZ _i =min (X _i ,Y _i ) andδ _i =1_{ x _i ≤Y _i}, whereY _i is a censoring variable being independent ofX _i . In this paper we investigate the strong consistency ofθ _n maximizing the modified likelihood function based on (Z _i ,δ _i , 1≤i≤n. The main result constitutes an extension of Wald’s theorem for complete data to censored data. Work partially supported by the “Deutsche Forschungsgemeinschaft”.     

Sulle operazioni associative tra variabili casuali

Si mostra che, sotto condizioni di regolarità, seo è un’operazione associativa tra variabili casuali reali e indipendenti, è definibile una trasformata integrale ξ delle loro funzioni di ripartizione con la proprietà: ξx ₀ Y (t)=ξx(t)·ξ _y (t). Si indicano alcune proprietà di tale trasformata e si tratta della possibilità di estendere a un’operazione associativa risultati noti per l’addizione tra variabili casuali. In particolare ci si occupa dell’« infinita divisibilità » fornendo condizioni perché una variabile casualeX ammetta la rappresentazioneX=X _1O X _2O … _o X _n per ognin naturale con leX _i indipendenti e identicamente ditribuite.     

On bootstrapping the mode in the nonparametric regression model with random design

In the nonparametric regression model with random design and based on i.i.d. pairs of observations (X _i, Y _i), where the regression function m is given by m(x)=?(Y _i|X _i=x), estimation of the location θ (mode) of a unique maximum of m by the location of a maximum of the Nadaraya-Watson kernel estimator for the curve m is considered. In order to obtain asymptotic confidence intervals for θ, the suitably normalized distribution of is bootstrapped in two ways: we present a paired bootstrap (PB) where resampling is done from the empirical distribution of the pairs of observations and a smoothed paired bootstrap (SPB) where the bootstrap variables are generated from a smooth bivariate density based on the pairs of observations. While the PB requires only relatively small computational effort when carried out in practice, it is shown to work only in the case of vanishing asymptotic bias, i.e. of “undersmoothing” when compared to optimal smoothing for mode estimation. On the other hand, the SPB, although causing more intricate computations, is able to capture the correct amount of bias if the pilot estimator for m oversmoothes. Received: May 2000     

Two-sample rank tests with truncated populations

LetX ₁,…,X _m andY ₁,…,Y _n be two independent samples from continuous distributionsF andG respectively. Using a Hoeffding (1951) type theorem, we obtain the distributions of the vector S=(S ₍₁₎,…,S _(n)), whereS _(j)=# (X _i ’s≤Y _(j)) andY _(j) is thej-th order statistic ofY sample, under three truncation models: (a)G is a left truncation ofF orG is a right truncation ofF, (b)F is a right truncation ofH andG is a left truncation ofH, whereH is some continuous distribution function, (c)G is a two tail truncation ofF. Exploiting the relation between S and the vectorR of the ranks of the order statistics of theY-sample in the pooled sample, we can obtain exact distributions of many rank tests. We use these to compare powers of the Hajek test (Hajek 1967), the Sidak Vondracek test (1957) and the Mann-Whitney-Wilcoxon test. We derive some order relations between the values of the probagility-functions under each model. Hence find that the tests based onS ₍₁₎ andS _(n) are the UMP rank tests for the alternative (a). We also find LMP rank tests under the alternatives (b) and (c).     

Locally minimax test of independence in elliptically symmetrical distributions with additional observations

Summary LetX=(X _ij )=(X ₁, ...,X _n )’,X’ _i =(X _i1, ...,X _ip )’,i=1,2, ...,n be a matrix having a multivariate elliptical distribution depending on a convex functionq with parameters, 0,σ. Let ϱ²=ϱ ₂ ^-2 be the squared multiple correlation coefficient between the first and the remainingp ₂+p ₃=p−1 components of eachX _i . We have considered here the problem of testingH ₀:ϱ²=0 against the alternativesH ₁:ϱ ₁ ^-2 =0, ϱ ₂ ^-2 >0 on the basis ofX andn ₁ additional observationsY ₁ (n ₁×1) on the first component,n ₂ observationsY ₂(n ₂×p ₂) on the followingp ₂ components andn ₃ additional observationsY ₃(n ₃×p ₃) on the lastp ₃ components and we have derived here the locally minimax test ofH ₀ againstH ₁ when ϱ ₂ ^-2 →0 for a givenq. This test, in general, depends on the choice ofq of the familyQ of elliptically symmetrical distributions and it is not optimality robust forQ.     

Regression analysis and dependence

In this paper we study the relationship between regression analysis and a multivariate dependency measure. If the general regression model Y=f() holds for some function f, where 1i₁< i₂<···i_m k, and X₁,...,X_k is a set of possible explanatory random variables for Y. Then there exists a dependency relation between the random variable Y and the random vector (). Using the dependency statistic defined below, we can detect such dependency even if the function f is not linear. We present several examples with real and simulated data to illustrate this assertion. We also present a way to select the appropriate subset among the random variables X₁,X₂,...,X_k, which better explain Y.     

Index

The Yule distribution is shown to have certain interesting properties in the area of regression analysis. In particular, it is shown that under certain conditions, a random variable Z will have linear regressions on another random variable X and on its observable part Y only when X has a Yule distribution. More generally, the regression on the observed part Y will be constant for a finite number of values of Y, say k, and linear otherwise, only when X has a Yule distribution with its first k frequencies truncated.     

The impact of public regulation on housing quality and on structural changes

This paper represents a theoretical investigation of profit maximizing behavior of a landlord under rent control. The situation envisioned is one that the landlord owns, free and clear, either two housing units ready to merge or one large one ready to convert into two small ones; in each case the units are already in place, such that capital costs are sunk and treated as bygones. Each unit has a technology of production of housing services with a fixed input of quantity of space (shelter) and a variable input of quality, that is affected by physical ageing (non-controllable) and by maintenance (controllable). At starting time t₀, we have a state of quality and a historical state of initial quality Q?_i, at the time the building was built, with the implicit constraint that Q_i(t) < Q?_i, for all t> t₀. The analysis addresses the general question of housing structural changes—conversion or merger—and how these changes are being accelerated under the threat of rent control.     

Bayesian estimation based on ranked set sample from Morgenstern type bivariate exponential distribution when ranking is imperfect

In this paper we consider Bayes estimation based on ranked set sample when ranking is imperfect, in which units are ranked based on measurements made on an easily and exactly measurable auxiliary variable X which is correlated with the study variable Y. Bayes estimators under squared error loss function and LINEX loss function for the mean of the study variate Y, when (X, Y) follows a Morgenstern type bivariate exponential distribution, are obtained based on both usual ranked set sample and extreme ranked set sample. Estimation procedures developed in this paper are illustrated using simulation studies and a real data.     

Identification and estimation of nonlinear dynamic panel data models with unobserved covariates

This paper considers nonparametric identification of nonlinear dynamic models for panel data with unobserved covariates. Including such unobserved covariates may control for both the individual-specific unobserved heterogeneity and the endogeneity of the explanatory variables. Without specifying the distribution of the initial condition with the unobserved variables, we show that the models are nonparametrically identified from two periods of the dependent variable _Yit$Y_{it}$ and three periods of the covariate _Xit$X_{it}$. The main identifying assumptions include high-level injectivity restrictions and require that the evolution of the observed covariates depends on the unobserved covariates but not on the lagged dependent variable. We also propose a sieve maximum likelihood estimator (MLE) and focus on two classes of nonlinear dynamic panel data models, i.e., dynamic discrete choice models and dynamic censored models. We present the asymptotic properties of the sieve MLE and investigate the finite sample properties of these sieve-based estimators through a Monte Carlo study. An intertemporal female labor force participation model is estimated as an empirical illustration using a sample from the Panel Study of Income Dynamics (PSID).     

Estimating covariation: Epps effect,microstructure noise

This paper is about how to estimate the integrated covariance _〈X,Y〉T${〈X,Y〉}_T$ of two assets over a fixed time horizon [0,T]$[0,T]$, when the observations of X$X$ and Y$Y$ are “contaminated” and when such noisy observations are at discrete, but not synchronized, times. We show that the usual previous-tick covariance estimator is biased, and the size of the bias is more pronounced for less liquid assets. This is an analytic characterization of the Epps effect. We also provide the optimal sampling frequency which balances the tradeoff between the bias and various sources of stochastic error terms, including nonsynchronous trading, microstructure noise, and time discretization. Finally, a two scales covariance estimator is provided which simultaneously cancels (to first order) the Epps effect and the effect of microstructure noise. The gain is demonstrated in data.     

Adaptive estimation of the dynamics of a discrete time stochastic volatility model

This paper is concerned with the discrete time stochastic volatility model _Yi=exp(_Xi/2)_ηi$Y_i=exp(X_i/2){\eta }_i$, _Xi+1=b(_Xi)+σ(_Xi)_ξi+1$X_{i+1}=b\left(X_i\right)+\sigma \left(X_i\right){\xi }_{i+1}$, where only (_Yi)$\left(Y_i\right)$ is observed. The model is rewritten as a particular hidden model: _Zi=_Xi+_εi$Z_i=X_i+{\epsilon }_i$, _Xi+1=b(_Xi)+σ(_Xi)_ξi+1$X_{i+1}=b\left(X_i\right)+\sigma \left(X_i\right){\xi }_{i+1}$, where (_ξi)$\left({\xi }_i\right)$ and (_εi)$\left({\epsilon }_i\right)$ are independent sequences of i.i.d. noise. Moreover, the sequences (_Xi)$\left(X_i\right)$ and (_εi)$\left({\epsilon }_i\right)$ are independent and the distribution of ε$\epsilon$ is known. Then, our aim is to estimate the functions b$b$ and σ²${\sigma }^2$ when only observations _Z1,…,_Zn$Z_1,\dots ,Z_n$ are available. We propose to estimate bf$bf$ and (b²+σ²)f$(b^2+{\sigma }^2)f$ and study the integrated mean square error of projection estimators of these functions on automatically selected projection spaces. By ratio strategy, estimators of b$b$ and σ²${\sigma }^2$ are then deduced. The mean square risk of the resulting estimators are studied and their rates are discussed. Lastly, simulation experiments are provided: constants in the penalty functions defining the estimators are calibrated and the quality of the estimators is checked on several examples.     

A bound on the expected overshoot for some concave boundaries

LetX ₁,X ₂,… be i.i.d. with finite meanμ>0,S _n =X ₁+…+X _n . Forf(n)=n ^β ,c>0 we consider the stopping timesT _c =inf{n:S _n >c+f(n)} with overshootR _c =S _{T c} −(c+f(T _c )). For 0<β<1 we give a bound for sup _c≥0 ER _c in the spirit of Lorden’s well-known inequality forf=0.     

Suffizienz bei verschiebungsinvarianten Schätzproblemen

Zusammenfassung Es wird gezeigt, daß beim Schätzen eines die Verteilung einer ZufallsgrößeX (mit Dichte) charakterisierenden Lageparameters verschiebungsinvariante FunktionenZ ₁=a ₁(X ₁,...,X _n ),...,Z _m =a _m (X ₁,...,X _n ) dern unabhängigen WiederholungenX ₁,...,X _n vonX genau dann suffizient sind, wenn für jede konvexe Schadensfunktion ein gleichmäßig bestes, nur vonZ ₁,...,Z _m abhängendes verschiebungsinvariantes Schätzverfahren existiert. Weiter wird bewiesen, daßX genau dann normalverteilt ist, wenn zu jeder konvexen Schadensfunktion ein existiert derart, daß ein gleichmäßig bestes verschiebungsinvariantes Schätzverfahren ist.

---

Summary LetX ₁,...,X _n be independent random variables with density functionf(x–) and unknown location parameter R ¹; furthermore leta _i (x ₁,...,x _n ),i=1,..., m, be functions which are invariant with respect to translations. ThenZ _i =a _i (X ₁,...,X _n ),i=1,...,m, are sufficient iff for every convex loss functions (.) there exists a functionh(z ₁,...,z _m ) such thath(Z ₁,...,Z _m ) is a best invariant estimate for the location parameter . Furthermore we show thatX ₁,...,X _n is a sample from a normal distribution if for every convex loss functions (.) there exists a constant such that is a best invariant estimate for .

     

Bayes sequential estimation for a particular exponential family of distributions under LINEX loss

The problem of sequentially estimating an unknown distribution parameter of a particular exponential family of distributions is considered under LINEX loss function for estimation error and a cost c > 0 for each of an i.i.d. sequence of potential observations X ₁, X ₂, . . . A Bayesian approach is adopted and conjugate prior distributions are assumed. Asymptotically pointwise optimal and asymptotically optimal procedures are derived.