theano have imcomplete gamma function and log-gamma implemented, you can probably write the cdf of gamma distribution in theano directly.
FYI, you might also find some inspiration of inferencing copula in this WIP notebook: https://github.com/junpenglao/Planet_Sakaar_Data_Science/blob/master/WIP/vector_transformation_copulas.ipynb