tests/testthat/test-mrfout.R

In mrf2d: Markov Random Field Models for Image Analysis

test_that("mrfout methods work", {
  # Print
  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1))
  prnt <- capture.output(print(a))
  expect_true(sum(grepl("Pseudolikelihood", prnt)) > 0)

  # Summary
  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "onepar")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")

  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "oneeach")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")

  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "absdif")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")

  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "dif")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")

  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "free")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")

  a <- fit_pl(Z_potts[1:30, 1:30], mrfi(1), family = "symmetric")
  prnt <- capture.output(summary(a))
  expect_true(sum(grepl("Interaction", prnt)) > 0)
  expect_is(plot(a), "ggplot")
})