[Rsiena-help] Accessing the simulated networks

[Ruth M. Ripley]

Dear Brian,

You need to ask for the simulations to be returned. Use returnDeps=TRUE in 
the call to siena07.

By the way, you can turn off estimation by using nsub=0 in the call to 
sienaModelCreate. If you do this, you do not need to fix the parameters.

Regards,

Ruth


On Tue, 8 Mar 2011, Brian Rubineau wrote:

> Dear RSiena experts,
> I'm trying to use the simulation engine in RSiena for fixed network models.
> That is, the output I'm interested in is the networks in the $sims portion of 
> an
> RSiena estimation object, not the estimation results itself. My problem is
> that I'm only getting "NULL" in the $sims list.
> I'm using R 2.12.0 and RSiena 1.0.11.119.
>
> Here's a simple example of what I'm trying to do, and the problem:
>
> library(statnet)  ## to get the random graph function rgraph
> two <- rgraph(20,m=2,tprob=0.2)
> two[1,,] <- 0  ## set the 1st graph to empty
> two[1,sample(1:20,1),sample(1:20,1)]<-1  ## make the 1st graph have only 1 
> random tie
> two.siena <- sienaNet(array(c(two[1,,],two[2,,]),dim=c(20,20,2)))
> two.siena.data <- sienaDataCreate(two.siena)
> two.siena.model <- sienaModelCreate(useStdInits = FALSE, projname = 
> 'two',n3=10)
> two.effects <- getEffects(two.siena.data)
> two.effects$fix <- two.effects$include ## to fix all the parameters
> two.1 <- siena07(two.siena.model,data=two.siena.data,effects=two.effects)
> two.1$sims
> [[1]]
> NULL
>
> [[2]]
> NULL
>
> [[3]]
> NULL
>
> [[4]]
> NULL
>
> [[5]]
> NULL
>
> The errors from the siena07 process are:
>> Error in solve.default(z$msfc) :
>>   Lapack routine dgesv: system is exactly singular
>> In addition: Warning message:
>> In max(abs(tstat)[!z$fixed & !z$BasicRateFunction & z$resist > 0.9]) :
>>   no non-missing arguments to max; returning -Inf
>
> I'm attaching the out file from this process.
>
> Any ideas as to what I'm doing wrong?
> Thanks so much!
> Best,
> Brian
>
>
>
> -- 
> Brian Rubineau
> Assistant Professor, Organizational Behavior        394 Ives Hall
> School of Industrial and Labor Relations         ph. 607.255.3048
> Cornell University                              fax. 607.255.2261
> Ithaca, NY 14853                            brubineau at cornell.edu
>
>