func_ci_single_Q_bias.m

function log_val = func_ci_single_Q_bias(n_pt, n_s, X, h, lambda, z, prior_mu, prior_nu, prior_mu_h, prior_nu_h, prior_a, prior_b)

log_val = 0;

for j = 1:n_s
    % prior z
    log_val = log_val + 0.5*log(prior_nu(j)) - 0.5*prior_nu(j)*(z(j) - prior_mu(j))^2;
    for i = 1:n_pt
        if ~isnan(X(i,j))
            % obj x
            log_val = log_val + 0.5*log(lambda(i)) - 0.5*lambda(i)*(X(i,j) - z(j) - h(i))^2;
        end
    end
end

for i = 1:n_pt
    % prior h
%     log_val = log_val + 0.5*log(prior_nu_h(i)) - 0.5*prior_nu_h(i)*(h(i) - prior_mu_h(i))^2;
    % prior lambda
    log_val = log_val + (prior_a(i)-1)*log(lambda(i)) - prior_b(i)*lambda(i); 
end