# NOT RUN {
#Note: number.samples = 50 below.
# However, please use a value of 1000 or higher
# View top few rows of goggles data set
# from Discovering Statistics Using R
set.seed(1)
head(album)
# Single block example
blk1 <- lm(sales ~ adverts + airplay, data=album)
apa.reg.boot.table(blk1)
apa.reg.boot.table(blk1,filename="exRegTable.doc")
# Two block example, more than two blocks can be used
blk1 <- lm(sales ~ adverts, data=album)
blk2 <- lm(sales ~ adverts + airplay + attract, data=album)
apa.reg.boot.table(blk1,blk2,filename="exRegBlocksTable.doc")
# Interaction product-term test with blocks
blk1 <- lm(sales ~ adverts + airplay, data=album)
blk2 <- lm(sales ~ adverts + airplay + I(adverts * airplay), data=album)
apa.reg.boot.table(blk1,blk2,filename="exInteraction1.doc")
# }
Run the code above in your browser using DataLab