\name{gbd} \alias{gbd} %- Also NEED an '\alias' for EACH other topic documented here. \title{A function to generate some binary regression data %% ~~function to do ... ~~ } \description{ %% ~~ A concise (1-5 lines) description of what the function does. ~~ } \usage{ gbd(n = 1000, p = 1, intercept = TRUE, betas = rep(0, p + intercept), mu = rep(0, p), Sigma = diag(1, p), seed = 1) } %- maybe also 'usage' for other objects documented here. \arguments{ \item{n}{ %% ~~Describe \code{n} here~~ } \item{p}{ %% ~~Describe \code{p} here~~ } \item{intercept}{ %% ~~Describe \code{intercept} here~~ } \item{betas}{ %% ~~Describe \code{betas} here~~ } \item{mu}{ %% ~~Describe \code{mu} here~~ } \item{Sigma}{ %% ~~Describe \code{Sigma} here~~ } \item{seed}{ %% ~~Describe \code{seed} here~~ } } \details{ %% ~~ If necessary, more details than the description above ~~ } \value{ %% ~Describe the value returned %% If it is a LIST, use %% \item{comp1 }{Description of 'comp1'} %% \item{comp2 }{Description of 'comp2'} %% ... } \references{ %% ~put references to the literature/web site here ~ } \author{ %% ~~who you are~~ } \note{ %% ~~further notes~~ } %% ~Make other sections like Warning with \section{Warning }{....} ~ \seealso{ %% ~~objects to See Also as \code{\link{help}}, ~~~ } \examples{ ##---- Should be DIRECTLY executable !! ---- ##-- ==> Define data, use random, ##-- or do help(data=index) for the standard data sets. ## The function is currently defined as function (n = 1000, p = 1, intercept = TRUE, betas = rep(0, p + intercept), mu = rep(0, p), Sigma = diag(1, p), seed = 1) { require(MASS) set.seed(seed) X <- mvrnorm(n, mu, Sigma) Xm <- X if (intercept) { Xm <- cbind(rep(1, n), X) } lp <- Xm \%*\% betas pv <- expit(lp) Y <- rbinom(n, 1, pv) return(list(Y = Y, X = X)) } } % Add one or more standard keywords, see file 'KEYWORDS' in the % R documentation directory. \keyword{ ~kwd1 } \keyword{ ~kwd2 }% __ONLY ONE__ keyword per line