context("Influence functions, robust standard errors")
fnorm <-
test_that("phreg.iid", {
data(TRACE)
out1 <- phreg(Surv(time,status==9)~vf+chf, data=TRACE)
expect_true(norm(estimate(out1)$vcov - solve(-out1$hessian))<1)
})
test_that("twinlm.iid", {
set.seed(1)
d <- twinsim(100,b1=c(1,-1),b2=c(),acde=c(1,1,0,1))
l <- twinlm(y ~ 1, data=d, DZ="DZ", zyg="zyg", id="id")
expect_true(norm(estimate(l)$vcov - vcov(l))<1)
})
test_that("biprobitid", {
set.seed(1)
d <- twinsim(100,b1=c(1,-1),b2=c(),acde=c(1,1,0,1))
l <- biprobit(y0 ~ x1, id="id", rho=~x1, data=d)
expect_true(norm(estimate(l)$vcov - vcov(l))<1)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.