Title

Title

Usage

model_null2(n, k, ...)

Arguments

n

TODO

k

TODO

...

Arguments passed on to model_mediator_block2

expected_degree

TODO

Value

TODO

Examples

set.seed(26)

mrdpg <- model_null2(n = 100, k = 5)

graph <- sample_tidygraph(mrdpg)
graph
#> # A tbl_graph: 100 nodes and 521 edges
#> #
#> # An undirected multigraph with 1 component
#> #
#> # Node Data: 100 × 6 (active)
#>    name  intercept   trt     C1     C2     y
#>    <chr>     <dbl> <dbl>  <dbl>  <dbl> <dbl>
#>  1 1             1 0.560 0.0687  2.69  2.84 
#>  2 2             1 1.89  0.261  -0.531 1.80 
#>  3 3             1 0.936 1.14    2.17  3.46 
#>  4 4             1 0.498 1.05    1.09  2.88 
#>  5 5             1 0.546 2.41    0.971 3.28 
#>  6 6             1 0.410 1.61    1.11  5.27 
#>  7 7             1 1.64  1.99   -0.485 0.486
#>  8 8             1 0.951 0.226   0.142 1.66 
#>  9 9             1 1.66  0.370   0.853 3.13 
#> 10 10            1 1.66  2.81   -1.41  0.985
#> # ℹ 90 more rows
#> #
#> # Edge Data: 521 × 2
#>    from    to
#>   <int> <int>
#> 1     4     6
#> 2     3    14
#> 3     4    12
#> # ℹ 518 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.138750  -0.256869   0.148496   0.058331  -0.069015
#> trt         0.009625   0.030977   0.012781   0.037112  -0.009563
#> C1          0.010491   0.009598  -0.005428   0.001196  -0.007176
#> C2          0.039830   0.085808  -0.003251   0.014024  -0.002797
#> 
o_fit
#> 
#> Call:
#> stats::lm(formula = formula, data = data)
#> 
#> Coefficients:
#> intercept        trt         C1         C2  US(A, 5)1  US(A, 5)2  US(A, 5)3  
#>   1.67245    0.03119    0.41527    0.53481    0.05254    1.00867   -0.11804  
#> US(A, 5)4  US(A, 5)5  
#>  -0.13847   -0.03136  
#>