added higher-level script for loglikelihood computation