# Create a simple 4D neuroimaging vector (10x10x10 volume with 20 timepoints)
space <- NeuroSpace(c(10,10,10,20), c(1,1,1))
vec <- NeuroVec(array(rnorm(10*10*10*20), c(10,10,10,20)), space)
# Extract time series using linear indices
ts1 <- series(vec, 1:10) # Get time series for first 10 voxels
# Extract time series using 3D coordinates
coords <- matrix(c(1,1,1, 2,2,2, 3,3,3), ncol=3, byrow=TRUE)
ts2 <- series(vec, coords) # Get time series for 3 specific voxel locations
# Extract single time series using x,y,z coordinates
ts3 <- series(vec, 5, 5, 5) # Get time series from middle voxel
Run the code above in your browser using DataLab