library(MethComp)
data( rainman )
str( rainman )
RM <- Meth( rainman, item=1, y=2:6 )
head( RM )
BA.est( RM, linked=FALSE )
library(lme4)
mf <- lmer( y ~ meth + item + (1|MI),
data = transform( RM, MI=interaction(meth,item) ) )
summary( mf )
mr <- lmer( y ~ (1|meth) + (1|item) + (1|MI),
data = transform( RM, MI=interaction(meth,item) ) )
summary( mr )
#
# Point swarms were generated by the following program
#
## Not run:
# set.seed(2) # Original
# npoints <- sample(4:30)*4
# nplots <- 10
# pdf(file="swarms.pdf", onefile=TRUE)
#
# s1 <- sample(npoints[1:nplots])
# print(s1)
# for (i in 1:nplots) {
# n <- s1[i]
# set.seed(n)
# x <- runif(n)
# y <- runif(n)
# plot(x,y, xlim=c(-.15, 1.15), ylim=c(-.15, 1.15), pch=20, axes=F,
# xlab="", ylab="")
# }
# s1 <- sample(npoints[1:nplots])
# print(s1)
# for (i in 1:nplots) {
# n <- s1[i]
# set.seed(n)
# x <- runif(n)
# y <- runif(n)
# plot(y,x, xlim=c(-.15, 1.15), ylim=c(-.15, 1.15), pch=20, axes=F,
# xlab="", ylab="")
# }
# s1 <- sample(npoints[1:nplots])
# print(s1)
# for (i in 1:nplots) {
# n <- s1[i]
# set.seed(n)
# x <- runif(n)
# y <- runif(n)
# plot(-x,y, xlim=c(-1.15, .15), ylim=c(-.15, 1.15), pch=20, axes=F,
# xlab="", ylab="")
# }
# dev.off()
# ## End(Not run)
Run the code above in your browser using DataLab