% retrieving one volume of data
niivol = niidata(:, :, :, 1581);

% setting one volume of data
niidata(:, :, :, 2711) = newvol;