joint likelihood function formula