R/testData2.R

In ReorderCluster: Reordering the Dendrogram According to the Class Labels

testData2<-function()
{
  num_gene=90
  num_exp=90
  num=10
  var=40
  p=0.01
  mat=runif(num_gene*num_exp)
  dim(mat)=c(num_gene,num_exp)
  class=numeric(num_gene)
  label=c(1,2,3)
  for(i in 1:6)
  {
    start=1+15*(i-1)
    end=start+15-1
    s=1+(i-1)*num
    e=s+var-1
    
    t=runif((end-start+1)*(e-s+1))
    index=which(t<=p)
    t=rep(-1,(end-start+1)*(e-s+1))
    pp=runif(length(index))
    for(j in 1:length(index))
    {
      if(pp[j]<0.5) 
        t[index[j]]=0
      else
        t[index[j]]=1
    }
    mat[start:end,s:e]=t
    if(i<=3)
    class[start:end]=label[i]
    else
    {
      if(i<=6)
    class[start:end]=label[i-3] 
      else
        class[start:end]=label[i-6] 
    }
  }
  
    
  image(mat)
  browser()
  mat=cbind(mat,class)
  colnames(mat)=c(paste("F",1:num_exp,sep=""),"class")
  mat
}