CTGEIM4：运行Gibbs采样器以使用边缘模型之前的无限混合模型对时程基因表达数据进行聚类

入门 library(CTGEIM4) 生成时程基因表达数据 library(mvtnorm) n=100;m=10 TT=seq(0,1,length=m) Y=mu=mu1=matrix(0,n,m) x1=rnorm(n);x2=rbinom(n,1,.1) b=rnorm(n) K=6 index1=1:round(n/K) index2=(round(n/K)+1):(round(2*n/K)) index3=(round(2*n/K)+1):round(3*n/K) index4=(round(3*n/K)+1):round(4*n/K) index5=(round(4*n/K)+1):round(5*n/K) index6=(round(5*n/K)+1):n cluster=rep(0,n) cluster[index1]=1 cluster[index2]=2 clust