<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Towards E ective and E cient Approximate Query Answering in Probabilistic DeLP?</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Mario A. Leiva</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Alejandro J. Garc a</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Gerardo I. Simari</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Depto. de Cs. e Ing. de la Comp., Universidad Nacional del Sur (UNS) Inst. de Cs. e Ing. de la Comp. (ICIC UNS-CONICET)</institution>
          ,
          <country country="AR">Argentina</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Rules</institution>
          ,
          <addr-line>Facts</addr-line>
        </aff>
      </contrib-group>
      <abstract>
        <p>This work presents an overview of a research project focused on R&amp;D for approximate query answering in probabilistic structured argumentation based on DeLP [2]. The ultimate goal is to develop a suite of algorithms for tackling the computational cost of this task and selection criteria for choosing the best one based on the analysis of available information. In the rest of this section, we brie y present the basics of DeLP3E knowledge bases.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>Analytical Model (AM)</p>
      <p>and
Presumptions
 ∪  ∪  ∪</p>
      <p>Arguments
Dialectical Process</p>
      <p>Annotation Function(AF)
 : Ω ∪ Θ ∪ Δ ∪ Φ →</p>
      <p>Environmental Model</p>
      <p>(EM)
Worlds       ...     (  )
 1
 2
 3
…</p>
      <p>T T T … T</p>
      <p>0.47
T T T … F 0.15
T T T … T 0.05
… … … … …</p>
      <p>…
 1024 F F F … F 0.003</p>
      <p>Probabilistic Model
is denoted with WEM . In the same way, the set of all subsets of elements of
PAM (subprograms) is called</p>
      <p>WAM . Intuitively, given PAM , every element of
(rules, facts or presumptions) only holds in certain worlds in the
, i.e., these elements are subject to probabilistic events. Each element
is thus associated with a formula over the set of all ground
elements of PEM (using conjunction, disjunction, and negation, as usual). The
notion of annotation function implements this association.</p>
      <p>In order to answer a query for a literal of interest, we need to compute the
probability interval with which it is warranted in the DeLP3E KB|for this, we
must sum the probability mass of the worlds that generate subprograms where
the queried literal is warranted (warranting scenarios, for the lower bound) and
the mass worlds whose generated subprogram warrants the complement of the
query (for the upper bound). The lower and upper bounds obtained in this
manner comprise the probability interval for the query. This is one of the main
sources of computational intractability, since we must answer the query either
present two approaches that we analyze to tackle this intractability.
for all the worlds in WEM or all the programs in WAM . In Section 2, we will</p>
      <p>Next, we present a simple example to illustrate the DeLP3E KB components
and show how the probability interval is calculated for a given query.
A Simple Example. Consider a small KB, with components PEM , PAM , and
af as shown in Figure 2.</p>
      <p>For the environmental model PEM , we use a very simple Bayesian network
that describes relationship between the nodes a, b and c (EM elements); we omit
a
b</p>
      <p>
        c
: ;
: 1 = l1
: 1 = l5 |&lt; l3
2 = l2 |&lt; l4
3 = l2 |&lt; l3
: 1 = l3 |&lt;
2 = l4 |&lt;
af(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) = true
af(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) = b _ c
af( 2) = true
af( 3) = true
af(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) = a
af( 2) = true
the detailed speci cation of the CPTs and instead provide directly the complete
probability distribution Pr over all possible worlds WEM . Thus, for example,
the probability associated with the world 5 where a is false and b and c are
true, is 0.30.
      </p>
      <p>The domain knowledge in the analytical model PAM is represented by the
DeLP program P = [ [ [ (see Figure 2 bottom left). In that program,
there is an empty set of strict rules = ;, only one fact = f 1g, two
presumptions = f 1; 2g and three defeasible rules = f 1; 2; 3g. The annotation
function af (Figure 2 bottom right) associates elements of P with probabilistic
events through formulas. As shown, in this example, almost every element hold
in every world (annotated as true), whereas for rule 1 the annotation means
that this presumption only holds in worlds which the formula a (probabilistic
event a) is satis ed, and the annotation for rule 1 means that this rule only
holds in worlds in which the formula b _ c, that relates the probabilistic events
b and c, is satis ed.</p>
      <p>To answer a query, consider for example the literal l2. To compute the
associated probability interval, we need to obtain the warranting scenarios for
both l2 (for the lower bound) and its complement l2 (for the upper bound).
It is easy to check that the condition for the former is that a does not hold,
which yields the set of worlds f 5; 6; 7; 8g; it is also easy to check that there
are no warranting scenarios for L2. Therefore, the lower bound is calculated by
summing P8</p>
      <p>i=5 Pr( i) = 0:7, while the upper bound is simply 1, so the answer
to the query is Pr( l2) 2 [0:70; 1].
 7</p>
      <p>5
 14 …  
 9
A well-known technique to address intractability in this kind of situation is
to sample a subset of the solution space in order to arrive at an approximate
answer. In this case, we can consider two types of sampling: world-based and
subprogram-based, as illustrated in Figure 3. In world-based sampling, a subset
of the possible EM worlds are chosen and, based on the annotation function,
di erent subprograms of PAM are obtained (as shown in Figure 1). So, each
PAM ( ) is a classical DeLP program [2] in which we can query for the status of
some literal. On the other hand, subprogram-based sampling divides the universe
of all possible subprograms into regions that represent programs that warrant the
query or its complement. Each of these programs can be generated by multiple
worlds through their annotations.</p>
      <p>In both cases, we also have those worlds or programs in which the status
of the query can be \undecided" or \unknown". The objective is to sum the
probability value corresponding to the green and red worlds (circles in the gure);
that is, those that are warranting scenarios for the query and its complement,
respectively. Given the incomplete nature of the process, some probability mass
will remain unexplored (the white elements in Figure 3).</p>
      <p>A Family of KB Metrics. In order to de ne a suite of sampling algorithms
to approximate the answer to a probabilistic query in the most e cient and
e ective way possible, we must consider what metrics we have available for the
input KB. Figure 4 shows a set of such possible metrics, organized with respect to
the component they apply to, and what attribute they are designed to measure.
For each one, we also analyze whether or not they can be tractably computed
(observable column) or approximated.</p>
      <p>For each component (or submodel of the KB), we can focus on two main
attributes: size and complexity. For PAM , #Rules and #Facts refers to the
total number of di erent rules and facts in the DeLP program, respectively.
Regarding the complexity component, the MDDL metric refers to the maximum
defeasible derivation length, i.e., the maximum number of defeasible rules present
in any argument constructed from the program. In a similar vein, h refers to
Component</p>
      <p>PAM
PEM
AF</p>
      <p>Attribute Example Metric</p>
      <p>Size ##RFaucletss</p>
      <p>MDDL</p>
      <sec id="sec-1-1">
        <title>Complexity h</title>
      </sec>
      <sec id="sec-1-2">
        <title>Size</title>
      </sec>
      <sec id="sec-1-3">
        <title>Complexity</title>
      </sec>
      <sec id="sec-1-4">
        <title>Size</title>
      </sec>
      <sec id="sec-1-5">
        <title>Complexity</title>
        <p>#RandomVars
#PGM Arcs
PGM TW
Entropy
%AF ann
AF Comp</p>
        <sec id="sec-1-5-1">
          <title>Observable?</title>
          <p>Yes
Yes
{
{
{
Yes
Yes
{
{
Yes
Yes</p>
        </sec>
        <sec id="sec-1-5-2">
          <title>Approx.?</title>
          <p>{
{
Yes
Yes
Yes
{
{
Yes</p>
        </sec>
      </sec>
      <sec id="sec-1-6">
        <title>Sometimes { { Fig. 4. Example metrics computable (exactly or approximately) for a DeLP3E KB.</title>
        <p>the maximum length of an argumentation line (a sequence of arguments where
each element of the sequence defeats its predecessor), and is the number of
dialectical trees (a tree structure formed by all argumentation lines arising from
one argument) for each literal in the program [2].</p>
        <p>On the other hand, for PEM , #RandomVars refers to the number of random
variables in the model (for instance, nodes in a BN), while #PGM Arcs refers
to the number of arcs in a probabilistic graphical model (PGM) [4]. The metric
PGM TW refers to the treewidth of a PGM (such as a BN), which intuitively
measures how \close to a tree" a graph is; approximations for this metric can be
computed [3], but the exact value cannot be tractably calculated for very large
models. The entropy metric is inherent to the joint probability distribution
function represented by PEM , and it is possible to approximate it in some particular
scenarios [1].</p>
        <p>Finally, for annotation functions (AF) an example metric for the size
attribute is %AF ann, which refers to the percentage of formulas of the program
that are annotated. For the complexity attribute, metric AF Comp refers to the
complexity of each annotation itself (that is, which operators can be used, if
nesting is allowed, etc.).</p>
        <p>Towards a Suite of Approximation Algorithms. Based on the information
obtained using these metrics, we can begin to explore the possible alternatives
for deriving sampling-based approximation algorithms. For example, consider a
scenario in which we have 10 elements in the PAM (#Rules + #Facts) and 50 in
the PEM (#RandomVars), and AF Comp is low (only one variable from PEM
is used in the annotations, without operators). Here, we can approximate the
probability of a query by sampling subprograms, since the number of
subprograms will be smaller in comparison with the number of worlds (210 &lt; 250),
and the annotations associated with subprograms are simple to evaluate
(conjunctions of variables or their negation). On the other hand, if the annotations
are complex (they use several variables, connectors and negations, for example
af (!1) = a _ b ! :c) then sampling by worlds is simpler (even if there are more
elements in PEM than rules in PAM ). This is because the task of obtaining the
subprogram is simple (substitute the random variable values according to the
sampled world and evaluating the annotation formula, then query for the literal
in the generated program) compared to computing the probability of a complex
expression in the EM.</p>
        <p>Note that the examples mentioned above expose a kind of asymmetry
between the PAM and the PEM |even though given a world it is only possible
to generate a single program, a program can be generated by a set of worlds
(since the formulas in the annotations can have more than one model). This can
be observed also in Figure 3. Consider a case in which the PEM is a BN with
high values of treewidth and entropy|for example, 5 (complex structure) and
14, respectively|meaning that the network is far from a tree structure and that
the underlying probability distribution has a high degree of uncertainty. This
leads to a slower, more complex, and less guided sampling process; therefore, an
algorithm to sample worlds randomly is not recommended, and in this case it
is better to decide in a previous step which worlds to sample (such as weighted
sampling given the BN), in order to optimize resources.</p>
        <p>These examples show the variety of possible approaches to approximate query
answering in probabilistic DeLP. Our ultimate goal is to develop a set of decision
criteria that allow to select the best algorithm for the job, contemplating the
inherent tradeo s between running time (including the cost of calculating any
necessary approximate metrics) and precision of the result obtained.
3</p>
      </sec>
    </sec>
    <sec id="sec-2">
      <title>Conclusions and Future Work</title>
      <p>In this work we presented a preliminary study of the aspects that must be
considered when developing a family of approximation algorithms for query answering
in the DeLP3E framework. Ongoing and future work consists of developing
approximation algorithms that can handle large instances based on the criteria
initially explored here, and implement them in the DAQAP platform [5]. In
order to do this e ectively, we are developing procedures for automatically
generating DeLP3E KBs in order to simulate di erent scenarios, running a suite of
approximation algorithms over them, and evaluating their performance. Other
approaches currently being evaluated are the application of machine learning
techniques to sample in a guided way during the approximation process.
Acknowledgments. This work was supported by funds provided by
Universidad Nacional del Sur (UNS) under grants PGI 24/N046 and PGI 24/ZN34),
Agencia Nacional de Promocion Cient ca y Tecnologica under grant
PICT-20180475 (PRH-2014-0007), and CONICET under grant PIP 11220170100871CO.</p>
      <sec id="sec-2-1">
        <title>2. Garc a, A.J., Simari, G.R.: Defeasible logic programming: An argumentative ap</title>
        <p>
          proach. TPLP 4(
          <xref ref-type="bibr" rid="ref1">1-2</xref>
          ), 95{138 (2004)
        </p>
      </sec>
      <sec id="sec-2-2">
        <title>3. Kloks, T.: Treewidth: computations and approximations, vol. 842. Springer Science</title>
        <p>&amp; Business Media (1994)</p>
      </sec>
      <sec id="sec-2-3">
        <title>4. Koller, D., Friedman, N.: Probabilistic graphical models: principles and techniques.</title>
        <p>MIT press (2009)
5. Leiva, M.A., Simari, G.I., Gottifredi, S., Garc a, A.J., Simari, G.R.: DAQAP:
defeasible argumentation query answering platform. In: Proc. FQAS 2019. LNCS, vol.
11529, pp. 126{138. Springer (2019)</p>
      </sec>
      <sec id="sec-2-4">
        <title>6. Nilsson, N.J.: Probabilistic logic. Arti cial intelligence 28(1), 71{87 (1986)</title>
      </sec>
      <sec id="sec-2-5">
        <title>7. Pearl, J.: Probabilistic reasoning in intelligent systems: networks of plausible infer</title>
        <p>ence. Elsevier (2014)</p>
      </sec>
      <sec id="sec-2-6">
        <title>8. Richardson, M., Domingos, P.: Markov logic networks. Machine learning 62(1-2),</title>
        <p>107{136 (2006)</p>
      </sec>
      <sec id="sec-2-7">
        <title>9. Shakarian, P., Simari, G.I., Moores, G., Paulo, D., Parsons, S., Falappa, M.A., Aleali,</title>
      </sec>
      <sec id="sec-2-8">
        <title>A.: Belief revision in structured probabilistic argumentation. Annals of Mathematics</title>
        <p>and Arti cial Intelligence 78(3-4), 259{301 (2016)</p>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Batu</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dasgupta</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kumar</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rubinfeld</surname>
            ,
            <given-names>R.:</given-names>
          </string-name>
          <article-title>The complexity of approximating the entropy</article-title>
          .
          <source>SIAM J. Comput</source>
          .
          <volume>35</volume>
          (
          <issue>1</issue>
          ),
          <volume>132</volume>
          {
          <fpage>150</fpage>
          (
          <year>2005</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>