func[?_?NumericQ, ?_?NumericQ, ?_?NumericQ, t_?NumericQ] :=
NIntegrate[Exp[t*y]*((E^(-y ?) (1 - E^(-y ?))^(-1 - ?) (-1 + (1 - E^(-y ?))^-?)^(-1 + ?) ? ? ?)/(1 + (-1 +
(1 - E^(-y ?))^-?)^?)^2), {y, 0, ?}, Method -> "AdaptiveQuasiMonteCarlo",
Exclusions -> {(1 + (-1 + (1 - E^(-y ?))^-?)^?)^2 == 0}]
func[2, 2, 1, 1]
(* 2.66667 *)
func[1/2, 2, 1/3, 1](* Probably integral is divergent !!! *)
(* Error-> Overflow, Indeterminate, or Infinity for all sampling points *)
Regards MI.