src/rat/sspl00.r

subroutine sspl00(x,y,w,n,p,knot,nk,method,tol,wp,match,nef,center,
		dfoff,dfmax,cost,lambda,df,cv,gcv,coef,s,lev,
		xrange,work,ier)
integer         n,p,nk,method,ier,nef, nefp1, n2,match(1),center
double precision x(n),y(n,p),w(n),knot(nk+4),tol,wp(p),
	 	dfoff,dfmax,cost,lambda,df,cv,gcv,coef(1),s(n,p),lev(nef),
		xrange(2),work(1)
#workspace must be (2*p+2)*nefp1 + (p+16)*nk + n +p

logical center2
if(center==0){center2=.false.}
else {center2=.true.}
call sspl0(x,y,w,n,p,knot,nk,method,tol,wp,match,nef,center2,
		dfoff,dfmax,cost,lambda,df,cv,gcv,coef,s,lev,
		xrange,work,ier)
return
end

Try the mda package in your browser

Any scripts or data that you put into this service are public.

mda documentation built on Nov. 17, 2017, 4:26 a.m.