search for: logitreg

...I am conducting logistic regression on data set on psi (ESP) ganzfeld trials. The response variable is binary (correct/incorrect), with a 25% guessing base rate. Dr. Ripley suggested that I modify a maximum likelihood fitting of a binomial logistic regression presented in MASS$ (p. 445): logitreg <- function(x, y, wt=rep(1, length(y)), intercept = T, start=rep(0,p), ...) { fmin <- function (beta, X, y, w){ p <- plogis(X %*% beta) -sum(2 * w * ifelse(y, log(p), log(1-p))) } gmin <- function (beta, X, y, w) { eta <- X %*% beta; p &...

...erly. I would expect the outputs of this function to be roughly equivalent to the results from object "p" as seen above, but they don't look similar at all. I'm not sure what this means. tBLparamsForOneSpecies <- function(theData) { x <- theData[ ,1] y <- theData[ ,2] logitReg <- glm(y ~ x + I(x^2), family=binomial) b0 <- logitReg$coefficients[1] b1 <- logitReg$coefficients[2] b2 <- logitReg$coefficients[3] opt <- (-b1)/(2*b2) tol <- 1 / sqrt(-2*b2) pmax<-1/(1+exp(b1^2/(4*b2)-b0)) theParams <- cbind(opt, tol, pmax) theParams } >theDat...

Dear R users, I am using the MICE package. Specifically, at some point in my code I have imp2=mice(PoptotalMICE,imputationMethod="logreg2") And R returns... iter imp variable 1 1 MICEYError in logitreg(xobs, yobs, intercept=F) : couldn't find function "ms" I have been looking for this ms function on the web, hoping it was just a matter of downloading a package. I didn't find anything. I have tried getAnywhere without success. What can it be? Where is the function ms?...

Hello, I am looking for a routine to do a logistic regression. In the book "Modern Applied Statistics with S-PLUS" a function is described which uses the 'ms' command. Is there an analogue for R, or an alternative approach that can accomplish the same thing? Thanks, John. -- ========================================== John Janmaat Department of Economics Acadia

Dear R Users/Experts, I am using a function called logitreg() originally described in MASS (the book 4th Ed.) by Venebles & Ripley, p445. I used the code as provided but made couple of changes to run a 'constrained' logistic regression, I set the method = "L-BFGS-B", set lower/upper values for the variables. Here is the function, lo...

Sorry, I deleted my old post. Pasting the new query below. Dear R Users/Experts, I am using a function called logitreg() originally described in MASS (the book 4th Ed.) by Venebles & Ripley, p445. I used the code as provided but made couple of changes to run a 'constrained' logistic regression, I set the method = "L-BFGS-B", set lower/upper values for the variables. Here is the function, l...

...using optim() give "L-BFGS-B" error, please help To: r-help at r-project.org Message-ID: <19733974.post at talk.nabble.com> Content-Type: text/plain; charset=us-ascii Sorry, I deleted my old post. Pasting the new query below. Dear R Users/Experts, I am using a function called logitreg() originally described in MASS (the book 4th Ed.) by Venebles & Ripley, p445. I used the code as provided but made couple of changes to run a 'constrained' logistic regression, I set the method = "L-BFGS-B", set lower/upper values for the variables. Here is the function, l...