Displaying 3 results from an estimated 3 matches for "logit_linkinv".
2025 Jan 08
1
binomial()$linkinv no longer accepts integer values
...were of interest I could 'git
bisect' to figure it out but maybe someone will save me the trouble ...)
Thoughts, insights? Should I post this to r-bugzilla? Or suck it up
and accept it as the new reality?
cheers
Ben Bolker
====
binomial()$linkinv is:
function (eta)
.Call(C_logit_linkinv, eta)
<environment: namespace:stats>
Here is the body of the function in C code:
https://github.com/r-devel/r-svn/blob/195ab0870a8131d01492f8f1d3a2ad514bc7e040/src/library/stats/src/family.c#L72C1-L89C2
SEXP logit_linkinv(SEXP eta)
{
SEXP ans = PROTECT(shallow_duplicate(eta));
in...
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