Extract lm.spike Residuals
Description
Get residuals from an lm.spike
object.
Usage
1 2 3 4 5 6 
Arguments
object 
An object of class 
burn 
The number of MCMC iterations in the object to be discarded as burnin. 
mean.only 
Logical. If 
... 
Unused, but present for compatibility with generic

Value
The posterior distribution (or posterior mean) of residuals from
the model object. If mean.only
is TRUE
then the return
value is the vector of residuals, otherwise the return value is a
matrix, with rows corresponding to MCMC iterations, and columns to
individual observations.
Author(s)
Steven L. Scott
See Also
lm.spike
SpikeSlabPrior
summary.lm.spike
plot.lm.spike
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16  niter < 1000
n < 100
p < 10
ngood < 3
x < cbind(1, matrix(rnorm(n * (p1)), nrow=n))
beta < rep(0, p)
good < sample(1:p, ngood)
beta[good] < rnorm(ngood)
sigma < 1
y < rnorm(n, x %*% beta, sigma)
model < lm.spike(y ~ x  1, niter=niter)
plot(model)
residuals(model)
residuals(model, mean.only = TRUE)
