# [Scipy-svn] r3276 - trunk/scipy/sandbox/models

scipy-svn@scip... scipy-svn@scip...
Tue Aug 28 23:19:15 CDT 2007

```Author: jonathan.taylor
Date: 2007-08-28 23:19:12 -0500 (Tue, 28 Aug 2007)
New Revision: 3276

Modified:
trunk/scipy/sandbox/models/bspline.py
Log:
references for bsplines

Modified: trunk/scipy/sandbox/models/bspline.py
===================================================================
--- trunk/scipy/sandbox/models/bspline.py	2007-08-29 04:05:33 UTC (rev 3275)
+++ trunk/scipy/sandbox/models/bspline.py	2007-08-29 04:19:12 UTC (rev 3276)
@@ -1,3 +1,20 @@
+'''
+Bspines and smoothing splines.
+
+General references:
+
+    Craven, P. and Wahba, G. (1978) "Smoothing noisy data with spline functions.
+    Estimating the correct degree of smoothing by
+    the method of generalized cross-validation."
+    Numerische Mathematik, 31(4), 377-403.
+
+    Hastie, Tibshirani and Friedman (2001). "The Elements of Statistical
+    Learning." Springer-Verlag. 536 pages.
+
+    Hutchison, M. and Hoog, F. "Smoothing noisy data with spline functions."
+    Numerische Mathematik, 47(1), 99-106.
+'''
+
import numpy as N
import numpy.linalg as L

@@ -259,7 +276,7 @@
x.shape = (1,)
x.shape = (N.product(_shape,axis=0),)
if i < self.tau.shape[0] - 1:
-            ## TODO: OWNDATA flags...
+           ## TODO: OWNDATA flags...
v = _bspline.evaluate(x, self.tau, self.m, d, i, i+1)
else:
return N.zeros(x.shape, N.float64)
@@ -499,7 +516,10 @@
"""
Generalized cross-validation score of current fit.

-        TODO: addin a reference to Wahba, and whoever else I used.
+        Craven, P. and Wahba, G.  "Smoothing noisy data with spline functions.
+        Estimating the correct degree of smoothing by
+        the method of generalized cross-validation."
+        Numerische Mathematik, 31(4), 377-403.
"""

norm_resid = (self.resid**2).sum()

```