add mask and label, get pvals in R