anova.ssr: Testing a Non-parametric Function Fitted via Smoothing...

View source: R/assist.R

anova.ssrR Documentation

Testing a Non-parametric Function Fitted via Smoothing Splines

Description

For smoothing spline models with a single smoothing parameter, test the hypothesis that the unknown funciton lies in the null space using the local most powerful (LMP) test and a GCV or GML test.

Usage

## S3 method for class 'ssr'
anova(object, simu.size=100, ...)

Arguments

object

an object of class "ssr" fitted with a single smooting parameter.

simu.size

an optional integer giving the number of simulations to calcualte p-values based on simulation. Default is 100.

...

other available arguments, currently unused.

Details

For Gaussian data with one smoothing parameter, test the hypothesis that the function is in the null space H_0, i.e. the parametric part of the fitted model is sufficient. Available are the LMP and GCV or GML methods. However, the p-values cannot be calculated analytically since the null distributions for these testing statistics under H_0 are unknown. Monte Carlo simulation is used to approximate the p-values for the LMP, and GCV (if spar="v") or GML (if spar="m") methods. Due to computation burden, the smoothing parameters are fixed at their estimate in the currect calculation.

When spar="m", an approximate p-value based on a mixture of two Chi-square distributions is also provided for the GML test, which tends to be conservationve (Pinherio and Bates, 2002).

Methods further developed in Liu and Wang (2004) and Liu, Meiring and Wang (2004) will be implemented in the future.

Value

a list including test values.

Author(s)

Chunlei Ke chunlei_ke@yahoo.com and Yuedong Wang yuedong@pstat.ucsb.edu

References

Cox, D. and Koh, E. (1989). A smoothing spline based test of model adequency in polynomial regression. Ann. Ins. Stat. Math. 41, 383-400.

Cox, D., Koh, E., Wahba, G. and Yandell, B.S. (1988). Testing the parameteric null model hypothesis in semi-parametric partial and generalized spline models. Ann. Statist. 16, 113-119.

Wahba, G. (1990). Spline Models for Observational Data. SIAM, Vol. 59.

Pinherio, J. C. and Bates, D. M. (2000) Mixed-effects Models in S and S-Plus. Springer.

Liu, A. and Wang, Y. (2004) Hypothesis Testing in Smoothing Spline Models. Journal of Statistical Computation and Simulation, to appear.

Liu, A., Meiring, W. and Wang, Y. (2004), Testing Generalized Linear Models Using Smoothing Spline Methods. Statistica Sinica, to appear,

See Also

ssr, print.anova.ssr

Examples

## Not run: 
data(acid)

# fit a partial thin-plate spline
temp <- ssr(ph~t1+x1+x2, rk=tp(t1), data=acid, spar="m")
anova(temp, 500)

## End(Not run)

assist documentation built on Aug. 22, 2023, 9:12 a.m.

Related to anova.ssr in assist...