Title: POLICY MESSAGE PASSING: A NEW ALGORITHM FOR PROBABILISTIC GRAPH INFERENCE Anonymous authors Paper under double-blind review

Abstract: A general graph-structured neural network architecture operates on graphs through two core components: (1) complex enough message functions; (2) a fixed information aggregation process. In this paper, we present the Policy Message Passing algorithm, which takes a probabilistic perspective and reformulates the whole information aggregation as stochastic sequential processes. The algorithm works on a much larger search space, utilizes reasoning history to perform inference, and is robust to noisy edges. We apply our algorithm to multiple complex graph reasoning and prediction tasks and show that our algorithm consistently outperforms state-of-the-art graph-structured models by a significant margin.

1. INTRODUCTION:

Not every path is created equal. Powerful sequential inference algorithms have been a core research topic across many tasks that involve partial information, large search spaces, or dynamics over time. An early algorithm example is dynamic programming which memorizes the local information and uses a fixed inference trace to acquire a global optimum. In the modern deep learning era, sequential memory architectures such as recurrent neural networks are used for sequential inference (e.g. language generation (Sutskever et al. (2014) )) to narrow down the search space. Neural-based iterative correction is developed in amortized inference for posterior approximation (Marino et al. (2018) ). Also, a large bulk of works are developed in reinforcement learning for sequential inference of the environment with partial observation, non-future aware states, and world dynamics.

Graph-structured models are another class of models that heavily rely on local observations and sequential inference. It contains nodes and message functions with partial information, the final output of model relies on a sequence of communication operations which transforms the nodes information dynamically over iterations. A general formulation of learning graph-structured models is maximizing the likelihood: max θ p(y|s 0 , x; F θ , P )

(1)

The objective describes the process where it takes in the information x over a graph, starts from initial states s 0 and maps to the target output y. This process involves two key components, message functions F and inference process P .

Message function designs focus on the choices of functions defined over the interaction between nodes or variables. Simple log-linear potential functions are commonly used in traditional probabilistic graphical models (Sutton et al. (2012) ). Kernelized potential functions (Lafferty et al. (2004) ), gaussian energy functions (Krähenbühl & Koltun (2011)) , or semantic-driven potential functions (Lan et al. (2011) ) are developed in varied cases to adapt to specific tasks. Deep graph neural networks generalizes the message functions to neural networks. There are a variety of function types developed such as graph convolution operations (Kipf & Welling (2016) ), attention-based functions (Veličković et al. (2017) ), gated functions (Li et al. (2015) ; Deng et al. (2016) ).

Inference procedures over graphs are mainly developed under the context of probabilistic graphical models, with the classic inference techniques such as belief propagation ( Pearl (1988) ), generalized belief propagation (Yedidia et al. (2001) ), and tree-reweighted message passing ( Wainwright et al. (2003) ). However, there has not been many research in developing powerful inference process in modern deep graph-structured models (Graph Neural Networks). Most graph-structured models still follow a fixed hand-crafted rule for performing inference over nodes.

In this paper, we take a different tack and instead propose to represent the inference process modeled by a neural network P φ :

EQUATION

The neural network P φ represents the distribution over all possible inference trajectories in the graph model. By taking the distribution perspective of message passing, the whole inference is reformulated as stochastic sequential processes in the deep graph-structured models. This formulation introduces several properties: Distribution nature -the possible inference trajectories in graph, besides synchronous or asynchronous, are in a vast space with varied possibilities, leading to different results; Consistent inference -maintaining a powerful neural network memorizing over time could lead to more effective and consistent inference over time; Uncertainty -inference itself will affect prediction, with a distribution, such uncertainty over the reasoning between nodes can be maintained.

The primary contributions of this paper are three folds. Firstly, we introduce a new perspective on message passing that generalizes the whole process to stochastic inference. Secondly, we propose a variational inference based framework for learning the inference models. Thirdly, We empirically demonstrate that having a powerful inference machine can enhance the prediction in graph-structured models, especially on difficult reasoning tasks, or graphs with noisy edges.

2. BACKGROUND: GRAPH NEURAL NETWORKS:

A large volume of Graph Neural Network (GNN) models have been developed recently (Defferrard et al. (2016) ; Li et al. (2015) ; Kipf & Welling (2016) ; Henaff et al. (2015) ; Duvenaud et al. ( 2015)) for varied applications and tasks. This line of methods typically share some common structures such as propagation networks and output functions. In our paper, without loss of generality, we follow a unified view of graph neural network proposed in a recent survey (Battaglia et al. (2018) ) and introduce the graph neural network architecture here.

In this formulation, there are two key components: propagation networks and output networks. In the propagation networks, the fundamental unit is a "graph-to-graph" module. This module takes a graph with states s t as input and performs information propagation (message passing) operations to modify the values to s t+1 . The graph states at each step t consists of sets of variables representing the attributes of nodes V t = {v} t and edges E t = {e} t , along with global properties u t . The graph-to-graph module updates the attributes of these nodes, edges, and global properties via a set of message passing, implemented as neural networks.

EQUATION

The message passing in the standard graph neural networks is a fixed aggregation process defined by the graph structure E. The output network is a general part of model that maps the information in the graph states to the target output y: p(y|s T ). In general, there is no restriction on the form of target y or the output mapping function, like the target can be a set of labels for nodes, edges, the whole graph or even sentences.

3. GRAPH MODELS WITH POLICY MESSAGE PASSING:

In this section we describe our stochastic graph inference algorithm under the formulation of graph neural networks. As the inference process under the algorithms is more of a reasoning procedure, we refer to it as reasoning in the later content. Our approach takes the initial state s 0 of a graph representing observations and latent variables in a domain, uses a set of parameterized reasoning agents to pass messages over this graph, and infers the output prediction y. Our model contains three main parts:

1. Graph s t , representing states of observed and latent variables. 2. A set of message functions {f k (•, •; θ k )}, parameterized by θ k , that determine the possible manner in which messages can be passed over edges in the graph.

3. An automatic reasoning agent A(•; φ), parameterized by φ, that determines which message function would be appropriate to use along an edge in the graph.

For every timestep t during reasoning, our model takes the current graph state s t and the previous reasoning history, then generates a distribution over actions for choosing which messages to pass (or pass no information). The new graph state s t+1 is obtained through performing message passing actions sampled from the distribution. The final prediction is thus inferred sequentially through the operation of the reasoning agents over the graph. Fig. 1 illustrates the algorithm.

3.1. MODEL DEFINITION:

Our inference procedure operates on graphs. We start with an initial graph state s 0 = (v 0 , e 0 , u 0 ), where v 0 , e 0 , and u 0 contain information about individual nodes, edges, and global information about the graph respectively. Nodes in this graph can represent observations or latent variables from our domain; we assume the set of nodes is known and fixed. Edges in this graph represent possible relationships between nodes, over which we will potentially pass messages during inference. The edges are assumed known and fixed, though our method will permit flexibility via its inference.

Given a graph state s, we can output y using a network f o (s; θ o ). Note that in general this formulation could have a variety of possible output functions, corresponding to different tasks. In the current instantiation we focus on a single such output function that is used for a single task. Given the initial graph state s 0 , we perform a set of reasoning steps that pass messages over this graph to obtain a refined graph representation. The goal of this message passing is to successively refine the graph state such that the graph state can be accurately decoded to match the correct output for the given input graph. The update functions we define will modify the node attributes v, though the formulation could be generalized to edge and global attributes as well.

Message Formulation: The message passing is decomposed into operations on edges in the graph.

For the i th node in the graph, we update its state by passing messages from all nodes in its neighbour set N (i).

v i = j∈N (i) f z ij (v i , v j ) (4)

The form of message to pass is a crucial component of our method. We select the function f z ij (•, •) to use in this message passing for each edge.

We define a common set of functions {f k (•, •)} 0≤k≤K that can be used across any edge in the graph.

To allow for no information to be passed along an edge, we define the function f 0 (•, •) = 0. Each of the other functions f k (•, •) is a non-linear function (e.g. MLP) with its own learnable parameters θ k .

Reasoning Agents: Alongside these functions, we define a reasoning agent that determines which function to use on a particular edge in the graph at the current time step of inference. We run one copy of this reasoning agent for each edge of the graph. The reasoning agent takes as input the current graph state and for its edge (i, j) outputs a distribution over the message to pass on that edge (computed by a network A). The reasoning agent maintains state h ij t (propagated via a network B) over iterations of inference so that it can potentially remember which types of message passing actions were taken at different time steps t. We also include a global state h u t = ij h ij t which can be used to share information across the agents operating on different edges.

The networks are updated and distributions computed using the following equations:

h ij t+1 = B(s t , i, j, h ij t ; θ B ) (5) π(z ij t+1 |s t ) ∝ A(s t , i, j, h ij t , h u t ; θ A ) + 1 (6)

where π(z ij t+1 |s t ) is a distribution over z ij t+1 , the message function to use for edge (i, j) on step t + 1 given the current graph state. Note that a uniform prior is added to message distribution for regularization, though in general other priors could be used for this purpose. The parameters φ = (θ A , θ B ) for the reasoning agent are shared; however, each edge maintains its own copy of these agents, via different states h ij t . Action sampling: Given these distributions π, we can now use them to perform inference by sampling actions from these distributions. In our implementation, the Gumbel softmax trick is used for differentiable sampling.

By executing a set of actions {z t } sampled from each of these distributions π, our model takes selected nodes and corresponding message functions to pass, and the graph state is updated to s t+1 after one reasoning/passing step. Over a sequence of time steps, this process corresponds to a reasoning trajectory τ = (z 1 , z 2 , . . . , z T ) that produces a sequence of graph states (s 0 , s 1 , . . . , s T ). At test time, we evaluate a learned model by sampling multiple such reasoning trajectory sequences and aggregating the results.

3.2. LEARNING OBJECTIVE:

We assume that we are provided with training data that specify an input graph structure along with its corresponding value to be inferred. For ease of exposition, we derive our learning objective for a single training data pair of graph s with target value y. For the learning objective, we adopt a probabilistic view on reasoning. As per the previous section, we can infer these values by running a particular reasoning process. However, rather than fixing a reasoning process, we could instead imagine inferring values y by marginalizing over all possible reasoning sequences τ .

EQUATION

The "prior" p(τ |s) could be either a learnable prior or a hand-designed prior. Belief propagation or other neural message passing frameworks define a delta function where the prior is deterministic and fixed. Here we instead adopt a learnable prior π for p(τ |s), where the reasoning process depends on the graph states. τ represents the latent reasoning sequence. Given a graph initial state and ground truth label y, we are interested in the posterior distribution over trajectories that describes an efficient and effective reasoning process. We adopt standard variational inference to approximate the posterior distribution and try to find such action spaces.

As per the previous section, we have a method for computing y given a reasoning sequence τ . However, as in variational inference applications, evaluation of the summation in Eq. 7 is challenging. Standard sampling approaches, i.e. choosing an inference procedure from an arbitrary proposal distribution, would be likely to lead to inference procedures ineffective for the given graph/inference task. Instead, we utilize our reasoning agent as a proposal distribution to propose inference procedures.

If the proposal distribution is close to the posterior distribution over reasoning actions, then we can conduct both effective and efficient inference to obtain the final desired output. Hence, we hope to minimize the following KL divergence:

KL(q(τ )||π * (τ |y, s)) (8)

where q is the proposal distribution over reasoning actions. π * is the true posterior distribution over reasoning sequences given the target inference values y and graph state s, which is unknown. We adopt the standard variational inference method to derive the following equations:

log p(y|s) -KL(q(τ |y, s)||π * (τ |y, s)) = E τ ∼q log p(y|τ , s) -KL(q(τ |s, y)||π(τ |s)) (9)

We can optimize the right hand side, which is also known as the evidence lowerbound (ELBO). The log likelihood and KL divergence in Eq. 8 are then implicitly optimized. We can also decide to define the learning objective as a product of probabilities over the T time steps in the reasoning process. In this case, the objective would be:

EQUATION

To summarize, we learn the parameters of the networks f o and {f k } which specify p(y|τ , s); and A and B which specify π(τ |s). This is accomplished by maximizing the ELBO in Eq. 10 over all graph-target value pairs in the training data.

Interpretation of Learning Objective: The learning objective consists in two terms. Likelihood term: This term provides main objective for the agent and model to maximize. The policy should be learned to perform accurate reasoning to fit the data and maximize the likelihood of the ground truth values for y. KL term: The policy we learn in proposal distribution can be used to guide the learnable prior agent toward good inference procedures. The KL divergence term in the objective function can be used to achieve this.

3.3. MIXED SAMPLING FROM SEQUENTIAL DISTRIBUTIONS:

As the whole inference over graphs is a sequential procedure, there is a common issue that will appear in the above formulation: the drifting problem when facing unfamiliar graph states in some inference step t. The proposal distribution in the variational inference formulation has access to ground truth data, which will lead to "good" graph states in most cases, while the prior agent will probably face some graph states that will lead to failure predictions and even drift away further. Similar issues have been found in sequential data generation (Lamb et al. (2016) ) or imitation learning (Ross et al. (2011) ). In our model, we adopt a easy fix by mix the sampling from both prior and proposal in the trajectory level, to have graph neural networks visit different states, and by firstly rollout using prior agent, then use the proposal agent to complete the rest trajectory. The objective is re-weighted using importance sampling weights when mixing the trajectories.

4. EXPERIMENTS:

We test our proposed model on two types of modalities of data: complex visual reasoning tasks and social network classification. Visual data has more information and more complex possible interactions or variations between nodes in the graph, while social network has larger scale of graphs.

4.1. VISUAL REASONING:

We introduce two visual reasoning tasks: (1) visual cue-based position reasoning and (2) large-scale image puzzle solving. To solve these visual based reasoning tasks, we require a model to have a strong ability in handling imperfect observations, performing inference, handling noisy edges, and flexibly making modifications to graph nodes.

Baseline methods::

We compare our proposed model with several popular graph neural network based variants. These models all more or less have the ability to learn and control the information flow, and have been demonstrated to have strong performance on many tasks or applications. (a) Illustration of the "Where am I" dataset and tasks. Each object is extracted from the original image along with its context. Between some patches there will be similar patterns appear. • Graph Neural Network (Battaglia et al. ( 2018)). We follow the graph neural network model described in the survey Battaglia et al. (2018) . The graph neural network contains three update functions: node, edge, global states. The aggregation process we used is mean pooling of messages.

• Attention-based Graph Neural Network (Veličković et al. (2017) ; Wang et al. (2018) ; Vaswani et al. (2017) ). An attention mechanism is used for modulating the information flow from different messages and can be viewed as a non-probabilistic single-step version of our model. We implement the attention by using sigmoid gates over different message functions.

• Recurrent Relational Network (Battaglia et al. ( 2018)). We follow the open source code and implement the recurrent relational network (RRN) model. The RRN model updates the node states through a Gated Recurrent Unit and adds dense supervision across all timesteps.

• Neural Relational Inference model (Kipf et al. (2018) ). The Neural Relational Inference model is the first proposed probabilistic model that captures the distribution over graph structures. We follow the same implementation and use four function types on each edge with entropy as an extra penalty in the objective function.

4.1.1 "Where am I": VISUAL POSITION REASONING Task design: We firstly design a synthetic dataset named "Where am I". The dataset is based on the Shapes Andreas et al. (2016) dataset, initially proposed in the Visual Question Answering (VQA) context. The dataset contains 30K images with shape objects randomly dropped on k × k grids. We assume that each shape object has a 3 × 3 context field of view and ensure that all those shape objects are connected within the visible context. The dataset is set up to have two grid settings, 9 × 9 and 12 × 12, which contain 25 objects and 36 objects respectively. For each image, the position for one object is known and the model needs to perform inference over objects by comparing the context to infer the correct position of other objects.

Model setup and problem difficulties: We build a graph to solve this problem. Each node in the graph represents one object, one of them contains a known position tuple (x, y), and the positions for the rest obejcts are unknown. The graph is assumed to be fully connected. We use graph networks with 50 hidden dimension on both node and edge, ELU as activation functions. Each model takes 7 and 8 steps of inference for 25 object and 36 object cases respectively. The difficulty of this problem is that for each object, there are 16 ways/configurations to overlap with another object. Each position configuration can contain N candidates. If N is larger than 1, there will be confusion, and models need to solve other objects first to provide more context information. The task in general needs model to perform iterative modification of nodes content and avoid false messages.

Experiment results: Table 2b shows the performance of each model measured by the classification accuracy for each object. Our proposed Policy Message Passing (PMP) model and its variants achieve a large performance boost over the baselines, demonstrating that there exists a strong need for this kind of dynamical reasoning in complex tasks. The models with fixed structure and rigid aggregation processes tend to not learn or generalize well. Attention based models slightly alleviate that problem, but are not able to perform consistent and strategic reasoning over objects. Experiment results: Tables 1 and 4a summarize the results of all models across varied task settings.

Our proposed models consistently outperform baseline models by a 10% -20% margin, even when baseline models have much larger capacity on the graphs. This is a further indication that having a powerful inference process is crucial to deep graph-structured models.

Analysis on number of inference steps and size of function set: We further examine the impact of inference steps and number of functions in function set on our models, shown in 4b. We found that empirically a function set of size 4 is enough for message diversity. For inference step analysis, we test our PMP model on 6 × 6 puzzles and found with 6 steps the model can already achieve a good performance. Similar behaviours are found for baseline models too. We only plot our main model in the figure. Comparison with state-of-the-art puzzle solving method: We compared our method with stateof-the-art deep learning based image puzzle solving on the CelebA dataset. The results are shown in Table 3b . Our model shows near perfect puzzle solving results on 2 × 2, 3 × 3, 4 × 4 and 5 × 5 settings.

4.2. SOCIAL NETWORK CLASSIFICATION WITH NOISY EDGES:

We further inspect on models' ability on handling noisy edges on large scale graphs. Current graph neural network benchmark datasets on social networks, such as Cora (Sen et al. (2008) ), heavily rely on well caliberated edge matrix as a strong prior information for inference. The burden on aggregation process are less, compare to more general cases where adjacency matrix information are often noisy or missing. In this experiment, we compare our model's performance with two state-ofthe-art graph-structured models: (1) Graph Convolutional Networks (GCN) and (2) Graph Attention Networks (GAT) on a standard benchmark Cora. Instead of using only the well designed adjacency matrix information, we also generate perturbation to the matrix by adding {%50, %100, %200} * |E| more noisy edges, |E| is the number of edges in the original graph. For implementation, we use the official release of the codes on two baseline models and default settings. In our proposed model, we follow the same hyper-parameters in graph networks and uses recurrent network agent with 32 hidden dimensions. Results summarized in table 2 show that our model is more robust to noisy edges in performing inference over large scale graphs. 

Models:



5. CONCLUSION:

Reasoning over graphs can be formulated as a learnable process in a Bayesian framework. We described a reasoning process that takes as input a graph and generates a distribution over messages that could be passed to propagate data over the graph. Both the functional form of the messages and the parameters that control the distribution over these messages are learned. Learning these parameters can be achieved via a variational approximation to the marginal distribution over provided graph label. Future possibilities include adapting our model to use a more general set of functions where some of them is not differentiable.

:



Figures and Tables:

z |v i ,v j ,h) Reasoning history (a) Single timestep action proposal.

Dynamic graph reasoning over multiple steps.

Figure 1: We design agents that learn to reason over graphs by sequentially performing actions for inference. (a) We learn a function set. Each function is a possible message along an edge. At each time step of inference, our reasoning agent examines the current graph states and produces a distribution over the function set to sample from. (b) A set of these agents, each one considering its own edge, is run to produce the full set of messages for a graph. Each agent maintains its own state and can communicate with the other agents via a shared memory.

Figure 2: Our "Where am I" dataset(left) and results compared with baseline methods measured by per object accuracy(right).

We compare our model's performance with state-of-the-art deep learning based puzzle solving model on Celeba dataset.

Figure 3: Image puzzle solving problem(left) and results compared with state-of-the-art method on Celeba dataset(right).

(a) Performance results on Visual Genome dataset, measured by per-patch classification accuracy and correlation score Kendall-tau. Results are shown across different puzzle sizes. (b) Analysis of our model on number of inference steps and number of functions.

Figure 4: Performance on Visual Genome(left) and model analysis plots(right).

4.1.2 IMAGE PUZZLE SOLVING ON BENCHMARKSTask design and datasets: Next we move on to large-scale benchmark datasets with real images containing complex visual information. We test our model on the Jigsaw Puzzle solving task for real images. Given an image with size H × W , we divide the image into d × d non-overlapping patches. The d 2 patches are un-ordered and need to be reassembled into the original image. We use two datasets to mainly test our model: (1) Visual Genome -The Visual Genome dataset contains 64346 images for training and 43903 images for testing. It is often used in image based scene graph generation tasks, indicating that there exists very rich interactions and relationships between visual elements in the image. We test the models on three different puzzle sizes: 3 × 3, 4 × 4, and 6 × 6 (2) COCO dataset -We use the COCO dataset (2014 version) that contains 82783 images for training and 40504 images for testing. We set up three tasks, 3 × 3, 4 × 4 and 6 × 6 for evaluation. To evaluate the models' performance, we use perpatch classification score over the d 2 classes and measure the Kendall-tau correlation score between predicted patch order and the ground-truth order. For all graph models, we use 50 hidden dimensions on both nodes and edges. Models will perform d + 2 steps of inference. For all baseline methods, we use 200 dimensions for nodes and messages, while still remain 50 dimensions for our main model. Performance results on COCO dataset, measured by per-patch classification accuracy and correlation score Kendall-tau. We use Kendall-tau×100 here for better comparison. Results are shown across different puzzle sizes.

Table content: <table><tr><td colspan="4">Evaluation metric and model implementation: Models Accuracy(%)</td><td/><td>Kendall-tau</td></tr><tr><td/><td>3x3</td><td>4x4</td><td>6x6</td><td>3x3</td><td>4x4</td><td>6x6</td></tr><tr><td>random</td><td colspan="2">11.11 6.25</td><td>2.78</td><td>-</td><td>-</td><td>-</td></tr><tr><td>GNN</td><td colspan="6">93.86 77.37 50.81 94.54 82.54 74.55</td></tr><tr><td colspan="7">GNN-attention 91.29 81.36 68.39 92.70 86.56 79.11</td></tr><tr><td>RRN</td><td colspan="6">89.91 78.14 62.52 85.33 84.46 76.63</td></tr><tr><td>NRI</td><td colspan="6">87.63 81.93 69.01 89.12 86.67 79.78</td></tr><tr><td>PMP</td><td colspan="6">98.31 97.96 90.53 99.27 98.67 96.01</td></tr></table>

± 1.3 23.53 ± 1.44 22.62 ± 1.84 GAT 58.83 ± 1.35 45.82 ± 2.55 34.53 ± 0.45 PMP 65.37 ± 1.36 61.42 ± 1.31 50.3 ± 1.13 The performance of node classification for Cora dataset under varied noisy edge settings. The experiment follows the standard transductive setup in Yang et al. (2016).

Table content: <table><tr><td/><td/><td>Noisy edges</td></tr><tr><td/><td>50%</td><td>100%</td><td>200%</td></tr><tr><td>GCN</td><td>26.55</td><td/></tr></table>

Peer review (iclr2020):
This paper introduces policy message passing: a graph neural network with multiple message types and an inference mechanism that assigns messages of a particular type to edges in a recurrent fashion. Experiments on visual reasoning tasks and on a noisy citation network classification task indicate competitive performance. There are many interesting ideas in this paper and I believe that the overall approach/idea has merit and is interesting for the community. In its current state, however, the paper has issues in terms of clarity of writing, technical and notational mistakes, unclear experimental settings and unclear implementation of baseline models, and needs a major revision, hence I recommend reject. This feedback is mostly meant as an encouragement to “polish” the ideas outlined in the paper and their presentation, as I believe that these can be quite impactful, if presented well and once all technical issues have been addressed. Discussion of the main idea outlined in this paper: ——————————————————————— The proposed dynamic message passing procedure (termed “reasoning agents” in the paper) can be understood as an extension of Neural Relational Inference (NRI; Kipf et al., ICML 2018), where a Gumbel softmax-based attention mechanism in parallel softly assigns a particular message function per edge in the graph. The main difference to NRI is that this is done sequentially over multiple message passing steps in the graph, whereas NRI (in the training phase) only performs one single such step (for all edges in the graph). Getting this to work well in practice is a challenging problem, as sequential sampling steps can accumulate biases in the approximate inference scheme. The authors propose to address this issue in a number of ways: the attention mechanism per edge is coupled over time steps via a local hidden state and coupled to all other edges via a global hidden state (this is novel), a vector of ones is added to the action probabilities and subsequently re-normalized to act as a regularizer (this appears to be a hack), they use a learned prior distribution instead of a fixed prior as in NRI, and lastly, the authors introduce a mixed sampling strategy which mixes samples from the (learned) prior distribution for edge types and the proposal distribution while re-weighing the objective with importance sampling weights — for this last contribution (mixing + reweighing) the description in the paper misses some details (how precisely is the mixing done over the course of training?). Clarity, technical correctness, experimental evaluation ——————————————————————————— * The paper contains a lot of spelling mistakes and grammatical errors, making it difficult to read. * Notational issues: Eq 3 is unclear (no node/edge indices provided, function arguments unspecific — all three functions have the same arguments); Eq. 8/9 switches between posteriors p, true posteriors pi^* (what’s the difference from p?) and priors pi without explanation— a derivation would be helpful. Eq. 10 only places a prior on the first time step — where does this come from? A derivation would be helpful. * “Belief propagation or other neural message passing frameworks define a delta function where the prior is deterministic and fixed” -> this is not the case for NRI * What temperature is used for the Gumbel softmax distribution? Does the forward pass use “hard” (straight-through trick) or soft activations? * Eq. 5: Why condition B on the indices i and j? This will likely not generalize to new edges * Eq. 6: What is the motivation for adding a vector of ones? Wouldn*t it be possible to achieve the desired regularization effect by increasing the softmax temperature? * How are the edge and global hidden states initialized? * The explanation for the mixed sampling procedure is unclear, please elaborate (potentially use an appendix) * What are the train/val/test splits, how are hyper parameters optimized in all experiments? * Why implement attention-based GNN baselines with sigmoid gates even though you cite Transformers and Graph Attention Networks (GAT) in this paragraph? I suspect that multi-head dot product attention (default for Transformers/GAT) instead of sigmoid gating would perform much better. * Please explain the architecture / experimental setting used in the GNN and RRN baselines in an appendix. Please also clarify how you use NRI (recurrent version, Markovian version etc.) * “Where am I” task: the model setup paragraph is unclear * Figure 3a: the example provided is not d x d (but 3x4 in terms of # puzzle pieces) * Ablation studies: please consider including experiments that analyze how the proposed model performs without global state, without the added prior (vector of ones), without the mixed sampling procedure, without recurrence and with different numbers of functions per edge. Other comments: * Please consider familiarizing yourself with the citep command in LaTeX for parenthetical citations to avoid sequences of closing parentheses.