Error running ModelHierarchicalEbertPrelec

[JingyiWangCatherine]

Dear Dr. Vincent,

Thank you for creating this useful toolbox and share for free.

I am running Hierarachical logK and EbertPrelec models for my data. Though logK can run without error, but the EbertPrelec model is giving trouble. I filtered the bad the participant out according logK results, but I still get the error for EbertPrelec model below:
matjags: Calculating statistics... Error using ^
Incorrect dimensions for raising a matrix to a power. Check that the matrix is
square and the power is a scalar. To perform elementwise matrix powers, use '.^'.

Error in kde (line 77)
t_cdf=(sqrt(pi)fN)^(-2/3);

Error in calcMode (line 15)
[bandwidth,density,xmesh,cdf] = kde(x,...

Error in computeStats (line 36)
stats.mode.(var_name)(n) = calcMode(
vec(var_samples(:,:,n)) );

Error in matjags/coda2matlab (line 355)
stats = computeStats(samples,opts.dic);

Error in matjags (line 117)
[samples, stats] = coda2matlab();

Error in sampleWithMatjags (line 18)
[samples, stats] = matjags(...

Error in Model/conductInference (line 116)
obj.coda = samplerFunction(...

Error in ModelHierarchicalEbertPrelec (line 12)
obj = obj.conductInference();

Error in run_me (line 107)
model = ModelHierarchicalEbertPrelec(...

Any insights could be very helpful.

Jingyi