score_test_statistic function

Calculate ML estimates