https://github.com/cran/fields
Tip revision: 6769ffc81115fbf0bf7d9c566cf7ac81be0049dc authored by Doug Nychka on 25 July 2005, 00:00:00 UTC
version 3.04
version 3.04
Tip revision: 6769ffc
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