data(testdata)
x=testdata[,9:248]
p=dim(x)[2]
for (i in 1:p){x[,i]=as.numeric(as.character(x[,i]))}
gr=testdata[,2]
gr=as.numeric(gr)
gr[which(gr<4)]=0
gr[which(gr==4)]=1
result=ZIBseq(data=x,outcome=gr)
Run the code above in your browser using DataLab