Title
Arguments
- n
TODO
- k
TODO
- ...
Arguments passed on to
model_mediator_block2expected_degreeTODO
Examples
set.seed(26)
mrdpg <- model_null(n = 100, k = 5)
graph <- sample_tidygraph(mrdpg)
graph
#> # A tbl_graph: 100 nodes and 64772 edges
#> #
#> # An undirected multigraph with 1 component
#> #
#> # Node Data: 100 × 7 (active)
#> name intercept trt C1 C2 C3 y
#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 1 1 0 0 0 0 0.243
#> 2 2 1 0 0 0 0 0.735
#> 3 3 1 0 0 0 0 1.87
#> 4 4 1 0 0 0 0 -1.19
#> 5 5 1 0 0 0 0 -1.60
#> 6 6 1 0 0 0 0 0.361
#> 7 7 1 0 0 0 0 0.900
#> 8 8 1 0 0 0 0 1.85
#> 9 9 1 0 0 0 0 0.651
#> 10 10 1 0 0 0 0 0.601
#> # ℹ 90 more rows
#> #
#> # Edge Data: 64,772 × 2
#> from to
#> <int> <int>
#> 1 2 8
#> 2 11 13
#> 3 1 4
#> # ℹ 64,769 more rows
m_fit <- nodelm(US(A, 5) ~ . - name - y - 1, graph = graph)
o_fit <- nodelm(y ~ . - name - 1 + US(A, 5), graph = graph)
m_fit
#>
#> Call:
#> stats::lm(formula = formula, data = data)
#>
#> Coefficients:
#> 1 2 3 4 5
#> intercept 0.6159 -0.6197 -0.6140 1.1294 -6.5077
#> trt 6.8738 2.2028 1.1629 -1.4438 6.6679
#> C1 0.3961 -0.8264 -6.0460 -2.1858 6.9156
#> C2 0.1789 -0.2456 -0.3578 5.6089 7.6635
#> C3 1.0436 -6.5795 2.2301 -1.7367 6.7858
#>
o_fit
#>
#> Call:
#> stats::lm(formula = formula, data = data)
#>
#> Coefficients:
#> intercept trt C1 C2 C3 US(A, 5)1 US(A, 5)2
#> 1.478602 -0.833679 0.642584 -0.686271 -0.299833 0.061360 -0.009302
#> US(A, 5)3 US(A, 5)4 US(A, 5)5
#> 0.100876 0.106714 0.140755
#>