Rcpp实现R与C++的无缝连接

         R的for循环非常慢，当看见这句话时：Thanks to the abstractions provided by Rcpp, the code implementing gibbs in C++ is nearly identical to the code you’d write in R, but runs 20 times faster.采用Rcpp实现R和C++的无缝连接是非常有必要的，至少从效率上讲。下面简单的介绍下我的初步成果：

         首先在R里面安装包Rcpp：install.packages('Rcpp')

         然后编写自己的c++文件，假设在Rstudio里面编写c/c++ script 如下：一个简单的循环输出循环变量程序

#include <Rcpp.h>//必须的头文件
using namespace Rcpp;

// [[Rcpp::export]]
void gibbs(int N, int thin) {
  for(int i = 0; i < N; i++) {
    for(int j = 0; j < thin; j++) {
      std::cout<<i<<" "<<j<<std::endl;//可以在前面using namespace std;
    }
  }
  //return(N);//可以定义返回值和函数返回值相同然
}

//下面是R脚本程序
/*** R
#格式非得这样:/*** R  
gibbs(3,1)
*/

     脚本输出：

0 0
1 0
2 0

      假设上面的程序保存为test.cpp，然后在R里先加载包library(Rcpp),载在R里运行命令：sourceCpp("test.cpp") 就会有输出结果了，注意这里假设test.cpp位于R的工作目录下，若不在R的工作目录下可以在test.cpp前加上路径,如：sourceCpp("../../../test.cpp")。

      若是在Rstudio里面的话还有个更简单的方法：在窗口内写c/c++ script后同时按:ctrl+shift+enter后会自动运行脚本。

      这里只是编写的c/c++程序为什么会有输出呢？这是因为test.cpp的最后几句:/*** R  后那几句是R脚本程序，这里相当于执行gibbs(3,1)函数，那个/*** R格式是固定的，改了后会报错~