I am fairly certain that this line is not kosher (and in any case it is a
confusing):
w = pd * 0.4 / (R::pnorm((qpd - sqrt(rsq) * (-0.42) * sgtemp), 0, 1, 1,
0));
try this instead:
w = qpd - sqrt(rsq) * (-0.42) * sgtemp
w = pd * 0.4 / (R::pnorm(w, 0, 1, 1, 0));
On Wed, Mar 13, 2013 at 11:42 P
On 14 March 2013 at 14:42, Aileen Lin wrote:
| My C code:
| //[[Rcpp::depends("Rcpp")]]
| #include
| #include
| using namespace Rcpp;
|
| //[[Rcpp::export]]
| NumericVector sigmutest(double pd, double rsq){
| double qpd = R::qnorm(pd, 0, 1, 1, 0);
| double sgtemp = 0.2;
| double sg = 0.3
My C code:
//[[Rcpp::depends("Rcpp")]]
#include
#include
using namespace Rcpp;
//[[Rcpp::export]]
NumericVector sigmutest(double pd, double rsq){
double qpd = R::qnorm(pd, 0, 1, 1, 0);
double sgtemp = 0.2;
double sg = 0.3;
double eor = 1;
double w = 0;
while (eor>=0.0001) {
