The problem is not with the BMA function, but rather that you have made the
common mistake of forgetting that you are smarter than the computer. The
variable Y below is a data frame (not a single numeric vector). Now a smart
person like you sees a single column in the data frame and understands that is
the y-variable, most intelligent people would do the same thing, but the
computer is not as smart as you, it sees a data frame (a form of a list, which
is what the error is talking about), knows that it wants a vector, but is not
smart enough to figure out the only column of the data frame is the vector to
use, so gives an error instead.
You can either read the data in as a vector, convert the data frame to a vector,
or use subscripting/indexing to pass only the vector to the function (or maybe a
4th alternative that I have not thought of yet).
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.snow at imail.org
801.408.8111
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of David S. Schwarz
> Sent: Thursday, May 14, 2009 12:55 PM
> To: r-help at r-project.org
> Subject: [R] Bayesian Model Averaging
>
> Hello R Group,
>
>
>
>
>
> Below is the code I submit:
>
>
>
> library("BMA")
>
> X <- read.table("C:/Documents and
> Settings/Administrator/Desktop/coding.txt",header=TRUE)
>
> Y <- read.table("C:/Documents and
> Settings/Administrator/Desktop/1DCS.txt",header=TRUE)
>
> IGout<- iBMA.glm(X, Y, glm.family= gaussian(), verbose = TRUE,
> thresProbne0
> = 5 )
>
> summary(IGout)
>
> IGout
>
>
>
> I get the following error message:
>
>
>
> Error in model.frame.default(formula = Y ~ 1, data > nastyHack_____x.df,
:
>
> invalid type (list) for variable 'Y'
>
> > summary(IGout)
>
> Error in summary(IGout) : object "IGout" not found
>
> > IGout
>
> Error: object "IGout" not found
>
>
>
> The data set X is 47 rows by 41 columns, all numeric. The vector Y is
> 47
> rows, all numeric.
>
>
>
> This closely follows an example in the help but I can't find the
> problem.
> Can anyone help?
>
>
>
> Dave Schwarz
>
>
>
>
>
> David S. Schwarz & Associates
>
> 4261 Americana Drive
>
> Suite 229
>
> Stow, OH 44224
>
> ph: 330-945-4733
>
> email: david.schwarz at goodstats.biz
>
> web site: www.goodstats.biz
>
>
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.