[Rsiena-help] Warnings "standard deviation is zero" in RSiena 1.1.-232 but not in RSiena 1.1-212

[Ruth Ripley]

My guess was based on an obvious difference between the versions. There are other differences which might have caused the problem. Don't know much about sienaRI but would be happy to investigate the fit if you send me the data privately. Am away from home and computer till 6th july, though.

Ruth

Sent from my iPad

> On 27 Jun 2014, at 16:55, Mark Manger <mark.manger at utoronto.ca> wrote:
> 
> Thanks for the quick reply. Tried that, but doesn’t make any difference — I still get 16 warnings and phase 3 doesn’t complete.
> 
> I wonder if I could produce the estimates with RSiena 1.1-212, and then update to 1.1-276 to use sienaRI()? Inconvenient, but better than nothing.
> 
> —Mark
> 
>> On Jun 26, 2014, at 5:26 PM, Ruth Ripley <ruth at stats.ox.ac.uk> wrote:
>> 
>> Dear Mark,
>> 
>> I am just guessing, but it could be that the message is coming from the
>> new Dolby option. Try setting Dolby off by dolby=FALSE in
>> sienaAlgorithmCreate (or sienaModelCreate).
>> 
>> Regards,
>> 
>> Ruth
>> 
>>> On 26/06/2014 22:01, Mark Manger wrote:
>>> Dear Tom, fellow RSiena users,
>>> 
>>> A small RSiena problem to resolve between World Cup matches:  I may have
>>> stumbled upon a bug/version conflict/problem of sorts with RSiena. I'm
>>> estimating model with two sets of one-mode networks.
>>> 
>>> The model can be estimated fine with R 2.15.3 and RSiena RSiena_1.1-212
>>> on my Mac running Mavericks, including using the parallel package. I've
>>> built RSiena_1.1-212 from source when downgrading to identify the
>>> potential problem, but it definitely worked back in April last year with
>>> the binary as well.
>>> 
>>> In the meantime I've upgraded to R 3.1.0 and RSiena 1.1-232. Alas, the
>>> identical model now estimates fine until Phase 3 Iteration 1, then stops
>>> with the following warnings (any number from 16 to 25 warnings are
>>> generated):
>>> 
>>> Warning messages:
>>> 1: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 2: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 3: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 4: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 5: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 6: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 7: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 8: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 9: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 10: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 11: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 12: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 13: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 14: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 15: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 16: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero
>>> 
>>> etc.
>>> 
>>> Updating RSiena to 1.1.276 changes nothing, i.e. I still get the
>>> warnings. I also get this when estimating the model without any
>>> covariates and three network effects.
>>> 
>>> First, I wonder if this should lead me to question the results obtained
>>> with Siena 1.1-212. If not, I'd be fine just using the old version. The
>>> only advantage of course would be that with 1.1-276 I could use sienaRI,
>>> which looks extremely promising. Second, what do these warning tell me?
>>> The .out file looks fine to me.
>>> 
>>> —Mark
>>> 
>>> 
>>> *Mark S. Manger* 
>>> Assistant Professor of Political Economy and Global Affairs
>>> 
>>> Munk School of Global Affairs & Department of Political Science |
>>> University of Toronto
>>> Observatory Site | 315 Bloor Street West | Room 212
>>> Toronto, ON   M5S 0A7
>>> Phone: 416-946-8927 | Fax: 416-946-8877
>>> 
>>> mark.manger at utoronto.ca <mailto:mark.manger at utoronto.ca>
>>> *www.munkschool.utoronto.ca/mga<http://www.munkschool.utoronto.ca/mga>*
>>> 
>>> JOIN THE GLOBAL CONVERSATION
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> _______________________________________________
>>> Rsiena-help mailing list
>>> Rsiena-help at lists.r-forge.r-project.org
>>> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help
>> _______________________________________________
>> Rsiena-help mailing list
>> Rsiena-help at lists.r-forge.r-project.org
>> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help
> 
> _______________________________________________
> Rsiena-help mailing list
> Rsiena-help at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/rsiena-help/attachments/20140627/b9129f34/attachment.html>