R/twoGaussians.R

In benmack/threshold: Histogram/density based threshold selection

#' @name twoGaussians
#' @title ...
#' @description ...
#' @param n ...
#' @param m.pos ...
#' @param m.neg ...
#' @param sd.pos ...
#' @param sd.neg ...
#' @param seed ...
#' @export 
twoGaussians <- function(n = c(1000, 5000), m.pos=2, m.neg=-2, 
                         sd.pos=1, sd.neg=1, seed=NULL) {
  if (!is.null(seed))
    set.seed(seed)
  pos <- rnorm(n[1], m.pos, sd.pos)
  neg <- rnorm(n[2], m.neg, sd.neg)
  x <- c(pos, neg)
  y <- rep(c(1, -1), n)
  return(data.frame(y=y, x=x))
}