Error when performing the Wilcoxon test on the peak counts about liger HOT 1 OPEN

# You can include as many other datasets to use here
peakDatasets <- c("atac") 
# Extract `ligerDataset` class objects requested
datasetObjList <- datasets(int.bmmc)[peakDatasets]
# Extract the peak matrices
# Make sure these datasets are of "atac" modality so it has the slots, otherwise it'll show you an error
rawPeakList <- lapply(datasetObjList, slot, "rawPeak")
peakLiger <- createLiger(rawPeakList)
peakLiger <- normalize(peakLiger)

# Transfer the cluster labels from the original object to the new temporary object for marker calling
peakLiger$cluster <- cellMeta(int.bmmc, "leiden_cluster", useDatasets = "atac")
defaultCluster(peakLiger) <- "cluster"

# It should work this time to run marker detection now
markerTable <-  runMarkerDEG(peakLiger)