# Load necessary libraries
library(testthat)
library(SUMO)  # Load your package
#output_sim <- OmixCraftHD(vector_features = c(2000,3000), sigmas_vector=c(4,5), n_samples=100, n_factors=5, num.factor='multiple', advanced_dist='mixed')
 output_obj <- OmixCraftHD(
   vector_features = c(4000,3000),
   n_samples=100,
   n_factors=2,
   signal.samples = NULL,
   signal.features.one = NULL,
   signal.features.two = NULL,
   snr = 2.5,
   num.factor='multiple',
   advanced_dist='mixed')
test_that("plot_simData runs without error", {
  # Test that the plot function runs without errors
  expect_silent(plot_simData(sim_object = output_obj, data = "merged", type = "heatmap"))
})