demo/logit.net.R

In IQSS/ZeligNetwork: Social Network Regressions for Zelig

library(ZeligNetwork)

data(friendship)

# fit model
z.out <- zelig(
               friends ~ advice + prestige + perpower, 
               model="logit.net", 
               data=friendship
               )
summary(z.out)

user.prompt()

# set explanatory variables
x.high <- setx(z.out, perpower = quantile(friendship$perpower, prob=0.75))
x.low <- setx(z.out, perpower = quantile(friendship$perpower, prob=0.25))

user.prompt()

# simulate quantities of interest
s.out <- sim(z.out, x = x.high)
summary(s.out)

user.prompt()

# plot quantities of interest
plot(s.out)