Usage
gdsidr(y, q, s, family, vmu="v", varht=NULL, limnla=c(-10, 3),
maxit=30, job=-1, tol1=0, tol2=0, prec=1e-06)
Arguments
y
a numerical vector representing the response, or a matrix of two columns for binomial data with the first column as the largest possible counts and the second column as the counts actually obsered.
q
a square matrix of the same order as the length of y, with elements equal to the reproducing kernel evaluated at the design points.
s
the design matrix of the null space $H_0$ of size (length-of-y,dim($H_0$)), with elements equal to the bases of $H_0$ evaluated at design points.
family
a string specifying the family of distribution. Families supported are "binary", "binomial", "poisson" and "gamma" for Bernoulli, binomial, poisson, and gamma distributions respectively. Canonical links are used except for Gamma family where a log link
vmu
a character string specifying a method for choosing the smoothing parameter. "v", "m" and "u" represent GCV, GML and UBR respectively. "u$\sim$", only used for non-Gaussian family, specifies UBR with estimated variance. Default is "v".
varht
needed only when vmu="u", which gives the fixed variance in calculation of the UBR function. Default is 1.0.
limnla
a vector of length 2, specifying a search range for the n times smoothing parameter on log10 scale. Default is (-10, 3).
maxit
maximum number of iterations allowed for the iteration in GRKPACK.
job
an integer representing the optimization method used to find the smoothing parameter.
The options are job=-1: golden-section search on (limnla(1), limnla(2));
job=0: golden-section search with interval specified automatically;
job >0: regular grid sea
tol1
the tolerance for elements of w's. Default is 0.0 which sets to square of machine precision.
tol2
tolerance for truncation used in `dsidr'. Default is 0.0 which sets to square of machine precision.
prec
precision requested for stopping the iteration. Default is $1e-06$.