similar to: aic() component in GLM-family objects

FWIW p. 206 of the White Book gives the following for names(binomial()): family, names, link, inverse, deriv, initialize, variance, deviance, weight. So $aic wasn't there In The Beginning. I haven't done any more archaeology to try to figure out when/by whom it was first introduced ... Section 6.3.3, on extending families, doesn't give any other relevant info. A patch for

Is it generally known/has it been previously discussed here that the $aic() component in GLM-family objects (e.g. results of binomial(), poisson(), etc.) does not as implemented actually return the AIC, but rather -2*log-likelihood + 2*(model_has_scale_parameter) ? Can anyone in this forum gauge how a documentation patch would be received? This behaviour does not seem to be documented in ?family

>>>>> Ben Bolker >>>>> on Sun, 3 Jun 2018 17:33:18 -0400 writes: > Is it generally known/has it been previously discussed here that the > $aic() component in GLM-family objects (e.g. results of binomial(), > poisson(), etc.) does not as implemented actually return the AIC, but > rather -2*log-likelihood + 2*(model_has_scale_parameter)

As of r87537, binomial()$linkinv no longer accepts integer arguments. binomial()$linkinv(1.0) ## 0.7310586 binomial()$linkinv(1L) Error in binomial()$linkinv(1L) : REAL() can only be applied to a 'numeric', not a 'integer' Since R is usually so permissive/sloppy with the distinction between integers and numeric/doubles, I'd argue that this is a bug ... (This came

On 1/8/25 17:26, Ben Bolker wrote: > Thanks, that makes sense. > > ? I guess if it never worked for integers (or hasn't worked in a long > time, at least) then it doesn't need to be fixed/changed ... Still you found out that the type check and use of REAL() in the function is in wrong order. Instead of "REAL() can only be applied to a 'numeric', not a

> On 9 Jan 2025, at 05:56, Tomas Kalibera <tomas.kalibera at gmail.com> wrote: > > > On 1/8/25 17:26, Ben Bolker wrote: >> Thanks, that makes sense. >> >> I guess if it never worked for integers (or hasn't worked in a long time, at least) then it doesn't need to be fixed/changed ... > > Still you found out that the type check and use of

Making data, especially data that have been the subject of published papers, widely available, can be a useful spinoff from the R project, another gift to the scientific community beyond the provision of computing and analytic tools. Nowadays, in a complete publication of a scientific result, there is every reason for the data to be part of that publication. The Gentleman and Lang 2004

Hi, I am a student in archaeology with some interest in statistics and R. Recently I've obtained the permission to distribute in the public domain a small dataset (named "spearheads" with 40 obs. of 14 variables) that was used in an introductory statistics book for archaeologists (published in 1994). I've rewritten most of the exercises of that book in R and made them

under "Details" (version 2015-11-29 r69717; not quite cutting-edge, but nothing has changed in src/library/stats/man/family.Rd in 5 months [sorry for using the Github mirror, but I prefer the interface ... <https://github.com/wch/r-source/blob/trunk/src/library/stats/man/family.Rd>]) it says: valid.mu: logical function. Returns ?TRUE? if a mean vector ?mu? is within the

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --27464147-1557463723-1147085467=:8118 Content-Type: TEXT/PLAIN; charset=iso-8859-1; format=flowed Content-Transfer-Encoding: 8BIT On Mon, 8 May 2006, x.sole at iconcologia.net wrote: > Full_Name: Xavier Sol? > Version: 2.3.0 > OS: Windows

Can anyone explain to me why the AIC values are so different when using glm.nb and glm with a negative.binomial family, from the MASS library? I'm using R 1.8.1 with Mac 0S 10.3.4. >library(MASS) > dfr <- data.frame(c=rnbinom(100,size=2,mu=rep(c(10,20,100,1000),rep(25,4))), + f=factor(rep(seq(1,4),rep(25,4)))) > AIC(nb1 <- glm.nb(c~f, data=dfr)) [1] 1047 >

Dear list I am relatively new to networking problems of this kind so apologies for the potentially simple question. I am trying to upgrade an existing network to one using Samba 3 to configure roaming XP profiles on a limited number of clients. I have re-written the smb.conf file to reflect what I think are the appropriate settings, and this passed testparm successfully, but I am unclear

On Wed, 8 Jan 2025 10:57:47 -0500 Ben Bolker <bbolker at gmail.com> wrote: > I haven't done the archaeology to figure out when this broke/exactly > what change in the R code base broke it: it happened within the last > month or so binomial() itself exhibits this property even in R-4.2.2 from more than two years ago: R -q -s -e 'getRversion(); binomial()$linkinv(1L)'

Thanks, that makes sense. I guess if it never worked for integers (or hasn't worked in a long time, at least) then it doesn't need to be fixed/changed ... cheers Ben On 2025-01-08 11:20 a.m., Ivan Krylov wrote: > On Wed, 8 Jan 2025 10:57:47 -0500 > Ben Bolker <bbolker at gmail.com> wrote: > >> I haven't done the archaeology to figure out when this

Dear all, I am trying to model number of samples from a given series. The series are modelled according Gamma function. In order to estimate the # samples, I use BIC/AIC with MLE (computed from dgamma function). Here is the code I have. __BEGIN__ mlogl <- function( x_func, theta_func, samp) { # computing log_likelihood return( - sum(dgamma(samp, shape = x_func, scale=theta_func, log

Hello List, I am modeling a binomial response (nest survival) and I want to incorporate a random effect, in this case site. I had previously been using glm with a custom link function, but my understanding is that lmer does not currently allow a custom link. Therefore, I was investigating if other procedures for mixed models will allow a custom link function. here is the custom link function:

Hi all, I've just installed 1.9.18p5 (although I've just got the mail about 1.9.18p6, so I'll be putting that on now) and am having a problem with browsing... [samba@suau1] ~> smbclient -L suau1 Added interface ip=134.225.16.5 bcast=134.225.19.255 nmask=255.255.252.0 Server time is Mon May 11 15:24:39 1998 Timezone is UTC+1.0 Password: Domain=[UNSEEN] OS=[Unix]

John, You have specified a model with E(y) = 1/eta where eta = X beta is the linear predictor and E(y) must be >0, since the family is Gamma and you have a lot of covariates in the model. glm now has to try to find a best linear predictor, but under the constraint that eta>0 for every single one of the observations (the log-likelihood involves a log(eta) term). The internal

I'm not sure who is responsible for quantile(), but I assume they read this list. Ivan Frohne and I have produced a revision of the quantile.default() function which enables the computation of alternative sample quantile definitions. The code and .Rd file are attached. This enables the user to produce quantiles that are equivalent to those in various statistics package. There is a type

I'm seeing various things fail when I try to next braces inside \eqn. This source \eqn{{\bf\beta}_j}{b(j)} is the vector produces this error ---------------------------------------------- [4] ! Missing $ inserted. <inserted text> $ l.258 \eqn{{\bf\beta}_j}{\bf\beta}_ j{{b(j)} is the vector of coefficients fo... I've inserted a