program _ggevaluator
version 12.0
	args lnf beta sigma p
	tempvar t d eps a
	qui gen double `a' = 1/`sigma'
	qui gen double `eps' = ln($ML_y1) - `beta'
	qui gen double `t' = ln(`a') + `a'*`p'*`eps' - exp(`a'*`eps')
	qui gen double `d' = lngamma(`p')
	qui replace `lnf'=  `t' - `d' - ln($ML_y1 )
end