test_bind.R
In simsem: SIMulated Structural Equation Modeling

library(simsem)
library(testthat)
## Needs more thorough tests, but this will work for now.
## Mostly just checks to see if the function works.


#context("bind - SimMatrix")
a <- matrix(0,3,3)
a[,1] <- NA
a[,2] <- "a1"
mat <- bind(free=a)
expect_true(class(mat) == "SimMatrix")
expect_true(all(mat@free[,2]=="a1"))
expect_true(all(is.na(mat@free)[,1]))
expect_true(all(is.nan(mat@popParam)))
expect_true(all(is.nan(mat@misspec)))

mat <- bind(free=a, popParam=.7, misspec=.01)
expect_true(class(mat) == "SimMatrix")
expect_true(all(mat@popParam[,1:2]==.7))
expect_true(all(mat@misspec[,3]==.01))

mat <- bind(free=a, popParam="runif(1,0,1)", misspec=.01)
expect_true(class(mat) == "SimMatrix")
expect_true(all(mat@popParam[,1:2]=="runif(1,0,1)"))
expect_true(all(mat@popParam[,3]==""))
expect_true(all(mat@misspec[,1:2]==""))

mat <- bind(free=a, popParam="runif(1,0,1)", misspec="runif(1,0,1)")
expect_true(class(mat) == "SimMatrix")
expect_true(all(mat@misspec[,3]=="runif(1,0,1)"))

pop <- matrix("runif(1,0,1)",3,3)
#expect_error(bind(a,pop))
  #Error - invalid expression
expect_error(bind(free=a, popParam="runif(1,0,1)", misspec="runif(1,0,1"))
expect_error(bind(free=a, popParam="runif(1,0,1", misspec="runif(1,0,1)"))


  ## Doesn't mean anything, but doesn't throw an error?
expect_true(class(bind(free=a, popParam=1)) == "SimMatrix")

  ## Error - different dimensions
expect_error(bind(free=a, popParam=matrix(0,2,2)))
expect_error(bind(free=a, misspec=matrix(0,2,2)))

#context("Symmetric Tests")
a <- matrix(NA,3,3)
a[upper.tri(a)] <- 0
expect_error(bind(free=a, symmetric=TRUE))

a[lower.tri(a)] <- 0
expect_true(class(bind(free=a, symmetric=TRUE)) == "SimMatrix")
expect_true(class(bind(free=a,popParam=.7,symmetric=TRUE)) == "SimMatrix")
expect_true(class(bind(free=a,misspec=.7,symmetric=TRUE)) == "SimMatrix")

a[upper.tri(a)] <- a[lower.tri(a)] <- "a1"
expect_true(class(bind(free=a, symmetric=TRUE)) == "SimMatrix")
expect_true(class(bind(free=a, popParam=.7, symmetric=TRUE)) == "SimMatrix")

a[lower.tri(a)] <- NA
expect_error(bind(a,.7,symmetric=TRUE))

a[lower.tri(a)] <- a[upper.tri(a)]
b <- matrix(0,3,3)
b[1,2] <- 7
expect_error(bind(a,b,symmetric=TRUE))
expect_error(bind(a,.7,b,symmetric=TRUE))



#context("Equality Constraints")
  a <- matrix(0,2,2)
  a[,1] <- 0
  a[,2] <- "a1"
expect_true(class(bind(free=a)) == "SimMatrix")

  b <- matrix(0,3,3)
  b[1:2,1] <- "b1"
  b[2:3,2] <- "b2"
expect_true(class(bind(free=b)) == "SimMatrix")

expect_true(class(bind(free=matrix(0,2,2))) == "SimMatrix")


###################################

#context("bind - SimVector")

a <- c(NA,NA,"a1","a1")
expect_true(class(bind(free=a)) == "SimVector")
expect_true(class(bind(free=a, popParam=.7, misspec=.01)) == "SimVector")
expect_true(class(bind(free=a, popParam="runif(1,0,1)", misspec=.01)) == "SimVector")
expect_true(class(bind(free=a, popParam="runif(1,0,1)", misspec="runif(1,0,1)")) == "SimVector")

  #Error - invalid expression
expect_error(bind(free=a, popParam="runif(1,0,1)", misspec="runif(1,0,1"))
expect_error(bind(free=a, popParam="runif(1,0,1", misspec="runif(1,0,1)"))

  ## Does throw an error in this case.
expect_error(class(bind(free=a, popParam="a1")) == "SimVector")

  ## Error - different dimensions
expect_error(bind(free=a, popParam=c(0,3,3)))
expect_error(bind(free=a, misspec=c(0,3,3)))

# context("Equality Constraints")

a <- c(0,0,"a1","a1")
expect_true(class(bind(free=a)) == "SimVector")

b <- c(0,0,"b1","b1",0,0,"b2","b2")
expect_true(class(bind(free=b)) == "SimVector")

expect_true(class(bind(free=c(0,2,2))) == "SimVector")

expect_error(bind(a,symmetric=TRUE))

# These are the kinds of tests that are actually useful
a <- bind(rep(NA, 3), 0, rep("runif(1,.01,.015)",3))
expect_equal(a@misspec[[3]],"runif(1,.01,.015)")
expect_equal(a@popParam[[3]],0)
expect_equal(a@free[[3]],NA)