[R] Bayesian regression stepwise function?
Charles C. Berry
cberry at tajo.ucsd.edu
Fri Oct 23 07:17:32 CEST 2009
On Thu, 22 Oct 2009, Ben Bolker wrote:
> Allan.Y wrote:
>> Hi everyone,
>> I am wondering if there exists a stepwise regression function for the
>> Bayesian regression model. I tried googling, but I couldn't find
>> anything. I know "step" function exists for regular stepwise regression,
>> but nothing for Bayes.
> Why? That seems so ... un-Bayesian ...
If 'fools rush in where angels fear to tread', then Bayesians 'jump' in
where frequentists fear to 'step'...
Seriously, there are Bayesian regression approaches that priorize the
model size (sometimes only implicitly by assigning a prior for the
inclusion of each candidate regressors). Then they 'jump' between models
of different sizes.
On CRAN, Package qtlbim (which is specialized to a particular genetics
problem) implements one such, I think.
Package bqtl does not implement the jumping approach, but does explore a
model space with differing numbers of regressors for the same (qtl)
Perhaps the closest to a general purpose 'stepwise flavored' Bayesian
regression is implemented in Package BMA, which IIRC borrows step() for
some of its work.
But CRAN now has more packages than my cortex has neurons, so there are
probably more packages that do something like this. Try
RSiteSearch("jump regression", restric='functions')
and start reading.
> View this message in context: http://www.nabble.com/Bayesian-regression-stepwise-function--tp26013725p26015081.html
> Sent from the R help mailing list archive at Nabble.com.
> R-help at r-project.org mailing list
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
Charles C. Berry (858) 534-2098
Dept of Family/Preventive Medicine
E mailto:cberry at tajo.ucsd.edu UC San Diego
http://famprevmed.ucsd.edu/faculty/cberry/ La Jolla, San Diego 92093-0901
More information about the R-help