I have a set of genes with expression data and p-values , but i want to know how to calculate Fold change(FC) values for these genes in R or MATLAB without assigning normal and tumor data in any function like mavolcanoplot. It would be better with any kind of statistical programming for calculating FC value for all.