https://github.com/cran/fields
Tip revision: a82f141190b29d59547a12b9a44a7e2d1af8fb60 authored by Doug Nychka on 05 January 2007, 21:27:53 UTC
version 3.3.1
version 3.3.1
Tip revision: a82f141
csstr.f
** finds value of h minimizing the generalized cross-validation
subroutine csstr(h,nobs,x,y,wght,c,offset,trace,vlam,work,ierr)
parameter(mxM=20000)
implicit double precision (a-h,o-z)
REAL*8 h,trace, vlam,c,offset
REAL*8 x(nobs),y(nobs),wght(nobs)
REAL*8 work(nobs),diag(mxM),dumm1(1),dumm2(1)
integer job(3),ideriv,ierr
data ideriv/0/
job(1)=3
job(2)=0
job(3)=0
diag(1)=c
diag(2)=offset
call css(h,nobs,x,y,wght,work,trace,diag,vlam,ndum,dumm1,dumm2,
- job,ideriv,ierr)
return
end