How would you write the R factor function to encode the batch and gender? [closed]