我用下面的程序产生了两组随机数,y0和y1,分别是5*20的矩阵形式。
行代表样品,列代表基因——第i行第j列的数表示第i个病人第j个基因的表达量
就现在想要比较y0和y1的基因是否有差异。请问改怎么做呢?
对于进来来看的前辈,小女子不胜感激。
B<-runif(20,min=0.1,max=1)<br />
BB<-rep(runif(1,min=0.1,max=1),each=20)<br />
BBB<-2*sqrt(B+BB)<br />
<br />
#creat delta0 and delta1<br />
<br />
delta0<-rep(0,each=20)<br />
delta1<-rep(0,each=20) <br />
for (j in 1:5){<br />
delta1[j]<-runif(1,min=0,max=BBB[j])<br />
j=j+1<br />
}<br />
<br />
<br />
#creat eta0 and eta1<br />
eta0<-rnorm(5,mean=0,sd=sqrt(BB[1]))<br />
eta1<-rnorm(5,mean=0,sd=sqrt(BB[1]))<br />
<br />
#creat epsilon0 and epsilon1<br />
epsilon0<-matrix(data=NA,nrow=5,ncol=20,byrow=T)<br />
epsilon1<-matrix(data=NA,nrow=5,ncol=20,byrow=T)<br />
for (j in 1:20 ) {<br />
for (i in 1:5) {<br />
epsilon0[i,j]=rnorm(1,mean=0,sd=sqrt(B[j]))<br />
epsilon1[i,j]=rnorm(1,mean=0,sd=sqrt(B[j]))<br />
i=i+1<br />
} <br />
j=j+1<br />
}<br />
<br />
#y = delta + eta + epsilon<br />
y0<-matrix(data=NA,nrow=5,ncol=20,byrow=T)<br />
y1<-matrix(data=NA,nrow=5,ncol=20,byrow=T)<br />
for (j in 1:20 ) {<br />
for (i in 1:5 ) {<br />
y0[i,j]=delta0[j]+eta0[i]+epsilon0[i,j]<br />
y1[i,j]=delta1[j]+eta1[i]+epsilon1[i,j]<br />
i=i+1<br />
} <br />
j=j+1<br />
}