Running examples in 'MCMC.qpcr-Ex.R' failed
The error most likely occurred in:

> ### Name: HPDsummary
> ### Title: Summarizes and plots results of mcmc.qpcr function series.
> ### Aliases: HPDsummary
> 
> ### ** Examples
> 
> 
> data(beckham.data)
> data(beckham.eff)
> 
> # analysing the first 5 genes 
> # (to try it with all 10 genes, change the line below to gcol=4:13)
> gcol=4:8 
> ccol=1:3 # columns containing experimental conditions
> 
> # recalculating into molecule counts, reformatting
> qs=cq2counts(data=beckham.data,genecols=gcol,
+ condcols=ccol,effic=beckham.eff,Cq1=37)
> 
> # creating a single factor, 'treatment.time', out of 'tr' and 'time'
> qs$treatment.time=as.factor(paste(qs$tr,qs$time,sep="."))
> 
> # fitting a naive model
> naive=mcmc.qpcr(
+ 	fixed="treatment.time",
+ 	data=qs,
+ 	nitt=3000,burnin=2000 # remove this line in actual analysis!
+ )
$PRIOR
$PRIOR$R
$PRIOR$R$V
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

$PRIOR$R$nu
[1] 4.002


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

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




$FIXED
[1] "count~0+gene++gene:treatment.time"

$RANDOM
[1] "~sample"


                       MCMC iteration = 0

 Acceptance ratio for liability set 1 = 0.000167

 Acceptance ratio for liability set 2 = 0.000425

 Acceptance ratio for liability set 3 = 0.000410

 Acceptance ratio for liability set 4 = 0.000436

 Acceptance ratio for liability set 5 = 0.000333

                       MCMC iteration = 1000

 Acceptance ratio for liability set 1 = 0.007310

 Acceptance ratio for liability set 2 = 0.026825

 Acceptance ratio for liability set 3 = 0.034256

 Acceptance ratio for liability set 4 = 0.119667

 Acceptance ratio for liability set 5 = 0.037381

                       MCMC iteration = 2000

 Acceptance ratio for liability set 1 = 0.005119

 Acceptance ratio for liability set 2 = 0.039875

 Acceptance ratio for liability set 3 = 0.048333

 Acceptance ratio for liability set 4 = 0.186205

 Acceptance ratio for liability set 5 = 0.056738

                       MCMC iteration = 3000

 Acceptance ratio for liability set 1 = 0.006119

 Acceptance ratio for liability set 2 = 0.046525

 Acceptance ratio for liability set 3 = 0.052897

 Acceptance ratio for liability set 4 = 0.202179

 Acceptance ratio for liability set 5 = 0.064643
> 
> #summary plot of inferred abundances
> # s1=HPDsummary(model=naive,data=qs)
> 
> #summary plot of fold-changes relative to the global control
> s0=HPDsummary(model=naive,data=qs,relative=TRUE)
Error in update_tibble_attrs(x, ...) : 
  object 'tibble_update_attrs' not found
Calls: HPDsummary ... recycle_columns -> new_tibble -> update_tibble_attrs
Execution halted