Results 1  10
of
28
Smoothing Spline ANOVA for Exponential Families, with Application to the Wisconsin Epidemiological Study of Diabetic Retinopathy
 ANN. STATIST
, 1995
"... Let y i ; i = 1; \Delta \Delta \Delta ; n be independent observations with the density of y i of the form h(y i ; f i ) = exp[y i f i \Gammab(f i )+c(y i )], where b and c are given functions and b is twice continuously differentiable and bounded away from 0. Let f i = f(t(i)), where t = (t 1 ; \De ..."
Abstract

Cited by 101 (46 self)
 Add to MetaCart
Let y i ; i = 1; \Delta \Delta \Delta ; n be independent observations with the density of y i of the form h(y i ; f i ) = exp[y i f i \Gammab(f i )+c(y i )], where b and c are given functions and b is twice continuously differentiable and bounded away from 0. Let f i = f(t(i)), where t = (t 1 ; \Delta \Delta \Delta ; t d ) 2 T (1)\Omega \Delta \Delta \Delta\Omega T (d) = T , the T (ff) are measureable spaces of rather general form, and f is an unknown function on T with some assumed `smoothness' properties. Given fy i ; t(i); i = 1; \Delta \Delta \Delta ; ng, it is desired to estimate f(t) for t in some region of interest contained in T . We develop the fitting of smoothing spline ANOVA models to this data of the form f(t) = C + P ff f ff (t ff ) + P ff!fi f fffi (t ff ; t fi ) + \Delta \Delta \Delta. The components of the decomposition satisfy side conditions which generalize the usual side conditions for parametric ANOVA. The estimate of f is obtained as the minimizer...
Finding Chaos in Noisy Systems
, 1991
"... In the past twenty years there has been much interest in the physical and biological sciences in nonlinear dynamical systems that appear to have random, unpredictable behavior. One important parameter of a dynamic system is the dominant Lyapunov exponent (LE). When the behavior of the system is comp ..."
Abstract

Cited by 70 (2 self)
 Add to MetaCart
In the past twenty years there has been much interest in the physical and biological sciences in nonlinear dynamical systems that appear to have random, unpredictable behavior. One important parameter of a dynamic system is the dominant Lyapunov exponent (LE). When the behavior of the system is compared for two similar initial conditions, this exponent is related to the rate at which the subsequent trajectories diverge. A bounded system with a positive LE is one operational definition of chaotic behavior. Most methods for determining the LE have assumed thousands of observations generated from carefully controlled physical experiments. Less attention has been given to estimating the LE for biological and economic systems that are subjected to random perturbations and observed over a limited amount of time. Using nonparametric regression techniques (Neural Networks and Thin Plate Splines) it is possible to consistently estimate the LE. The properties of these methods have been studied using simulated data and are applied to a biological time series: marten fur returns for the Hudson Bay Company (18201900). Based on a nonparametric analysis there is little evidence for lowdimensional chaos in these data. Although these methods appear to work well for systems perturbed by small amounts of noise, finding chaos in a system with a significant stochastic component may be difficult.
Bayesian Smoothing and Regression Splines for Measurement Error Problems
 Journal of the American Statistical Association
, 2001
"... In the presence of covariate measurement error, estimating a regression function nonparametrically is extremely dicult, the problem being related to deconvolution. Various frequentist approaches exist for this problem, but to date there has been no Bayesian treatment. In this paper we describe Bayes ..."
Abstract

Cited by 40 (8 self)
 Add to MetaCart
In the presence of covariate measurement error, estimating a regression function nonparametrically is extremely dicult, the problem being related to deconvolution. Various frequentist approaches exist for this problem, but to date there has been no Bayesian treatment. In this paper we describe Bayesian approaches to modeling a exible regression function when the predictor variable is measured with error. The regression function is modeled with smoothing splines and regression P{splines. Two methods are described for exploration of the posterior. The rst is called iterative conditional modes (ICM) and is only partially Bayesian. ICM uses a componentwise maximization routine to nd the mode of the posterior. It also serves to create starting values for the second method, which is fully Bayesian and uses Markov chain Monte Carlo techniques to generate observations from the joint posterior distribution. Using the MCMC approach has the advantage that interval estimates that directly model and adjust for the measurement error are easily calculated. We provide simulations with several nonlinear regression functions and provide an illustrative example. Our simulations indicate that the frequentist mean squared error properties of the fully Bayesian method are better than those of ICM and also of previously proposed frequentist methods, at least in the examples we have studied. KEY WORDS: Bayesian methods; Eciency; Errors in variables; Functional method; Generalized linear models; Kernel regression; Measurement error; Nonparametric regression; P{splines; Regression Splines; SIMEX; Smoothing Splines; Structural modeling. Short title. Nonparametric Regression with Measurement Error Author Aliations Scott M. Berry (Email: scott@berryconsultants.com) is Statistical Scientist,...
Bootstrap Confidence Intervals for Smoothing Splines and their Comparison to Bayesian `Confidence Intervals'
 J. Statist. Comput. Simulation
, 1994
"... We construct bootstrap confidence intervals for smoothing spline and smoothing spline ANOVA estimates based on Gaussian data, and penalized likelihood smoothing spline estimates based on data from exponential families. Several variations of bootstrap confidence intervals are considered and compared. ..."
Abstract

Cited by 22 (6 self)
 Add to MetaCart
We construct bootstrap confidence intervals for smoothing spline and smoothing spline ANOVA estimates based on Gaussian data, and penalized likelihood smoothing spline estimates based on data from exponential families. Several variations of bootstrap confidence intervals are considered and compared. We find that the commonly used bootstrap percentile intervals are inferior to the T intervals and to intervals based on bootstrap estimation of mean squared errors. The best variations of the bootstrap confidence intervals behave similar to the well known Bayesian confidence intervals. These bootstrap confidence intervals have an average coverage probability across the function being estimated, as opposed to a pointwise property. Keywords: BAYESIAN CONFIDENCE INTERVALS, BOOTSTRAP CONFIDENCE INTERVALS, PENALIZED LOG LIKELIHOOD ESTIMATES, SMOOTHING SPLINES, SMOOTHING SPLINE ANOVA'S. 1 Introduction Smoothing splines and smoothing spline ANOVAs (SS ANOVAs) have been used successfully in a bro...
Behavior near zero of the distribution of GCV smoothing parameter estimates for splines
 Statistics and Probability Letters
, 1993
"... It has been noticed by several authors that there is a small but nonzero probability that the GCV estimate 2 of the smoothing parameter in spline and related smoothing problems will he extremely small, leading to gross undersmoothing. We obtain an upper bound to the probability that the GCV functio ..."
Abstract

Cited by 13 (7 self)
 Add to MetaCart
It has been noticed by several authors that there is a small but nonzero probability that the GCV estimate 2 of the smoothing parameter in spline and related smoothing problems will he extremely small, leading to gross undersmoothing. We obtain an upper bound to the probability that the GCV function, whose minimizer provides,~, has a (possibly local) minimum at 0. This upper bound goes to 0 exponentially fast as the sample size gets large. For the mediumto smallsample case we study this probability both by Monte Carlo evaluation of a formula for the exact probability that the GCV function has a minimum at 0 as well as by replicated calculations of ~..
Confidence Intervals for Nonparametric Curve Estimates Based on Local Smoothing
 J. Am. Stat. Assoc
, 1998
"... Numerous nonparametric regression methods exist which yield consistent estimators of function curves. Often one is also interested in constructing confidence intervals for the unknown function. Pointwise confidence intervals are available using globally crossvalidated smoothing spline (GCV) estim ..."
Abstract

Cited by 11 (0 self)
 Add to MetaCart
Numerous nonparametric regression methods exist which yield consistent estimators of function curves. Often one is also interested in constructing confidence intervals for the unknown function. Pointwise confidence intervals are available using globally crossvalidated smoothing spline (GCV) estimation. When the function estimate is based on a single global smoothing parameter the resulting confidence intervals may hold their desired confidence level 1 \Gamma ff on average but because bias in nonparametric estimation is not uniform, they do not hold the desired level uniformly at all design points. To deal with this problem, a new smoothing spline estimator is developed which uses a local crossvalidation (LCV) criterion to determine a separate smoothing parameter for each design point. The local smoothing parameters are then used to compute the point estimators of the regression curve and the corresponding pointwise confidence intervals. Incorporation of local information th...
Testing generalized linear models using smoothing spline methods, Statistica Sinica 15
, 2005
"... Abstract: This article considers testing the hypothesis of Generalized Linear Models (GLM) versus general smoothing spline models for data from exponential families. The tests developed are based on the connection between smoothing spline models and Bayesian models (Gu (1992)). They are extensions o ..."
Abstract

Cited by 4 (2 self)
 Add to MetaCart
Abstract: This article considers testing the hypothesis of Generalized Linear Models (GLM) versus general smoothing spline models for data from exponential families. The tests developed are based on the connection between smoothing spline models and Bayesian models (Gu (1992)). They are extensions of the locally most powerful (LMP) test of Cox, Koh, Wahba and Yandell (1988), the generalized maximum likelihood ratio (GML) test and the generalized cross validation (GCV) test of Wahba (1990) for Gaussian data. Null distribution approximations are considered and simulations are done to evaluate these approximations. Simulations show that the LMP and GML tests are more powerful for low frequency functions while the GCV test is more powerful for high frequency functions, which is also true for Gaussian data (Liu and Wang (2004)). The tests are applied to data from the Wisconsin Epidemiology Study of Diabetic Retinopathy, the results of which conrm and provide more denite analysis than those of previous studies. The good performances of the tests make them useful tools for diagnosis of GLM. Key words and phrases: Diagnosis, generalized cross validation, generalized maximum likelihood, locally most powerful test, hypothesis test, reproducing kernel hilbert space. 1.
Backfitting in smoothing spline ANOVA, with application to historical global temperature data
, 1996
"... In the attempt to estimate the temperature history of the earth using the surface observations, various biases can exist. An important source of bias is the incompleteness of sampling over both time and space. There have been a few methods proposed to deal with this problem. Although they can correc ..."
Abstract

Cited by 3 (2 self)
 Add to MetaCart
In the attempt to estimate the temperature history of the earth using the surface observations, various biases can exist. An important source of bias is the incompleteness of sampling over both time and space. There have been a few methods proposed to deal with this problem. Although they can correct some biases resulting from incomplete sampling, they have ignored some other significant biases. In this dissertation, a smoothing spline ANOVA approach which is a multivariate function estimation method is proposed to deal simultaneously with various biases resulting from incomplete sampling. Besides that, an advantage of this method is that we can get various components of the estimated temperature history with a limited amount of information stored. This method can also be used for detecting erroneous observations in the data base. The method is illustrated through an example of modeling winter surface air temperature as a function of year and location. Extension to more complicated mod...
Quantitative Study of Smoothing SplineANOVA Based Fingerprint Methods for Attribution of Global Warming
, 1999
"... A fingerprintbased method for climate change detection and attribution with some novel features is proposed. The method is based on a functional ANOVA (ANalysis Of VAriance) decomposition of a time and space signal, further decomposed into global timetrend and timetrend anomaly as a function of s ..."
Abstract

Cited by 3 (1 self)
 Add to MetaCart
A fingerprintbased method for climate change detection and attribution with some novel features is proposed. The method is based on a functional ANOVA (ANalysis Of VAriance) decomposition of a time and space signal, further decomposed into global timetrend and timetrend anomaly as a function of space. The method estimates the signal as a component of forced minus background climate model output, and then uses a partial spline model to estimate and test for the existence of signal in historical data. The method is based on the classical detection of signal in noise, however there are several features apparently novel to the fingerprint literature, in particular, the analysis takes place directly in observation space, anomalies are tted directly and there is possibility for estimating certain parameters of covariance models for the historical data as part of the analysis. Simulation studies using climate model runs from GFDL and NCAR and historical data for NH Winter average...
Parameter estimation in a coupled system of nonlinear sizestructured populations. Mathematical biosciences and engineering
 MBE
"... A leastsquares technique is developed for identifying unknown parameters in a coupled system of nonlinear sizestructured populations. Convergence results for the parameter estimation technique are established. Ample numerical simulations and statistical evidence are provided to demonstrate the fea ..."
Abstract

Cited by 3 (0 self)
 Add to MetaCart
(Show Context)
A leastsquares technique is developed for identifying unknown parameters in a coupled system of nonlinear sizestructured populations. Convergence results for the parameter estimation technique are established. Ample numerical simulations and statistical evidence are provided to demonstrate the feasibility of this approach. 1 1