<html 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=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:"Courier New";}
span.gnkrckgcasb
{mso-style-name:gnkrckgcasb;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=EN-US link=blue vlink="#954F72"><div class=WordSection1><p class=MsoNormal>@Dirk: I didn’t want to overwhelm those on the list with tons of code. Below is the full source for the function. What else would you need to see? </p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>After clearing my environment variables, I am now getting the following error:</p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Error in logit(30, 0, 100) : </p><p class=MsoNormal> function 'enterRNGScope' not provided by package 'Rcpp'</p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I do have some roxygen documentation errors in my package? Might this have an effect in producing the error?</p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here is the function:</p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-family:"Courier New"'>// [[Rcpp::export]]<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'>Rcpp::NumericVector logit(Rcpp::NumericVector x, <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> Rcpp::NumericVector xMin,<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> Rcpp::NumericVector xMax, <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> double adjustment = 0.3) {<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> int N = x.length();<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> int nMin = xMin.length();<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> int nMax = xMax.length();<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> Rcpp::NumericVector logOgive(N,0.0);<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(!(nMin == 1 || nMin == N)) <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> return logOgive;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(!(nMax == 1 || nMax == N)) <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> return logOgive;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(xMin.length() != xMax.length()) <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> return logOgive;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> int justOne;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(xMin.length() == 1) <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> justOne = 1; <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> else <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> justOne = 0;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> NumericVector upperBound = xMax - adjustment;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> NumericVector lowerBound = xMin + adjustment;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> for(int i = 0; i < N; i++) {<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(x[i] < 0) continue;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(justOne == 1) {<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(x[i] <= xMin[0]) x[i] = lowerBound[i];<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(x[i] >= xMax[0]) x[i] = upperBound[i];<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> logOgive[i] = log((x[i] - xMin[0])/(xMax[0] - x[i]));<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> } else {<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(x[i] <= xMin[i]) x[i] = lowerBound[i];<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> if(x[i] >= xMax[i]) x[i] = upperBound[i];<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> logOgive[i] = log((x[i] - xMin[i])/(xMax[i] - x[i]));<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> }<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> }<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'> return logOgive;<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Courier New"'>}<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Barth</p><p class=MsoNormal><o:p> </o:p></p></div></body></html>