Nothing
context("lmWithCov raw results")
test_that("Correct raw results using lmWithCov", {
nam<-c("SES","IQ","nAch","GPA")
cor.m<-matrix(c(
1 ,.3 ,.41 ,.33,
.3 ,1 , .16 , .57 ,
.41 , .16, 1 , .50 ,
.33 , .57, .50 , 1
),4,4,byrow=T,
dimnames=list(nam,nam)
)
lwith <- lmWithCov(GPA ~ SES + IQ + nAch, cor.m)
rawR <- daRawResults(lwith)
base.fits <-
c(
"1" = 0.0000000,
SES = 0.1089000,
"IQ" = 0.3249000,
nAch = 0.2500000,
"SES+IQ" = 0.3526813,
"SES+nAch" = 0.2687823,
"IQ+nAch" = 0.4964080,
"SES+IQ+nAch" = 0.4964735
)
expect_equal(rawR$base.fits[, 1], base.fits, tolerance = 0.0001)
})
test_that("Correct raw results using lmWithCov with manual terms", {
nam<-c("SES","IQ","nAch","GPA")
cor.m<-matrix(c(
1 ,.3 ,.41 ,.33,
.3 ,1 , .16 , .57 ,
.41 , .16, 1 , .50 ,
.33 , .57, .50 , 1
),4,4,byrow=T,
dimnames=list(nam,nam)
)
lwith <- lmWithCov(GPA ~ SES + IQ + nAch, cor.m)
rawR <- daRawResults(lwith,terms=c("SES+IQ","nAch"))
base.fits <-
c(
"1" = 0.0000000,
"SES+IQ" = 0.3526813,
"nAch" = 0.2500000,
"SES+IQ+nAch" = 0.4964735
)
expect_equal(rawR$base.fits[, 1], base.fits, tolerance = 0.0001)
expect_equal(rawR$level, c(0,1,1,2), tolerance = 0.0001)
})
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.