Nothing
context("daSubmodels tests")
x1<-rnorm(1000)
x2<-rnorm(1000)
x3<-rnorm(1000)
y1<-x1+x2+x3+rnorm(1000)
d.f<-data.frame(xa=x1,xb=x2,xc=x3,y=y1)
test_that("Correct submodels for standard data.frame and lm", {
pred.matrix<-matrix(
c(0,0,0,
1,0,0,
0,1,0,
0,0,1,
1,1,0,
1,0,1,
0,1,1,
1,1,1),8,3,byrow=T)
predictors<-c("xa","xb","xc")
level<-c(0,1,1,1,2,2,2,3)
response<-c("y")
constants<-c("")
lm.1<-lm(y~xa+xb+xc,d.f)
ds<-daSubmodels(lm.1)
expect_that(as.numeric(ds$pred.matrix),equals(as.numeric(pred.matrix)))
expect_that(ds$predictors,equals(predictors))
expect_that(ds$response,equals(response))
expect_that(ds$level,equals(level))
})
test_that("Correct submodels for standard data.frame, lm and manual terms", {
lm.1<-lm(y~xa+xb+xc, d.f)
pred.matrix<-matrix(
c(0,0,
1,0,
0,1,
1,1),4,2,byrow=T)
predictors<-c("xa+xb","xc")
level<-c(0,1,1,2)
response<-c("y")
constants<-c("")
ds<-daSubmodels(lm.1,terms=c("xa+xb","xc"))
expect_that(as.numeric(ds$pred.matrix),equals(as.numeric(pred.matrix)))
expect_that(ds$predictors,equals(predictors))
expect_that(ds$response,equals(response))
expect_that(ds$level,equals(level))
})
test_that("Correct submodels for standard data.frame and lmer", {
require(lme4)
g<-factor(gl(50,20))
y1<-x1+x2+x3+rnorm(1000)
d.f<-data.frame(xa=x1,xb=x2,xc=x3,y=y1)
pred.matrix<-matrix(
c(0,0,0,
1,0,0,
0,1,0,
0,0,1,
1,1,0,
1,0,1,
0,1,1,
1,1,1),8,3,byrow=T)
predictors<-c("xa","xb","xc")
level<-c(0,1,1,1,2,2,2,3)
response<-c("y")
constants<-c("( 1 | g )")
lm.1<-lmer(y~xa+xb+xc+(1|g),d.f)
ds<-daSubmodels(lm.1)
expect_that(as.numeric(ds$pred.matrix),equals(as.numeric(pred.matrix)))
expect_that(ds$predictors,equals(predictors))
expect_that(ds$response,equals(response))
expect_that(ds$constants,equals(constants))
expect_that(ds$level,equals(level))
})
test_that("Correct submodels for standard data.frame with constants", {
x4<-rnorm(1000)
y1<-x1+x2+x3+x4+rnorm(1000)
d.f<-data.frame(xa=x1,xb=x2,xc=x3,xd=x4,y=y1)
pred.matrix<-matrix(
c(0,0,0,
1,0,0,
0,1,0,
0,0,1,
1,1,0,
1,0,1,
0,1,1,
1,1,1),8,3,byrow=T)
predictors<-c("xa","xb","xc")
response<-c("y")
constants<-c("xd")
level<-c(0,1,1,1,2,2,2,3)
lm.1<-lm(y~xa+xb+xc+xd,d.f)
ds<-daSubmodels(lm.1,constants="xd")
expect_that(as.numeric(ds$pred.matrix),equals(as.numeric(pred.matrix)))
expect_that(ds$predictors,equals(predictors))
expect_that(ds$response,equals(response))
expect_that(ds$constants,equals(constants))
expect_that(ds$level,equals(level))
})
test_that("should raise an error if a constant is not available in predictors" ,{
y1<-x1+x2+x3+rnorm(1000)
d.f<-data.frame(xa=x1,xb=x2,xc=x3,y=y1)
lm.1<-lm(y~xa+xb+xc,d.f)
expect_error(daSubmodels(lm.1,constants="xe"),"all constants should be in predictors")
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.