Title

Title

Usage

model_uninformative(n, k, ...)

Arguments

n

TODO

k

TODO

...

Arguments passed on to model_mediator_uninformative

dim_w

TODO

Value

TODO

Examples

set.seed(26)

mrdpg <- model_uninformative(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   0.149
#>  2 2             1 1.89  0.261  -0.531 -0.468
#>  3 3             1 0.936 1.14    2.17   1.01 
#>  4 4             1 0.498 1.05    1.09   0.344
#>  5 5             1 0.546 2.41    0.971  0.873
#>  6 6             1 0.410 1.61    1.11   2.97 
#>  7 7             1 1.64  1.99   -0.485 -1.31 
#>  8 8             1 0.951 0.226   0.142 -0.274
#>  9 9             1 1.66  0.370   0.853  1.52 
#> 10 10            1 1.66  2.81   -1.41  -0.359
#> # ℹ 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.4048     0.2453     0.4385     0.5516     3.4429    -1.6113    -2.8710  
#> US(A, 5)4  US(A, 5)5  
#>   -2.1474     2.1055  
#> 

mrdpg$nde
#>       trt 
#> 0.2273358 
coef(o_fit)
#>  intercept        trt         C1         C2  US(A, 5)1  US(A, 5)2  US(A, 5)3 
#>  1.4048210  0.2452682  0.4385177  0.5516146  3.4429378 -1.6113163 -2.8710301 
#>  US(A, 5)4  US(A, 5)5 
#> -2.1473541  2.1054510 

coef(m_fit)
#>                     1            2            3           4            5
#> intercept 0.138749563 -0.256869263  0.148495852 0.058331167 -0.069015136
#> trt       0.009624984  0.030976510  0.012780530 0.037111693 -0.009562659
#> C1        0.010490521  0.009597515 -0.005427650 0.001196029 -0.007176028
#> C2        0.039830301  0.085808405 -0.003251126 0.014023819 -0.002797388
coef(o_fit)
#>  intercept        trt         C1         C2  US(A, 5)1  US(A, 5)2  US(A, 5)3 
#>  1.4048210  0.2452682  0.4385177  0.5516146  3.4429378 -1.6113163 -2.8710301 
#>  US(A, 5)4  US(A, 5)5 
#> -2.1473541  2.1054510