PFsample3.R