-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathupdatebeta_RW.h
More file actions
15 lines (14 loc) · 811 Bytes
/
updatebeta_RW.h
File metadata and controls
15 lines (14 loc) · 811 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
//update beta with Random walk proposal
#include <Rcpp.h>
#include <gsl_rng.h>
#include <gsl_randist.h>
//#include <math.h>
#include <stdlib.h>
//#include <time.h>
#include <R_ext/Utils.h>
using namespace std;
void updatebeta_RW(vector<double>& xmust, vector<double>& xmuut, vector<int>& xn_s, vector<int>& xn_u, int xI, int xK, int xM,
int K1, int t1, double xp_var, double sqrt_var1, double sqrt_var2, int xtt, vector<int>& xgammat,
Rcpp::IntegerMatrix& xd, Rcpp::NumericMatrix& xybar_s, Rcpp::NumericMatrix& xybar_u, Rcpp::NumericMatrix& pt1, Rcpp::NumericMatrix& pt2,
Rcpp::NumericMatrix& pt3, double xlambda, Rcpp::NumericVector& beta, double xalpha, double xsig_beta1,
double xbeta1, Rcpp::IntegerVector& xAbeta);