Determination of the number of components in the PARAFAC model with a nonnegative tensor structure: a simulated EEG data study