Calculates ERDS values in the specified frequency band. This function calculates ERDS values in the specified frequency band and returns, in addition to the ERDS values, the power in the reference interval and a confidence interval. Usage: [erds, refp, cl, cu] = getErds(s, h, triallen, f, trls, class, rmartif, ref, submean); Input parameters: s ... Input signal (as obtained by sload) h ... Header structure (as obtained by sload) triallen ... Length of one trial (s) f ... Frequency band (Hz) Optional input parameters: trls ... Trials class ... Class rmartif ... Remove artifacts ref ... Reference interval (s) submean ... Subtract mean signal alpha ... Significance level Output parameters: erds ... ERDS values of each channel refp ... Reference power of each channel cl ... Lower confidence level cu ... Upper confidence level