入门
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=(
