<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Dear Mark,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The error message is produced by a comparison of a variance with 0, and this of course may differ between systems depending on machine accuracy. I attach a
file in which this comparison is replaced by a comparison with e-12 (scientific notation for 10^{-12}).<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">You can use this by putting this file in your working directory and then, after<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">library(RSienaTest)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">giving the command<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">insertSource("phase3.r", package='RSienaTest')<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Does the error then still occur?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Best wishes,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Tom<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">================================================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Tom A.B. Snijders<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Professor of Statistics in the Social Sciences<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Department of Politics and Department of Statistics<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Nuffield College<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">University of Oxford<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">tel. +44-01865-278599<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> rsiena-help-bounces@lists.r-forge.r-project.org [mailto:rsiena-help-bounces@lists.r-forge.r-project.org]
<b>On Behalf Of </b>Mark Manger<br>
<b>Sent:</b> 27 June 2014 16:55<br>
<b>To:</b> rsiena-help@lists.r-forge.r-project.org<br>
<b>Subject:</b> Re: [Rsiena-help] Warnings "standard deviation is zero" in RSiena 1.1.-232 but not in RSiena 1.1-212<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Thanks for the quick reply. Tried that, but doesn’t make any difference — I still get 16 warnings and phase 3 doesn’t complete.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">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.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">—Mark<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">On Jun 26, 2014, at 5:26 PM, Ruth Ripley <<a href="mailto:ruth@stats.ox.ac.uk">ruth@stats.ox.ac.uk</a>> wrote:<o:p></o:p></p>
</div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<div>
<p class="MsoNormal"><span style="font-size:13.5pt">Dear Mark,<br>
<br>
I am just guessing, but it could be that the message is coming from the<br>
new Dolby option. Try setting Dolby off by dolby=FALSE in<br>
sienaAlgorithmCreate (or sienaModelCreate).<br>
<br>
Regards,<br>
<br>
Ruth<br>
<br>
On 26/06/2014 22:01, Mark Manger wrote:<br>
<br>
<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:13.5pt"><span style="font-size:13.5pt">Dear Tom, fellow RSiena users,<br>
<br>
A small RSiena problem to resolve between World Cup matches: I may have<br>
stumbled upon a bug/version conflict/problem of sorts with RSiena. I'm<br>
estimating model with two sets of one-mode networks.<br>
<br>
The model can be estimated fine with R 2.15.3 and RSiena RSiena_1.1-212<br>
on my Mac running Mavericks, including using the parallel package. I've<br>
built RSiena_1.1-212 from source when downgrading to identify the<br>
potential problem, but it definitely worked back in April last year with<br>
the binary as well.<br>
<br>
In the meantime I've upgraded to R 3.1.0 and RSiena 1.1-232. Alas, the<br>
identical model now estimates fine until Phase 3 Iteration 1, then stops<br>
with the following warnings (any number from 16 to 25 warnings are<br>
generated):<br>
<br>
Warning messages:<br>
1: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
2: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
3: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
4: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
5: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
6: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
7: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
8: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
9: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
10: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
11: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
12: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
13: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
14: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
15: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
16: In cor(z$sf[, i], scores[, i]) : the standard deviation is zero<br>
<br>
etc.<br>
<br>
Updating RSiena to 1.1.276 changes nothing, i.e. I still get the<br>
warnings. I also get this when estimating the model without any<br>
covariates and three network effects.<br>
<br>
First, I wonder if this should lead me to question the results obtained<br>
with Siena 1.1-212. If not, I'd be fine just using the old version. The<br>
only advantage of course would be that with 1.1-276 I could use sienaRI,<br>
which looks extremely promising. Second, what do these warning tell me?<br>
The .out file looks fine to me.<br>
<br>
—Mark<br>
<br>
<br>
*Mark S. Manger*<span class="apple-converted-space"> </span><br>
Assistant Professor of Political Economy and Global Affairs<br>
<br>
Munk School of Global Affairs & Department of Political Science |<br>
University of Toronto<br>
Observatory Site | 315 Bloor Street West | Room 212<br>
Toronto, ON M5S 0A7<br>
Phone: 416-946-8927 | Fax: 416-946-8877<br>
<br>
<a href="mailto:mark.manger@utoronto.ca">mark.manger@utoronto.ca</a><span class="apple-converted-space"> </span><<a href="mailto:mark.manger@utoronto.ca">mailto:mark.manger@utoronto.ca</a>><br>
*<a href="http://www.munkschool.utoronto.ca/mga">www.munkschool.utoronto.ca/mga</a><<a href="http://www.munkschool.utoronto.ca/mga">http://www.munkschool.utoronto.ca/mga</a>>*<br>
<br>
JOIN THE GLOBAL CONVERSATION<br>
<br>
<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Rsiena-help mailing list<br>
<a href="mailto:Rsiena-help@lists.r-forge.r-project.org">Rsiena-help@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help</a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:13.5pt">_______________________________________________<br>
Rsiena-help mailing list<br>
<a href="mailto:Rsiena-help@lists.r-forge.r-project.org">Rsiena-help@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rsiena-help</a><o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>