analyses_focals.R