Add buhlmann function to core.R