By linear-response U approach, It can be found using the bare and self-consistent response functions from linear regression of the relationship between the occupations that result from a range of rigid potential shifts, alpha.
http://hjkgrp.mit.edu/content/calculating-hubbard-u