# Create a simple 5D dataset (10x10x10 spatial, 5 trials, 3 features)
dims <- c(10, 10, 10)
space <- NeuroSpace(c(dims, 5, 3))
# Create a sparse mask (20% of voxels)
mask_data <- array(runif(prod(dims)) < 0.2, dims)
mask <- LogicalNeuroVol(mask_data, NeuroSpace(dims))
# Generate random data for active voxels
n_voxels <- sum(mask_data)
data <- array(rnorm(3 * 5 * n_voxels), dim = c(3, 5, n_voxels)) # [features x trials x voxels]
# Create NeuroHyperVec object
hvec <- NeuroHyperVec(data, space, mask)
# Access operations
# Get data for specific voxel across all trials/features
series(hvec, 5, 5, 5)
# Extract a 3D volume for specific trial and feature
hvec[,,,2,1]
Run the code above in your browser using DataLab