\name{est.bes}
\alias{est.bes}
%- Also NEED an '\alias' for EACH other topic documented here.
\title{A function to estimate betas in a logistic regression
%%  ~~function to do ... ~~
}
\description{
%%  ~~ A concise (1-5 lines) description of what the function does. ~~
}
\usage{
est.bes(Y, X, intercept = TRUE, start.vals = rep(0, ncol(X) + as.numeric(intercept)))
}
%- maybe also 'usage' for other objects documented here.
\arguments{
  \item{Y}{
%%     ~~Describe \code{Y} here~~
}
  \item{X}{
%%     ~~Describe \code{X} here~~
}
  \item{intercept}{
%%     ~~Describe \code{intercept} here~~
}
  \item{start.vals}{
%%     ~~Describe \code{start.vals} 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 (Y, X, intercept = TRUE, start.vals = rep(0, ncol(X) + 
    as.numeric(intercept))) 
{
    x <- start.vals
    opt <- optim(x, fn = lr.loglike, gr = NULL, Y = Y, X = X, 
        intercept = intercept, method = "BFGS", control = list(fnscale = -1), 
        hessian = T)
    beta.hat <- opt$par
    beta.cov <- solve(opt$hessian)/length(Y)
    if (intercept) {
        Xm <- cbind(rep(1, length(Y)), X)
    }
    else {
        Xm <- X
    }
    fitted.lp <- Xm \%*\% beta.hat
    fitted.Y <- expit(fitted.lp)
    conv <- TRUE
    if (opt$convergence != 0) {
        conv <- FALSE
        print("Optim did not converge")
    }
    return(list(ests = beta.hat, varcov = beta.cov, fitted.lp = fitted.lp, 
        fitted.Y = fitted.Y, conv = conv))
  }
}
% Add one or more standard keywords, see file 'KEYWORDS' in the
% R documentation directory.
\keyword{ ~kwd1 }
\keyword{ ~kwd2 }% __ONLY ONE__ keyword per line