<div dir="ltr"><div>Dear Sir</div><div><br></div><div><pre tabindex="0" class="gmail-GD15MCFCEUB" id="gmail-rstudio_console_output" style="font-family:"Lucida Console",monospace;font-size:16px;outline:currentcolor none medium;border:medium none;word-break:break-all;margin:0px;white-space:pre-wrap;line-height:1.2;color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><span><span class="gmail-GD15MCFCOTB gmail-ace_constant" style="color:rgb(197,6,11)"></span></span></pre></div><div>Here is my MWE. <br></div><div><br></div><div>I am still getting the follwoing error.<br></div><div><br></div><div>
Error in dmvnormC(x, mean, sigma, 0) : 
  Cannot convert object to a function: [type=NULL; target=CLOSXP, SPECIALSXP, or BUILTINSXP].

<table class="gmail-GD15MCFCNTB gmail-ace_text-layer gmail-ace_line gmail-GD15MCFCAW" style="font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-numeric:inherit;font-variant-east-asian:inherit;font-weight:400;font-stretch:inherit;font-size:12pt;line-height:1.2;font-family:"Lucida Console",monospace;padding-left:6px;padding-bottom:8px;border:medium none;outline:currentcolor none medium;white-space:pre-wrap;color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;width:468.546px" cellspacing="0" cellpadding="0"><tbody><tr><td style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" align="left"><table style="width:468.788px" cellspacing="0" cellpadding="0"><tbody><tr><td rowspan="1" style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" width="1" height="" align="left"><br>Could you please guide what is going wrong?<br>
<div>y <- matrix(1:10, nrow=5)</div><div>mean <- c(4, 5)</div><div>sigma <- matrix(c(0.01, 0.009, 0.009, 1.2), nrow =2, byrow = TRUE)</div><div>dmvnormC(y, mean, sigma, 0) # Function is below.</div>
#include <RcppArmadillo.h><br>#include <mvtnormAPI.h><br>using namespace Rcpp;<br>using namespace RcppArmadillo;<br>using namespace arma;<br><br>//[[Rcpp::depends(RcppArmadillo)]]<br>//[[Rcpp::depends(mvtnorm)]]<br>//[[Rcpp::export]]<br><br>arma::vec  dmvnormC(const arma::mat& x, const arma::vec& mean, const arma::mat& sigma, int give_log)<br>{<br>  Environment pkg = Environment::namespace_env("mvtnorm");<br>  Function dmvnorm = pkg["mvtnorm"];<br>  int N = x.n_rows;<br>  arma::rowvec f(N);<br>  for(int k=0; k<N; k++)<br>    {<br>      f(k) = Rcpp::as<double>(dmvnorm(x.row(k), mean,  sigma, give_log));<br>    }<br>  return f;<br>}<br><br><br>Thank you <br></td></tr></tbody></table></td></tr></tbody></table><table class="gmail-GD15MCFCNTB gmail-ace_text-layer gmail-ace_line gmail-GD15MCFCAW" style="font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-numeric:inherit;font-variant-east-asian:inherit;font-weight:400;font-stretch:inherit;font-size:12pt;line-height:1.2;font-family:"Lucida Console",monospace;padding-left:6px;padding-bottom:8px;border:medium none;outline:currentcolor none medium;white-space:pre-wrap;color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;width:468.546px" cellspacing="0" cellpadding="0"><tbody><tr><td style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" align="left"><table style="width:468.788px" cellspacing="0" cellpadding="0"><tbody><tr><td rowspan="1" style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" width="1" height="" align="left"><br></td></tr></tbody></table></td></tr></tbody></table></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Dec 20, 2019 at 12:19 AM Ralf Stubner <<a href="mailto:ralf.stubner@gmail.com">ralf.stubner@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 19, 2019 at 8:01 AM Shaami <<a href="mailto:nzshaam@gmail.com" target="_blank">nzshaam@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">  Environment pkg = Environment::namespace_env("mvtnorm");<br><div>  Function dmvnorm = pkg["mvtnorm"];</div></div></blockquote><div><br></div><div>There is no mvtnorm function in the mvtnorm package. You probably want to dmvnorm function instead.</div><div><br></div><div>BTW, the mvtnorm package provides a C API for the dmvnorm function, c.f. <a href="https://stackoverflow.com/a/51294364/8416610" target="_blank">https://stackoverflow.com/a/51294364/8416610</a></div><div><br></div><div>cheerio</div><div>ralf<br></div></div></div>
</blockquote></div>