crp_gibbs: Markov chain Monte Carlo methods for CRP clustering

In CRPClustering: Bayesian Nonparametric Chinese Restaurant Process Clustering with Entropy

Description

Markov chain Monte Carlo methods for CRP clustering

Usage

crp_gibbs(data, mu = c(0, 0), sigma_table = 14, alpha = 0.3, ro_0 = 0.1,
  burn_in = 40, iteration = 200)

Arguments

data	: a matrix of data for clustering. Row is each data_i and column is dimensions of each data_i.
mu	: a vector of center points of data. If data is 3 dimensions, a vector of 3 elements like c(2,4,2).
sigma_table	: a numeric of CRP variance.
alpha	: a numeric of a CRP concentrate rate.
ro_0	: a numeric of a CRP mu change rate.
burn_in	: an iteration integer of burn in.
iteration	: an iteration integer.

Value

z_result : an array expresses cluster numbers for each data_i.

Examples

data <- matrix(c(1.8,1.9,2.1,2.5,5.6,5.2,6,6.1), 4, 2)
z_result <- crp_gibbs(
                      data,
                      mu=c(0,0),
                      sigma_table=14,
                      alpha=0.3,
                      ro_0=0.1,
                      burn_in=10,
                      iteration=100
                     )

Example output

[1] "simulation start."
[1] "In preparation :  1"
[1] "In preparation :  2"
[1] "In preparation :  3"
[1] "In preparation :  4"
[1] "Iteration start."
[1] "iteration :  2 , cluster numbers in the serch spaace :  1"
[1] "iteration :  3 , cluster numbers in the serch spaace :  1"
[1] "iteration :  4 , cluster numbers in the serch spaace :  1"
[1] "iteration :  5 , cluster numbers in the serch spaace :  1"
[1] "iteration :  6 , cluster numbers in the serch spaace :  1"
[1] "iteration :  7 , cluster numbers in the serch spaace :  1"
[1] "iteration :  8 , cluster numbers in the serch spaace :  1"
[1] "iteration :  9 , cluster numbers in the serch spaace :  1"
[1] "iteration :  10 , cluster numbers in the serch spaace :  1"
[1] "iteration :  11 , cluster numbers in the serch spaace :  1"
[1] "iteration :  12 , cluster numbers in the serch spaace :  1"
[1] "iteration :  13 , cluster numbers in the serch spaace :  1"
[1] "iteration :  14 , cluster numbers in the serch spaace :  1"
[1] "iteration :  15 , cluster numbers in the serch spaace :  1"
[1] "iteration :  16 , cluster numbers in the serch spaace :  1"
[1] "iteration :  17 , cluster numbers in the serch spaace :  1"
[1] "iteration :  18 , cluster numbers in the serch spaace :  1"
[1] "iteration :  19 , cluster numbers in the serch spaace :  1"
[1] "iteration :  20 , cluster numbers in the serch spaace :  1"
[1] "iteration :  21 , cluster numbers in the serch spaace :  1"
[1] "iteration :  22 , cluster numbers in the serch spaace :  1"
[1] "iteration :  23 , cluster numbers in the serch spaace :  1"
[1] "iteration :  24 , cluster numbers in the serch spaace :  1"
[1] "iteration :  25 , cluster numbers in the serch spaace :  1"
[1] "iteration :  26 , cluster numbers in the serch spaace :  1"
[1] "iteration :  27 , cluster numbers in the serch spaace :  1"
[1] "iteration :  28 , cluster numbers in the serch spaace :  1"
[1] "iteration :  29 , cluster numbers in the serch spaace :  1"
[1] "iteration :  30 , cluster numbers in the serch spaace :  1"
[1] "iteration :  31 , cluster numbers in the serch spaace :  1"
[1] "iteration :  32 , cluster numbers in the serch spaace :  1"
[1] "iteration :  33 , cluster numbers in the serch spaace :  1"
[1] "iteration :  34 , cluster numbers in the serch spaace :  1"
[1] "iteration :  35 , cluster numbers in the serch spaace :  1"
[1] "iteration :  36 , cluster numbers in the serch spaace :  1"
[1] "iteration :  37 , cluster numbers in the serch spaace :  1"
[1] "iteration :  38 , cluster numbers in the serch spaace :  1"
[1] "iteration :  39 , cluster numbers in the serch spaace :  1"
[1] "iteration :  40 , cluster numbers in the serch spaace :  1"
[1] "iteration :  41 , cluster numbers in the serch spaace :  1"
[1] "iteration :  42 , cluster numbers in the serch spaace :  1"
[1] "iteration :  43 , cluster numbers in the serch spaace :  1"
[1] "iteration :  44 , cluster numbers in the serch spaace :  2"
[1] "iteration :  45 , cluster numbers in the serch spaace :  1"
[1] "iteration :  46 , cluster numbers in the serch spaace :  1"
[1] "iteration :  47 , cluster numbers in the serch spaace :  1"
[1] "iteration :  48 , cluster numbers in the serch spaace :  1"
[1] "iteration :  49 , cluster numbers in the serch spaace :  1"
[1] "iteration :  50 , cluster numbers in the serch spaace :  1"
[1] "iteration :  51 , cluster numbers in the serch spaace :  1"
[1] "iteration :  52 , cluster numbers in the serch spaace :  1"
[1] "iteration :  53 , cluster numbers in the serch spaace :  1"
[1] "iteration :  54 , cluster numbers in the serch spaace :  1"
[1] "iteration :  55 , cluster numbers in the serch spaace :  1"
[1] "iteration :  56 , cluster numbers in the serch spaace :  1"
[1] "iteration :  57 , cluster numbers in the serch spaace :  1"
[1] "iteration :  58 , cluster numbers in the serch spaace :  1"
[1] "iteration :  59 , cluster numbers in the serch spaace :  1"
[1] "iteration :  60 , cluster numbers in the serch spaace :  1"
[1] "iteration :  61 , cluster numbers in the serch spaace :  2"
[1] "iteration :  62 , cluster numbers in the serch spaace :  1"
[1] "iteration :  63 , cluster numbers in the serch spaace :  1"
[1] "iteration :  64 , cluster numbers in the serch spaace :  1"
[1] "iteration :  65 , cluster numbers in the serch spaace :  1"
[1] "iteration :  66 , cluster numbers in the serch spaace :  1"
[1] "iteration :  67 , cluster numbers in the serch spaace :  1"
[1] "iteration :  68 , cluster numbers in the serch spaace :  1"
[1] "iteration :  69 , cluster numbers in the serch spaace :  1"
[1] "iteration :  70 , cluster numbers in the serch spaace :  1"
[1] "iteration :  71 , cluster numbers in the serch spaace :  1"
[1] "iteration :  72 , cluster numbers in the serch spaace :  1"
[1] "iteration :  73 , cluster numbers in the serch spaace :  1"
[1] "iteration :  74 , cluster numbers in the serch spaace :  1"
[1] "iteration :  75 , cluster numbers in the serch spaace :  1"
[1] "iteration :  76 , cluster numbers in the serch spaace :  1"
[1] "iteration :  77 , cluster numbers in the serch spaace :  1"
[1] "iteration :  78 , cluster numbers in the serch spaace :  1"
[1] "iteration :  79 , cluster numbers in the serch spaace :  1"
[1] "iteration :  80 , cluster numbers in the serch spaace :  1"
[1] "iteration :  81 , cluster numbers in the serch spaace :  1"
[1] "iteration :  82 , cluster numbers in the serch spaace :  1"
[1] "iteration :  83 , cluster numbers in the serch spaace :  1"
[1] "iteration :  84 , cluster numbers in the serch spaace :  1"
[1] "iteration :  85 , cluster numbers in the serch spaace :  1"
[1] "iteration :  86 , cluster numbers in the serch spaace :  1"
[1] "iteration :  87 , cluster numbers in the serch spaace :  1"
[1] "iteration :  88 , cluster numbers in the serch spaace :  1"
[1] "iteration :  89 , cluster numbers in the serch spaace :  1"
[1] "iteration :  90 , cluster numbers in the serch spaace :  1"
[1] "iteration :  91 , cluster numbers in the serch spaace :  1"
[1] "iteration :  92 , cluster numbers in the serch spaace :  1"
[1] "iteration :  93 , cluster numbers in the serch spaace :  1"
[1] "iteration :  94 , cluster numbers in the serch spaace :  1"
[1] "iteration :  95 , cluster numbers in the serch spaace :  1"
[1] "iteration :  96 , cluster numbers in the serch spaace :  1"
[1] "iteration :  97 , cluster numbers in the serch spaace :  1"
[1] "iteration :  98 , cluster numbers in the serch spaace :  1"
[1] "iteration :  99 , cluster numbers in the serch spaace :  1"
[1] "iteration :  100 , cluster numbers in the serch spaace :  1"
[1] "Center coordinates of Clusters"
  cluster   mu_k.1   mu_k.2
1       1 2.123107 6.369863
[1] "The Number of Total Clusters  :  1"
[1] "Entropy All :  0"
[1] "z_result has numbers of clusters."

CRPClustering documentation built on May 1, 2019, 7:11 p.m.