Nothing
require(MASS)
library(MASS)
data(UScrime)
#UScrime[,-2] = log(UScrime[,-2])
crime.bic = bas.lm(log(y) ~ log(M) + So + log(Ed) + log(Po1) + log(Po2)
+ log(LF) + log(M.F) + log(Pop) + log(NW) +
log(U1) + log(U2) + log(GDP) + log(Ineq) + log(Prob)+
log(Time),
data=UScrime, n.models=2^15, prior="BIC",
modelprior=beta.binomial(1,1),
initprobs= "eplogp")
summary(crime.bic)
plot(crime.bic)
image(crime.bic, subset=-1)
# takes a while to run:
# crime.coef = coefficients(crime.bic)
# crime.coef
# par(mfrow=c(3,2))
# plot(crime.coef, ask=FALSE)
# see update
#crime.aic = update(crime.bic, newprior="AIC")
#crime.zs = update(crime.bic, newprior="ZS-null")
#crime.EBG = EB.global.bma(crime.bic)
# same as update(crime.bic, newprior="EB-global")
#image(crime.EBG, subset=-1)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.