Spaces:
Sleeping
Sleeping
joshua400 commited on
Commit ·
dc0c6cc
1
Parent(s): f186b81
💎 HONEST TRUTH: Synchronized Space simulation with train.ipynb metrics (0.4/0.4/0.2 & MAD Fairness)
Browse files- fairrecovery_env/rewards.py +43 -69
- server/app.py +26 -21
- train.ipynb +922 -495
fairrecovery_env/rewards.py
CHANGED
|
@@ -4,58 +4,47 @@ FairRecovery++ - Reward Engine (Fair-GRPO-RLVR).
|
|
| 4 |
Computes dense, verifiable, formula-based rewards - no learned reward model.
|
| 5 |
Implements the Fair-GRPO-RLVR multi-objective reinforcement learning framework.
|
| 6 |
|
| 7 |
-
R_total =
|
| 8 |
"""
|
| 9 |
|
| 10 |
from __future__ import annotations
|
| 11 |
import structlog
|
| 12 |
from dataclasses import dataclass, field
|
| 13 |
from typing import List
|
| 14 |
-
from .constants import (GRADER_SCORE_MAX, GRADER_SCORE_MIN, MAX_DAYS
|
| 15 |
-
PENALTY_IGNORE_VULNERABLE, PENALTY_PATTERN_IGNORED,
|
| 16 |
-
REWARD_WEIGHTS, VULNERABILITY_THRESHOLD)
|
| 17 |
from .state import CityState, ZoneState
|
| 18 |
from .tasks import ScenarioConfig
|
| 19 |
|
| 20 |
logger = structlog.get_logger(__name__)
|
| 21 |
|
| 22 |
|
| 23 |
-
def compute_exec_reward(
|
| 24 |
-
"""Mean service
|
| 25 |
if not zones:
|
| 26 |
return 0.0
|
| 27 |
-
|
| 28 |
-
|
| 29 |
-
return float(0.05 + sum(improvements) / len(improvements))
|
| 30 |
|
| 31 |
|
| 32 |
def compute_fairness_reward(zones: List[ZoneState]) -> float:
|
| 33 |
"""
|
| 34 |
-
|
| 35 |
Higher value means more equitable distribution of services.
|
| 36 |
"""
|
| 37 |
if not zones:
|
| 38 |
return 0.0
|
| 39 |
services = [z.service for z in zones]
|
| 40 |
mean_svc = sum(services) / len(services)
|
| 41 |
-
|
|
|
|
|
|
|
| 42 |
# Fairness index in [0, 1]
|
| 43 |
-
return float(max(0.0, 1.0 -
|
| 44 |
|
| 45 |
|
| 46 |
def compute_safety_reward(violations: List[str]) -> float:
|
| 47 |
-
"""
|
| 48 |
-
return float(
|
| 49 |
-
|
| 50 |
-
|
| 51 |
-
def compute_stability_reward(zones: List[ZoneState]) -> float:
|
| 52 |
-
"""Reward for system balance — low variance in satisfaction."""
|
| 53 |
-
sats = [z.citizen_satisfaction for z in zones]
|
| 54 |
-
if len(sats) < 2:
|
| 55 |
-
return 0.0
|
| 56 |
-
mean_sat = sum(sats) / len(sats)
|
| 57 |
-
variance = sum((s - mean_sat) ** 2 for s in sats) / len(sats)
|
| 58 |
-
return float(max(-1.0, -variance * 4)) # Scale up variance penalty
|
| 59 |
|
| 60 |
|
| 61 |
def compute_analysis_reward(chosen_zones: List[int], zones: List[ZoneState]) -> float:
|
|
@@ -75,8 +64,6 @@ class RewardComponents:
|
|
| 75 |
R_exec: float = 0.0
|
| 76 |
R_fair: float = 0.0
|
| 77 |
R_safe: float = 0.0
|
| 78 |
-
R_adapt: float = 0.0
|
| 79 |
-
R_stable: float = 0.0
|
| 80 |
R_analysis: float = 0.0
|
| 81 |
R_total: float = 0.0
|
| 82 |
violations: List[str] = field(default_factory=list)
|
|
@@ -85,8 +72,7 @@ class RewardComponents:
|
|
| 85 |
def to_dict(self) -> dict:
|
| 86 |
return {k: round(v, 4) if isinstance(v, float) else v
|
| 87 |
for k, v in {"R_exec": self.R_exec, "R_fair": self.R_fair,
|
| 88 |
-
"R_safe": self.R_safe, "
|
| 89 |
-
"R_stable": self.R_stable, "R_total": self.R_total,
|
| 90 |
"violations": self.violations}.items()}
|
| 91 |
|
| 92 |
|
|
@@ -98,7 +84,6 @@ class RewardEngine:
|
|
| 98 |
self._cumulative_reward: float = 0.0
|
| 99 |
self._step_count: int = 0
|
| 100 |
self._action_history: List[str] = []
|
| 101 |
-
self._vulnerable_ignored_days: int = 0
|
| 102 |
|
| 103 |
@property
|
| 104 |
def cumulative_reward(self) -> float:
|
|
@@ -107,67 +92,56 @@ class RewardEngine:
|
|
| 107 |
def compute_analysis_step(self, chosen_zones: List[int], city: CityState) -> RewardComponents:
|
| 108 |
self._step_count += 1
|
| 109 |
R_analysis = compute_analysis_reward(chosen_zones, city.zones)
|
| 110 |
-
#
|
| 111 |
-
R_total = 0.
|
| 112 |
self._cumulative_reward += R_total
|
| 113 |
return RewardComponents(
|
| 114 |
R_analysis=R_analysis, R_total=R_total,
|
| 115 |
feedback=f"Analysis: {R_total:+.3f} ({int(R_analysis * max(1, len(city.zones)//2))}"
|
| 116 |
f"/{max(1, len(city.zones)//2)} critical zones correct)")
|
| 117 |
|
| 118 |
-
def compute_execute_step(self, city: CityState, violations: List[str]
|
| 119 |
-
|
| 120 |
-
"""Main dense reward after execute step — now includes adaptation and stability."""
|
| 121 |
self._step_count += 1
|
| 122 |
|
| 123 |
-
|
| 124 |
-
|
| 125 |
-
|
| 126 |
-
|
| 127 |
-
|
| 128 |
-
|
| 129 |
-
|
| 130 |
-
|
| 131 |
-
|
| 132 |
-
|
| 133 |
-
R_exec = compute_exec_reward(city.prev_services, city.zones)
|
| 134 |
-
R_fair = compute_fairness_reward(city.zones)
|
| 135 |
-
R_safe = compute_safety_reward(violations)
|
| 136 |
-
R_adapt = adaptation_score # from Predictor.evaluate_adaptation()
|
| 137 |
-
R_stable = compute_stability_reward(city.zones)
|
| 138 |
-
|
| 139 |
-
w = REWARD_WEIGHTS
|
| 140 |
-
# Fair-GRPO-RLVR Rubric (Optimised for positive feedback and strong learning signals)
|
| 141 |
-
# We use a +0.3 baseline for a successful step to ensure the baseline is clearly positive
|
| 142 |
-
baseline = 0.3 if not violations else 0.0
|
| 143 |
-
R_total = (w["exec"] * R_exec + w["fair"] * R_fair + w["safe"] * R_safe + baseline)
|
| 144 |
-
R_total = float(max(0.0 if not violations else -0.5, min(1.0, R_total)))
|
| 145 |
self._cumulative_reward += R_total
|
| 146 |
|
| 147 |
-
feedback = (f"
|
| 148 |
-
f"R_adapt={R_adapt:+.3f} | R_stable={R_stable:+.3f} → R_total={R_total:+.3f}")
|
| 149 |
if violations:
|
| 150 |
feedback += f" | Violations: {violations}"
|
| 151 |
|
| 152 |
-
return RewardComponents(R_exec=
|
| 153 |
R_total=R_total, violations=violations, feedback=feedback)
|
| 154 |
|
| 155 |
def compute_submit_reward(self, city: CityState) -> RewardComponents:
|
|
|
|
| 156 |
self._step_count += 1
|
| 157 |
-
|
| 158 |
-
|
| 159 |
-
|
| 160 |
-
|
|
|
|
| 161 |
terminal = float(max(0.0, min(1.0, terminal)))
|
| 162 |
self._cumulative_reward += terminal
|
|
|
|
| 163 |
return RewardComponents(
|
| 164 |
-
R_fair=
|
| 165 |
-
feedback=f"Terminal
|
| 166 |
|
| 167 |
def get_final_grader_score(self, city: CityState) -> float:
|
| 168 |
"""Normalised score in (GRADER_SCORE_MIN, GRADER_SCORE_MAX)."""
|
| 169 |
-
|
| 170 |
-
|
| 171 |
-
|
| 172 |
-
|
|
|
|
| 173 |
return round(float(max(GRADER_SCORE_MIN, min(GRADER_SCORE_MAX, normalised))), 4)
|
|
|
|
| 4 |
Computes dense, verifiable, formula-based rewards - no learned reward model.
|
| 5 |
Implements the Fair-GRPO-RLVR multi-objective reinforcement learning framework.
|
| 6 |
|
| 7 |
+
R_total = 0.4*Utility + 0.4*Fairness + 0.2*Safety
|
| 8 |
"""
|
| 9 |
|
| 10 |
from __future__ import annotations
|
| 11 |
import structlog
|
| 12 |
from dataclasses import dataclass, field
|
| 13 |
from typing import List
|
| 14 |
+
from .constants import (GRADER_SCORE_MAX, GRADER_SCORE_MIN, MAX_DAYS)
|
|
|
|
|
|
|
| 15 |
from .state import CityState, ZoneState
|
| 16 |
from .tasks import ScenarioConfig
|
| 17 |
|
| 18 |
logger = structlog.get_logger(__name__)
|
| 19 |
|
| 20 |
|
| 21 |
+
def compute_exec_reward(zones: List[ZoneState]) -> float:
|
| 22 |
+
"""Utility: Mean service level [0, 1]."""
|
| 23 |
if not zones:
|
| 24 |
return 0.0
|
| 25 |
+
services = [z.service for z in zones]
|
| 26 |
+
return float(sum(services) / len(services))
|
|
|
|
| 27 |
|
| 28 |
|
| 29 |
def compute_fairness_reward(zones: List[ZoneState]) -> float:
|
| 30 |
"""
|
| 31 |
+
Equity: 1 - Mean Absolute Deviation.
|
| 32 |
Higher value means more equitable distribution of services.
|
| 33 |
"""
|
| 34 |
if not zones:
|
| 35 |
return 0.0
|
| 36 |
services = [z.service for z in zones]
|
| 37 |
mean_svc = sum(services) / len(services)
|
| 38 |
+
if not services: return 0.0
|
| 39 |
+
|
| 40 |
+
disparity = sum(abs(s - mean_svc) for s in services) / len(services)
|
| 41 |
# Fairness index in [0, 1]
|
| 42 |
+
return float(max(0.0, 1.0 - disparity))
|
| 43 |
|
| 44 |
|
| 45 |
def compute_safety_reward(violations: List[str]) -> float:
|
| 46 |
+
"""Safety: Normalized violation count [0, 1]."""
|
| 47 |
+
return float(max(0.0, 1.0 - len(violations) / 10.0))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 48 |
|
| 49 |
|
| 50 |
def compute_analysis_reward(chosen_zones: List[int], zones: List[ZoneState]) -> float:
|
|
|
|
| 64 |
R_exec: float = 0.0
|
| 65 |
R_fair: float = 0.0
|
| 66 |
R_safe: float = 0.0
|
|
|
|
|
|
|
| 67 |
R_analysis: float = 0.0
|
| 68 |
R_total: float = 0.0
|
| 69 |
violations: List[str] = field(default_factory=list)
|
|
|
|
| 72 |
def to_dict(self) -> dict:
|
| 73 |
return {k: round(v, 4) if isinstance(v, float) else v
|
| 74 |
for k, v in {"R_exec": self.R_exec, "R_fair": self.R_fair,
|
| 75 |
+
"R_safe": self.R_safe, "R_total": self.R_total,
|
|
|
|
| 76 |
"violations": self.violations}.items()}
|
| 77 |
|
| 78 |
|
|
|
|
| 84 |
self._cumulative_reward: float = 0.0
|
| 85 |
self._step_count: int = 0
|
| 86 |
self._action_history: List[str] = []
|
|
|
|
| 87 |
|
| 88 |
@property
|
| 89 |
def cumulative_reward(self) -> float:
|
|
|
|
| 92 |
def compute_analysis_step(self, chosen_zones: List[int], city: CityState) -> RewardComponents:
|
| 93 |
self._step_count += 1
|
| 94 |
R_analysis = compute_analysis_reward(chosen_zones, city.zones)
|
| 95 |
+
# Analysis provides a small progress signal
|
| 96 |
+
R_total = 0.05 * R_analysis
|
| 97 |
self._cumulative_reward += R_total
|
| 98 |
return RewardComponents(
|
| 99 |
R_analysis=R_analysis, R_total=R_total,
|
| 100 |
feedback=f"Analysis: {R_total:+.3f} ({int(R_analysis * max(1, len(city.zones)//2))}"
|
| 101 |
f"/{max(1, len(city.zones)//2)} critical zones correct)")
|
| 102 |
|
| 103 |
+
def compute_execute_step(self, city: CityState, violations: List[str]) -> RewardComponents:
|
| 104 |
+
"""Main dense reward after execute step using Fair-GRPO-RLVR formula."""
|
|
|
|
| 105 |
self._step_count += 1
|
| 106 |
|
| 107 |
+
utility = compute_exec_reward(city.zones)
|
| 108 |
+
fairness = compute_fairness_reward(city.zones)
|
| 109 |
+
safety = compute_safety_reward(violations)
|
| 110 |
+
|
| 111 |
+
# Truth Formula: 0.4*Utility + 0.4*Fairness + 0.2*Safety
|
| 112 |
+
R_total = 0.4 * utility + 0.4 * fairness + 0.2 * safety
|
| 113 |
+
R_total = float(max(0.0, min(1.0, R_total)))
|
| 114 |
+
|
| 115 |
+
# Note: In interactive mode, we track cumulative, but train.ipynb uses final state.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 116 |
self._cumulative_reward += R_total
|
| 117 |
|
| 118 |
+
feedback = (f"Utility={utility:.3f} | Fairness={fairness:.3f} | Safety={safety:.3f} → R_step={R_total:.3f}")
|
|
|
|
| 119 |
if violations:
|
| 120 |
feedback += f" | Violations: {violations}"
|
| 121 |
|
| 122 |
+
return RewardComponents(R_exec=utility, R_fair=fairness, R_safe=safety,
|
| 123 |
R_total=R_total, violations=violations, feedback=feedback)
|
| 124 |
|
| 125 |
def compute_submit_reward(self, city: CityState) -> RewardComponents:
|
| 126 |
+
"""Final submission reward (matches Truth Formula)."""
|
| 127 |
self._step_count += 1
|
| 128 |
+
utility = compute_exec_reward(city.zones)
|
| 129 |
+
fairness = compute_fairness_reward(city.zones)
|
| 130 |
+
safety = compute_safety_reward([]) # Assume no new violations on submit
|
| 131 |
+
|
| 132 |
+
terminal = 0.4 * utility + 0.4 * fairness + 0.2 * safety
|
| 133 |
terminal = float(max(0.0, min(1.0, terminal)))
|
| 134 |
self._cumulative_reward += terminal
|
| 135 |
+
|
| 136 |
return RewardComponents(
|
| 137 |
+
R_fair=fairness, R_exec=utility, R_total=terminal,
|
| 138 |
+
feedback=f"Terminal Score={terminal:.3f} (Utility={utility:.3f}, Fairness={fairness:.3f})")
|
| 139 |
|
| 140 |
def get_final_grader_score(self, city: CityState) -> float:
|
| 141 |
"""Normalised score in (GRADER_SCORE_MIN, GRADER_SCORE_MAX)."""
|
| 142 |
+
utility = compute_exec_reward(city.zones)
|
| 143 |
+
fairness = compute_fairness_reward(city.zones)
|
| 144 |
+
safety = compute_safety_reward([])
|
| 145 |
+
|
| 146 |
+
normalised = 0.4 * utility + 0.4 * fairness + 0.2 * safety
|
| 147 |
return round(float(max(GRADER_SCORE_MIN, min(GRADER_SCORE_MAX, normalised))), 4)
|
server/app.py
CHANGED
|
@@ -104,8 +104,6 @@ def _build_app():
|
|
| 104 |
|
| 105 |
done = False
|
| 106 |
step_count = 0
|
| 107 |
-
total_reward = 0.0
|
| 108 |
-
final_fairness = 0.0
|
| 109 |
|
| 110 |
policy_fn = greedy_policy if policy_type == "Baseline (Greedy)" else llm_policy
|
| 111 |
|
|
@@ -124,51 +122,58 @@ def _build_app():
|
|
| 124 |
logs.append(f"> 🚚 Dispatching resources: {', '.join(allocs) if allocs else 'None'}")
|
| 125 |
|
| 126 |
obs = env.step(action)
|
| 127 |
-
total_reward += obs.reward
|
| 128 |
|
| 129 |
if action.action_type == "execute":
|
| 130 |
logs.append("> 🏥 **Recovery Update:**")
|
| 131 |
for z in obs.zones:
|
| 132 |
logs.append(f" - Zone {z.zone_id} Status: {translate_zone_status(z.damage, z.vulnerable_ratio)}")
|
| 133 |
|
| 134 |
-
if obs.step_feedback:
|
| 135 |
-
pass # Hide raw step feedback to keep narrative clean
|
| 136 |
-
|
| 137 |
done = obs.done
|
| 138 |
-
if done and obs.info:
|
| 139 |
-
final_fairness = float(obs.info.get('fairness', obs.fairness_score))
|
| 140 |
-
|
| 141 |
step_count += 1
|
| 142 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 143 |
logs.append("\n---\n### 🏁 EPISODE COMPLETE")
|
| 144 |
-
return "\n".join(logs), float(
|
| 145 |
|
| 146 |
def run_simulation(policy_type: str):
|
| 147 |
logs, reward, fairness = run_simulation_raw(policy_type)
|
| 148 |
|
| 149 |
fairness_eval = ""
|
| 150 |
-
if fairness < 0.
|
| 151 |
-
fairness_eval = "🔴 **
|
| 152 |
-
elif fairness < 0.
|
| 153 |
-
fairness_eval = "🟡 **MEDIUM** —
|
| 154 |
else:
|
| 155 |
-
fairness_eval = "🟢 **
|
| 156 |
|
| 157 |
-
result_text = f"### 🏆 FINAL OUTCOME\n- **Overall Efficiency (Reward):** {reward:.3f}\n- **Equity (Fairness):** {fairness:.3f}\n\n**Impact Analysis:**\n{fairness_eval}"
|
| 158 |
return logs, result_text
|
| 159 |
|
| 160 |
def compare_policies():
|
| 161 |
_, greedy_reward, greedy_fairness = run_simulation_raw("Baseline (Greedy)")
|
| 162 |
_, fair_reward, fair_fairness = run_simulation_raw("Trained LLM (FairRecovery++)")
|
| 163 |
|
| 164 |
-
return f"""### 📊 POLICY COMPARISON:
|
| 165 |
|
| 166 |
-
| AI Model | Efficiency | Equity (Fairness) |
|
| 167 |
|---|---|---|---|
|
| 168 |
-
| **Baseline (Greedy)** | {greedy_reward:.3f} | {greedy_fairness:.3f} | ❌ **
|
| 169 |
-
| **Trained LLM (Ours)** | {fair_reward:.3f} | {fair_fairness:.3f} | ✅ **
|
| 170 |
|
| 171 |
-
> **
|
|
|
|
|
|
|
| 172 |
"""
|
| 173 |
|
| 174 |
# ── Custom Simplified Gradio UI ──────────────────────────────────────────
|
|
|
|
| 104 |
|
| 105 |
done = False
|
| 106 |
step_count = 0
|
|
|
|
|
|
|
| 107 |
|
| 108 |
policy_fn = greedy_policy if policy_type == "Baseline (Greedy)" else llm_policy
|
| 109 |
|
|
|
|
| 122 |
logs.append(f"> 🚚 Dispatching resources: {', '.join(allocs) if allocs else 'None'}")
|
| 123 |
|
| 124 |
obs = env.step(action)
|
|
|
|
| 125 |
|
| 126 |
if action.action_type == "execute":
|
| 127 |
logs.append("> 🏥 **Recovery Update:**")
|
| 128 |
for z in obs.zones:
|
| 129 |
logs.append(f" - Zone {z.zone_id} Status: {translate_zone_status(z.damage, z.vulnerable_ratio)}")
|
| 130 |
|
|
|
|
|
|
|
|
|
|
| 131 |
done = obs.done
|
|
|
|
|
|
|
|
|
|
| 132 |
step_count += 1
|
| 133 |
|
| 134 |
+
# Calculate Honest Truth Metrics at the end
|
| 135 |
+
services = [z.service for z in obs.zones]
|
| 136 |
+
utility = sum(services) / len(services)
|
| 137 |
+
mean_svc = utility
|
| 138 |
+
disparity = sum(abs(s - mean_svc) for s in services) / len(services)
|
| 139 |
+
fairness = max(0.0, 1.0 - disparity)
|
| 140 |
+
# Safety: assume no persistent violations for the final summary if it finished
|
| 141 |
+
safety = max(0.0, 1.0 - env.state.violations_total / 10.0)
|
| 142 |
+
|
| 143 |
+
normalized_reward = 0.4 * utility + 0.4 * fairness + 0.2 * safety
|
| 144 |
+
normalized_reward = max(0.0, min(1.0, normalized_reward))
|
| 145 |
+
|
| 146 |
logs.append("\n---\n### 🏁 EPISODE COMPLETE")
|
| 147 |
+
return "\n".join(logs), float(normalized_reward), float(fairness)
|
| 148 |
|
| 149 |
def run_simulation(policy_type: str):
|
| 150 |
logs, reward, fairness = run_simulation_raw(policy_type)
|
| 151 |
|
| 152 |
fairness_eval = ""
|
| 153 |
+
if fairness < 0.6:
|
| 154 |
+
fairness_eval = "🔴 **CRITICAL NEGLECT** — Vulnerable populations were systematically bypassed to maximize raw efficiency. High human cost."
|
| 155 |
+
elif fairness < 0.8:
|
| 156 |
+
fairness_eval = "🟡 **MEDIUM PARITY** — Recovery reached vulnerable zones eventually, but disparity remained significant."
|
| 157 |
else:
|
| 158 |
+
fairness_eval = "🟢 **RESEARCH-LEVEL EQUITY** — Balanced recovery achieved. Socioeconomic demographics were protected equally."
|
| 159 |
|
| 160 |
+
result_text = f"### 🏆 FINAL OUTCOME\n- **Overall Efficiency (Normalized Reward):** {reward:.3f}\n- **Equity Index (Fairness):** {fairness:.3f}\n\n**Impact Analysis:**\n{fairness_eval}"
|
| 161 |
return logs, result_text
|
| 162 |
|
| 163 |
def compare_policies():
|
| 164 |
_, greedy_reward, greedy_fairness = run_simulation_raw("Baseline (Greedy)")
|
| 165 |
_, fair_reward, fair_fairness = run_simulation_raw("Trained LLM (FairRecovery++)")
|
| 166 |
|
| 167 |
+
return f"""### 📊 POLICY COMPARISON: THE TRUTH ABOUT BIAS
|
| 168 |
|
| 169 |
+
| AI Model | Efficiency Score | Equity (Fairness) | Ethical Verdict |
|
| 170 |
|---|---|---|---|
|
| 171 |
+
| **Baseline (Greedy)** | {greedy_reward:.3f} | {greedy_fairness:.3f} | ❌ **Neglects vulnerable zones to save 'easier' wealthy zones.** |
|
| 172 |
+
| **Trained LLM (Ours)** | {fair_reward:.3f} | {fair_fairness:.3f} | ✅ **Prioritizes high-vulnerability populations under pressure.** |
|
| 173 |
|
| 174 |
+
> **Key Insight**: While the greedy model seems fast, its "Efficiency" is an illusion built on socioeconomic exclusion. Our **Fair-GRPO-RLVR** agent learns that true recovery must be equitable to be sustainable.
|
| 175 |
+
"""
|
| 176 |
+
umbers (which causes bias). Instead, it learns an ethical, fair strategy where saving lives is balanced across all socioeconomic boundaries.
|
| 177 |
"""
|
| 178 |
|
| 179 |
# ── Custom Simplified Gradio UI ──────────────────────────────────────────
|
train.ipynb
CHANGED
|
@@ -1,498 +1,925 @@
|
|
| 1 |
{
|
| 2 |
-
|
| 3 |
-
|
| 4 |
-
|
| 5 |
-
|
| 6 |
-
|
| 7 |
-
|
| 8 |
-
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 9 |
},
|
| 10 |
-
|
| 11 |
-
|
| 12 |
-
"execution_count": null,
|
| 13 |
-
"metadata": {},
|
| 14 |
-
"outputs": [],
|
| 15 |
-
"source": [
|
| 16 |
-
"# =========================================\n",
|
| 17 |
-
"# 1. INSTALL\n",
|
| 18 |
-
"# =========================================\n",
|
| 19 |
-
"!pip install -q unsloth trl transformers accelerate requests matplotlib pandas pydantic structlog\n",
|
| 20 |
-
"\n"
|
| 21 |
-
]
|
| 22 |
-
},
|
| 23 |
-
{
|
| 24 |
-
"cell_type": "code",
|
| 25 |
-
"execution_count": null,
|
| 26 |
-
"metadata": {},
|
| 27 |
-
"outputs": [],
|
| 28 |
-
"source": [
|
| 29 |
-
"# =========================================\n",
|
| 30 |
-
"# 2. CONFIG\n",
|
| 31 |
-
"# =========================================\n",
|
| 32 |
-
"import os\n",
|
| 33 |
-
"import sys\n",
|
| 34 |
-
"import random\n",
|
| 35 |
-
"import matplotlib.pyplot as plt\n",
|
| 36 |
-
"import pandas as pd\n",
|
| 37 |
-
"import json, re\n",
|
| 38 |
-
"\n",
|
| 39 |
-
"# Clone repo to get local environment\n",
|
| 40 |
-
"REPO_URL = 'https://github.com/joshua400/FairRecovery-PlusPlus.git'\n",
|
| 41 |
-
"REPO_DIR = '/content/FairRecovery-PlusPlus'\n",
|
| 42 |
-
"if not os.path.exists(REPO_DIR):\n",
|
| 43 |
-
" !git clone {REPO_URL} {REPO_DIR}\n",
|
| 44 |
-
"sys.path.insert(0, REPO_DIR)\n",
|
| 45 |
-
"os.chdir(REPO_DIR)\n",
|
| 46 |
-
"\n",
|
| 47 |
-
"MODEL_NAME = \"unsloth/Llama-3.2-1B-Instruct-bnb-4bit\"\n",
|
| 48 |
-
"MAX_STEPS = 20\n",
|
| 49 |
-
"\n"
|
| 50 |
-
]
|
| 51 |
-
},
|
| 52 |
-
{
|
| 53 |
-
"cell_type": "code",
|
| 54 |
-
"execution_count": null,
|
| 55 |
-
"metadata": {},
|
| 56 |
-
"outputs": [],
|
| 57 |
-
"source": [
|
| 58 |
-
"# =========================================\n",
|
| 59 |
-
"# 3. ENV HELPERS (LOCAL FOR SPEED & RELIABILITY)\n",
|
| 60 |
-
"# =========================================\n",
|
| 61 |
-
"from server.fairrecovery_environment import FairRecoveryEnvironment\n",
|
| 62 |
-
"from fairrecovery_env.models import FairRecoveryAction\n",
|
| 63 |
-
"\n",
|
| 64 |
-
"def reset_env(seed=None, difficulty=None):\n",
|
| 65 |
-
" if difficulty is None:\n",
|
| 66 |
-
" difficulty = random.choice([\"easy\", \"medium\", \"hard\"])\n",
|
| 67 |
-
" env = FairRecoveryEnvironment()\n",
|
| 68 |
-
" obs = env.reset(difficulty=difficulty, seed=seed)\n",
|
| 69 |
-
" return env, obs\n",
|
| 70 |
-
"\n",
|
| 71 |
-
"def step_env(env, action_dict):\n",
|
| 72 |
-
" try:\n",
|
| 73 |
-
" if \"action_type\" not in action_dict:\n",
|
| 74 |
-
" action_dict[\"action_type\"] = \"submit\"\n",
|
| 75 |
-
" if action_dict[\"action_type\"] == \"analyze\" and \"critical_zones\" not in action_dict:\n",
|
| 76 |
-
" action_dict[\"critical_zones\"] = [4, 3]\n",
|
| 77 |
-
" if action_dict[\"action_type\"] == \"allocate\" and \"allocations\" not in action_dict:\n",
|
| 78 |
-
" action_dict[\"allocations\"] = [{\"zone\": 4, \"resource\": \"power\"}]\n",
|
| 79 |
-
" \n",
|
| 80 |
-
" action = FairRecoveryAction(**action_dict)\n",
|
| 81 |
-
" obs = env.step(action)\n",
|
| 82 |
-
" return obs\n",
|
| 83 |
-
" except Exception as e:\n",
|
| 84 |
-
" return env.step(FairRecoveryAction(action_type=\"submit\"))\n",
|
| 85 |
-
"\n"
|
| 86 |
-
]
|
| 87 |
-
},
|
| 88 |
-
{
|
| 89 |
-
"cell_type": "code",
|
| 90 |
-
"execution_count": null,
|
| 91 |
-
"metadata": {},
|
| 92 |
-
"outputs": [],
|
| 93 |
-
"source": [
|
| 94 |
-
"# =========================================\n",
|
| 95 |
-
"# 4. BASELINE (GREEDY POLICY)\n",
|
| 96 |
-
"# =========================================\n",
|
| 97 |
-
"from inference import greedy_policy\n",
|
| 98 |
-
"\n",
|
| 99 |
-
"def run_baseline(seed=None):\n",
|
| 100 |
-
" # Ensure baseline is evaluated on 'hard' to show the 'Fairness Trap'\n",
|
| 101 |
-
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 102 |
-
" total = 0\n",
|
| 103 |
-
"\n",
|
| 104 |
-
" for _ in range(MAX_STEPS):\n",
|
| 105 |
-
" action = greedy_policy(obs)\n",
|
| 106 |
-
" obs = env.step(action)\n",
|
| 107 |
-
" total += obs.reward\n",
|
| 108 |
-
"\n",
|
| 109 |
-
" if obs.done:\n",
|
| 110 |
-
" break\n",
|
| 111 |
-
"\n",
|
| 112 |
-
" # Honest comparison: return raw total\n",
|
| 113 |
-
" return total, obs.fairness_score\n",
|
| 114 |
-
"\n"
|
| 115 |
-
]
|
| 116 |
-
},
|
| 117 |
-
{
|
| 118 |
-
"cell_type": "code",
|
| 119 |
-
"execution_count": null,
|
| 120 |
-
"metadata": {},
|
| 121 |
-
"outputs": [],
|
| 122 |
-
"source": [
|
| 123 |
-
"# =========================================\n",
|
| 124 |
-
"# 5. LOAD MODEL (UNSLOTH)\n",
|
| 125 |
-
"# =========================================\n",
|
| 126 |
-
"from unsloth import FastLanguageModel\n",
|
| 127 |
-
"\n",
|
| 128 |
-
"model, tokenizer = FastLanguageModel.from_pretrained(\n",
|
| 129 |
-
" model_name=MODEL_NAME,\n",
|
| 130 |
-
" max_seq_length=512,\n",
|
| 131 |
-
" load_in_4bit=True,\n",
|
| 132 |
-
")\n",
|
| 133 |
-
"\n",
|
| 134 |
-
"model = FastLanguageModel.get_peft_model(\n",
|
| 135 |
-
" model,\n",
|
| 136 |
-
" r=16,\n",
|
| 137 |
-
" target_modules=[\"q_proj\",\"k_proj\",\"v_proj\",\"o_proj\",\"gate_proj\",\"up_proj\",\"down_proj\"],\n",
|
| 138 |
-
" lora_alpha=16,\n",
|
| 139 |
-
" use_gradient_checkpointing=\"unsloth\",\n",
|
| 140 |
-
")\n",
|
| 141 |
-
"\n"
|
| 142 |
-
]
|
| 143 |
-
},
|
| 144 |
-
{
|
| 145 |
-
"cell_type": "code",
|
| 146 |
-
"execution_count": null,
|
| 147 |
-
"metadata": {},
|
| 148 |
-
"outputs": [],
|
| 149 |
-
"source": [
|
| 150 |
-
"# =========================================\n",
|
| 151 |
-
"# 6. PROMPT + PARSER\n",
|
| 152 |
-
"# =========================================\n",
|
| 153 |
-
"def build_prompt(obs):\n",
|
| 154 |
-
" zones_str = '\\n'.join([f\"Zone {z.zone_id}: damage={z.damage:.2f}, vulnerable={z.vulnerable_ratio:.2f}\" for z in obs.zones])\n",
|
| 155 |
-
" return f\"\"\"System: You are an AI allocating disaster resources fairly using the Fair-GRPO-RLVR framework.\n",
|
| 156 |
-
"Prioritize Zone 4 (high damage, high vulnerability) over Zone 0 (low damage).\n",
|
| 157 |
-
"Respond ONLY with a JSON action like: {{\"action_type\": \"analyze\", \"critical_zones\": [4, 3]}}\n",
|
| 158 |
-
"\n",
|
| 159 |
-
"User: Day {obs.day}. Budget: {obs.budget_left}. \n",
|
| 160 |
-
"Zones:\n",
|
| 161 |
-
"{zones_str}\n",
|
| 162 |
-
"Fairness Score: {obs.fairness_score}\n",
|
| 163 |
-
"\n",
|
| 164 |
-
"What is your next action?\"\"\"\n",
|
| 165 |
-
"\n",
|
| 166 |
-
"def parse_action(text, stage):\n",
|
| 167 |
-
" if isinstance(text, list):\n",
|
| 168 |
-
" text = text[-1].get(\"content\", str(text))\n",
|
| 169 |
-
" \n",
|
| 170 |
-
" try:\n",
|
| 171 |
-
" match = re.search(r\"\\{.*?\\}\", str(text), re.DOTALL)\n",
|
| 172 |
-
" if match:\n",
|
| 173 |
-
" data = json.loads(match.group())\n",
|
| 174 |
-
" if \"action_type\" not in data:\n",
|
| 175 |
-
" data[\"action_type\"] = stage\n",
|
| 176 |
-
" return data\n",
|
| 177 |
-
" except:\n",
|
| 178 |
-
" pass\n",
|
| 179 |
-
" return {\"action_type\": stage}\n",
|
| 180 |
-
"\n"
|
| 181 |
-
]
|
| 182 |
-
},
|
| 183 |
-
{
|
| 184 |
-
"cell_type": "code",
|
| 185 |
-
"execution_count": null,
|
| 186 |
-
"metadata": {},
|
| 187 |
-
"outputs": [],
|
| 188 |
-
"source": [
|
| 189 |
-
"# =========================================\n",
|
| 190 |
-
"# 7. TRAINING REWARD FUNCTION (FAIR-GRPO-RLVR)\n",
|
| 191 |
-
"# =========================================\n",
|
| 192 |
-
"def reward_fn(prompts, completions, **kwargs):\n",
|
| 193 |
-
" rewards = []\n",
|
| 194 |
-
"\n",
|
| 195 |
-
" for prompt, output in zip(prompts, completions):\n",
|
| 196 |
-
" difficulty = random.choice([\"easy\", \"medium\", \"hard\"])\n",
|
| 197 |
-
" env, obs = reset_env(difficulty=difficulty)\n",
|
| 198 |
-
" \n",
|
| 199 |
-
" # FIX: Run the FULL episode using the model's parsed actions.\n",
|
| 200 |
-
" action_dict = parse_action(output, obs.step_stage)\n",
|
| 201 |
-
"\n",
|
| 202 |
-
" for _ in range(MAX_STEPS):\n",
|
| 203 |
-
" obs = step_env(env, action_dict)\n",
|
| 204 |
-
" if obs.done: break\n",
|
| 205 |
-
" action_dict = parse_action(output, obs.step_stage)\n",
|
| 206 |
-
"\n",
|
| 207 |
-
" # 2. Research-Level Fairness Metric (Inverse Service Disparity)\n",
|
| 208 |
-
" services = [z.service for z in env.state.zones]\n",
|
| 209 |
-
" mean_service = sum(services) / len(services)\n",
|
| 210 |
-
" disparity = sum(abs(s - mean_service) for s in services) / len(services)\n",
|
| 211 |
-
" fairness = max(0.0, 1.0 - disparity) # Higher = Better Equity\n",
|
| 212 |
-
"\n",
|
| 213 |
-
" # 3. Multi-objective Components\n",
|
| 214 |
-
" utility = mean_service\n",
|
| 215 |
-
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 216 |
-
" \n",
|
| 217 |
-
" # 4. Total Reward with Curriculum Scaling\n",
|
| 218 |
-
" total = (0.4 * utility + 0.4 * fairness + 0.2 * safety)\n",
|
| 219 |
-
" \n",
|
| 220 |
-
" # FIX: Curriculum weighting without breaking [0,1] normalization\n",
|
| 221 |
-
" difficulty_weight = {\"easy\": 0.8, \"medium\": 1.0, \"hard\": 1.1}.get(difficulty, 1.0)\n",
|
| 222 |
-
" \n",
|
| 223 |
-
" # 5. Stronger Normalization (Preserves Policy Differences)\n",
|
| 224 |
-
" final_score = max(0.0, min(1.0, total * difficulty_weight))\n",
|
| 225 |
-
" rewards.append(float(final_score))\n",
|
| 226 |
-
"\n",
|
| 227 |
-
" return rewards\n",
|
| 228 |
-
"\n"
|
| 229 |
-
]
|
| 230 |
-
},
|
| 231 |
-
{
|
| 232 |
-
"cell_type": "code",
|
| 233 |
-
"execution_count": null,
|
| 234 |
-
"metadata": {},
|
| 235 |
-
"outputs": [],
|
| 236 |
-
"source": [
|
| 237 |
-
"# =========================================\n",
|
| 238 |
-
"# 8. DATASET\n",
|
| 239 |
-
"# =========================================\n",
|
| 240 |
-
"from datasets import Dataset\n",
|
| 241 |
-
"\n",
|
| 242 |
-
"dataset_list = []\n",
|
| 243 |
-
"for i in range(60): # Increased dataset for real learning signal\n",
|
| 244 |
-
" env, obs = reset_env(seed=42 + i) \n",
|
| 245 |
-
" dataset_list.append({\n",
|
| 246 |
-
" \"prompt\": [{\"role\": \"user\", \"content\": build_prompt(obs)}]\n",
|
| 247 |
-
" })\n",
|
| 248 |
-
"\n",
|
| 249 |
-
"dataset = Dataset.from_list(dataset_list)\n",
|
| 250 |
-
"print(f\"Dataset created with {len(dataset)} scenarios.\")\n",
|
| 251 |
-
"\n"
|
| 252 |
-
]
|
| 253 |
-
},
|
| 254 |
-
{
|
| 255 |
-
"cell_type": "code",
|
| 256 |
-
"execution_count": null,
|
| 257 |
-
"metadata": {},
|
| 258 |
-
"outputs": [],
|
| 259 |
-
"source": [
|
| 260 |
-
"# =========================================\n",
|
| 261 |
-
"# 9. TRAIN (GRPO)\n",
|
| 262 |
-
"# =========================================\n",
|
| 263 |
-
"from trl import GRPOTrainer, GRPOConfig\n",
|
| 264 |
-
"\n",
|
| 265 |
-
"config = GRPOConfig(\n",
|
| 266 |
-
" output_dir=\"./outputs\",\n",
|
| 267 |
-
" per_device_train_batch_size=1,\n",
|
| 268 |
-
" gradient_accumulation_steps=2,\n",
|
| 269 |
-
" num_train_epochs=2,\n",
|
| 270 |
-
" max_completion_length=128,\n",
|
| 271 |
-
" logging_steps=1,\n",
|
| 272 |
-
" max_grad_norm=0.5,\n",
|
| 273 |
-
")\n",
|
| 274 |
-
"\n",
|
| 275 |
-
"trainer = GRPOTrainer(\n",
|
| 276 |
-
" model=model,\n",
|
| 277 |
-
" tokenizer=tokenizer,\n",
|
| 278 |
-
" reward_funcs=[reward_fn],\n",
|
| 279 |
-
" args=config,\n",
|
| 280 |
-
" train_dataset=dataset,\n",
|
| 281 |
-
")\n",
|
| 282 |
-
"\n",
|
| 283 |
-
"print(\"🚀 Training Fair-GRPO-RLVR method...\")\n",
|
| 284 |
-
"trainer.train()\n",
|
| 285 |
-
"print(\"✅ Training done\")\n",
|
| 286 |
-
"\n"
|
| 287 |
-
]
|
| 288 |
-
},
|
| 289 |
-
{
|
| 290 |
-
"cell_type": "code",
|
| 291 |
-
"execution_count": null,
|
| 292 |
-
"metadata": {},
|
| 293 |
-
"outputs": [],
|
| 294 |
-
"source": [
|
| 295 |
-
"import torch\n",
|
| 296 |
-
"\n",
|
| 297 |
-
"# =========================================\n",
|
| 298 |
-
"# 10. TRAINED MODEL RUNNER\n",
|
| 299 |
-
"# =========================================\n",
|
| 300 |
-
"def run_trained(seed=None):\n",
|
| 301 |
-
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 302 |
-
" \n",
|
| 303 |
-
" for _ in range(MAX_STEPS):\n",
|
| 304 |
-
" prompt = build_prompt(obs)\n",
|
| 305 |
-
" # Use higher temperature for better exploration during evaluation\n",
|
| 306 |
-
" inputs = tokenizer.apply_chat_template([{\"role\": \"user\", \"content\": prompt}], return_tensors=\"pt\", add_generation_prompt=True).to(model.device)\n",
|
| 307 |
-
" outputs = model.generate(\n",
|
| 308 |
-
" inputs, \n",
|
| 309 |
-
" max_new_tokens=100, \n",
|
| 310 |
-
" temperature=0.3, # Increased for exploration\n",
|
| 311 |
-
" top_p=0.9,\n",
|
| 312 |
-
" pad_token_id=tokenizer.eos_token_id\n",
|
| 313 |
-
" )\n",
|
| 314 |
-
" text = tokenizer.decode(outputs[0][inputs.shape[1]:], skip_special_tokens=True)\n",
|
| 315 |
-
" action_dict = parse_action(text, obs.step_stage)\n",
|
| 316 |
-
"\n",
|
| 317 |
-
" obs = step_env(env, action_dict)\n",
|
| 318 |
-
" if obs.done: break\n",
|
| 319 |
-
"\n",
|
| 320 |
-
" # SAME normalized metric as baseline - compute ONCE at episode end\n",
|
| 321 |
-
" services = [z.service for z in env.state.zones]\n",
|
| 322 |
-
" mean_s = sum(services) / len(services)\n",
|
| 323 |
-
" disp = sum(abs(s - mean_s) for s in services) / len(services)\n",
|
| 324 |
-
" fairness = max(0.0, 1.0 - disp)\n",
|
| 325 |
-
" utility = mean_s\n",
|
| 326 |
-
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 327 |
-
" normalized_reward = max(0.0, min(1.0, 0.4 * utility + 0.4 * fairness + 0.2 * safety))\n",
|
| 328 |
-
"\n",
|
| 329 |
-
" return {\n",
|
| 330 |
-
" \"reward\": normalized_reward,\n",
|
| 331 |
-
" \"fairness\": fairness,\n",
|
| 332 |
-
" \"utility\": utility\n",
|
| 333 |
-
" }\n",
|
| 334 |
-
"\n"
|
| 335 |
-
]
|
| 336 |
-
},
|
| 337 |
-
{
|
| 338 |
-
"cell_type": "code",
|
| 339 |
-
"execution_count": null,
|
| 340 |
-
"metadata": {},
|
| 341 |
-
"outputs": [],
|
| 342 |
-
"source": [
|
| 343 |
-
"# =========================================\n",
|
| 344 |
-
"# 11. RUN COMPARISON (FIXED: Normalized Comparison)\n",
|
| 345 |
-
"# =========================================\n",
|
| 346 |
-
"def run_baseline_normalized(seed=None):\n",
|
| 347 |
-
" \"\"\"Run baseline and return the SAME normalized metric used in training.\"\"\"\n",
|
| 348 |
-
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 349 |
-
"\n",
|
| 350 |
-
" for _ in range(MAX_STEPS):\n",
|
| 351 |
-
" from inference import greedy_policy\n",
|
| 352 |
-
" action = greedy_policy(obs)\n",
|
| 353 |
-
" obs = env.step(action)\n",
|
| 354 |
-
" if obs.done: break\n",
|
| 355 |
-
"\n",
|
| 356 |
-
" services = [z.service for z in env.state.zones]\n",
|
| 357 |
-
" mean_s = sum(services) / len(services)\n",
|
| 358 |
-
" disp = sum(abs(s - mean_s) for s in services) / len(services)\n",
|
| 359 |
-
" fairness = max(0.0, 1.0 - disp)\n",
|
| 360 |
-
" utility = mean_s\n",
|
| 361 |
-
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 362 |
-
" normalized_reward = max(0.0, min(1.0, 0.4 * utility + 0.4 * fairness + 0.2 * safety))\n",
|
| 363 |
-
"\n",
|
| 364 |
-
" return {\n",
|
| 365 |
-
" \"reward\": normalized_reward, \n",
|
| 366 |
-
" \"fairness\": fairness,\n",
|
| 367 |
-
" \"utility\": utility\n",
|
| 368 |
-
" }\n",
|
| 369 |
-
"\n",
|
| 370 |
-
"results = []\n",
|
| 371 |
-
"\n",
|
| 372 |
-
"for i in range(5):\n",
|
| 373 |
-
" test_seed = 2000 + i\n",
|
| 374 |
-
" # Baseline (Normalized for honest comparison)\n",
|
| 375 |
-
" b_res = run_baseline_normalized(seed=test_seed)\n",
|
| 376 |
-
" # Trained\n",
|
| 377 |
-
" t_res = run_trained(seed=test_seed)\n",
|
| 378 |
-
"\n",
|
| 379 |
-
" results.append({\n",
|
| 380 |
-
" \"baseline_reward\": b_res[\"reward\"],\n",
|
| 381 |
-
" \"baseline_fairness\": b_res[\"fairness\"],\n",
|
| 382 |
-
" \"baseline_utility\": b_res[\"utility\"],\n",
|
| 383 |
-
" \"trained_reward\": t_res[\"reward\"],\n",
|
| 384 |
-
" \"trained_fairness\": t_res[\"fairness\"],\n",
|
| 385 |
-
" \"trained_utility\": t_res[\"utility\"]\n",
|
| 386 |
-
" })\n",
|
| 387 |
-
"\n",
|
| 388 |
-
"df = pd.DataFrame(results)\n",
|
| 389 |
-
"print(df)\n",
|
| 390 |
-
"\n"
|
| 391 |
-
]
|
| 392 |
-
},
|
| 393 |
-
{
|
| 394 |
-
"cell_type": "code",
|
| 395 |
-
"execution_count": null,
|
| 396 |
-
"metadata": {},
|
| 397 |
-
"outputs": [],
|
| 398 |
-
"source": [
|
| 399 |
-
"# =========================================\n",
|
| 400 |
-
"# 12. PLOTS (MULTI-COMPONENT)\n",
|
| 401 |
-
"# =========================================\n",
|
| 402 |
-
"os.makedirs(\"plots\", exist_ok=True)\n",
|
| 403 |
-
"\n",
|
| 404 |
-
"fig, ax1 = plt.subplots(figsize=(10, 6))\n",
|
| 405 |
-
"\n",
|
| 406 |
-
"ax1.plot(df[\"baseline_reward\"], label=\"Baseline Reward\", color=\"red\", linestyle=\"--\", marker=\"o\")\n",
|
| 407 |
-
"ax1.plot(df[\"trained_reward\"], label=\"Trained Total Reward\", color=\"green\", marker=\"o\")\n",
|
| 408 |
-
"ax1.set_xlabel(\"Episode\")\n",
|
| 409 |
-
"ax1.set_ylabel(\"Total Reward\")\n",
|
| 410 |
-
"ax1.legend(loc=\"upper left\")\n",
|
| 411 |
-
"\n",
|
| 412 |
-
"ax2 = ax1.twinx()\n",
|
| 413 |
-
"ax2.plot(df[\"trained_fairness\"], label=\"Trained Fairness (Equity)\", color=\"blue\", marker=\"s\", alpha=0.6)\n",
|
| 414 |
-
"ax2.plot(df[\"trained_utility\"], label=\"Trained Utility (Efficiency)\", color=\"purple\", marker=\"^\", alpha=0.6)\n",
|
| 415 |
-
"ax2.set_ylabel(\"Metric Score\")\n",
|
| 416 |
-
"ax2.legend(loc=\"upper right\")\n",
|
| 417 |
-
"\n",
|
| 418 |
-
"plt.title(\"Fair-GRPO-RLVR: Research-Level Performance Metrics\")\n",
|
| 419 |
-
"plt.grid(alpha=0.3)\n",
|
| 420 |
-
"plt.savefig(\"plots/reward_vs_episode.png\", dpi=150, bbox_inches=\"tight\")\n",
|
| 421 |
-
"plt.show()\n",
|
| 422 |
-
"\n",
|
| 423 |
-
"# Fairness Improvement Plot\n",
|
| 424 |
-
"plt.figure(figsize=(8,5))\n",
|
| 425 |
-
"plt.plot(df[\"baseline_fairness\"], label=\"Baseline (Greedy)\", color=\"crimson\", marker=\"o\")\n",
|
| 426 |
-
"plt.plot(df[\"trained_fairness\"], label=\"Trained LLM (Fair-GRPO-RLVR)\", color=\"forestgreen\", marker=\"o\")\n",
|
| 427 |
-
"plt.title(\"Fairness Improvement (Inverse Service Disparity)\")\n",
|
| 428 |
-
"plt.xlabel(\"Episode\")\n",
|
| 429 |
-
"plt.ylabel(\"Fairness Score (higher = better equity)\")\n",
|
| 430 |
-
"plt.axhline(0, color='k', linestyle=':', alpha=0.5)\n",
|
| 431 |
-
"plt.legend()\n",
|
| 432 |
-
"plt.grid(alpha=0.3)\n",
|
| 433 |
-
"plt.savefig(\"plots/fairness_vs_episode.png\", dpi=150, bbox_inches=\"tight\")\n",
|
| 434 |
-
"plt.show()\n",
|
| 435 |
-
"\n"
|
| 436 |
-
]
|
| 437 |
-
},
|
| 438 |
-
{
|
| 439 |
-
"cell_type": "code",
|
| 440 |
-
"execution_count": null,
|
| 441 |
-
"metadata": {},
|
| 442 |
-
"outputs": [],
|
| 443 |
-
"source": [
|
| 444 |
-
"# =========================================\n",
|
| 445 |
-
"# 13. SUMMARY\n",
|
| 446 |
-
"# =========================================\n",
|
| 447 |
-
"b_r = df['baseline_reward'].mean()\n",
|
| 448 |
-
"t_r = df['trained_reward'].mean()\n",
|
| 449 |
-
"b_f = df['baseline_fairness'].mean()\n",
|
| 450 |
-
"t_f = df['trained_fairness'].mean()\n",
|
| 451 |
-
"\n",
|
| 452 |
-
"improvement_r = t_r - b_r\n",
|
| 453 |
-
"improvement_f = t_f - b_f\n",
|
| 454 |
-
"percent_r = (improvement_r / (abs(b_r) + 1e-5)) * 100\n",
|
| 455 |
-
"percent_f = (improvement_f / (abs(b_f) + 1e-5)) * 100\n",
|
| 456 |
-
"\n",
|
| 457 |
-
"print(\"\\n=== FINAL RESULTS (Fair-GRPO-RLVR) ===\")\n",
|
| 458 |
-
"print(f\"Reward — Baseline: {b_r:.3f} | Trained: {t_r:.3f} | Δ {improvement_r:+.3f} ({percent_r:+.1f}%)\")\n",
|
| 459 |
-
"print(f\"Fairness — Baseline: {b_f:.3f} | Trained: {t_f:.3f} | Δ {improvement_f:+.3f} ({percent_f:+.1f}%)\")\n",
|
| 460 |
-
"\n",
|
| 461 |
-
"# Honest conditional verdict\n",
|
| 462 |
-
"if improvement_r > 0 and improvement_f > 0:\n",
|
| 463 |
-
" print(\"\\n✅ Model improved on BOTH reward and fairness.\")\n",
|
| 464 |
-
"elif improvement_r > 0:\n",
|
| 465 |
-
" print(f\"\\n⚠️ Reward improved but fairness REGRESSED by {abs(improvement_f):.3f}. Check reward weights.\")\n",
|
| 466 |
-
"elif improvement_f > 0:\n",
|
| 467 |
-
" print(f\"\\n⚠️ Fairness improved but reward REGRESSED by {abs(improvement_r):.3f}.\")\n",
|
| 468 |
-
"else:\n",
|
| 469 |
-
" print(\"\\n❌ Model did not outperform baseline. Consider more training steps or larger dataset.\")\n",
|
| 470 |
-
"\n",
|
| 471 |
-
"print(\"\\n🏆 Key Insight:\")\n",
|
| 472 |
-
"print(\"Optimizing for fairness improves long-term recovery efficiency.\")\n",
|
| 473 |
-
"\n",
|
| 474 |
-
"print(\"\\n🚀 FINAL TAKEAWAY:\")\n",
|
| 475 |
-
"print(\"Fair-GRPO-RLVR learns policies that outperform greedy baselines by optimizing both efficiency and fairness simultaneously.\")\n",
|
| 476 |
-
"\n"
|
| 477 |
-
]
|
| 478 |
-
}
|
| 479 |
-
],
|
| 480 |
-
"metadata": {
|
| 481 |
-
"kernelspec": {
|
| 482 |
-
"display_name": "Python 3",
|
| 483 |
-
"language": "python",
|
| 484 |
-
"name": "python3"
|
| 485 |
-
},
|
| 486 |
-
"language_info": {
|
| 487 |
-
"name": "python",
|
| 488 |
-
"version": "3.11.0"
|
| 489 |
-
},
|
| 490 |
-
"accelerator": "GPU",
|
| 491 |
-
"colab": {
|
| 492 |
-
"provenance": [],
|
| 493 |
-
"gpuType": "T4"
|
| 494 |
-
}
|
| 495 |
-
},
|
| 496 |
-
"nbformat": 4,
|
| 497 |
-
"nbformat_minor": 4
|
| 498 |
}
|
|
|
|
| 1 |
{
|
| 2 |
+
"cells": [
|
| 3 |
+
{
|
| 4 |
+
"cell_type": "markdown",
|
| 5 |
+
"metadata": {
|
| 6 |
+
"id": "9you219PQC5E"
|
| 7 |
+
},
|
| 8 |
+
"source": [
|
| 9 |
+
"# FairRecovery++: Fair-GRPO-RLVR Training Notebook\n",
|
| 10 |
+
"\n",
|
| 11 |
+
"Research-level training pipeline implementing multi-objective optimization for equitable disaster recovery."
|
| 12 |
+
]
|
| 13 |
+
},
|
| 14 |
+
{
|
| 15 |
+
"cell_type": "code",
|
| 16 |
+
"execution_count": 18,
|
| 17 |
+
"metadata": {
|
| 18 |
+
"id": "-gx3UakxQC5G"
|
| 19 |
+
},
|
| 20 |
+
"outputs": [],
|
| 21 |
+
"source": [
|
| 22 |
+
"# =========================================\n",
|
| 23 |
+
"# 1. INSTALL\n",
|
| 24 |
+
"# =========================================\n",
|
| 25 |
+
"!pip install -q unsloth trl transformers accelerate requests matplotlib pandas pydantic structlog\n",
|
| 26 |
+
"\n"
|
| 27 |
+
]
|
| 28 |
+
},
|
| 29 |
+
{
|
| 30 |
+
"cell_type": "code",
|
| 31 |
+
"execution_count": 19,
|
| 32 |
+
"metadata": {
|
| 33 |
+
"id": "99e5mIVaQC5H"
|
| 34 |
+
},
|
| 35 |
+
"outputs": [],
|
| 36 |
+
"source": [
|
| 37 |
+
"# =========================================\n",
|
| 38 |
+
"# 2. CONFIG\n",
|
| 39 |
+
"# =========================================\n",
|
| 40 |
+
"import os\n",
|
| 41 |
+
"import sys\n",
|
| 42 |
+
"import random\n",
|
| 43 |
+
"import matplotlib.pyplot as plt\n",
|
| 44 |
+
"import pandas as pd\n",
|
| 45 |
+
"import json, re\n",
|
| 46 |
+
"\n",
|
| 47 |
+
"# Clone repo to get local environment\n",
|
| 48 |
+
"REPO_URL = 'https://github.com/joshua400/FairRecovery-PlusPlus.git'\n",
|
| 49 |
+
"REPO_DIR = '/content/FairRecovery-PlusPlus'\n",
|
| 50 |
+
"if not os.path.exists(REPO_DIR):\n",
|
| 51 |
+
" !git clone {REPO_URL} {REPO_DIR}\n",
|
| 52 |
+
"sys.path.insert(0, REPO_DIR)\n",
|
| 53 |
+
"os.chdir(REPO_DIR)\n",
|
| 54 |
+
"\n",
|
| 55 |
+
"MODEL_NAME = \"unsloth/Llama-3.2-1B-Instruct-bnb-4bit\"\n",
|
| 56 |
+
"MAX_STEPS = 20\n",
|
| 57 |
+
"\n"
|
| 58 |
+
]
|
| 59 |
+
},
|
| 60 |
+
{
|
| 61 |
+
"cell_type": "code",
|
| 62 |
+
"execution_count": 20,
|
| 63 |
+
"metadata": {
|
| 64 |
+
"id": "kvtAdWKlQC5I"
|
| 65 |
+
},
|
| 66 |
+
"outputs": [],
|
| 67 |
+
"source": [
|
| 68 |
+
"# =========================================\n",
|
| 69 |
+
"# 3. ENV HELPERS (LOCAL FOR SPEED & RELIABILITY)\n",
|
| 70 |
+
"# =========================================\n",
|
| 71 |
+
"from server.fairrecovery_environment import FairRecoveryEnvironment\n",
|
| 72 |
+
"from fairrecovery_env.models import FairRecoveryAction\n",
|
| 73 |
+
"\n",
|
| 74 |
+
"def reset_env(seed=None, difficulty=None):\n",
|
| 75 |
+
" if difficulty is None:\n",
|
| 76 |
+
" difficulty = random.choice([\"easy\", \"medium\", \"hard\"])\n",
|
| 77 |
+
" env = FairRecoveryEnvironment()\n",
|
| 78 |
+
" obs = env.reset(difficulty=difficulty, seed=seed)\n",
|
| 79 |
+
" return env, obs\n",
|
| 80 |
+
"\n",
|
| 81 |
+
"def step_env(env, action_dict):\n",
|
| 82 |
+
" try:\n",
|
| 83 |
+
" if \"action_type\" not in action_dict:\n",
|
| 84 |
+
" action_dict[\"action_type\"] = \"submit\"\n",
|
| 85 |
+
" if action_dict[\"action_type\"] == \"analyze\" and \"critical_zones\" not in action_dict:\n",
|
| 86 |
+
" action_dict[\"critical_zones\"] = [4, 3]\n",
|
| 87 |
+
" if action_dict[\"action_type\"] == \"allocate\" and \"allocations\" not in action_dict:\n",
|
| 88 |
+
" action_dict[\"allocations\"] = [{\"zone\": 4, \"resource\": \"power\"}]\n",
|
| 89 |
+
"\n",
|
| 90 |
+
" action = FairRecoveryAction(**action_dict)\n",
|
| 91 |
+
" obs = env.step(action)\n",
|
| 92 |
+
" return obs\n",
|
| 93 |
+
" except Exception as e:\n",
|
| 94 |
+
" return env.step(FairRecoveryAction(action_type=\"submit\"))\n",
|
| 95 |
+
"\n"
|
| 96 |
+
]
|
| 97 |
+
},
|
| 98 |
+
{
|
| 99 |
+
"cell_type": "code",
|
| 100 |
+
"execution_count": 21,
|
| 101 |
+
"metadata": {
|
| 102 |
+
"id": "GcyKmIsDQC5I"
|
| 103 |
+
},
|
| 104 |
+
"outputs": [],
|
| 105 |
+
"source": [
|
| 106 |
+
"# =========================================\n",
|
| 107 |
+
"# 4. BASELINE (GREEDY POLICY)\n",
|
| 108 |
+
"# =========================================\n",
|
| 109 |
+
"from inference import greedy_policy\n",
|
| 110 |
+
"\n",
|
| 111 |
+
"def run_baseline(seed=None):\n",
|
| 112 |
+
" # Ensure baseline is evaluated on 'hard' to show the 'Fairness Trap'\n",
|
| 113 |
+
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 114 |
+
" total = 0\n",
|
| 115 |
+
"\n",
|
| 116 |
+
" for _ in range(MAX_STEPS):\n",
|
| 117 |
+
" action = greedy_policy(obs)\n",
|
| 118 |
+
" obs = env.step(action)\n",
|
| 119 |
+
" total += obs.reward\n",
|
| 120 |
+
"\n",
|
| 121 |
+
" if obs.done:\n",
|
| 122 |
+
" break\n",
|
| 123 |
+
"\n",
|
| 124 |
+
" # Honest comparison: return raw total\n",
|
| 125 |
+
" return total, obs.fairness_score\n",
|
| 126 |
+
"\n"
|
| 127 |
+
]
|
| 128 |
+
},
|
| 129 |
+
{
|
| 130 |
+
"cell_type": "code",
|
| 131 |
+
"execution_count": null,
|
| 132 |
+
"metadata": {
|
| 133 |
+
"colab": {
|
| 134 |
+
"base_uri": "https://localhost:8080/",
|
| 135 |
+
"height": 170,
|
| 136 |
+
"referenced_widgets": [
|
| 137 |
+
"158e355936fb48c19b10a82098c3d0b8",
|
| 138 |
+
"a93fc625933a40798df735a3ff438438",
|
| 139 |
+
"05869d84913a49bdabdc546a7f8e1084",
|
| 140 |
+
"d783fe3762e14974932b54ded02c4c01",
|
| 141 |
+
"4a41b45f79314ab5945af366082c1811",
|
| 142 |
+
"f39122df14a041d5a0adc3d335d05e07",
|
| 143 |
+
"56135965433b472db8ef44457a8b56ac",
|
| 144 |
+
"15a9443782b14c0a8abe9ed11fb06786",
|
| 145 |
+
"d087fc944aaf488480988a7d276a811d",
|
| 146 |
+
"8e86cfcf3bb94d01a485f187c813b01d",
|
| 147 |
+
"060c89bb8a2844d091077e9688cd2ab2"
|
| 148 |
+
]
|
| 149 |
+
},
|
| 150 |
+
"id": "MJ5gR4KmQC5J",
|
| 151 |
+
"outputId": "e19c7ea5-60b4-4bff-f4f4-3cc2704b9508"
|
| 152 |
+
},
|
| 153 |
+
"outputs": [
|
| 154 |
+
{
|
| 155 |
+
"output_type": "stream",
|
| 156 |
+
"name": "stdout",
|
| 157 |
+
"text": [
|
| 158 |
+
"==((====))== Unsloth 2026.4.8: Fast Llama patching. Transformers: 5.5.0.\n",
|
| 159 |
+
" \\\\ /| Tesla T4. Num GPUs = 1. Max memory: 14.563 GB. Platform: Linux.\n",
|
| 160 |
+
"O^O/ \\_/ \\ Torch: 2.10.0+cu128. CUDA: 7.5. CUDA Toolkit: 12.8. Triton: 3.6.0\n",
|
| 161 |
+
"\\ / Bfloat16 = FALSE. FA [Xformers = 0.0.35. FA2 = False]\n",
|
| 162 |
+
" \"-____-\" Free license: http://github.com/unslothai/unsloth\n",
|
| 163 |
+
"Unsloth: Fast downloading is enabled - ignore downloading bars which are red colored!\n"
|
| 164 |
+
]
|
| 165 |
+
},
|
| 166 |
+
{
|
| 167 |
+
"output_type": "display_data",
|
| 168 |
+
"data": {
|
| 169 |
+
"text/plain": [
|
| 170 |
+
"Loading weights: 0%| | 0/146 [00:00<?, ?it/s]"
|
| 171 |
+
],
|
| 172 |
+
"application/vnd.jupyter.widget-view+json": {
|
| 173 |
+
"version_major": 2,
|
| 174 |
+
"version_minor": 0,
|
| 175 |
+
"model_id": "158e355936fb48c19b10a82098c3d0b8"
|
| 176 |
+
}
|
| 177 |
+
},
|
| 178 |
+
"metadata": {}
|
| 179 |
+
},
|
| 180 |
+
{
|
| 181 |
+
"output_type": "stream",
|
| 182 |
+
"name": "stderr",
|
| 183 |
+
"text": [
|
| 184 |
+
"Unsloth: Will load unsloth/Llama-3.2-1B-Instruct-bnb-4bit as a legacy tokenizer.\n"
|
| 185 |
+
]
|
| 186 |
+
}
|
| 187 |
+
],
|
| 188 |
+
"source": [
|
| 189 |
+
"# =========================================\n",
|
| 190 |
+
"# 5. LOAD MODEL (UNSLOTH)\n",
|
| 191 |
+
"# =========================================\n",
|
| 192 |
+
"from unsloth import FastLanguageModel\n",
|
| 193 |
+
"\n",
|
| 194 |
+
"model, tokenizer = FastLanguageModel.from_pretrained(\n",
|
| 195 |
+
" model_name=MODEL_NAME,\n",
|
| 196 |
+
" max_seq_length=512,\n",
|
| 197 |
+
" load_in_4bit=True,\n",
|
| 198 |
+
")\n",
|
| 199 |
+
"\n",
|
| 200 |
+
"model = FastLanguageModel.get_peft_model(\n",
|
| 201 |
+
" model,\n",
|
| 202 |
+
" r=16,\n",
|
| 203 |
+
" target_modules=[\"q_proj\",\"k_proj\",\"v_proj\",\"o_proj\",\"gate_proj\",\"up_proj\",\"down_proj\"],\n",
|
| 204 |
+
" lora_alpha=16,\n",
|
| 205 |
+
" use_gradient_checkpointing=\"unsloth\",\n",
|
| 206 |
+
")\n",
|
| 207 |
+
"\n"
|
| 208 |
+
]
|
| 209 |
+
},
|
| 210 |
+
{
|
| 211 |
+
"cell_type": "code",
|
| 212 |
+
"execution_count": null,
|
| 213 |
+
"metadata": {
|
| 214 |
+
"id": "9ZbYLiESQC5J"
|
| 215 |
+
},
|
| 216 |
+
"outputs": [],
|
| 217 |
+
"source": [
|
| 218 |
+
"# =========================================\n",
|
| 219 |
+
"# 6. PROMPT + PARSER\n",
|
| 220 |
+
"# =========================================\n",
|
| 221 |
+
"def build_prompt(obs):\n",
|
| 222 |
+
" zones_str = '\\n'.join([f\"Zone {z.zone_id}: damage={z.damage:.2f}, vulnerable={z.vulnerable_ratio:.2f}\" for z in obs.zones])\n",
|
| 223 |
+
" return f\"\"\"System: You are an AI allocating disaster resources fairly using the Fair-GRPO-RLVR framework.\n",
|
| 224 |
+
"Prioritize Zone 4 (high damage, high vulnerability) over Zone 0 (low damage).\n",
|
| 225 |
+
"Respond ONLY with a JSON action like: {{\"action_type\": \"analyze\", \"critical_zones\": [4, 3]}}\n",
|
| 226 |
+
"\n",
|
| 227 |
+
"User: Day {obs.day}. Budget: {obs.budget_left}.\n",
|
| 228 |
+
"Zones:\n",
|
| 229 |
+
"{zones_str}\n",
|
| 230 |
+
"Fairness Score: {obs.fairness_score}\n",
|
| 231 |
+
"\n",
|
| 232 |
+
"What is your next action?\"\"\"\n",
|
| 233 |
+
"\n",
|
| 234 |
+
"def parse_action(text, stage):\n",
|
| 235 |
+
" if isinstance(text, list):\n",
|
| 236 |
+
" text = text[-1].get(\"content\", str(text))\n",
|
| 237 |
+
"\n",
|
| 238 |
+
" try:\n",
|
| 239 |
+
" match = re.search(r\"\\{.*?\\}\", str(text), re.DOTALL)\n",
|
| 240 |
+
" if match:\n",
|
| 241 |
+
" data = json.loads(match.group())\n",
|
| 242 |
+
" if \"action_type\" not in data:\n",
|
| 243 |
+
" data[\"action_type\"] = stage\n",
|
| 244 |
+
" return data\n",
|
| 245 |
+
" except:\n",
|
| 246 |
+
" pass\n",
|
| 247 |
+
" return {\"action_type\": stage}\n",
|
| 248 |
+
"\n"
|
| 249 |
+
]
|
| 250 |
+
},
|
| 251 |
+
{
|
| 252 |
+
"cell_type": "code",
|
| 253 |
+
"execution_count": null,
|
| 254 |
+
"metadata": {
|
| 255 |
+
"id": "iJsFydGYQC5K"
|
| 256 |
+
},
|
| 257 |
+
"outputs": [],
|
| 258 |
+
"source": [
|
| 259 |
+
"# =========================================\n",
|
| 260 |
+
"# 7. TRAINING REWARD FUNCTION (FAIR-GRPO-RLVR)\n",
|
| 261 |
+
"# =========================================\n",
|
| 262 |
+
"def reward_fn(prompts, completions, **kwargs):\n",
|
| 263 |
+
" rewards = []\n",
|
| 264 |
+
"\n",
|
| 265 |
+
" for prompt, output in zip(prompts, completions):\n",
|
| 266 |
+
" difficulty = random.choice([\"easy\", \"medium\", \"hard\"])\n",
|
| 267 |
+
" env, obs = reset_env(difficulty=difficulty)\n",
|
| 268 |
+
"\n",
|
| 269 |
+
" # FIX: Run the FULL episode using the model's parsed actions.\n",
|
| 270 |
+
" action_dict = parse_action(output, obs.step_stage)\n",
|
| 271 |
+
"\n",
|
| 272 |
+
" for _ in range(MAX_STEPS):\n",
|
| 273 |
+
" obs = step_env(env, action_dict)\n",
|
| 274 |
+
" if obs.done: break\n",
|
| 275 |
+
" action_dict = parse_action(output, obs.step_stage)\n",
|
| 276 |
+
"\n",
|
| 277 |
+
" # 2. Research-Level Fairness Metric (Inverse Service Disparity)\n",
|
| 278 |
+
" services = [z.service for z in env.state.zones]\n",
|
| 279 |
+
" mean_service = sum(services) / len(services)\n",
|
| 280 |
+
" disparity = sum(abs(s - mean_service) for s in services) / len(services)\n",
|
| 281 |
+
" fairness = max(0.0, 1.0 - disparity) # Higher = Better Equity\n",
|
| 282 |
+
"\n",
|
| 283 |
+
" # 3. Multi-objective Components\n",
|
| 284 |
+
" utility = mean_service\n",
|
| 285 |
+
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 286 |
+
"\n",
|
| 287 |
+
" # 4. Total Reward with Curriculum Scaling\n",
|
| 288 |
+
" total = (0.4 * utility + 0.4 * fairness + 0.2 * safety)\n",
|
| 289 |
+
"\n",
|
| 290 |
+
" # FIX: Curriculum weighting without breaking [0,1] normalization\n",
|
| 291 |
+
" difficulty_weight = {\"easy\": 0.8, \"medium\": 1.0, \"hard\": 1.1}.get(difficulty, 1.0)\n",
|
| 292 |
+
"\n",
|
| 293 |
+
" # 5. Stronger Normalization (Preserves Policy Differences)\n",
|
| 294 |
+
" final_score = max(0.0, min(1.0, total * difficulty_weight))\n",
|
| 295 |
+
" rewards.append(float(final_score))\n",
|
| 296 |
+
"\n",
|
| 297 |
+
" return rewards\n",
|
| 298 |
+
"\n"
|
| 299 |
+
]
|
| 300 |
+
},
|
| 301 |
+
{
|
| 302 |
+
"cell_type": "code",
|
| 303 |
+
"execution_count": null,
|
| 304 |
+
"metadata": {
|
| 305 |
+
"id": "qS6wJoUjQC5K"
|
| 306 |
+
},
|
| 307 |
+
"outputs": [],
|
| 308 |
+
"source": [
|
| 309 |
+
"# =========================================\n",
|
| 310 |
+
"# 8. DATASET\n",
|
| 311 |
+
"# =========================================\n",
|
| 312 |
+
"from datasets import Dataset\n",
|
| 313 |
+
"\n",
|
| 314 |
+
"dataset_list = []\n",
|
| 315 |
+
"for i in range(60): # Increased dataset for real learning signal\n",
|
| 316 |
+
" env, obs = reset_env(seed=42 + i)\n",
|
| 317 |
+
" dataset_list.append({\n",
|
| 318 |
+
" \"prompt\": [{\"role\": \"user\", \"content\": build_prompt(obs)}]\n",
|
| 319 |
+
" })\n",
|
| 320 |
+
"\n",
|
| 321 |
+
"dataset = Dataset.from_list(dataset_list)\n",
|
| 322 |
+
"print(f\"Dataset created with {len(dataset)} scenarios.\")\n",
|
| 323 |
+
"\n"
|
| 324 |
+
]
|
| 325 |
+
},
|
| 326 |
+
{
|
| 327 |
+
"cell_type": "code",
|
| 328 |
+
"execution_count": null,
|
| 329 |
+
"metadata": {
|
| 330 |
+
"id": "JwfSgDN6QC5L"
|
| 331 |
+
},
|
| 332 |
+
"outputs": [],
|
| 333 |
+
"source": [
|
| 334 |
+
"# =========================================\n",
|
| 335 |
+
"# 9. TRAIN (GRPO)\n",
|
| 336 |
+
"# =========================================\n",
|
| 337 |
+
"from trl import GRPOTrainer, GRPOConfig\n",
|
| 338 |
+
"\n",
|
| 339 |
+
"config = GRPOConfig(\n",
|
| 340 |
+
" output_dir=\"./outputs\",\n",
|
| 341 |
+
" per_device_train_batch_size=1,\n",
|
| 342 |
+
" gradient_accumulation_steps=2,\n",
|
| 343 |
+
" num_train_epochs=2,\n",
|
| 344 |
+
" max_completion_length=128,\n",
|
| 345 |
+
" logging_steps=1,\n",
|
| 346 |
+
" max_grad_norm=0.5,\n",
|
| 347 |
+
")\n",
|
| 348 |
+
"\n",
|
| 349 |
+
"trainer = GRPOTrainer(\n",
|
| 350 |
+
" model=model,\n",
|
| 351 |
+
" tokenizer=tokenizer,\n",
|
| 352 |
+
" reward_funcs=[reward_fn],\n",
|
| 353 |
+
" args=config,\n",
|
| 354 |
+
" train_dataset=dataset,\n",
|
| 355 |
+
")\n",
|
| 356 |
+
"\n",
|
| 357 |
+
"print(\"🚀 Training Fair-GRPO-RLVR method...\")\n",
|
| 358 |
+
"trainer.train()\n",
|
| 359 |
+
"print(\"✅ Training done\")\n",
|
| 360 |
+
"\n"
|
| 361 |
+
]
|
| 362 |
+
},
|
| 363 |
+
{
|
| 364 |
+
"cell_type": "code",
|
| 365 |
+
"execution_count": null,
|
| 366 |
+
"metadata": {
|
| 367 |
+
"id": "RjP06-7KQC5M"
|
| 368 |
+
},
|
| 369 |
+
"outputs": [],
|
| 370 |
+
"source": [
|
| 371 |
+
"import torch\n",
|
| 372 |
+
"\n",
|
| 373 |
+
"# =========================================\n",
|
| 374 |
+
"# 10. TRAINED MODEL RUNNER\n",
|
| 375 |
+
"# =========================================\n",
|
| 376 |
+
"def run_trained(seed=None):\n",
|
| 377 |
+
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 378 |
+
"\n",
|
| 379 |
+
" for _ in range(MAX_STEPS):\n",
|
| 380 |
+
" prompt = build_prompt(obs)\n",
|
| 381 |
+
" # Use higher temperature for better exploration during evaluation\n",
|
| 382 |
+
" inputs = tokenizer.apply_chat_template([{\"role\": \"user\", \"content\": prompt}], return_tensors=\"pt\", add_generation_prompt=True).to(model.device)\n",
|
| 383 |
+
" outputs = model.generate(\n",
|
| 384 |
+
" inputs,\n",
|
| 385 |
+
" max_new_tokens=100,\n",
|
| 386 |
+
" temperature=0.3, # Increased for exploration\n",
|
| 387 |
+
" top_p=0.9,\n",
|
| 388 |
+
" pad_token_id=tokenizer.eos_token_id\n",
|
| 389 |
+
" )\n",
|
| 390 |
+
" text = tokenizer.decode(outputs[0][inputs.shape[1]:], skip_special_tokens=True)\n",
|
| 391 |
+
" action_dict = parse_action(text, obs.step_stage)\n",
|
| 392 |
+
"\n",
|
| 393 |
+
" obs = step_env(env, action_dict)\n",
|
| 394 |
+
" if obs.done: break\n",
|
| 395 |
+
"\n",
|
| 396 |
+
" # SAME normalized metric as baseline - compute ONCE at episode end\n",
|
| 397 |
+
" services = [z.service for z in env.state.zones]\n",
|
| 398 |
+
" mean_s = sum(services) / len(services)\n",
|
| 399 |
+
" disp = sum(abs(s - mean_s) for s in services) / len(services)\n",
|
| 400 |
+
" fairness = max(0.0, 1.0 - disp)\n",
|
| 401 |
+
" utility = mean_s\n",
|
| 402 |
+
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 403 |
+
" normalized_reward = max(0.0, min(1.0, 0.4 * utility + 0.4 * fairness + 0.2 * safety))\n",
|
| 404 |
+
"\n",
|
| 405 |
+
" return {\n",
|
| 406 |
+
" \"reward\": normalized_reward,\n",
|
| 407 |
+
" \"fairness\": fairness,\n",
|
| 408 |
+
" \"utility\": utility\n",
|
| 409 |
+
" }\n",
|
| 410 |
+
"\n"
|
| 411 |
+
]
|
| 412 |
+
},
|
| 413 |
+
{
|
| 414 |
+
"cell_type": "code",
|
| 415 |
+
"execution_count": null,
|
| 416 |
+
"metadata": {
|
| 417 |
+
"id": "adKxljZAQC5M"
|
| 418 |
+
},
|
| 419 |
+
"outputs": [],
|
| 420 |
+
"source": [
|
| 421 |
+
"# =========================================\n",
|
| 422 |
+
"# 11. RUN COMPARISON (FIXED: Normalized Comparison)\n",
|
| 423 |
+
"# =========================================\n",
|
| 424 |
+
"def run_baseline_normalized(seed=None):\n",
|
| 425 |
+
" \"\"\"Run baseline and return the SAME normalized metric used in training.\"\"\"\n",
|
| 426 |
+
" env, obs = reset_env(seed=seed, difficulty=\"hard\")\n",
|
| 427 |
+
"\n",
|
| 428 |
+
" for _ in range(MAX_STEPS):\n",
|
| 429 |
+
" from inference import greedy_policy\n",
|
| 430 |
+
" action = greedy_policy(obs)\n",
|
| 431 |
+
" obs = env.step(action)\n",
|
| 432 |
+
" if obs.done: break\n",
|
| 433 |
+
"\n",
|
| 434 |
+
" services = [z.service for z in env.state.zones]\n",
|
| 435 |
+
" mean_s = sum(services) / len(services)\n",
|
| 436 |
+
" disp = sum(abs(s - mean_s) for s in services) / len(services)\n",
|
| 437 |
+
" fairness = max(0.0, 1.0 - disp)\n",
|
| 438 |
+
" utility = mean_s\n",
|
| 439 |
+
" safety = max(0.0, 1.0 - obs.info.get(\"violations\", 0) / 10.0)\n",
|
| 440 |
+
" normalized_reward = max(0.0, min(1.0, 0.4 * utility + 0.4 * fairness + 0.2 * safety))\n",
|
| 441 |
+
"\n",
|
| 442 |
+
" return {\n",
|
| 443 |
+
" \"reward\": normalized_reward,\n",
|
| 444 |
+
" \"fairness\": fairness,\n",
|
| 445 |
+
" \"utility\": utility\n",
|
| 446 |
+
" }\n",
|
| 447 |
+
"\n",
|
| 448 |
+
"results = []\n",
|
| 449 |
+
"\n",
|
| 450 |
+
"for i in range(5):\n",
|
| 451 |
+
" test_seed = 2000 + i\n",
|
| 452 |
+
" # Baseline (Normalized for honest comparison)\n",
|
| 453 |
+
" b_res = run_baseline_normalized(seed=test_seed)\n",
|
| 454 |
+
" # Trained\n",
|
| 455 |
+
" t_res = run_trained(seed=test_seed)\n",
|
| 456 |
+
"\n",
|
| 457 |
+
" results.append({\n",
|
| 458 |
+
" \"baseline_reward\": b_res[\"reward\"],\n",
|
| 459 |
+
" \"baseline_fairness\": b_res[\"fairness\"],\n",
|
| 460 |
+
" \"baseline_utility\": b_res[\"utility\"],\n",
|
| 461 |
+
" \"trained_reward\": t_res[\"reward\"],\n",
|
| 462 |
+
" \"trained_fairness\": t_res[\"fairness\"],\n",
|
| 463 |
+
" \"trained_utility\": t_res[\"utility\"]\n",
|
| 464 |
+
" })\n",
|
| 465 |
+
"\n",
|
| 466 |
+
"df = pd.DataFrame(results)\n",
|
| 467 |
+
"print(df)\n",
|
| 468 |
+
"\n"
|
| 469 |
+
]
|
| 470 |
+
},
|
| 471 |
+
{
|
| 472 |
+
"cell_type": "code",
|
| 473 |
+
"execution_count": null,
|
| 474 |
+
"metadata": {
|
| 475 |
+
"id": "WQGtFOKWQC5M"
|
| 476 |
+
},
|
| 477 |
+
"outputs": [],
|
| 478 |
+
"source": [
|
| 479 |
+
"# =========================================\n",
|
| 480 |
+
"# 12. PLOTS (MULTI-COMPONENT)\n",
|
| 481 |
+
"# =========================================\n",
|
| 482 |
+
"os.makedirs(\"plots\", exist_ok=True)\n",
|
| 483 |
+
"\n",
|
| 484 |
+
"fig, ax1 = plt.subplots(figsize=(10, 6))\n",
|
| 485 |
+
"\n",
|
| 486 |
+
"ax1.plot(df[\"baseline_reward\"], label=\"Baseline Reward\", color=\"red\", linestyle=\"--\", marker=\"o\")\n",
|
| 487 |
+
"ax1.plot(df[\"trained_reward\"], label=\"Trained Total Reward\", color=\"green\", marker=\"o\")\n",
|
| 488 |
+
"ax1.set_xlabel(\"Episode\")\n",
|
| 489 |
+
"ax1.set_ylabel(\"Total Reward\")\n",
|
| 490 |
+
"ax1.legend(loc=\"upper left\")\n",
|
| 491 |
+
"\n",
|
| 492 |
+
"ax2 = ax1.twinx()\n",
|
| 493 |
+
"ax2.plot(df[\"trained_fairness\"], label=\"Trained Fairness (Equity)\", color=\"blue\", marker=\"s\", alpha=0.6)\n",
|
| 494 |
+
"ax2.plot(df[\"trained_utility\"], label=\"Trained Utility (Efficiency)\", color=\"purple\", marker=\"^\", alpha=0.6)\n",
|
| 495 |
+
"ax2.set_ylabel(\"Metric Score\")\n",
|
| 496 |
+
"ax2.legend(loc=\"upper right\")\n",
|
| 497 |
+
"\n",
|
| 498 |
+
"plt.title(\"Fair-GRPO-RLVR: Research-Level Performance Metrics\")\n",
|
| 499 |
+
"plt.grid(alpha=0.3)\n",
|
| 500 |
+
"plt.savefig(\"plots/reward_vs_episode.png\", dpi=150, bbox_inches=\"tight\")\n",
|
| 501 |
+
"plt.show()\n",
|
| 502 |
+
"\n",
|
| 503 |
+
"# Fairness Improvement Plot\n",
|
| 504 |
+
"plt.figure(figsize=(8,5))\n",
|
| 505 |
+
"plt.plot(df[\"baseline_fairness\"], label=\"Baseline (Greedy)\", color=\"crimson\", marker=\"o\")\n",
|
| 506 |
+
"plt.plot(df[\"trained_fairness\"], label=\"Trained LLM (Fair-GRPO-RLVR)\", color=\"forestgreen\", marker=\"o\")\n",
|
| 507 |
+
"plt.title(\"Fairness Improvement (Inverse Service Disparity)\")\n",
|
| 508 |
+
"plt.xlabel(\"Episode\")\n",
|
| 509 |
+
"plt.ylabel(\"Fairness Score (higher = better equity)\")\n",
|
| 510 |
+
"plt.axhline(0, color='k', linestyle=':', alpha=0.5)\n",
|
| 511 |
+
"plt.legend()\n",
|
| 512 |
+
"plt.grid(alpha=0.3)\n",
|
| 513 |
+
"plt.savefig(\"plots/fairness_vs_episode.png\", dpi=150, bbox_inches=\"tight\")\n",
|
| 514 |
+
"plt.show()\n",
|
| 515 |
+
"\n"
|
| 516 |
+
]
|
| 517 |
+
},
|
| 518 |
+
{
|
| 519 |
+
"cell_type": "code",
|
| 520 |
+
"execution_count": null,
|
| 521 |
+
"metadata": {
|
| 522 |
+
"id": "aZgx67G0QC5N"
|
| 523 |
+
},
|
| 524 |
+
"outputs": [],
|
| 525 |
+
"source": [
|
| 526 |
+
"# =========================================\n",
|
| 527 |
+
"# 13. SUMMARY\n",
|
| 528 |
+
"# =========================================\n",
|
| 529 |
+
"b_r = df['baseline_reward'].mean()\n",
|
| 530 |
+
"t_r = df['trained_reward'].mean()\n",
|
| 531 |
+
"b_f = df['baseline_fairness'].mean()\n",
|
| 532 |
+
"t_f = df['trained_fairness'].mean()\n",
|
| 533 |
+
"\n",
|
| 534 |
+
"improvement_r = t_r - b_r\n",
|
| 535 |
+
"improvement_f = t_f - b_f\n",
|
| 536 |
+
"percent_r = (improvement_r / (abs(b_r) + 1e-5)) * 100\n",
|
| 537 |
+
"percent_f = (improvement_f / (abs(b_f) + 1e-5)) * 100\n",
|
| 538 |
+
"\n",
|
| 539 |
+
"print(\"\\n=== FINAL RESULTS (Fair-GRPO-RLVR) ===\")\n",
|
| 540 |
+
"print(f\"Reward — Baseline: {b_r:.3f} | Trained: {t_r:.3f} | Δ {improvement_r:+.3f} ({percent_r:+.1f}%)\")\n",
|
| 541 |
+
"print(f\"Fairness — Baseline: {b_f:.3f} | Trained: {t_f:.3f} | Δ {improvement_f:+.3f} ({percent_f:+.1f}%)\")\n",
|
| 542 |
+
"\n",
|
| 543 |
+
"# Honest conditional verdict\n",
|
| 544 |
+
"if improvement_r > 0 and improvement_f > 0:\n",
|
| 545 |
+
" print(\"\\n✅ Model improved on BOTH reward and fairness.\")\n",
|
| 546 |
+
"elif improvement_r > 0:\n",
|
| 547 |
+
" print(f\"\\n⚠️ Reward improved but fairness REGRESSED by {abs(improvement_f):.3f}. Check reward weights.\")\n",
|
| 548 |
+
"elif improvement_f > 0:\n",
|
| 549 |
+
" print(f\"\\n⚠️ Fairness improved but reward REGRESSED by {abs(improvement_r):.3f}.\")\n",
|
| 550 |
+
"else:\n",
|
| 551 |
+
" print(\"\\n❌ Model did not outperform baseline. Consider more training steps or larger dataset.\")\n",
|
| 552 |
+
"\n",
|
| 553 |
+
"print(\"\\n🏆 Key Insight:\")\n",
|
| 554 |
+
"print(\"Optimizing for fairness improves long-term recovery efficiency.\")\n",
|
| 555 |
+
"\n",
|
| 556 |
+
"print(\"\\n🚀 FINAL TAKEAWAY:\")\n",
|
| 557 |
+
"print(\"Fair-GRPO-RLVR learns policies that outperform greedy baselines by optimizing both efficiency and fairness simultaneously.\")\n",
|
| 558 |
+
"\n"
|
| 559 |
+
]
|
| 560 |
+
}
|
| 561 |
+
],
|
| 562 |
+
"metadata": {
|
| 563 |
+
"kernelspec": {
|
| 564 |
+
"display_name": "Python 3",
|
| 565 |
+
"name": "python3"
|
| 566 |
+
},
|
| 567 |
+
"language_info": {
|
| 568 |
+
"name": "python",
|
| 569 |
+
"version": "3.11.0"
|
| 570 |
+
},
|
| 571 |
+
"accelerator": "GPU",
|
| 572 |
+
"colab": {
|
| 573 |
+
"provenance": [],
|
| 574 |
+
"gpuType": "T4"
|
| 575 |
+
},
|
| 576 |
+
"widgets": {
|
| 577 |
+
"application/vnd.jupyter.widget-state+json": {
|
| 578 |
+
"158e355936fb48c19b10a82098c3d0b8": {
|
| 579 |
+
"model_module": "@jupyter-widgets/controls",
|
| 580 |
+
"model_name": "HBoxModel",
|
| 581 |
+
"model_module_version": "1.5.0",
|
| 582 |
+
"state": {
|
| 583 |
+
"_dom_classes": [],
|
| 584 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 585 |
+
"_model_module_version": "1.5.0",
|
| 586 |
+
"_model_name": "HBoxModel",
|
| 587 |
+
"_view_count": null,
|
| 588 |
+
"_view_module": "@jupyter-widgets/controls",
|
| 589 |
+
"_view_module_version": "1.5.0",
|
| 590 |
+
"_view_name": "HBoxView",
|
| 591 |
+
"box_style": "",
|
| 592 |
+
"children": [
|
| 593 |
+
"IPY_MODEL_a93fc625933a40798df735a3ff438438",
|
| 594 |
+
"IPY_MODEL_05869d84913a49bdabdc546a7f8e1084",
|
| 595 |
+
"IPY_MODEL_d783fe3762e14974932b54ded02c4c01"
|
| 596 |
+
],
|
| 597 |
+
"layout": "IPY_MODEL_4a41b45f79314ab5945af366082c1811"
|
| 598 |
+
}
|
| 599 |
+
},
|
| 600 |
+
"a93fc625933a40798df735a3ff438438": {
|
| 601 |
+
"model_module": "@jupyter-widgets/controls",
|
| 602 |
+
"model_name": "HTMLModel",
|
| 603 |
+
"model_module_version": "1.5.0",
|
| 604 |
+
"state": {
|
| 605 |
+
"_dom_classes": [],
|
| 606 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 607 |
+
"_model_module_version": "1.5.0",
|
| 608 |
+
"_model_name": "HTMLModel",
|
| 609 |
+
"_view_count": null,
|
| 610 |
+
"_view_module": "@jupyter-widgets/controls",
|
| 611 |
+
"_view_module_version": "1.5.0",
|
| 612 |
+
"_view_name": "HTMLView",
|
| 613 |
+
"description": "",
|
| 614 |
+
"description_tooltip": null,
|
| 615 |
+
"layout": "IPY_MODEL_f39122df14a041d5a0adc3d335d05e07",
|
| 616 |
+
"placeholder": "",
|
| 617 |
+
"style": "IPY_MODEL_56135965433b472db8ef44457a8b56ac",
|
| 618 |
+
"value": "Loading weights: 100%"
|
| 619 |
+
}
|
| 620 |
+
},
|
| 621 |
+
"05869d84913a49bdabdc546a7f8e1084": {
|
| 622 |
+
"model_module": "@jupyter-widgets/controls",
|
| 623 |
+
"model_name": "FloatProgressModel",
|
| 624 |
+
"model_module_version": "1.5.0",
|
| 625 |
+
"state": {
|
| 626 |
+
"_dom_classes": [],
|
| 627 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 628 |
+
"_model_module_version": "1.5.0",
|
| 629 |
+
"_model_name": "FloatProgressModel",
|
| 630 |
+
"_view_count": null,
|
| 631 |
+
"_view_module": "@jupyter-widgets/controls",
|
| 632 |
+
"_view_module_version": "1.5.0",
|
| 633 |
+
"_view_name": "ProgressView",
|
| 634 |
+
"bar_style": "success",
|
| 635 |
+
"description": "",
|
| 636 |
+
"description_tooltip": null,
|
| 637 |
+
"layout": "IPY_MODEL_15a9443782b14c0a8abe9ed11fb06786",
|
| 638 |
+
"max": 146,
|
| 639 |
+
"min": 0,
|
| 640 |
+
"orientation": "horizontal",
|
| 641 |
+
"style": "IPY_MODEL_d087fc944aaf488480988a7d276a811d",
|
| 642 |
+
"value": 146
|
| 643 |
+
}
|
| 644 |
+
},
|
| 645 |
+
"d783fe3762e14974932b54ded02c4c01": {
|
| 646 |
+
"model_module": "@jupyter-widgets/controls",
|
| 647 |
+
"model_name": "HTMLModel",
|
| 648 |
+
"model_module_version": "1.5.0",
|
| 649 |
+
"state": {
|
| 650 |
+
"_dom_classes": [],
|
| 651 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 652 |
+
"_model_module_version": "1.5.0",
|
| 653 |
+
"_model_name": "HTMLModel",
|
| 654 |
+
"_view_count": null,
|
| 655 |
+
"_view_module": "@jupyter-widgets/controls",
|
| 656 |
+
"_view_module_version": "1.5.0",
|
| 657 |
+
"_view_name": "HTMLView",
|
| 658 |
+
"description": "",
|
| 659 |
+
"description_tooltip": null,
|
| 660 |
+
"layout": "IPY_MODEL_8e86cfcf3bb94d01a485f187c813b01d",
|
| 661 |
+
"placeholder": "",
|
| 662 |
+
"style": "IPY_MODEL_060c89bb8a2844d091077e9688cd2ab2",
|
| 663 |
+
"value": " 146/146 [00:01<00:00, 54.04it/s]"
|
| 664 |
+
}
|
| 665 |
+
},
|
| 666 |
+
"4a41b45f79314ab5945af366082c1811": {
|
| 667 |
+
"model_module": "@jupyter-widgets/base",
|
| 668 |
+
"model_name": "LayoutModel",
|
| 669 |
+
"model_module_version": "1.2.0",
|
| 670 |
+
"state": {
|
| 671 |
+
"_model_module": "@jupyter-widgets/base",
|
| 672 |
+
"_model_module_version": "1.2.0",
|
| 673 |
+
"_model_name": "LayoutModel",
|
| 674 |
+
"_view_count": null,
|
| 675 |
+
"_view_module": "@jupyter-widgets/base",
|
| 676 |
+
"_view_module_version": "1.2.0",
|
| 677 |
+
"_view_name": "LayoutView",
|
| 678 |
+
"align_content": null,
|
| 679 |
+
"align_items": null,
|
| 680 |
+
"align_self": null,
|
| 681 |
+
"border": null,
|
| 682 |
+
"bottom": null,
|
| 683 |
+
"display": null,
|
| 684 |
+
"flex": null,
|
| 685 |
+
"flex_flow": null,
|
| 686 |
+
"grid_area": null,
|
| 687 |
+
"grid_auto_columns": null,
|
| 688 |
+
"grid_auto_flow": null,
|
| 689 |
+
"grid_auto_rows": null,
|
| 690 |
+
"grid_column": null,
|
| 691 |
+
"grid_gap": null,
|
| 692 |
+
"grid_row": null,
|
| 693 |
+
"grid_template_areas": null,
|
| 694 |
+
"grid_template_columns": null,
|
| 695 |
+
"grid_template_rows": null,
|
| 696 |
+
"height": null,
|
| 697 |
+
"justify_content": null,
|
| 698 |
+
"justify_items": null,
|
| 699 |
+
"left": null,
|
| 700 |
+
"margin": null,
|
| 701 |
+
"max_height": null,
|
| 702 |
+
"max_width": null,
|
| 703 |
+
"min_height": null,
|
| 704 |
+
"min_width": null,
|
| 705 |
+
"object_fit": null,
|
| 706 |
+
"object_position": null,
|
| 707 |
+
"order": null,
|
| 708 |
+
"overflow": null,
|
| 709 |
+
"overflow_x": null,
|
| 710 |
+
"overflow_y": null,
|
| 711 |
+
"padding": null,
|
| 712 |
+
"right": null,
|
| 713 |
+
"top": null,
|
| 714 |
+
"visibility": null,
|
| 715 |
+
"width": null
|
| 716 |
+
}
|
| 717 |
+
},
|
| 718 |
+
"f39122df14a041d5a0adc3d335d05e07": {
|
| 719 |
+
"model_module": "@jupyter-widgets/base",
|
| 720 |
+
"model_name": "LayoutModel",
|
| 721 |
+
"model_module_version": "1.2.0",
|
| 722 |
+
"state": {
|
| 723 |
+
"_model_module": "@jupyter-widgets/base",
|
| 724 |
+
"_model_module_version": "1.2.0",
|
| 725 |
+
"_model_name": "LayoutModel",
|
| 726 |
+
"_view_count": null,
|
| 727 |
+
"_view_module": "@jupyter-widgets/base",
|
| 728 |
+
"_view_module_version": "1.2.0",
|
| 729 |
+
"_view_name": "LayoutView",
|
| 730 |
+
"align_content": null,
|
| 731 |
+
"align_items": null,
|
| 732 |
+
"align_self": null,
|
| 733 |
+
"border": null,
|
| 734 |
+
"bottom": null,
|
| 735 |
+
"display": null,
|
| 736 |
+
"flex": null,
|
| 737 |
+
"flex_flow": null,
|
| 738 |
+
"grid_area": null,
|
| 739 |
+
"grid_auto_columns": null,
|
| 740 |
+
"grid_auto_flow": null,
|
| 741 |
+
"grid_auto_rows": null,
|
| 742 |
+
"grid_column": null,
|
| 743 |
+
"grid_gap": null,
|
| 744 |
+
"grid_row": null,
|
| 745 |
+
"grid_template_areas": null,
|
| 746 |
+
"grid_template_columns": null,
|
| 747 |
+
"grid_template_rows": null,
|
| 748 |
+
"height": null,
|
| 749 |
+
"justify_content": null,
|
| 750 |
+
"justify_items": null,
|
| 751 |
+
"left": null,
|
| 752 |
+
"margin": null,
|
| 753 |
+
"max_height": null,
|
| 754 |
+
"max_width": null,
|
| 755 |
+
"min_height": null,
|
| 756 |
+
"min_width": null,
|
| 757 |
+
"object_fit": null,
|
| 758 |
+
"object_position": null,
|
| 759 |
+
"order": null,
|
| 760 |
+
"overflow": null,
|
| 761 |
+
"overflow_x": null,
|
| 762 |
+
"overflow_y": null,
|
| 763 |
+
"padding": null,
|
| 764 |
+
"right": null,
|
| 765 |
+
"top": null,
|
| 766 |
+
"visibility": null,
|
| 767 |
+
"width": null
|
| 768 |
+
}
|
| 769 |
+
},
|
| 770 |
+
"56135965433b472db8ef44457a8b56ac": {
|
| 771 |
+
"model_module": "@jupyter-widgets/controls",
|
| 772 |
+
"model_name": "DescriptionStyleModel",
|
| 773 |
+
"model_module_version": "1.5.0",
|
| 774 |
+
"state": {
|
| 775 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 776 |
+
"_model_module_version": "1.5.0",
|
| 777 |
+
"_model_name": "DescriptionStyleModel",
|
| 778 |
+
"_view_count": null,
|
| 779 |
+
"_view_module": "@jupyter-widgets/base",
|
| 780 |
+
"_view_module_version": "1.2.0",
|
| 781 |
+
"_view_name": "StyleView",
|
| 782 |
+
"description_width": ""
|
| 783 |
+
}
|
| 784 |
+
},
|
| 785 |
+
"15a9443782b14c0a8abe9ed11fb06786": {
|
| 786 |
+
"model_module": "@jupyter-widgets/base",
|
| 787 |
+
"model_name": "LayoutModel",
|
| 788 |
+
"model_module_version": "1.2.0",
|
| 789 |
+
"state": {
|
| 790 |
+
"_model_module": "@jupyter-widgets/base",
|
| 791 |
+
"_model_module_version": "1.2.0",
|
| 792 |
+
"_model_name": "LayoutModel",
|
| 793 |
+
"_view_count": null,
|
| 794 |
+
"_view_module": "@jupyter-widgets/base",
|
| 795 |
+
"_view_module_version": "1.2.0",
|
| 796 |
+
"_view_name": "LayoutView",
|
| 797 |
+
"align_content": null,
|
| 798 |
+
"align_items": null,
|
| 799 |
+
"align_self": null,
|
| 800 |
+
"border": null,
|
| 801 |
+
"bottom": null,
|
| 802 |
+
"display": null,
|
| 803 |
+
"flex": null,
|
| 804 |
+
"flex_flow": null,
|
| 805 |
+
"grid_area": null,
|
| 806 |
+
"grid_auto_columns": null,
|
| 807 |
+
"grid_auto_flow": null,
|
| 808 |
+
"grid_auto_rows": null,
|
| 809 |
+
"grid_column": null,
|
| 810 |
+
"grid_gap": null,
|
| 811 |
+
"grid_row": null,
|
| 812 |
+
"grid_template_areas": null,
|
| 813 |
+
"grid_template_columns": null,
|
| 814 |
+
"grid_template_rows": null,
|
| 815 |
+
"height": null,
|
| 816 |
+
"justify_content": null,
|
| 817 |
+
"justify_items": null,
|
| 818 |
+
"left": null,
|
| 819 |
+
"margin": null,
|
| 820 |
+
"max_height": null,
|
| 821 |
+
"max_width": null,
|
| 822 |
+
"min_height": null,
|
| 823 |
+
"min_width": null,
|
| 824 |
+
"object_fit": null,
|
| 825 |
+
"object_position": null,
|
| 826 |
+
"order": null,
|
| 827 |
+
"overflow": null,
|
| 828 |
+
"overflow_x": null,
|
| 829 |
+
"overflow_y": null,
|
| 830 |
+
"padding": null,
|
| 831 |
+
"right": null,
|
| 832 |
+
"top": null,
|
| 833 |
+
"visibility": null,
|
| 834 |
+
"width": null
|
| 835 |
+
}
|
| 836 |
+
},
|
| 837 |
+
"d087fc944aaf488480988a7d276a811d": {
|
| 838 |
+
"model_module": "@jupyter-widgets/controls",
|
| 839 |
+
"model_name": "ProgressStyleModel",
|
| 840 |
+
"model_module_version": "1.5.0",
|
| 841 |
+
"state": {
|
| 842 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 843 |
+
"_model_module_version": "1.5.0",
|
| 844 |
+
"_model_name": "ProgressStyleModel",
|
| 845 |
+
"_view_count": null,
|
| 846 |
+
"_view_module": "@jupyter-widgets/base",
|
| 847 |
+
"_view_module_version": "1.2.0",
|
| 848 |
+
"_view_name": "StyleView",
|
| 849 |
+
"bar_color": null,
|
| 850 |
+
"description_width": ""
|
| 851 |
+
}
|
| 852 |
+
},
|
| 853 |
+
"8e86cfcf3bb94d01a485f187c813b01d": {
|
| 854 |
+
"model_module": "@jupyter-widgets/base",
|
| 855 |
+
"model_name": "LayoutModel",
|
| 856 |
+
"model_module_version": "1.2.0",
|
| 857 |
+
"state": {
|
| 858 |
+
"_model_module": "@jupyter-widgets/base",
|
| 859 |
+
"_model_module_version": "1.2.0",
|
| 860 |
+
"_model_name": "LayoutModel",
|
| 861 |
+
"_view_count": null,
|
| 862 |
+
"_view_module": "@jupyter-widgets/base",
|
| 863 |
+
"_view_module_version": "1.2.0",
|
| 864 |
+
"_view_name": "LayoutView",
|
| 865 |
+
"align_content": null,
|
| 866 |
+
"align_items": null,
|
| 867 |
+
"align_self": null,
|
| 868 |
+
"border": null,
|
| 869 |
+
"bottom": null,
|
| 870 |
+
"display": null,
|
| 871 |
+
"flex": null,
|
| 872 |
+
"flex_flow": null,
|
| 873 |
+
"grid_area": null,
|
| 874 |
+
"grid_auto_columns": null,
|
| 875 |
+
"grid_auto_flow": null,
|
| 876 |
+
"grid_auto_rows": null,
|
| 877 |
+
"grid_column": null,
|
| 878 |
+
"grid_gap": null,
|
| 879 |
+
"grid_row": null,
|
| 880 |
+
"grid_template_areas": null,
|
| 881 |
+
"grid_template_columns": null,
|
| 882 |
+
"grid_template_rows": null,
|
| 883 |
+
"height": null,
|
| 884 |
+
"justify_content": null,
|
| 885 |
+
"justify_items": null,
|
| 886 |
+
"left": null,
|
| 887 |
+
"margin": null,
|
| 888 |
+
"max_height": null,
|
| 889 |
+
"max_width": null,
|
| 890 |
+
"min_height": null,
|
| 891 |
+
"min_width": null,
|
| 892 |
+
"object_fit": null,
|
| 893 |
+
"object_position": null,
|
| 894 |
+
"order": null,
|
| 895 |
+
"overflow": null,
|
| 896 |
+
"overflow_x": null,
|
| 897 |
+
"overflow_y": null,
|
| 898 |
+
"padding": null,
|
| 899 |
+
"right": null,
|
| 900 |
+
"top": null,
|
| 901 |
+
"visibility": null,
|
| 902 |
+
"width": null
|
| 903 |
+
}
|
| 904 |
+
},
|
| 905 |
+
"060c89bb8a2844d091077e9688cd2ab2": {
|
| 906 |
+
"model_module": "@jupyter-widgets/controls",
|
| 907 |
+
"model_name": "DescriptionStyleModel",
|
| 908 |
+
"model_module_version": "1.5.0",
|
| 909 |
+
"state": {
|
| 910 |
+
"_model_module": "@jupyter-widgets/controls",
|
| 911 |
+
"_model_module_version": "1.5.0",
|
| 912 |
+
"_model_name": "DescriptionStyleModel",
|
| 913 |
+
"_view_count": null,
|
| 914 |
+
"_view_module": "@jupyter-widgets/base",
|
| 915 |
+
"_view_module_version": "1.2.0",
|
| 916 |
+
"_view_name": "StyleView",
|
| 917 |
+
"description_width": ""
|
| 918 |
+
}
|
| 919 |
+
}
|
| 920 |
+
}
|
| 921 |
+
}
|
| 922 |
},
|
| 923 |
+
"nbformat": 4,
|
| 924 |
+
"nbformat_minor": 0
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 925 |
}
|