"sim.data" <-
function(N, prop, n, mu.A.mu=0, mu.A.sd=.5, mu.B=0, sd.A=.45, sd.B=.45)
  {
    nA <- as.integer(prop*N)
    nB <- N-nA
    mu.A <- rnorm(nA, mu.A.mu, mu.A.sd)
    MBd.M <- matrix(nrow=N,ncol=n)
    for(i in 1:n) {
      MBd.M[1:nA,i] <- rnorm(nA,mu.A,sd.A)
      MBd.M[(nA+1):N,i] <- rnorm(nB,mu.B,sd.B)
    }
    MBd <- new("MBdataDE")
    MBd <- addData(MBd,N=N,n=n,M=MBd.M)
    MBd <- setProp(MBd,prop)

    return(MBd)
  }