htna-named wrapper for Nestimate::certainty(). The closed-form
counterpart of bootstrap_htna(): rather than resampling, it derives
analytic confidence intervals and an inference verdict for each edge
weight, which is fast enough to run where a full bootstrap would be
costly.
Arguments
- x
An htna network from
build_htna()(or other object accepted byNestimate::certainty()).- ...
Additional arguments passed to
Nestimate::certainty(), such asprior,ci_level,inference,consistency_range, andedge_threshold. See that function for the full, current argument list.
Value
An object of class net_certainty (also inheriting
net_bootstrap). See Nestimate::certainty() for the full
component list and the corresponding plot() method.
Details
Works on htna networks and grouped htna networks directly.
Suffixed _htna to sit alongside bootstrap_htna() and
reliability_htna() in the htna confirmatory-analysis family.
See also
bootstrap_htna() for the resampling analogue,
reliability_htna(), casedrop_reliability_htna().
Examples
# \donttest{
data(human_ai)
net <- build_htna(human_ai, actor_type = "actor_type")
#> Warning: A network with one long sequence is not recommended and can't be validated using bootstrap and other confirmatory testings.
#> Metadata aggregated per session: ties resolved by first occurrence in 'session_date' (1 sessions), 'cluster' (42 sessions), 'actor_type' (24 sessions)
certainty_htna(net)
#> Edge Mean 95% CI p
#> -----------------------------------------------
#> Delegate → Plan 0.600 [0.543, 0.656] ***
#> Refine → Specify 0.410 [0.377, 0.445] ***
#> Check → Execute 0.410 [0.383, 0.437] ***
#> Ask → Plan 0.408 [0.388, 0.428] ***
#> Repair → Execute 0.384 [0.326, 0.444] **
#> ... and 63 more certain edges
#>
#> Certainty (Dirichlet) [Transition Network (relative) | directed]
#> Prior : Dirichlet(0.50) | Nodes : 12
#> Edges : 68 certain / 135 total
#> CI : 95% | Inference: stability | CR [0.75, 1.25]
# }
