Factored Approches for MDP & RL

1. Factored Approches for MDP & RL (Some Slides taken from Alan Fern’s course)

2. Factored MDP/RL Representations Advantages Specification: is far easier Inference: Novel lifted versions of the Value and Policy iterations possible Bellman backup directly in terms of ADDs Policy gradient approach where you do direct search in the policy space Learning: Generalization possibilities Q-learning etc. will now directly update the factored representations (e.g. weights of the features) Thus giving implicit generalization Approaches such as FF-HOP can recognize and reuse common substructure • States made of features • Boolean vs. Continuous • Actions modify the features (probabilistically) • Representations include Probabilistic STRIPS, 2-Time-slice Dynamic Bayes Nets etc. • Reward and Value functions • Representations include ADDs, linear weighted sums of features etc.

3. Problems with transition systems • Transition systems are a great conceptual tool to understand the differences between the various planning problems • …However direct manipulation of transition systems tends to be too cumbersome • The size of the explicit graph corresponding to a transition system is often very large • The remedy is to provide “compact” representations for transition systems • Start by explicating the structure of the “states” • e.g. states specified in terms of state variables • Represent actions not as incidence matrices but rather functions specified directly in terms of the state variables • An action will work in any state where some state variables have certain values. When it works, it will change the values of certain (other) state variables

4. State Variable Models • World is made up of states which are defined in terms of state variables • Can be boolean (or multi-ary or continuous) • States are complete assignments over state variables • So, k boolean state variables can represent how many states? • Actions change the values of the state variables • Applicability conditions of actions are also specified in terms of partial assignments over state variables

5. Init: Ontable(A),Ontable(B), Clear(A), Clear(B), hand-empty Goal: ~clear(B), hand-empty Blocks world State variables: Ontable(x) On(x,y) Clear(x) hand-empty holding(x) Initial state: Complete specification of T/F values to state variables --By convention, variables with F values are omitted STRIPS ASSUMPTION: If an action changes a state variable, this must be explicitly mentioned in its effects Goal state: A partial specification of the desired state variable/value combinations --desired values can be both positive and negative Pickup(x) Prec: hand-empty,clear(x),ontable(x) eff: holding(x),~ontable(x),~hand-empty,~Clear(x) Putdown(x) Prec: holding(x) eff: Ontable(x), hand-empty,clear(x),~holding(x) Unstack(x,y) Prec: on(x,y),hand-empty,cl(x) eff: holding(x),~clear(x),clear(y),~hand-empty Stack(x,y) Prec: holding(x), clear(y) eff: on(x,y), ~cl(y), ~holding(x), hand-empty All the actions here have only positive preconditions; but this is not necessary

6. Why is STRIPS representation compact?(than explicit transition systems) • In explicit transition systems actions are represented as state-to-state transitions where in each action will be represented by an incidence matrix of size |S|x|S| • In state-variable model, actions are represented only in terms of state variables whose values they care about, and whose value they affect. • Consider a state space of 1024 states. It can be represented by log21024=10 state variables. If an action needs variable v1 to be true and makes v7 to be false, it can be represented by just 2 bits (instead of a 1024x1024 matrix) • Of course, if the action has a complicated mapping from states to states, in the worst case the action rep will be just as large • The assumption being made here is that the actions will have effects on a small number of state variables. First order Sit. Calc Rel/ Prop STRIPS rep Atomic Transition rep

7. Factored Representations fo MDPs: Actions • Actions can be represented directly in terms of their effects on the individual state variables (fluents). The CPTs of the BNs can be represented compactly too! • Write a Bayes Network relating the value of fluents at the state before and after the action • Bayes networks representing fluents at different time points are called “Dynamic Bayes Networks” • We look at 2TBN (2-time-slice dynamic bayes nets) • Go further by using STRIPS assumption • Fluents not affected by the action are not represented explicitly in the model • Called Probabilistic STRIPS Operator (PSO) model

8. Action CLK

11. Factored Representations: Reward, Value and Policy Functions • Reward functions can be represented in factored form too. Possible representations include • Decision trees (made up of fluents) • ADDs (Algebraic decision diagrams) • Value functions are like reward functions (so they too can be represented similarly) • Bellman update can then be done directly using factored representations..

13. SPUDDs use of ADDs

14. Direct manipulation of ADDs in SPUDD

15. Use heuristic search (and reachability information) LAO*, RTDP Use execution and/or Simulation “Actual Execution” Reinforcement learning (Main motivation for RL is to “learn” the model) “Simulation” –simulate the given model to sample possible futures Policy rollout, hindsight optimization etc. Use “factored” representations Factored representations for Actions, Reward Functions, Values and Policies Directly manipulating factored representations during the Bellman update Ideas for Efficient Algorithms..

16. Probabilistic Planning --The competition (IPPC) --The Action language.. PPDDL was based on PSO A new standard RDDL is based on 2-TBN

18. Not ergodic

21. Reducing Heuristic Computation Cost by exploiting factored representations • The heuristics computed for a state might give us an idea about the heuristic value of other “similar” states • Similarity is possible to determine in terms of the state structure • Exploit overlapping structure of heuristics for different states • E.g. SAG idea for McLUG • E.g. Triangle tables idea for plans (c.f. Kolobov)

22. A Plan is a Terrible Thing to Waste • Suppose we have a plan • s0—a0—s1—a1—s2—a2—s3…an—sG • We realized that this tells us not just the estimated value of s0, but also of s1,s2…sn • So we don’t need to compute the heuristic for them again • Is that all? • If we have states and actions in factored representation, then we can explain exactly what aspects of si are relevant for the plan’s success. • The “explanation” is a proof of correctness of the plan • Can be based on regression (if the plan is a sequence) or causal proof (if the plan is a partially ordered one. • The explanation will typically be just a subset of the literals making up the state • That means actually, the plan suffix from si may actually be relevant in many more states that are consistent with that explanation

23. Triangle Table Memoization • Use triangle tables / memoization C B A A B C If the above problem is solved, then we don’t need to call FF again for the below: B A A B

24. Explanation-based Generalization (of Successes and Failures) • Suppose we have a plan P that solves a problem [S, G]. • We can first find out what aspects of S does this plan actually depend on • Explain (prove) the correctness of the plan, and see which parts of S actually contribute to this proof • Now you can memoize this plan for just that subset of S

25. Relaxations for Stochastic Planning • Determinizations can also be used as a basis for heuristics to initialize the V for value iteration [mGPT; GOTH etc] • Heuristics come from relaxation • We can relax along two separate dimensions: • Relax –ve interactions • Consider +ve interactions alone using relaxed planning graphs • Relax uncertainty • Consider determinizations • Or a combination of both!

26. --Factored TD and Q-learning --Policy search (has to be factored..)

27. Large State Spaces • When a problem has a large state space we can not longer represent the V or Q functions as explicit tables • Even if we had enough memory • Never enough training data! • Learning takes too long • What to do?? [Slides from Alan Fern]

28. Function Approximation • Never enough training data! • Must generalize what is learned from one situation to other “similar” new situations • Idea: • Instead of using large table to represent V or Q, use a parameterized function • The number of parameters should be small compared to number of states (generally exponentially fewer parameters) • Learn parameters from experience • When we update the parameters based on observations in one state, then our V or Q estimate will also change for other similar states • I.e. the parameterization facilitates generalization of experience

29. Linear Function Approximation • Define a set of state features f1(s), …, fn(s) • The features are used as our representation of states • States with similar feature values will be considered to be similar • A common approximation is to represent V(s) as a weighted sum of the features (i.e. a linear approximation) • The approximation accuracy is fundamentally limited by the information provided by the features • Can we always define features that allow for a perfect linear approximation? • Yes. Assign each state an indicator feature. (I.e. i’th feature is 1 iff i’th state is present and i represents value of i’th state) • Of course this requires far to many features and gives no generalization.

30. 10 Example • Consider grid problem with no obstacles, deterministic actions U/D/L/R (49 states) • Features for state s=(x,y): f1(s)=x, f2(s)=y (just 2 features) • V(s) = 0 + 1 x + 2 y • Is there a good linear approximation? • Yes. • 0 =10, 1 = -1, 2 = -1 • (note upper right is origin) • V(s) = 10 - x - ysubtracts Manhattan dist.from goal reward 6 0 0 6

31. 10 But What If We Change Reward … • V(s) = 0 + 1 x + 2 y • Is there a good linear approximation? • No. 0 0

32. 10 But What If We Change Reward … • V(s) = 0 + 1 x + 2 y • Is there a good linear approximation? • No. 0 0

33. 10 But What If… + 3 z • Include new feature z • z= |3-x| + |3-y| • z is dist. to goal location • Does this allow a good linear approx? • 0 =10, 1 = 2 = 0, 0 = -1 • V(s) = 0 + 1 x + 2 y 0 3 0 3 Feature Engineering….

34. Linear Function Approximation • Define a set of features f1(s), …, fn(s) • The features are used as our representation of states • States with similar feature values will be treated similarly • More complex functions require more complex features • Our goal is to learn good parameter values (i.e. feature weights) that approximate the value function well • How can we do this? • Use TD-based RL and somehow update parameters based on each experience.

35. TD-based RL for Linear Approximators • Start with initial parameter values • Take action according to an explore/exploit policy(should converge to greedy policy, i.e. GLIE) • Update estimated model • Perform TD update for each parameter • Goto 2 What is a “TD update” for a parameter?

36. Aside: Gradient Descent • Given a function f(1,…, n) of n real values =(1,…, n) suppose we want to minimize f with respect to  • A common approach to doing this is gradient descent • The gradient of f at point , denoted by  f(), is an n-dimensional vector that points in the direction where f increases most steeply at point  • Vector calculus tells us that  f() is just a vector of partial derivativeswhere can decrease f by moving in negative gradient direction This will be used Again with Graphical Model Learning

37. Aside: Gradient Descent for Squared Error • Suppose that we have a sequence of states and target values for each state • E.g. produced by the TD-based RL loop • Our goal is to minimize the sum of squared errors between our estimated function and each target value: • After seeing j’th state the gradient descent rule tells us that we can decrease error by updating parameters by: squared error of example j target value for j’th state our estimated valuefor j’th state learning rate

38. Aside: continued depends on form of approximator • For a linear approximation function: • Thus the update becomes: • For linear functions this update is guaranteed to converge to best approximation for suitable learning rate schedule

39. Use the TD prediction based on the next state s’ • this is the same as previous TD method only with approximation TD-based RL for Linear Approximators • Start with initial parameter values • Take action according to an explore/exploit policy(should converge to greedy policy, i.e. GLIE) Transition from s to s’ • Update estimated model • Perform TD update for each parameter • Goto 2 What should we use for “target value” v(s)? Note that we are generalizing w.r.t. possibly faulty data.. (the neighbor’s value may not be correct yet..)

40. TD-based RL for Linear Approximators • Start with initial parameter values • Take action according to an explore/exploit policy(should converge to greedy policy, i.e. GLIE) • Update estimated model • Perform TD update for each parameter • Goto 2 • Step 2 requires a model to select greedy action • For applications such as Backgammon it is easy to get a simulation-based model • For others it is difficult to get a good model • But we can do the same thing for model-free Q-learning

41. Q-learning with Linear Approximators Features are a function of states and actions. • Start with initial parameter values • Take action a according to an explore/exploit policy(should converge to greedy policy, i.e. GLIE) transitioning from s to s’ • Perform TD update for each parameter • Goto 2 • For both Q and V, these algorithms converge to the closest linear approximation to optimal Q or V.

46. Policy Gradient Ascent • Let () be the expected value of policy . • () is just the expected discounted total reward for a trajectory of . • For simplicity assume each trajectory starts at a single initial state. • Our objective is to find a  that maximizes () • Policy gradient ascent tells us to iteratively update parameters via: • Problem: ()is generally very complex and it is rare that we can compute a closed form for the gradient of (). • We will instead estimate the gradient based on experience

47. Gradient Estimation • Concern: Computing or estimating the gradient of discontinuous functions can be problematic. • For our example parametric policy is () continuous? • No. • There are values of  where arbitrarily small changes, cause the policy to change. • Since different policies can have different values this means that changing  can cause discontinuous jump of ().

48. Example: Discontinous () • Consider a problem with initial state s and two actions a1 and a2 • a1 leads to a very large terminal reward R1 • a2 leads to a very small terminal reward R2 • Fixing 2 to a constant we can plot the ranking assigned to each action by Q and the corresponding value () Discontinuity in () when ordering of a1 and a2 change R1 () R2 1 1

49. Probabilistic Policies • We would like to avoid policies that drastically change with small parameter changes, leading to discontinuities • A probabilistic policy  takes a state as input and returns a distribution over actions • Given a state s (s,a) returns the probability that  selects action a in s • Note that () is still well defined for probabilistic policies • Now uncertainty of trajectories comes from environment and policy • Importantly if (s,a) is continuous relative to changing  then () is also continuous relative to changing  • A common form for probabilistic policies is the softmax function or Boltzmann exploration function Aka Mixed Policy (not needed for Optimality…)

50. Empirical Gradient Estimation • Our first approach to estimating  () is to simply compute empirical gradient estimates • Recall that  = (1,…, n) and so we can compute the gradient by empirically estimating each partial derivative • So for small  we can estimate the partial derivatives by • This requires estimating n+1 values: