search for: linktemp

...o find any other published examples using this variance function. I'd really like to verify that my code above is working before applying it to real data sets. Can anybody help? Thanks, Henric - - - - - quasi2 <- function (link = "identity", variance = "constant") { linktemp <- substitute(link) if (is.expression(linktemp) || is.call(linktemp)) linktemp <- link else if (!is.character(linktemp)) linktemp <- deparse(linktemp) if (is.character(linktemp)) stats <- make.link(linktemp) else stats <- linktemp v...

I am resubmitting this to r-bugs, since Thomas Lumley indicates that it might be an error: On Wed, 5 Jul 2000, Thomas Lumley wrote: > On Wed, 5 Jul 2000, halvorsen wrote: > > > Hola! > > > > I have the following simple function: > > > > > testcar > > function(pow){ > > ob <- glm(Pound~CG+Age+Vage,data=car,weights=No, > >

Dear R users I have been looking for functions that can deal with overdispersed poisson models. Some (one) of the observations are negative. According to actuarial literature (England & Verall, Stochastic Claims Reserving in General Insurance , Institute of Actiuaries 2002) this can be handled through the use of quasi likelihoods instead of normal likelihoods. The presence of negatives is not

This is not so much a bug as an infelicity in the code that can easily be fixed. The initialize expression in the quasi family function is, (uniformly for all links and all variance functions): initialize <- expression({ n <- rep.int(1, nobs) mustart <- y + 0.1 * (y == 0) }) This is inappropriate (and often fails) for variance function "mu(1-mu)".

Hi, I would like apply different link functions from Aranda-Ordaz (1981) family to large binary dataset (n = 2000). The existing links in glm for binomial data (logit, probit, cloglog) are not adequate for my data, and I need to test some other transformations. Is it possible to do this in R? And how? Thank you for your help, /Sharon