Displaying 6 results from an estimated 6 matches for "nconvergence".
Did you mean:
convergence
2005 Feb 20
1
logistic regression and 3PL model
...;, ...)
names(fit$par) <- dn
cat("\nCoefficients:\n"); print(fit$par)
# R: use fit$value and fit$convergence
cat("\nResidual Deviance:", format(fit$objective), "\n")
cat("\nConvergence message:", fit$message, "\n")
invisible(fit)
}
I have been using "lower=.25" in the call for logitreg:
options(contrasts = c("contr.treatment", "contr.poly"))
X <- model.matrix(longhit ~ age*gender, data=data)[,-1]
logitre...
2002 Feb 22
1
Logit / 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
2008 Sep 28
0
constrained logistic regression: Error in optim() with method = "L-BFGS-B"
...list(trace = 6), method = "L-BFGS-B", lower = c(-Inf, 0.05), upper = c(Inf,
Inf), ...)
names(fit$par) <- dn
cat("\nCoefficients:\n"); print(fit$par)
cat("\nResidual Deviance:", format(fit$value), "\n")
if(fit$convergence > 0) cat("\nConvergence code:", fit$convergence,
"\n")
return(list(fitpar = fit$par))
invisible(fit)
}
And here is my data,
# ----------------------
# Dose 0 1 emp.prop
# ----------------------
# 1 3 0 0.000
# 2 1 0 0.000
# 3.3 4 0 0.000
# 5.016 3 0 0.000
#...
2008 Sep 29
0
Logistic Regression using optim() give "L-BFGS-B" error, please help
...ntrol = list(trace =
6), method = "L-BFGS-B", lower = c(-Inf, 0.05), upper = c(Inf, Inf), ...)
names(fit$par) <- dn
cat("\nCoefficients:\n"); print(fit$par)
cat("\nResidual Deviance:", format(fit$value), "\n")
if(fit$convergence > 0) cat("\nConvergence code:", fit$convergence, "\n")
return(list(fitpar = fit$par))
invisible(fit)
}
And here is my data,
# ----------------------
# Dose 0 1 emp.prop
# ----------------------
# 1 3 0 0.000
# 2 1 0 0.000
# 3.3 4 0 0.000
# 5.016 3 0 0.000
# 7.0224 4 0 0.000
# 9.3398 2 0 0.000...
2006 Apr 16
3
second try; writing user-defined GLM link function
I apologize for my earlier posting that, unbeknownst to me before,
apparently was not in the correct format for this list. Hopefully this
attempt will go through, and no-one will hold the newbie mistake
against me.
I could really use some help in writing a new glm link function in
order to run an analysis of daily nest survival rates. I've struggled
with this for weeks now, and can at least
2008 Oct 15
0
R-help Digest, Vol 67, Issue 31
...ntrol = list(trace =
6), method = "L-BFGS-B", lower = c(-Inf, 0.05), upper = c(Inf, Inf), ...)
names(fit$par) <- dn
cat("\nCoefficients:\n"); print(fit$par)
cat("\nResidual Deviance:", format(fit$value), "\n")
if(fit$convergence > 0) cat("\nConvergence code:", fit$convergence, "\n")
return(list(fitpar = fit$par))
invisible(fit)
}
And here is my data,
# ----------------------
# Dose 0 1 emp.prop
# ----------------------
# 1 3 0 0.000
# 2 1 0 0.000
# 3.3 4 0 0.000
# 5.016 3 0 0.000
# 7.0224 4 0 0.000
# 9.3398 2 0 0.000...