BiRepeated: Data from a correlated bivariate beta distribution with... In zoib: Bayesian Inference for Beta Regression and Zero-or-One Inflated Beta Regression

Description

BiRepeated is a simulated data set with bi-Beta variables y1 and y2 on 200 independent cases. Both y1 and y2 are repeatedly measured at a given set of covariate values x = (0.1, 0.2, 0.3, 0.4, 0.5, 0.6). Joint zoib modeling is applied to examine the effect of x on y1 and y2 simultaneously.

Usage

 1 data(BiRepeated)

Format

A data frame with 200 independent cases, from which 6 measurements are taken on 2 response variables.

id

id of the 200 cases.

y1

value of one beta variable (6 measurements per case) ranged from 0 to 1.

y2

value of the other beta variables (6 measurements per case) ranged from 0 to 1.

x

numerical; explanatory variable.

Examples

 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 ## Not run: library(zoib) data("BiRepeated", package = "zoib") eg2 <- zoib(y1|y2 ~ x|1|x, data= BiRepeated, random=1,n.response=2, EUID= BiRepeated\$id, joint=TRUE,zero.inflation = FALSE, one.inflation = FALSE, prior.Sigma = "VC.unif", n.iter=7000,n.thin=25,n.burn=2000) coeff <- eg2\$coeff summary(coeff) ### check convergence traceplot(coeff); autocorr.plot(coeff); check.psrf(coeff) ### plot posterior mean of y vs. observed y to check on goodness of fit. n= nrow(BiRepeated) ypred1 = rbind(eg2\$ypred[][,1:n],eg2\$ypred[][,1:n]) ypred2 = rbind(eg2\$ypred[][,(n+1):(2*n)],eg2\$ypred[][,(n+1):(2*n)]) post.mean1 = apply(ypred1,2,mean); post.mean2 = apply(ypred2,2,mean); plot(BiRepeated\$y1, post.mean1, xlim=c(0,1),ylim=c(0,1), col='green2', pch=2,xlab='Observed y', ylab='Predicted y', main="") points(BiRepeated\$y2,post.mean2,col='purple') abline(0,1,col='red') legend(0.1,0.9,col=c('green2','purple'),c("y1","y2"),pch=c(2,1)) ## End(Not run)

zoib documentation built on May 29, 2021, 5:06 p.m.