Update imputation.R

Add bpca method, used as default