def validation_step(self, batch: rlt.PreprocessedRankingInput,
batch_idx: int):
seq2slate_net = self.seq2slate_net
assert seq2slate_net.training is False
logged_slate_rank_prob = torch.exp(
seq2slate_net(batch, mode=Seq2SlateMode.PER_SEQ_LOG_PROB_MODE).
log_probs.detach().flatten().cpu())
ranked_slate_output = seq2slate_net(batch,
Seq2SlateMode.RANK_MODE,
greedy=True)
ranked_slate_rank_prob = ranked_slate_output.ranked_per_seq_probs.cpu()
self.reporter.log(
logged_slate_rank_probs=logged_slate_rank_prob,
ranked_slate_rank_probs=ranked_slate_rank_prob,
)
if not self.calc_cpe:
return
edp_g = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
batch,
eval_greedy=True,
)
edp_ng = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
batch,
eval_greedy=False,
)
return edp_g, edp_ng
def test_seq2slate_eval_data_page(self):
"""
Create 3 slate ranking logs and evaluate using Direct Method, Inverse
Propensity Scores, and Doubly Robust.
The logs are as follows:
state: [1, 0, 0], [0, 1, 0], [0, 0, 1]
indices in logged slates: [3, 2], [3, 2], [3, 2]
model output indices: [2, 3], [3, 2], [2, 3]
logged reward: 4, 5, 7
logged propensities: 0.2, 0.5, 0.4
predicted rewards on logged slates: 2, 4, 6
predicted rewards on model outputted slates: 1, 4, 5
predicted propensities: 0.4, 0.3, 0.7
When eval_greedy=True:
Direct Method uses the predicted rewards on model outputted slates.
Thus the result is expected to be (1 + 4 + 5) / 3
Inverse Propensity Scores would scale the reward by 1.0 / logged propensities
whenever the model output slate matches with the logged slate.
Since only the second log matches with the model output, the IPS result
is expected to be 5 / 0.5 / 3
Doubly Robust is the sum of the direct method result and propensity-scaled
reward difference; the latter is defined as:
1.0 / logged_propensities * (logged reward - predicted reward on logged slate)
* Indicator(model slate == logged slate)
Since only the second logged slate matches with the model outputted slate,
the DR result is expected to be (1 + 4 + 5) / 3 + 1.0 / 0.5 * (5 - 4) / 3
When eval_greedy=False:
Only Inverse Propensity Scores would be accurate. Because it would be too
expensive to compute all possible slates' propensities and predicted rewards
for Direct Method.
The expected IPS = (0.4 / 0.2 * 4 + 0.3 / 0.5 * 5 + 0.7 / 0.4 * 7) / 3
"""
batch_size = 3
state_dim = 3
src_seq_len = 2
tgt_seq_len = 2
candidate_dim = 2
reward_net = FakeSeq2SlateRewardNetwork()
seq2slate_net = FakeSeq2SlateTransformerNet()
src_seq = torch.eye(candidate_dim).repeat(batch_size, 1, 1)
tgt_out_idx = torch.LongTensor([[3, 2], [3, 2], [3, 2]])
tgt_out_seq = src_seq[
torch.arange(batch_size).repeat_interleave(tgt_seq_len),
tgt_out_idx.flatten() - 2, ].reshape(batch_size, tgt_seq_len,
candidate_dim)
ptb = rlt.PreprocessedTrainingBatch(
training_input=rlt.PreprocessedRankingInput(
state=rlt.FeatureData(float_features=torch.eye(state_dim)),
src_seq=rlt.FeatureData(float_features=src_seq),
tgt_out_seq=rlt.FeatureData(float_features=tgt_out_seq),
src_src_mask=torch.ones(batch_size, src_seq_len, src_seq_len),
tgt_out_idx=tgt_out_idx,
tgt_out_probs=torch.tensor([0.2, 0.5, 0.4]),
slate_reward=torch.tensor([4.0, 5.0, 7.0]),
),
extras=rlt.ExtraData(
sequence_number=torch.tensor([0, 0, 0]),
mdp_id=np.array(["0", "1", "2"]),
),
)
edp = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net, reward_net, ptb.training_input, eval_greedy=True)
logger.info(
"---------- Start evaluating eval_greedy=True -----------------")
doubly_robust_estimator = OPEstimatorAdapter(DoublyRobustEstimator())
dm_estimator = OPEstimatorAdapter(DMEstimator())
ips_estimator = OPEstimatorAdapter(IPSEstimator())
switch_estimator = OPEstimatorAdapter(SwitchEstimator())
switch_dr_estimator = OPEstimatorAdapter(SwitchDREstimator())
doubly_robust = doubly_robust_estimator.estimate(edp)
inverse_propensity = ips_estimator.estimate(edp)
direct_method = dm_estimator.estimate(edp)
# Verify that Switch with low exponent is equivalent to IPS
switch_ips = switch_estimator.estimate(edp, exp_base=1)
# Verify that Switch with no candidates is equivalent to DM
switch_dm = switch_estimator.estimate(edp, candidates=0)
# Verify that SwitchDR with low exponent is equivalent to DR
switch_dr_dr = switch_dr_estimator.estimate(edp, exp_base=1)
# Verify that SwitchDR with no candidates is equivalent to DM
switch_dr_dm = switch_dr_estimator.estimate(edp, candidates=0)
logger.info(f"{direct_method}, {inverse_propensity}, {doubly_robust}")
avg_logged_reward = (4 + 5 + 7) / 3
self.assertAlmostEqual(direct_method.raw, (1 + 4 + 5) / 3, delta=1e-6)
self.assertAlmostEqual(direct_method.normalized,
direct_method.raw / avg_logged_reward,
delta=1e-6)
self.assertAlmostEqual(inverse_propensity.raw, 5 / 0.5 / 3, delta=1e-6)
self.assertAlmostEqual(
inverse_propensity.normalized,
inverse_propensity.raw / avg_logged_reward,
delta=1e-6,
)
self.assertAlmostEqual(doubly_robust.raw,
direct_method.raw + 1 / 0.5 * (5 - 4) / 3,
delta=1e-6)
self.assertAlmostEqual(doubly_robust.normalized,
doubly_robust.raw / avg_logged_reward,
delta=1e-6)
self.assertAlmostEqual(switch_ips.raw,
inverse_propensity.raw,
delta=1e-6)
self.assertAlmostEqual(switch_dm.raw, direct_method.raw, delta=1e-6)
self.assertAlmostEqual(switch_dr_dr.raw, doubly_robust.raw, delta=1e-6)
self.assertAlmostEqual(switch_dr_dm.raw, direct_method.raw, delta=1e-6)
logger.info(
"---------- Finish evaluating eval_greedy=True -----------------")
logger.info(
"---------- Start evaluating eval_greedy=False -----------------")
edp = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net, reward_net, ptb.training_input, eval_greedy=False)
doubly_robust_estimator = OPEstimatorAdapter(DoublyRobustEstimator())
dm_estimator = OPEstimatorAdapter(DMEstimator())
ips_estimator = OPEstimatorAdapter(IPSEstimator())
doubly_robust = doubly_robust_estimator.estimate(edp)
inverse_propensity = ips_estimator.estimate(edp)
direct_method = dm_estimator.estimate(edp)
self.assertAlmostEqual(
inverse_propensity.raw,
(0.4 / 0.2 * 4 + 0.3 / 0.5 * 5 + 0.7 / 0.4 * 7) / 3,
delta=1e-6,
)
self.assertAlmostEqual(
inverse_propensity.normalized,
inverse_propensity.raw / avg_logged_reward,
delta=1e-6,
)
logger.info(
"---------- Finish evaluating eval_greedy=False -----------------")
def evaluate(self, eval_tdp: PreprocessedTrainingBatch) -> None:
seq2slate_net = self.trainer.seq2slate_net
seq2slate_net_prev_mode = seq2slate_net.training
seq2slate_net.eval()
logged_slate_rank_prob = torch.exp(
seq2slate_net(eval_tdp.training_input,
mode=Seq2SlateMode.PER_SEQ_LOG_PROB_MODE).log_probs.
detach().flatten().cpu())
eval_baseline_loss = torch.tensor([0.0]).reshape(1)
if self.trainer.baseline_net:
baseline_net = self.trainer.baseline_net
# pyre-fixme[16]: `Optional` has no attribute `training`.
baseline_net_prev_mode = baseline_net.training
# pyre-fixme[16]: `Optional` has no attribute `eval`.
baseline_net.eval()
# pyre-fixme[29]: `Optional[reagent.models.seq2slate.BaselineNet]` is
# not a function.
b = baseline_net(eval_tdp.training_input).detach()
eval_baseline_loss = (F.mse_loss(
b, eval_tdp.training_input.slate_reward).cpu().reshape(1))
# pyre-fixme[16]: `Optional` has no attribute `train`.
baseline_net.train(baseline_net_prev_mode)
else:
b = torch.zeros_like(eval_tdp.training_input.slate_reward)
eval_advantage = (
# pyre-fixme[16]: `Optional` has no attribute `__sub__`.
(eval_tdp.training_input.slate_reward - b).flatten().cpu())
ranked_slate_output = seq2slate_net(eval_tdp.training_input,
Seq2SlateMode.RANK_MODE,
greedy=True)
ranked_slate_rank_prob = torch.prod(
torch.gather(
ranked_slate_output.ranked_tgt_out_probs,
2,
ranked_slate_output.ranked_tgt_out_idx.unsqueeze(-1),
).squeeze(),
-1,
).cpu()
seq2slate_net.train(seq2slate_net_prev_mode)
if not self.calc_cpe:
return
edp_g = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
# pyre-fixme[6]: Expected `Module` for 2nd param but got
# `Optional[nn.Module]`.
self.reward_network,
eval_tdp.training_input,
eval_greedy=True,
)
if self.eval_data_pages_g is None:
self.eval_data_pages_g = edp_g
else:
# pyre-fixme[16]: `Optional` has no attribute `append`.
self.eval_data_pages_g = self.eval_data_pages_g.append(edp_g)
edp_ng = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
# pyre-fixme[6]: Expected `Module` for 2nd param but got
# `Optional[nn.Module]`.
self.reward_network,
eval_tdp.training_input,
eval_greedy=False,
)
if self.eval_data_pages_ng is None:
self.eval_data_pages_ng = edp_ng
else:
self.eval_data_pages_ng = self.eval_data_pages_ng.append(edp_ng)
# pyre-fixme[16]: `RankingPolicyGradientEvaluator` has no attribute
# `notify_observers`.
self.notify_observers(
eval_baseline_loss=eval_baseline_loss,
eval_advantages=eval_advantage,
logged_slate_rank_probs=logged_slate_rank_prob,
ranked_slate_rank_probs=ranked_slate_rank_prob,
)
def evaluate(self, eval_tdp: PreprocessedTrainingBatch) -> None:
seq2slate_net = self.trainer.seq2slate_net
seq2slate_net_prev_mode = seq2slate_net.training
seq2slate_net.eval()
logged_slate_log_prob = (seq2slate_net(
eval_tdp.training_input, mode=Seq2SlateMode.PER_SEQ_LOG_PROB_MODE).
log_probs.detach().flatten().cpu().numpy())
if self.trainer.baseline_net:
baseline_net = self.trainer.baseline_net
# pyre-fixme[16]: `Optional` has no attribute `training`.
baseline_net_prev_mode = baseline_net.training
# pyre-fixme[16]: `Optional` has no attribute `eval`.
baseline_net.eval()
# pyre-fixme[29]: `Optional[reagent.models.seq2slate.BaselineNet]` is
# not a function.
b = baseline_net(eval_tdp.training_input).detach()
self.baseline_loss.append(
F.mse_loss(b, eval_tdp.training_input.slate_reward).item())
# pyre-fixme[16]: `Optional` has no attribute `train`.
baseline_net.train(baseline_net_prev_mode)
else:
b = torch.zeros_like(eval_tdp.training_input.slate_reward)
self.baseline_loss.append(0.0)
advantage = (eval_tdp.training_input.slate_reward -
b).flatten().cpu().numpy()
self.advantages.append(advantage)
self.logged_slate_log_probs.append(logged_slate_log_prob)
ranked_slate_output = seq2slate_net(eval_tdp.training_input,
Seq2SlateMode.RANK_MODE,
greedy=True)
ranked_slate_prob = (torch.prod(
torch.gather(
ranked_slate_output.ranked_tgt_out_probs,
2,
ranked_slate_output.ranked_tgt_out_idx.unsqueeze(-1),
).squeeze(),
-1,
).cpu().numpy())
self.ranked_slate_probs.append(ranked_slate_prob)
seq2slate_net.train(seq2slate_net_prev_mode)
if not self.calc_cpe:
return
edp_g = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
# pyre-fixme[6]: Expected `Module` for 2nd param but got
# `Optional[nn.Module]`.
self.reward_network,
eval_tdp.training_input,
eval_greedy=True,
)
if self.eval_data_pages_g is None:
self.eval_data_pages_g = edp_g
else:
# pyre-fixme[16]: `Optional` has no attribute `append`.
self.eval_data_pages_g = self.eval_data_pages_g.append(edp_g)
edp_ng = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
# pyre-fixme[6]: Expected `Module` for 2nd param but got
# `Optional[nn.Module]`.
self.reward_network,
eval_tdp.training_input,
eval_greedy=False,
)
if self.eval_data_pages_ng is None:
self.eval_data_pages_ng = edp_ng
else:
self.eval_data_pages_ng = self.eval_data_pages_ng.append(edp_ng)
def evaluate(self, eval_tdp: PreprocessedRankingInput) -> None:
seq2slate_net = self.trainer.seq2slate_net
seq2slate_net_prev_mode = seq2slate_net.training
seq2slate_net.eval()
logged_slate_rank_prob = torch.exp(
seq2slate_net(eval_tdp, mode=Seq2SlateMode.PER_SEQ_LOG_PROB_MODE).
log_probs.detach().flatten().cpu())
eval_baseline_loss = torch.tensor([0.0]).reshape(1)
if self.trainer.baseline_net:
baseline_net = self.trainer.baseline_net
# pyre-fixme[16]: `Optional` has no attribute `training`.
baseline_net_prev_mode = baseline_net.training
# pyre-fixme[16]: `Optional` has no attribute `eval`.
baseline_net.eval()
# pyre-fixme[29]: `Optional[reagent.models.seq2slate.BaselineNet]` is
# not a function.
b = baseline_net(eval_tdp).detach()
eval_baseline_loss = F.mse_loss(
b, eval_tdp.slate_reward).cpu().reshape(1)
# pyre-fixme[16]: `Optional` has no attribute `train`.
baseline_net.train(baseline_net_prev_mode)
else:
b = torch.zeros_like(eval_tdp.slate_reward)
eval_advantage = (
# pyre-fixme[58]: `-` is not supported for operand types
# `Optional[torch.Tensor]` and `Any`.
(eval_tdp.slate_reward - b).flatten().cpu())
ranked_slate_output = seq2slate_net(eval_tdp,
Seq2SlateMode.RANK_MODE,
greedy=True)
ranked_slate_rank_prob = ranked_slate_output.ranked_per_seq_probs.cpu()
seq2slate_net.train(seq2slate_net_prev_mode)
if not self.calc_cpe:
return
edp_g = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
eval_tdp,
eval_greedy=True,
)
if self.eval_data_pages_g is None:
self.eval_data_pages_g = edp_g
else:
# pyre-fixme[16]: `Optional` has no attribute `append`.
self.eval_data_pages_g = self.eval_data_pages_g.append(edp_g)
edp_ng = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
eval_tdp,
eval_greedy=False,
)
if self.eval_data_pages_ng is None:
self.eval_data_pages_ng = edp_ng
else:
self.eval_data_pages_ng = self.eval_data_pages_ng.append(edp_ng)
# pyre-fixme[16]: `RankingPolicyGradientEvaluator` has no attribute
# `notify_observers`.
self.notify_observers(
eval_baseline_loss=eval_baseline_loss,
eval_advantages=eval_advantage,
logged_slate_rank_probs=logged_slate_rank_prob,
ranked_slate_rank_probs=ranked_slate_rank_prob,
)
def evaluate(self, eval_tdp: PreprocessedTrainingBatch) -> None:
seq2slate_net = self.trainer.seq2slate_net
seq2slate_net_prev_mode = seq2slate_net.training
seq2slate_net.eval()
logged_slate_log_prob = (seq2slate_net(
eval_tdp.training_input, mode=Seq2SlateMode.PER_SEQ_LOG_PROB_MODE).
log_probs.detach().flatten().cpu().numpy())
if self.trainer.baseline_net:
baseline_net = self.trainer.baseline_net
baseline_net_prev_mode = baseline_net.training
baseline_net.eval()
b = baseline_net(eval_tdp.training_input).detach()
self.baseline_loss.append(
F.mse_loss(b, eval_tdp.training_input.slate_reward).item())
baseline_net.train(baseline_net_prev_mode)
else:
b = torch.zeros_like(eval_tdp.training_input.slate_reward)
self.baseline_loss.append(0.0)
advantage = (eval_tdp.training_input.slate_reward -
b).flatten().cpu().numpy()
self.advantages.append(advantage)
self.logged_slate_log_probs.append(logged_slate_log_prob)
ranked_slate_output = seq2slate_net(eval_tdp.training_input,
Seq2SlateMode.RANK_MODE,
greedy=True)
ranked_slate_prob = (torch.prod(
torch.gather(
ranked_slate_output.ranked_tgt_out_probs,
2,
ranked_slate_output.ranked_tgt_out_idx.unsqueeze(-1),
).squeeze(),
-1,
).cpu().numpy())
self.ranked_slate_probs.append(ranked_slate_prob)
seq2slate_net.train(seq2slate_net_prev_mode)
if not self.calc_cpe:
return
edp_g = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
eval_tdp.training_input,
eval_greedy=True,
)
if self.eval_data_pages_g is None:
self.eval_data_pages_g = edp_g
else:
self.eval_data_pages_g = self.eval_data_pages_g.append(edp_g)
edp_ng = EvaluationDataPage.create_from_tensors_seq2slate(
seq2slate_net,
self.reward_network,
eval_tdp.training_input,
eval_greedy=False,
)
if self.eval_data_pages_ng is None:
self.eval_data_pages_ng = edp_ng
else:
self.eval_data_pages_ng = self.eval_data_pages_ng.append(edp_ng)
