vkt <- (p[1]/(1 - (1 - p[1]) * fft(fj)))^r[1] * (p[2]/(1 - (1 - p[2]) * fft(fi)))^r[2] ...