Running examples in ‘MCMC.OTU-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: MCMC.OTU-package
> ### Title: Bayesian analysis of multivariate counts data in DNA
> ###   metabarcoding and ecology.
> ### Aliases: MCMC.OTU-package MCMC.OTU
> 
> ### ** Examples
> 
> 
> # Symbiodinium sp diversity in two coral species at two reefs (banks)
> data(green.data)
> 
> # removing outliers
> goods=purgeOutliers(
+ 	data=green.data,
+ 	count.columns=c(4:length(green.data[1,])),
+ 	zero.cut=0.25 # remove this line for real analysis
+ )
[1] "samples with counts below z-score -2.5 :"
[1] "EFAV153"
[1] "zscores:"
  EFAV153 
-3.884085 
[1] "OTUs passing frequency cutoff  0.001 : 10"
[1] "OTUs with counts in 0.25 of samples:"

FALSE  TRUE 
    3     7 
> 
> # stacking the data table
> gs=otuStack(
+ 	data=goods,
+ 	count.columns=c(4:length(goods[1,])),
+ 	condition.columns=c(1:3)
+ 	)
> 
> # fitting the model
> mm=mcmc.otu(
+ 	fixed="bank+species+bank:species",
+ 	data=gs,
+ 	nitt=3000,burnin=2000 # remove this line for real analysis!
+ 	)
$PRIOR
$PRIOR$R
$PRIOR$R$V
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]    1    0    0    0    0    0    0    0
[2,]    0    1    0    0    0    0    0    0
[3,]    0    0    1    0    0    0    0    0
[4,]    0    0    0    1    0    0    0    0
[5,]    0    0    0    0    1    0    0    0
[6,]    0    0    0    0    0    1    0    0
[7,]    0    0    0    0    0    0    1    0
[8,]    0    0    0    0    0    0    0    1

$PRIOR$R$nu
[1] 7.02


$PRIOR$G
$PRIOR$G$G1
$PRIOR$G$G1$V
[1] 1

$PRIOR$G$G1$nu
[1] 0




$FIXED
[1] "count~otu+bank+species+bank:species+otu:bank+otu:species+otu:bank:species"

$RANDOM
[1] "~sample"


                       MCMC iteration = 0

 Acceptance ratio for liability set 1 = 0.000123

 Acceptance ratio for liability set 2 = 0.000298

 Acceptance ratio for liability set 3 = 0.000175

 Acceptance ratio for liability set 4 = 0.000246

 Acceptance ratio for liability set 5 = 0.000333

 Acceptance ratio for liability set 6 = 0.000596

 Acceptance ratio for liability set 7 = 0.000333

 Acceptance ratio for liability set 8 = 0.000719

                       MCMC iteration = 1000

 Acceptance ratio for liability set 1 = 0.207895

 Acceptance ratio for liability set 2 = 0.432649

 Acceptance ratio for liability set 3 = 0.219211

 Acceptance ratio for liability set 4 = 0.434351

 Acceptance ratio for liability set 5 = 0.345684

 Acceptance ratio for liability set 6 = 0.424140

 Acceptance ratio for liability set 7 = 0.367088

 Acceptance ratio for liability set 8 = 0.514702

                       MCMC iteration = 2000

 Acceptance ratio for liability set 1 = 0.295263

 Acceptance ratio for liability set 2 = 0.425228

 Acceptance ratio for liability set 3 = 0.304000

 Acceptance ratio for liability set 4 = 0.457351

 Acceptance ratio for liability set 5 = 0.377842

 Acceptance ratio for liability set 6 = 0.422105

 Acceptance ratio for liability set 7 = 0.373526

 Acceptance ratio for liability set 8 = 0.518614

                       MCMC iteration = 3000

 Acceptance ratio for liability set 1 = 0.326211

 Acceptance ratio for liability set 2 = 0.384561

 Acceptance ratio for liability set 3 = 0.303105

 Acceptance ratio for liability set 4 = 0.493246

 Acceptance ratio for liability set 5 = 0.347789

 Acceptance ratio for liability set 6 = 0.433667

 Acceptance ratio for liability set 7 = 0.363596

 Acceptance ratio for liability set 8 = 0.477526
> 
> # selecting the OTUs that were modeled reliably
> acpass=otuByAutocorr(mm,gs)
> 
> # calculating effect sizes and p-values:
> ss=OTUsummary(mm,gs,summ.plot=FALSE)
> 
> # correcting for mutliple comparisons (FDR)
> ss=padjustOTU(ss)
> 
> # getting significatly changing OTUs (FDR<0.05)
> sigs=signifOTU(ss)
> 
> # plotting them
> ss2=OTUsummary(mm,gs,otus=sigs)
Error in Summary.unit(c(4.4, 4.4), na.rm = FALSE) : 
  VECTOR_ELT() can only be applied to a 'list', not a 'double'
Calls: OTUsummary ... render_strips -> build_strip -> ggstrip -> Summary.unit
Execution halted