predict.shim: Make predictions from a shim object
In sahirbhatnagar/shim: High Dimensional Interactions with an Environment Variable

this function only works for tuning parameter values defined by the shim_multiple_faster function. The interpolation feature is not working yet

Get coefficients from a "shim" object

## S3 method for class 'shim'
predict(object, newx, s = NULL, type = c("link", "response",
  "coefficients", "nonzero", "class"))

## S3 method for class 'shim'
coef(object, s = NULL)

`object`	Fitted shim model object
`s`	index of tuning parameter. Must be a character and an element of "s1","s2",...."s100", where "s100" is the index of the last pair of tuning parameters. Default is `NULL`
`type`	Type of prediction required. Type "link" gives the fitted values. Type "response" for "gaussian" type "response" is equivalent to type "link". Type "coefficients" computes the coefficients at the requested values for s. Type "nonzero" returns a list of the indices of the nonzero coefficients for each value of s.