Consider a discrete-time insurance risk model. Within period i, the net insurance loss is denoted by a real-valued random variable X
i
. The insurer makes both risk-free and risky investments, leading to an overall stochastic discount factor Y
i
from time i to time i − 1. Assume that (X
i
, Y
i
), i ∈ N, form a sequence of independent and identically distributed random pairs following a common bivariate Farlie-Gumbel-Morgenstern distribution with marginal distribution functions F and G. When F is subexponential and G fulfills some constraints in order for the product convolution of F and G to be subexponential too, we derive a general asymptotic formula for the finite-time ruin probability. Then, for special cases in which F belongs to the Fréchet or Weibull maximum domain of attraction, we improve this general formula to be transparent.