一些代码在SIBRE的捆绑演示数据集上执行此操作 .
在这个例子中,我们尝试使用ggplot2创建后椭圆的多个样本的一些图 .
library(SIBER)
library(ggplot2)
library(dplyr)
library(ellipse)
使基本SIBER模型适合与包捆绑在一起的示例数据 .
# load in the included demonstration dataset
data("demo.siber.data")
#
# create the siber object
siber.example
# Calculate summary statistics for each group: TA, SEA and SEAc
group.ML
# options for running jags
parms
parms$n.iter
parms$n.burnin
parms$n.thin
parms$n.chains
# define the priors
priors
priors$R
priors$k
priors$tau.mu
# fit the ellipses which uses an Inverse Wishart prior
# on the covariance matrix Sigma, and a vague normal prior on the
# means. Fitting is via the JAGS method.
ellipses.posterior
# The pos