A bivariate distribution whose marginal are gamma and beta prime distribution is introduced. The distribution is derived and the generation of such bivariate sample is shown. Extension of the results are given in the multivariate case under a joint independent component analysis method. Simulated applications are given and they show consistency of our approach. Estimation procedures for the bivariate case are provided.