r>
    <tr><td class="paramdir">[in]</td><td class="paramname">priora,priorb</td><td>: prior parameters (todo: need to be made dimension-specific) for uniform, it is the range, for normal it is the moments </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">xdata</td><td>: x-values of data, nx-by-xdim </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">ydata</td><td>: y-values of data, nx-by-neach </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">xg