Dear all,

I could find the formula explaining how to calculate McFadden R-squared in the gretl user's guide.
I tried to replicate it in Probit model but I fail to get the same value. Is anything wrong with my calculation as follows.
Thanks

Yi-Nung Yang

<script>
nulldata 10
set seed 89675430 
series u=normal()
series y=(u>0)
series x = uniform()

probit y const x
scalar lnL=$lnl

probit y const
scalar lnL0=$lnl

# McFadden R-squared as in Greene's Econometric Analysis
scalar McR2= 1-lnL/lnL0

</script>