[Rsiena-help] Multiple networks

[Ruth M. Ripley]

Dear Paola,

You should be able to use something like
mynet1 <- sienaNet(array(c(s501, s502), dim=c(50, 50, 2)))
mynet2 <- sienaNet(array(c(s502, s503), dim=c(50, 50, 2)))
mydata <- sienaDataCreate(mynet1, mynet2)
myeff <- getEffects(mydata)
## for the first named network
myeff <- includeEffects(myeff, transTrip, balance)
## for other networks
myeff <- includeEffects(myeff, transTrip, balance, name='mynet2')

Effects that include two networks will need interaction1 set equal to the 
second network of the pair, with name specified if necessary (the default 
is the first network that appears in the effects object).

Please send me the code if you still cannot get it to work.

Regards,

Ruth


On Thu, 22 Jul 2010, paola zappa wrote:

> Dear Ruth,
> I'm trying to model the evolution of multiple one-mode networks (so far just
> two: RTRw and RTRs) and I'm not sure about how to specify for which network
> I want to estimate the effects I selected.
> I read in RSiena manual (page 84) it should be done by including 'name= ...'
> in the includeEffects function, but I'm keeping on receiving error messages!
>
>
> Thank you in advance.
> Regards,
> Paola
>