"blup" procedure added to scan1 and scan1.threshold