Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import numpy as np
- from rpy2.robjects.packages import importr
- import rpy2.robjects as robjects
- r = robjects.r #allows access to r object with r.
- bcp = importr('bcp') #import bayesian change point package in python
- values = bcp.bcp( r.c( r.rnorm(50) , r.rnorm(50,5,1), r.rnorm(50) ) ) #use bcp function on vector
- posterior_means = np.array(values[5]).flatten()
- posterior_probability = np.array(values[7]).flatten()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement