[Rsiena-help] Accessing the simulated networks

[Brian Rubineau]

Perfect! I'm off and running, thanks!
Brian

On 3/8/2011 5:35 PM, Ruth M. Ripley wrote:
> 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
>>
>>


-- 
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