statnet / ergm.multi

Fit, Simulate and Diagnose Exponential-Family Models for Multiple or Multilayer Networks
Other
14 stars 1 forks source link

In N() operator, provide a way to specify different network-level predictors for different terms. #3

Closed krivit closed 2 years ago

krivit commented 3 years ago

Motivating example: we want GWESP with decay common to all networks in the sample but coefficient varying, say, as a function of net size. Specifying N(~gwesp(), ~log(n)) will have net size effect for both sample and decay. How do we want to specify the UI?

krivit commented 3 years ago

UI already allows this to be done: offset(N(~gwesp(), ~log(n)), 4) with offset.coef=0 will "lock" the slope for the decay parameter at 0. TODO: Document this somewhere.

cmarcum commented 3 years ago

would the interpretation of the statistic effectively be a grand mean over all networks and all layers? I can add the documentation: I've been meaning to do that to clean up some typos anyway

On Tue, Jul 20, 2021, 8:22 AM Pavel N. Krivitsky @.***> wrote:

UI already allows this to be done: offset(N(~gwesp(), ~log(n)), 4) with offset.coef=0 will "lock" the slope for the decay parameter at 0. TODO: Document this somewhere.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <statnet/issue-mover#27>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAI377F7ZL32AZ5K7SZ25STTYVS7NANCNFSM5AQVXTFA .

krivit commented 3 years ago

I guess it could be viewed as a weighted mean, but generally it's a pooled MLE (e.g., https://pubmed.ncbi.nlm.nih.gov/32547745/). That said, I would suggest holding off on updating the documentation, since it's going to be ported to the new documentation system that should let you do it in Rmarkdown.

krivit commented 3 years ago

And thanks for offering!

krivit commented 2 years ago

@cmarcum , the documentation has been ported, in case you've been waiting.