CountDataSet1 <- function (n, p, K, param, sdsignal, drate) 

{

    if (n < 4 * K) 

        stop("We require n to be at least 4*K.")

    q0 <- rexp(p, rate = 1/25)

    isDE <-  runif(p)< drate    #runif(p)

    classk <- matrix(NA, nrow = K, ncol = p)

    for (k in 1:K) {

        lfc <- rnorm(p, sd = sdsignal)

        classk[k, ] <- ifelse(isDE, q0 * exp(lfc), q0)

    }

    truesf <- runif(n) * 2 + 0.2

    truesfte <- runif(n) * 2 + 0.2

    conds <- sample(c(rep(1:K, 4), sample(1:K, n - 4 * K, replace = TRUE)))

    condste <- sample(c(rep(1:K, 4), sample(1:K, n - 4 * K, replace = TRUE)))

    x <- xte <- matrix(NA, nrow = n, ncol = p)

    for (i in 1:n) {

        for (k in 1:K) {

            if (conds[i] == k) 

                x[i, ] <- rnbinom(p, mu = truesf[i] * classk[k, 

                  ], size = param)

            if (condste[i] == k) 

                xte[i, ] <- rnbinom(p, mu = truesfte[i] * classk[k, 

                  ], size = param)

        }

    }

    rm <- apply(x, 2, sum) == 0

    return(list(x = x[, !rm], xte = xte[, !rm], y = conds, yte = condste, 

        truesf = truesf, truesfte = truesfte, dkg=classk))

}