# define Rcpp file
code1 <- "
// array A
Rcpp::NumericMatrix AA(A);
// Rcpp::IntegerVector dimAA(dimA);
int nrows = AA.nrow();
int ncolumns = AA.ncol();
NumericMatrix Alogis(nrows,ncolumns) ;
// compute logistic distribution
for (int ii=0; ii<nrows; ++ii){
NumericVector h1=AA.row(ii) ;
NumericVector res = plogis( h1 ) ;
for (int jj=0;jj<ncolumns;++jj){
Alogis(ii,jj) = res[jj] ;
}
}
return( wrap(Alogis) );
"
# compile Rcpp code
calc1 <- cxxfunction( signature( A= "matrix"), code1, plugin = "Rcpp", verbose=TRUE )
cxxfunction.copy( cppfct=calc1, name="calclogis" )
Run the code above in your browser using DataLab