View source: R/read.bayesx.output.R
| read.bayesx.output | R Documentation |
This function automatically reads in BayesX estimation output which is stored in an output directory.
read.bayesx.output(dir, model.name = NULL)
dir |
a character string, specifies the directory file where BayesX output is stored. |
model.name |
a character string, specifies the base name of the model that should be read in,
also see the examples. If not supplied |
The function searches for model term objects in the specified directory, which are then stored in
a list. Each model term object will be of class xx.bayesx, so the generic functions
described in plot.bayesx may be applied for visualizing the results. In addition
summary statistics of the models may be printed to the R console with
summary.bayesx.
read.bayesx.output typically returns a list of class "bayesx" with the first element
containing a list with the following objects:
formula |
the STAR formula used, |
bayesx.setup |
an object of class |
bayesx.prg |
a character containing the |
bayesx.run |
details on processing with |
call |
the original function call, |
fitted.values |
the fitted values of the estimated model, |
residuals |
the residuals of the estimated model, |
effects |
a |
fixed.effects |
a |
variance |
estimation results for the variance parameter of the model, |
smooth.hyp |
a |
model.fit |
list containing additional information to be supplied to
|
Nikolaus Umlauf, Thomas Kneib, Stefan Lang, Achim Zeileis.
summary.bayesx, plot.bayesx, samples.
## Not run: ## load example data from
## package example folder
dir <- file.path(find.package("R2BayesX"), "/examples/ex01")
b <- read.bayesx.output(dir)
## some model summaries
print(b)
summary(b)
## now plot estimated effects
plot(b)
## 2nd example
dir <- file.path(find.package("R2BayesX"), "/examples/ex02")
list.files(dir)
## dir contains of 2 different
## base names
## 01 only one nonparametric effect
b <- read.bayesx.output(dir, model.name = "nonparametric")
plot(b)
## 02 only one bivariate
## nonparametric effect
b <- read.bayesx.output(dir, model.name = "surface")
plot(b)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.