Skip to contents

Title

Usage

model_mediator_informative(
  n,
  k = 5,
  ztheta_0 = NULL,
  ztheta_t = NULL,
  ztheta_c = NULL,
  ztheta_tc = NULL,
  expected_degree = NULL
)

Arguments

n

TODO

k

TODO

ztheta_0

TODO

ztheta_t

TODO

ztheta_c

TODO

ztheta_tc

TODO

expected_degree

If specified, the desired expected degree of the graph. Specifying expected_degree simply rescales S to achieve this. Defaults to NULL. Do not specify both expected_degree and expected_density at the same time.

Value

TODO

Examples


inf <- model_mediator_informative(n = 100, k = 5)

graph <- sample_tidygraph(inf)
graph
#> # A tbl_graph: 100 nodes and 6772 edges
#> #
#> # An undirected multigraph with 1 component
#> #
#> # Node Data: 100 × 6 (active)
#>    name     C1    C2    C3    C4    C5
#>    <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
#>  1 1      2.95     0     0     0     0
#>  2 2      2.94     0     0     0     0
#>  3 3      2.66     0     0     0     0
#>  4 4      2.65     0     0     0     0
#>  5 5      2.48     0     0     0     0
#>  6 6      2.38     0     0     0     0
#>  7 7      2.34     0     0     0     0
#>  8 8      2.03     0     0     0     0
#>  9 9      1.97     0     0     0     0
#> 10 10     1.97     0     0     0     0
#> # ℹ 90 more rows
#> #
#> # Edge Data: 6,772 × 2
#>    from    to
#>   <int> <int>
#> 1    10    10
#> 2     3     5
#> 3    18    21
#> # ℹ 6,769 more rows

nodelm(US(A, 5) ~ . - name - 1, graph = graph)
#> 
#> Call:
#> stats::lm(formula = formula, data = data)
#> 
#> Coefficients:
#>     1          2          3          4          5        
#> C1  -1.357054   0.102438   0.025337   0.020521  -0.028681
#> C2  -0.425916  -0.278725  -0.537582   0.154448  -0.043429
#> C3  -0.064309  -0.309419  -0.473340   0.119659   0.608472
#> C4  -0.076867  -0.952825   0.567663   0.126887  -0.002168
#> C5  -0.060119  -0.293400  -0.188738  -1.237245  -0.032641
#> 

coef(inf)
#>            US1         US2         US3          US4         US5
#> C1 -1.36611137  0.09520364 -0.04821554  0.030836413 -0.01943658
#> C2 -0.45311261 -0.09629745  0.63843018  0.065601836 -0.02784641
#> C3 -0.05173586 -0.12960543  0.54207556 -0.002474788  0.62139094
#> C4 -0.08621943 -1.14471268 -0.21373428  0.058457948  0.01055509
#> C5 -0.05313442 -0.09609422  0.09965366 -1.253199397 -0.08662016

# should be very close to zero
inf$zX - inf$zC_true
#> 100 x 5 sparse Matrix of class "dgCMatrix"
#>              
#> 4   0 . . . .
#> 9   0 . . . .
#> 2   0 . . . .
#> 1   0 . . . .
#> 18  0 . . . .
#> 23  0 . . . .
#> 5   0 . . . .
#> 19  0 . . . .
#> 21  0 . . . .
#> 25  0 . . . .
#> 20  0 . . . .
#> 22  0 . . . .
#> 10  0 . . . .
#> 17  0 . . . .
#> 15  0 . . . .
#> 6   0 . . . .
#> 11  0 . . . .
#> 12  0 . . . .
#> 24  0 . . . .
#> 8   0 . . . .
#> 16  0 . . . .
#> 3   0 . . . .
#> 7   0 . . . .
#> 14  0 . . . .
#> 13  0 . . . .
#> 30  . 0 . . .
#> 27  . 0 . . .
#> 28  . 0 . . .
#> 42  . 0 . . .
#> 37  . 0 . . .
#> 32  . 0 . . .
#> 35  . 0 . . .
#> 33  . 0 . . .
#> 34  . 0 . . .
#> 26  . 0 . . .
#> 39  . 0 . . .
#> 41  . 0 . . .
#> 31  . 0 . . .
#> 38  . 0 . . .
#> 40  . 0 . . .
#> 29  . 0 . . .
#> 36  . 0 . . .
#> 59  . . 0 . .
#> 53  . . 0 . .
#> 60  . . 0 . .
#> 51  . . 0 . .
#> 48  . . 0 . .
#> 49  . . 0 . .
#> 43  . . 0 . .
#> 58  . . 0 . .
#> 54  . . 0 . .
#> 47  . . 0 . .
#> 50  . . 0 . .
#> 46  . . 0 . .
#> 44  . . 0 . .
#> 55  . . 0 . .
#> 57  . . 0 . .
#> 52  . . 0 . .
#> 56  . . 0 . .
#> 45  . . 0 . .
#> 81  . . . 0 .
#> 72  . . . 0 .
#> 80  . . . 0 .
#> 61  . . . 0 .
#> 71  . . . 0 .
#> 62  . . . 0 .
#> 76  . . . 0 .
#> 64  . . . 0 .
#> 73  . . . 0 .
#> 68  . . . 0 .
#> 63  . . . 0 .
#> 70  . . . 0 .
#> 79  . . . 0 .
#> 75  . . . 0 .
#> 65  . . . 0 .
#> 74  . . . 0 .
#> 66  . . . 0 .
#> 69  . . . 0 .
#> 67  . . . 0 .
#> 77  . . . 0 .
#> 78  . . . 0 .
#> 82  . . . 0 .
#> 96  . . . . 0
#> 86  . . . . 0
#> 100 . . . . 0
#> 95  . . . . 0
#> 83  . . . . 0
#> 99  . . . . 0
#> 98  . . . . 0
#> 87  . . . . 0
#> 93  . . . . 0
#> 85  . . . . 0
#> 91  . . . . 0
#> 97  . . . . 0
#> 89  . . . . 0
#> 92  . . . . 0
#> 84  . . . . 0
#> 90  . . . . 0
#> 88  . . . . 0
#> 94  . . . . 0

inf$zX - inf$zC_obs
#> 100 x 5 sparse Matrix of class "dgCMatrix"
#>                                                    
#> 4   0.00000000  .          .          .           .
#> 9   0.00000000  .          .          .           .
#> 2   0.00000000  .          .          .           .
#> 1   0.00000000  .          .          .           .
#> 18  0.00000000  .          .          .           .
#> 23  0.01206701  .          .          .           .
#> 5   0.03899010  .          .          .           .
#> 19  0.30915533  .          .          .           .
#> 21  0.30172164  .          .          .           .
#> 25  0.23911078  .          .          .           .
#> 20  0.23019508  .          .          .           .
#> 22  0.13220641  .          .          .           .
#> 10  0.21675651  .          .          .           .
#> 17  0.35017114  .          .          .           .
#> 15  0.41757253  .          .          .           .
#> 6   0.36425015  .          .          .           .
#> 11  0.64981798  .          .          .           .
#> 12  0.64903406  .          .          .           .
#> 24  0.59789837  .          .          .           .
#> 8   1.62283678 -2.9593289  .          .           .
#> 16  1.54970751 -2.9499810  .          .           .
#> 3   1.54209579 -2.8148440  .          .           .
#> 7   1.25622645 -2.6173262  .          .           .
#> 14  1.16616357 -2.3921666  .          .           .
#> 13  1.14226919 -2.2747296  .          .           .
#> 30  .           0.7529381  .          .           .
#> 27  .           0.8134400  .          .           .
#> 28  .           0.7765931  .          .           .
#> 42  .           0.7766368  .          .           .
#> 37  .           0.7846188  .          .           .
#> 32  .           0.8771586  .          .           .
#> 35  .           1.5126161  .          .           .
#> 33  .           2.2789057 -2.9128068  .           .
#> 34  .           1.9850142 -2.7309134  .           .
#> 26  .           1.9542766 -2.6482965  .           .
#> 39  .           1.9376598 -2.6382565  .           .
#> 41  .           1.9095708 -2.5389675  .           .
#> 31  .           1.8636777 -2.3209790  .           .
#> 38  .           1.6624749 -2.2789057  .           .
#> 40  .           1.5366793 -1.9956923  .           .
#> 29  .           1.0263513 -1.9850142  .           .
#> 36  .           1.0108907 -1.9376598  .           .
#> 59  .           .          0.9775275  .           .
#> 53  .           .          0.7051769  .           .
#> 60  .           .          0.4767255  .           .
#> 51  .           .          0.4220860  .           .
#> 48  .           .          0.1658179  .           .
#> 49  .           .          0.1315501  .           .
#> 43  .           .          0.1802310  .           .
#> 58  .           .          0.2364181  .           .
#> 54  .           .          0.2573972  .           .
#> 47  .           .          0.2670500  .           .
#> 50  .           .          0.2443710  .           .
#> 46  .           .          0.1073722  .           .
#> 44  .           .          0.1927989  .           .
#> 55  .           .          0.2925164  .           .
#> 57  .           .          0.1822232  .           .
#> 52  .           .          0.1433857  .           .
#> 56  .           .          0.1225647  .           .
#> 45  .           .          1.0292029 -2.976612888 .
#> 81  .           .          .          0.088833918 .
#> 72  .           .          .          0.215995147 .
#> 80  .           .          .          0.000904772 .
#> 61  .           .          .          0.064599419 .
#> 71  .           .          .          0.027025355 .
#> 62  .           .          .          0.192957972 .
#> 76  .           .          .          0.000000000 .
#> 64  .           .          .          0.000000000 .
#> 73  .           .          .          0.000000000 .
#> 68  .           .          .          0.000000000 .
#> 63  .           .          .          0.000000000 .
#> 70  .           .          .          0.000000000 .
#> 79  .           .          .          0.000000000 .
#> 75  .           .          .          0.000000000 .
#> 65  .           .          .          0.000000000 .
#> 74  .           .          .          0.000000000 .
#> 66  .           .          .          0.000000000 .
#> 69  .           .          .          0.000000000 .
#> 67  .           .          .          0.000000000 .
#> 77  .           .          .          0.000000000 .
#> 78  .           .          .          0.000000000 .
#> 82  .           .          .          0.000000000 .
#> 96  .           .          .          .           0
#> 86  .           .          .          .           0
#> 100 .           .          .          .           0
#> 95  .           .          .          .           0
#> 83  .           .          .          .           0
#> 99  .           .          .          .           0
#> 98  .           .          .          .           0
#> 87  .           .          .          .           0
#> 93  .           .          .          .           0
#> 85  .           .          .          .           0
#> 91  .           .          .          .           0
#> 97  .           .          .          .           0
#> 89  .           .          .          .           0
#> 92  .           .          .          .           0
#> 84  .           .          .          .           0
#> 90  .           .          .          .           0
#> 88  .           .          .          .           0
#> 94  .           .          .          .           0
inf$W - inf$zC_obs
#> 100 x 5 sparse Matrix of class "dgCMatrix"
#>     C1 C2 C3 C4 C5
#> 4    0  .  .  .  .
#> 9    0  .  .  .  .
#> 2    0  .  .  .  .
#> 1    0  .  .  .  .
#> 18   0  .  .  .  .
#> 5    0  .  .  .  .
#> 19   0  .  .  .  .
#> 10   0  .  .  .  .
#> 17   0  .  .  .  .
#> 15   0  .  .  .  .
#> 6    0  .  .  .  .
#> 11   0  .  .  .  .
#> 12   0  .  .  .  .
#> 8    0  .  .  .  .
#> 16   0  .  .  .  .
#> 3    0  .  .  .  .
#> 7    0  .  .  .  .
#> 14   0  .  .  .  .
#> 13   0  .  .  .  .
#> 30   .  0  .  .  .
#> 27   .  0  .  .  .
#> 28   .  0  .  .  .
#> 32   .  0  .  .  .
#> 23   .  0  .  .  .
#> 21   .  0  .  .  .
#> 25   .  0  .  .  .
#> 20   .  0  .  .  .
#> 22   .  0  .  .  .
#> 26   .  0  .  .  .
#> 31   .  0  .  .  .
#> 24   .  0  .  .  .
#> 29   .  0  .  .  .
#> 59   .  .  0  .  .
#> 42   .  .  0  .  .
#> 37   .  .  0  .  .
#> 53   .  .  0  .  .
#> 35   .  .  0  .  .
#> 51   .  .  0  .  .
#> 33   .  .  0  .  .
#> 48   .  .  0  .  .
#> 34   .  .  0  .  .
#> 39   .  .  0  .  .
#> 49   .  .  0  .  .
#> 43   .  .  0  .  .
#> 41   .  .  0  .  .
#> 58   .  .  0  .  .
#> 54   .  .  0  .  .
#> 47   .  .  0  .  .
#> 50   .  .  0  .  .
#> 38   .  .  0  .  .
#> 46   .  .  0  .  .
#> 40   .  .  0  .  .
#> 44   .  .  0  .  .
#> 55   .  .  0  .  .
#> 57   .  .  0  .  .
#> 52   .  .  0  .  .
#> 56   .  .  0  .  .
#> 45   .  .  0  .  .
#> 36   .  .  0  .  .
#> 81   .  .  .  0  .
#> 72   .  .  .  0  .
#> 80   .  .  .  0  .
#> 61   .  .  .  0  .
#> 71   .  .  .  0  .
#> 62   .  .  .  0  .
#> 60   .  .  .  0  .
#> 76   .  .  .  0  .
#> 64   .  .  .  0  .
#> 73   .  .  .  0  .
#> 68   .  .  .  0  .
#> 63   .  .  .  0  .
#> 70   .  .  .  0  .
#> 79   .  .  .  0  .
#> 75   .  .  .  0  .
#> 65   .  .  .  0  .
#> 74   .  .  .  0  .
#> 66   .  .  .  0  .
#> 69   .  .  .  0  .
#> 67   .  .  .  0  .
#> 77   .  .  .  0  .
#> 78   .  .  .  0  .
#> 82   .  .  .  0  .
#> 96   .  .  .  .  0
#> 86   .  .  .  .  0
#> 100  .  .  .  .  0
#> 95   .  .  .  .  0
#> 83   .  .  .  .  0
#> 99   .  .  .  .  0
#> 98   .  .  .  .  0
#> 87   .  .  .  .  0
#> 93   .  .  .  .  0
#> 85   .  .  .  .  0
#> 91   .  .  .  .  0
#> 97   .  .  .  .  0
#> 89   .  .  .  .  0
#> 92   .  .  .  .  0
#> 84   .  .  .  .  0
#> 90   .  .  .  .  0
#> 88   .  .  .  .  0
#> 94   .  .  .  .  0