mstdesign <- "
    B1 =~ c(i1, i2, i3, i4, i5)
    B2 =~ c(i6, i7, i8, i9, i10)
    B3 =~ c(i11, i12, i13, i14, i15)
    B4 =~ c(i16, i17, i18, i19, i20)
    B5 =~ c(i21, i22, i23, i24, i25)
    B6 =~ c(i26, i27, i28, i29, i30)

    # define branches
    b1 := B4(0,2) + B2(0,2) + B1(0,5)
    b2 := B4(0,2) + B2(3,5) + B3(0,5)
    b3 := B4(3,5) + B5(0,2) + B3(0,5)
    b4 := B4(3,5) + B5(3,5) + B6(0,5)
  "

items <- seq(-2,2,length.out=30)
persons <- 1000

dat <- tmt:::sim.rm(persons,items,1111)
datna <- dat
datna[sample(seq_len(length(datna)),50,replace = FALSE)] <- NA
datrm_1 <- tmt_rm(dat, optimization="optim")


names(items) <- c(paste0("i",1:30))
dat_mst <- tmt_sim(mstdesign = mstdesign,
			items = items,
			persons = persons,
			seed = 1111)
datrm_1p <- tmt_rm(dat, optimization="optim")
datrm_2p <- tmt_rm(dat_mst$data,mstdesign=mstdesign, optimization="optim")

datrm_3p <- tmt_rm(dat, sum0 = FALSE)
datrm_4p <- tmt_rm(dat_mst$data,mstdesign=mstdesign, sum0 = FALSE)

datrm_5p <- tmt_rm(dat, se=FALSE, optimization="optim")
datrm_6p <- tmt_rm(dat_mst$data, se=FALSE,mstdesign=mstdesign, optimization="optim")


datlrt_1p <- tmt_lrtest(datrm_1p, optimization="optim")
datlrt_2p <- tmt_lrtest(datrm_2p, optimization="optim")

datrm_1p.out <- capture.output(summary(datrm_1p))
datrm_2p.out <- capture.output(summary(datrm_2p))
datrm_3p.out <- capture.output(summary(datrm_3p))
datrm_4p.out <- capture.output(summary(datrm_4p))
datrm_5p.out <- capture.output(summary(datrm_5p))
datrm_6p.out <- capture.output(summary(datrm_6p))

datlrt_1p.out <- capture.output(summary(datlrt_1p))
datlrt_2p.out <- capture.output(summary(datlrt_2p))

# -----------------------------------------------------------------
context("test-summary")
# -----------------------------------------------------------------
test_that("print functions raschmodel nmst", {
  expect_that(datrm_1p.out[6], equals(c("Results of Rasch model (nmst) estimation: ")))
  expect_that(datrm_1p.out[8], equals(c("Difficulty parameters: ")))
  expect_that(datrm_3p.out[6], equals(c("Results of Rasch model (nmst) estimation: ")))
  expect_that(datrm_3p.out[8], equals(c("Difficulty parameters: ")))
  expect_that(datrm_5p.out[8], equals(c("Difficulty parameters: ")))
  expect_that(grep("Std. Error", datrm_4p.out[11]), equals(1))
})

test_that("print functions raschmodel mst", {
  expect_that(datrm_2p.out[6], equals(c("Results of Rasch model (mst) estimation: ")))
  expect_that(datrm_2p.out[8], equals(c("Difficulty parameters: ")))
  expect_that(datrm_4p.out[6], equals(c("Results of Rasch model (mst) estimation: ")))
  expect_that(datrm_6p.out[9], equals(c("Difficulty parameters: ")))
  expect_that(grep("Std. Error", datrm_2p.out[11]), equals(1))
})


# -----------------------------------------------------------------
context("test-print functions lrtest")
# -----------------------------------------------------------------
test_that("print functions lrtest", {
  expect_that(datlrt_1p.out[2], 
  	equals(c("Likelihood ratio test (Andersen):")))
  expect_that(datlrt_2p.out[2], 
  	equals(c("Likelihood ratio test (Andersen) for multistage designs:")))
})