code: def _init_individual(subjs, voxels, TRs):
          
          return [np.zeros((voxels[i], TRs)) for i in range(subjs)]
docs: Initializes the individual components `S_i` to empty (all zeros).
  
          Parameters
          ----------
  
          subjs : int
              The number of subjects.
  
          voxels : list of int
              A list with the number of voxels per subject.
  
          TRs : int
              The number of timepoints in the data.
  
          Returns
          -------
  
          S : list of 2D array, element i has shape=[voxels_i, timepoints]
              The individual component :math:`S_i` for each subject initialized
              to zero.