<span class="Apple-style-span" style="font-family: tahoma, sans-serif; border-collapse: collapse; "><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; ">
<span class="Apple-style-span" style="font-size: small; ">To the Maintainers of the package fBasics,</span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Arial; min-height: 15px; ">
<span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif"><br></font></span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Arial; ">
<span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif">I thought you might be interested in the attached code written by Robert B. Gramacy and me which. It is inspired by the rgig function in ghyp which is essentially identical to the .rgigjd and .rgigjd1 functions in fBasics. </font></span><span style="font: normal normal normal 13px/normal Arial; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif">The attached code is a merging of a C-version of the code in rgig and some code from <a href="http://www.cran.r-project.org/web/packages/monomvn/index.html">monomvn</a> (</font></span></span><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif">a package for estimation of multivariate normal and student-t data of arbitrary dimension where the pattern of missing data is monotone).</font></span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 16px/normal Times; min-height: 19px; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif"><br>
</font></span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">I am sending this in the spirit of </span><span style="font: normal normal normal 13px/normal Arial; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif">the open source community because I believe that </font></span></span><span class="Apple-style-span" style="font-size: small; ">the C-version, being much faster, could contribute to your package (see file time.R for a comparison between the R and C codes). Some details about the code follow. The R function uniroot has been replaced by a modified version of ZEROIN from the C math library (algorithm from G.Forsythe, M.Malcolm, C.Moler). Inverse gaussian variates (case lambda=-0.5) are generated by the function rinvgauss (algorithm from Michael,Schucany,Haas as given in Gentle's book on page 193) essentially in the same way as in ghyp. </span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 16px/normal Times; min-height: 19px; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif"><br>
</font></span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(51, 51, 51); line-height: 16px; "><span class="Apple-style-span" style="font-size: small; ">Should you have any questions, please do not hesitate to ask.</span></span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 16px/normal Times; min-height: 19px; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif"><br>
</font></span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">Best regards,</span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 16px/normal Times; min-height: 19px; "><span class="Apple-style-span" style="font-size: small; "><font class="Apple-style-span" face="tahoma, sans-serif"><br>
</font></span></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">Ester Pantaleo</span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 16px/normal Times; min-height: 19px; "><font class="Apple-style-span" face="tahoma, sans-serif"><span class="Apple-style-span" style="font-size: small; "><br>
</span></font></p><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">Dipartimento di Fisica</span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">Universita' di Bari</span></p>
<p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Tahoma; "><span class="Apple-style-span" style="font-size: small; ">Italy</span></p></span>