cV = 'Robs';
cV = strvcat(cV,'R_Fobs');
cV = strvcat(cV,'NEER');

cS = 'E_R';
cS = strvcat(cS,'E_R_F');
cS = strvcat(cS,'E_NEER');

mcValue = zeros(size(cV,1),1);
mcValue(1,1) = oo_.SmoothedVariables.Robs(end);
mcValue(2,1) = oo_.SmoothedVariables.R_Fobs(end);

cL = 8;
H  = 40; 
B  = 10000;
ci = .80;

imcforecast('mode',cV,cS,cL,H,mcValue,B,ci);