# splsmu.r          smoothing cubic splines
 #
  all <- read.table("eppright.dat")
  x <- all[,1]
  y <- all[,2]
  smuspl <- smooth.spline(x,y,spar=0)
  smuspl$lambda
  smuspl$lev
  smuspl <- smooth.spline(x,y,spar=.654278)
  smuspl$lambda
  smuspl$lev
  cbind(x,y,smuspl$y)
 # clean up & go
  rm(list=ls())
  q()