


generalized cross-validation. Usage: g=gcv(varargin) Author: Catherine Loader.


0001 function g=gcv(varargin) 0002 % generalized cross-validation. 0003 % 0004 % Usage: g=gcv(varargin) 0005 % 0006 % 0007 % Author: Catherine Loader. 0008 0009 fit = locfit(varargin{:}); 0010 rs = rsum(fit); 0011 0012 df0 = rs(1); 0013 df1 = rs(2); 0014 llk = rs(3); 0015 n = size(fit{1}{1},1); 0016 gc = -2*n*llk/((n-df0)^2); 0017 g = [llk df0 df1 gc]; 0018 0019 return;