suppressMessages({
  # Unit test for Monte Carlo nodes creation
  test_that("create_mcnodes works", {
    # Test that function fails without mctable parameter
    expect_error(create_mcnodes(data = imports_data))

    # Create nodes with both required parameters
    create_mcnodes(data = imports_data, mctable = imports_mctable)

    # Verify dimensions of h_prev match expected values
    expect_equal(dim(h_prev), c(ndvar(), 1, nrow(imports_data)))

    # Compare automatic vs manual node creation
    automatic_node <- round(mean(extractvar(h_prev)), 2)
    manual_node <- round(mean(mcstoc(runif, min = imports_data$h_prev_min[1], max = imports_data$h_prev_max[1])), 2)

    # Check if both methods produce same result
    expect_equal(manual_node, automatic_node)

    # Test node creation after setting mctable
    set_mctable(imports_mctable)
    expect_no_error(create_mcnodes(data = imports_data))
    reset_mctable()
  })
  
  test_that("create_mcnodes handles out-of-order columns correctly", {
    # Create test data with columns in alphabetical order (different from rpert parameter order)
    # rpert expects: min, mode, max
    # Alphabetical order: max, min, mode
    test_data <- data.frame(
      n_animals_max = c(100, 120),
      n_animals_min = c(50, 60),
      n_animals_mode = c(75, 90)
    )
    
    # Create mctable for rpert distribution
    test_mctable <- data.frame(
      mcnode = "n_animals",
      description = "Number of animals",
      mc_func = "rpert",
      from_variable = NA,
      transformation = NA,
      sensi_analysis = FALSE
    )
    
    # Create environment for testing
    test_env <- new.env()
    
    # Create the mcnode - should not error
    expect_no_error(
      create_mcnodes(data = test_data, mctable = test_mctable, envir = test_env)
    )
    
    # Verify the node was created
    expect_true(exists("n_animals", envir = test_env))
    
    # Verify the created node has correct dimensions
    expect_equal(dim(test_env$n_animals), c(ndvar(), 1, nrow(test_data)))
    
    # Verify the mcnode does not contain NAs
    expect_false(any(is.na(test_env$n_animals)))
  })
})