<div>My C code:</div>
<div>//[[Rcpp::depends("Rcpp")]]<br>#include <Rcpp.h><br>#include <iostream><br>using namespace Rcpp;</div>
<div> </div>
<div>//[[Rcpp::export]]<br>NumericVector sigmutest(double pd, double rsq){<br>  double qpd = R::qnorm(pd, 0, 1, 1, 0);<br>  double sgtemp = 0.2;<br>  double sg = 0.3;<br>  double eor = 1;<br>  <br>  double w = 0;<br>  while (eor>=0.0001) {  <br>
        sg = sgtemp;<br>        w = pd * 0.4 / (R::pnorm((qpd - sqrt(rsq) * (-0.42) * sgtemp), 0, 1, 1, 0));<br>        sgtemp = (-0.5) * w + 0.4;<br>        std::cout << "sg " << sg << std::endl;<br>
        std::cout << "sgtemp " << sgtemp << std::endl;<br>        eor = abs(sg - sgtemp);<br>        <br>        std::cout << "error " << eor << std::endl;<br>        <br>
      } <br>  NumericVector out(3);<br>  out(0) = sg;<br>  out(1) = sgtemp;<br>  out(2) = eor;<br>  return out;  <br>}<br clear="all"><br>My R code:</div>
<div><span style="TEXT-TRANSFORM:none;BACKGROUND-COLOR:rgb(225,226,229);TEXT-INDENT:0px;BORDER-COLLAPSE:separate;FONT:13px/15px 'Lucida Console';WHITE-SPACE:pre-wrap;LETTER-SPACING:normal;COLOR:rgb(0,0,0);WORD-SPACING:0px" class="Apple-style-span"><pre style="BORDER-BOTTOM-STYLE:none;LINE-HEIGHT:1.2;BORDER-RIGHT-STYLE:none;MARGIN:0px;OUTLINE-STYLE:none;FONT-FAMILY:'Lucida Console';WHITE-SPACE:pre-wrap!important;BORDER-TOP-STYLE:none;FONT-SIZE:10pt!important;BORDER-LEFT-STYLE:none" class="GNVMTOMCABB" tabindex="0">
<span style="COLOR:blue" class="GNVMTOMCHAB ace_keyword">> Rcpp::sourceCpp('src/sbi.cpp')
</span><span style="WHITE-SPACE:pre;COLOR:blue" class="GNVMTOMCDBB ace_keyword">> </span><span style="COLOR:blue" class="GNVMTOMCHAB ace_keyword">x <- sigmutest(0.0002327279, 0.1025499338)
</span>sg 0.2
sgtemp 0.219135
error 0</pre></span></div>
<div> </div>
<div>Does anyone know what is going on? Thanks.</div>
<div>-- <br>Aileen L.</div>
<div><br>
<div>View my Linkedin profile: <a style="PADDING-BOTTOM:0px;LINE-HEIGHT:15px;BORDER-RIGHT-WIDTH:0px;BACKGROUND-COLOR:rgb(255,255,255);MARGIN:0px;OUTLINE-STYLE:none;PADDING-LEFT:0px;PADDING-RIGHT:0px;FONT-FAMILY:Arial,Helvetica,'Nimbus Sans L',sans-serif;BORDER-TOP-WIDTH:0px;BORDER-BOTTOM-WIDTH:0px;COLOR:rgb(0,102,153);FONT-SIZE:13px;VERTICAL-ALIGN:baseline;BORDER-LEFT-WIDTH:0px;PADDING-TOP:0px" title="View public profile" href="http://au.linkedin.com/in/aileen2" name="SafeHtmlFilter_SafeHtmlFilter_webProfileURL" target="_blank">http://au.linkedin.com/in/aileen2</a></div>

<div><br></div>
<div>
<h2 style="TEXT-ALIGN:left;PADDING-BOTTOM:0px;BACKGROUND-COLOR:rgb(255,255,255);MARGIN:0px;MIN-HEIGHT:35px;PADDING-LEFT:0px;PADDING-RIGHT:2px;FONT-FAMILY:Geneva,Arial,Helvetica,sans-serif;FONT-WEIGHT:normal;PADDING-TOP:2px">
<font color="#993399">Being happy doesn't mean you're perfect. It just means you've decided to look beyond the imperfectio</font><font color="#cc33cc">ns- <a style="PADDING-BOTTOM:0px;MARGIN:0px;PADDING-LEFT:0px;PADDING-RIGHT:0px;TEXT-DECORATION:none;PADDING-TOP:0px" href="http://www.boardofwisdom.com/default.asp?topic=1010&search=K%2EB+Indiana+%28age+14%29" target="_blank">K.B Indiana (age 14)</a></font></h2>
</div></div>