how to calculate log2 fold change in r