Ejemplo n.º 1
0
// [[Rcpp::export]]
void runit_islocked( Environment env ){
    env.assign( "x1", 1 ) ;
    env.assign( "x2", 10.0 ) ;
    env.assign( "x3", std::string( "foobar" ) ) ;
    env.assign( "x4", "foobar" ) ;
    std::vector< std::string > aa(2) ; aa[0] = "foo" ; aa[1] = "bar" ;
    env.assign( "x5", aa ) ;
}
Ejemplo n.º 2
0
// [[Rcpp::export]]
Environment huffman_encode(NumericVector probabilities) {
  CharacterVector items = probabilities.names();
  Environment env = Environment::empty_env().new_child(true);
  
  visitNodes(internal_huffman_encode(probabilities), "", [&items, &env] (Node *n, string s) {
    if (n->index >= 0) {
      env.assign((const string) items[n->index], s);
    }
    delete n;
  });
  return env;
}
Ejemplo n.º 3
0
// [[Rcpp::export]]
NumericVector evalKernelCpp2(StringMatrix Integrand2, 
                             ExpressionVector Integrand2expr, 
                             Environment myenvd1, Environment myenvd2, LogicalVector CondIntensity,
                             StringVector NameCountingVar,
                             StringVector Namecovariates,
                             LogicalVector ExistdN, LogicalVector ExistdX,
                             List gridTime, IntegerVector dimCol,
                             StringVector NameCol, StringVector JumpTimeName){
  int FinT = gridTime.size();
  int lengtcov= Namecovariates.size();
  NumericVector GlobalKernel(FinT);
  for(int t=0;t<FinT;t++){
    double IntegralKernel = 0;
    if(ExistdN[0]){
      String ttime = myenvd1.get("t.time");
      NumericVector PosInTimedN = gridTime[t];
      myenvd1.assign(ttime, PosInTimedN[0]);
    }
    if(CondIntensity[0]){
      StringVector PosListCountingVariableC = myenvd1.get("PosListCountingVariable");
      int dimCountingVar = NameCountingVar.size();
      for(int i=0;i<dimCountingVar;i++){
        String dummyStringListCount = as<std::string>(PosListCountingVariableC[i]);
        List DummyListCount = myenvd1.get(dummyStringListCount);
        String dummyStringCountN = as<std::string>(NameCountingVar[i]);
        myenvd1.assign(dummyStringCountN,DummyListCount[t]);
      }
    }
    if(lengtcov>0){
      StringVector PosListCovariatesC = myenvd1.get("PosListCovariates");
      for(int h=0;h<lengtcov;h++){
        String dummyStringListcovariates = as<std::string>(PosListCovariatesC[h]);
        List DummyListcovariates = myenvd1.get(dummyStringListcovariates);
        String dummyStringcovariatesN = as<std::string>(Namecovariates[h]);
        myenvd1.assign(dummyStringcovariatesN,DummyListcovariates[t]);
      }
    }
    
    for(int i=0;i<dimCol[0];i++){  
      if(ExistdN[0]){
        StringVector namedJumpTimeX1 = myenvd1.get("namedJumpTimeX"); 
        StringVector namedX1 = myenvd1.get("namedX");
        for(int j=0;j<namedJumpTimeX1.size();j++){
          if(JumpTimeName[i]==namedJumpTimeX1[j]){
            String dummyString = as<std::string>(JumpTimeName[i]);
            List JumpTimeDiff = myenvd1.get(dummyString);
            String vartime =myenvd1.get("var.time");
            myenvd1.assign(vartime,JumpTimeDiff[t]);
            String dumdN = NameCol[i];
            List dNAll = myenvd1.get(dumdN);
            NumericVector dN = dNAll[t];
            if(dN.size()==0){
              dN[0] = 0;
            }
            Language eval_call( "eval", Integrand2expr[i], myenvd1);
            SEXP X = Rf_eval( eval_call, myenvd1);
            NumericVector XX(X);
            NumericVector dumY = na_omit(XX);
            double dumSumKer = 0;
            for(int l=0;l<dumY.size();l++){
              dumSumKer = dumSumKer + dumY[l]*dN[l];
            }
            IntegralKernel = IntegralKernel+dumSumKer;            
          }
        }
      }
      if(ExistdX[0]){
        StringVector namedJumpTimeX2 = myenvd2.get("namedJumpTimeX");
        for(int j=0;j<namedJumpTimeX2.size();j++){
          if(JumpTimeName[i]==namedJumpTimeX2[j]){
            String ttimedX = myenvd2.get("t.time");
            NumericVector PosInTimedX = gridTime[t];
            myenvd2.assign(ttimedX,PosInTimedX[0]);
            String dummyString2 = as<std::string>(JumpTimeName[i]);
            NumericVector JumpTime2 = myenvd2.get(dummyString2);
            String vartime2 =myenvd2.get("var.time");
            IntegerVector idx = seq(0,t);
            myenvd2.assign(vartime2,JumpTime2[idx]);
            Language eval_call( "eval", Integrand2expr[i], myenvd2);
            SEXP X2 = Rf_eval( eval_call, myenvd2);
            NumericVector XX2(X2);
            
            String dummyString3 = as<std::string>(NameCol[i]);
            NumericVector dXdummy = myenvd2.get(dummyString3);
            NumericVector dumY1 = XX2[idx]*dXdummy[idx];
            NumericVector dumY2 = na_omit(dumY1);
            IntegralKernel = IntegralKernel+sum(dumY2);
            
          }
        }
      }
    }
    GlobalKernel[t] = IntegralKernel;
  }
  return GlobalKernel;
}
Ejemplo n.º 4
0
// [[Rcpp::export]]
NumericVector evalKernelCpp(StringMatrix Integrand2, 
                             ExpressionVector Integrand2expr, 
                             Environment myenvd1, Environment myenvd2,
                             LogicalVector ExistdN, LogicalVector ExistdX,
                             NumericVector gridTime, IntegerVector dimCol,
                             StringVector NameCol, StringVector JumpTimeName){
  int FinT = gridTime.size();
  NumericVector GlobalKernel(FinT);
  for(int t=0;t<FinT;t++){
      if(ExistdN[0]){
        String ttime = myenvd1.get("t.time");
        myenvd1.assign(ttime, gridTime[t]);
      }
      if(ExistdX[0]){
        String ttimedX = myenvd2.get("t.time");
        myenvd2.assign(ttimedX,gridTime[t]);
      }
      double IntegralKernel = 0;
      for(int i=0;i<dimCol[0];i++){
        if(ExistdN[0]){
          StringVector namedJumpTimeX1 = myenvd1.get("namedJumpTimeX"); 
          StringVector namedX1 = myenvd1.get("namedX");
          //  printf(namedJumpTimeX1[0]);
          for(int j=0;j<namedJumpTimeX1.size();j++){
            if(JumpTimeName[i]==namedJumpTimeX1[j]){
              String dummyString = as<std::string>(JumpTimeName[i]);
              NumericVector JumpTimeDiff = myenvd1.get(dummyString);
              String vartime =myenvd1.get("var.time");
              myenvd1.assign(vartime,JumpTimeDiff);
            }
          }
          for(int hh=0;hh<NameCol.size();hh++){
            if(namedX1[i]==NameCol[hh]){
              Language eval_call( "eval", Integrand2expr[hh], myenvd1);
              SEXP X = Rf_eval( eval_call, myenvd1);
              NumericVector XX(X);
              NumericVector dumY = na_omit(XX);
              IntegralKernel = IntegralKernel+sum(XX);
            }
          }
/*
          if(any(cond)){
          assign(my.envd1$var.time,my.envd1[[my.envd1$namedJumpTimeX[cond]]],envir=my.envd1)
          condpos <- my.envd1$namedX %in% NameCol[i]  
          if(any(condpos)){
          IntegralKernelDum<- sum(eval(Integrand2expr[condpos], envir=my.envd1))
          IntegralKernel<-IntegralKernel+IntegralKernelDum
          }
          }
          */ 
        }
        if(ExistdX[0]){
          StringVector namedJumpTimeX2 = myenvd2.get("namedJumpTimeX");
          StringVector namedX2 = myenvd2.get("namedX");
          for(int j=0;j<namedJumpTimeX2.size();j++){
            if(JumpTimeName[i]==namedJumpTimeX2[j]){
              String dummyString2 = as<std::string>(JumpTimeName[i]);
              NumericVector JumpTimeDiff2 = myenvd2.get(dummyString2);
              String vartime2 =myenvd2.get("var.time");
              myenvd2.assign(vartime2,JumpTimeDiff2);
            }
          }
          for(int hh=0;hh<NameCol.size();hh++){
            if(namedX2[i]==NameCol[hh]){
              Language eval_call( "eval", Integrand2expr[hh], myenvd2);
              SEXP X2 = Rf_eval( eval_call, myenvd2);
              NumericVector XX2(X2);
              NumericVector dumY2 = na_omit(XX2);
              IntegralKernel = IntegralKernel+sum(dumY2);
            }
          }
          
        }
/*
        if(ExistdX){  
          cond <- my.envd2$namedJumpTimeX %in% JumpTimeName[i]
          if(any(cond)){
            assign(my.envd2$var.time,my.envd2[[my.envd2$namedJumpTimeX[cond]]],envir=my.envd2)
            condpos <- my.envd2$namedX %in% NameCol[i]
            if(any(condpos)){
              IntegralKernelDum<- sum(eval(Integrand2expr[condpos], envir=my.envd2))
              IntegralKernel<-IntegralKernel+IntegralKernelDum
            }
          }
        }
*/
    }
      GlobalKernel[t] = IntegralKernel;
  }
  return GlobalKernel;
}
Ejemplo n.º 5
0
// [[Rcpp::export]]
void runit_assign( Environment env, std::string st, SEXP object ){
    env.assign(st, object) ;
}