I am trying to perform anti-correlation analysis between unmatched miRNA-mRNA expression data. In the first step I am doing limma differential expression for the extraction of P-Values. In the second step, I am trying to perform inverse chi square between an up-regulated miRNA versus down-regulated differential expressed genes.

I have found R package geoR which have following arguments can be used.

dinvchisq(x, df, scale, log = FALSE)
rinvchisq(n, df, scale = 1/df)

Lets say if I have

miRNA_P_Value = 4.242297e-27

mRNA_P_Valule = 1.119698e-23

How can I find/combine these values?

Thanks in advance.