Displaying 2 results from an estimated 2 matches for "logit_linkinv".
2012 Aug 10
0
error applying user-defined link function to lmer
...# .m goes global
if (.m < 2)
stop(".m must be an integer > 1")
linkfun <- function(mu) {
mu <- pmax(mu, 1/.m + .Machine$double.eps)
qlogis((.m * mu - 1)/(.m - 1))
}
linkinv <- function(eta) {
1/.m + (.m - 1)/.m * .Call("logit_linkinv", eta, PACKAGE = "stats")
}
mu.eta <- function(eta) ((.m - 1)/.m) * .Call("logit_mu_eta",
eta, PACKAGE = "stats")
valideta <- function(eta) TRUE
link <- paste("mafc.logit(", .m, ")", sep = "")
stru...
2009 Jun 30
1
fitting in logistic model
I would like to know how R computes the probability of an event in a
logistic model (P(y=1)) from the score s, linear combination of x and
beta.
I noticed that there are differences (small, less than e-16) between the
fitting values automatically computed in the glm procedure by R, and the
values "manually" computed by me applying the reverse formula
p=e^s/(1+e^s); moreover I noticed