Title
Arguments
- n
TODO
- k
TODO
- ...
Arguments passed on to
model_mediator_block2expected_degreeTODO
Examples
set.seed(26)
mrdpg <- model_block(n = 100, k = 5)
graph <- sample_tidygraph(mrdpg)
graph
#> # A tbl_graph: 100 nodes and 64849 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 -18.7
#> 2 2 1 0 0 0 0 -16.8
#> 3 3 1 0 0 0 0 -16.8
#> 4 4 1 0 0 0 0 -17.4
#> 5 5 1 0 0 0 0 -16.0
#> 6 6 1 0 0 0 0 -12.8
#> 7 7 1 0 0 0 0 -15.5
#> 8 8 1 0 0 0 0 -12.3
#> 9 9 1 0 0 0 0 -11.1
#> 10 10 1 0 0 0 0 -8.66
#> # ℹ 90 more rows
#> #
#> # Edge Data: 64,849 × 2
#> from to
#> <int> <int>
#> 1 3 15
#> 2 2 2
#> 3 5 13
#> # ℹ 64,846 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.6634 -0.6109 -0.7455 1.1750 -6.5163
#> trt 6.8154 2.2159 1.2954 -1.4991 6.6958
#> C1 0.3209 -0.7589 -5.9138 -2.2826 7.0145
#> C2 0.1433 -0.2490 -0.2667 5.5889 7.7282
#> C3 1.0306 -6.5937 2.3031 -1.7866 6.7774
#>
o_fit
#>
#> Call:
#> stats::lm(formula = formula, data = data)
#>
#> Coefficients:
#> intercept trt C1 C2 C3 US(A, 5)1 US(A, 5)2
#> 2.5317 -1.2148 0.0885 -1.1350 -0.2397 0.8030 1.5277
#> US(A, 5)3 US(A, 5)4 US(A, 5)5
#> 1.1536 -1.3185 1.5893
#>