similar to: plotting profile likelihood curves

Hi, I have below GLM fit clotting <- data.frame( u = c(5,10,15,20,30,40,60,80,100), lot1 = c(118,58,42,35,27,25,21,19,18), lot2 = c(69,35,26,21,18,16,13,12,12)) summary(glm(lot1 ~ log(u), data = clotting, family = gaussian)) Is there any direct function to extract estimate of Error standard deviation?

... but do note: glm(lot1 ~ log(u), data = clotting, family = gaussian) is a plain old *linear model*, which is of course a specific type of glm, but not one that requires the machinery of glm() to fit. That is, the above is exactly the same as: lm(lot1 ~ log(u), data = clotting) and gives exactly the same sigma() ! (and I would therefore hazard the guess that the poster may misunderstand

Hello, silly question I suppose, but somehow I can't manage to extract the probabilities from a glm.predict() result: > str(res) Named num [1:9] 0.00814 0.01877 0.025 0.02941 0.03563 ... - attr(*, "names")= chr [1:9] "1" "2" "3" "4" ... I got from: # A Gamma example, from McCullagh & Nelder (1989, pp. 300-2) clotting <-

Dear all, I'm trying to use the "optim" function to replicate the results from the "glm" using an example from the help page of "glm", but I could not get the "optim" function to work. Would you please point out where I did wrong? Thanks a lot. The following is the code: # Step 1: fit the glm clotting <- data.frame( u =

vcov(). ? On Tue, Dec 24, 2024, 8:45 AM Christofer Bogaso <bogaso.christofer at gmail.com> wrote: > Hi, > > I have below GLM fit > > clotting <- data.frame( > u = c(5,10,15,20,30,40,60,80,100), > lot1 = c(118,58,42,35,27,25,21,19,18), > lot2 = c(69,35,26,21,18,16,13,12,12)) > summary(glm(lot1 ~ log(u), data = clotting, family = gaussian)) > >

I think vcov() gives estimates of VCV for coefficients. I want estimate of SD for residuals On Tue, Dec 24, 2024 at 7:24?PM Ben Bolker <bbolker at gmail.com> wrote: > > vcov(). ? > > > On Tue, Dec 24, 2024, 8:45 AM Christofer Bogaso <bogaso.christofer at gmail.com> wrote: >> >> Hi, >> >> I have below GLM fit >> >> clotting <-

?deviance ?anova Bert On Tue, Dec 24, 2024 at 6:22?AM Christofer Bogaso <bogaso.christofer at gmail.com> wrote: > > I think vcov() gives estimates of VCV for coefficients. > > I want estimate of SD for residuals > > On Tue, Dec 24, 2024 at 7:24?PM Ben Bolker <bbolker at gmail.com> wrote: > > > > vcov(). ? > > > > > > On Tue, Dec 24,

> d = data.frame(gender=rep(c('f','m'), 5), pos=rep(c('worker', 'manager', 'speaker', 'sales', 'investor'), 2), lot1=rnorm(10), lot2=rnorm(10)) > d gender pos lot1 lot2 1 f worker 1.1035316 0.8710510 2 m manager -0.4824027 -0.2595865 3 f speaker 0.8933589 -0.5966119 4 m sales

Hi R users, I want to use dglm Package. I run the examples and it give me an error: Error en dglm(lot1 ~ log(u), ~1, data = clotting, family = Gamma) : no se pudo encontrar la funci?n "as.family" dglm can't find "as.family" function why ? Thank you for your help

Hi, I have a data set (weight) that does not follow the Gaussian (Normal) distribution. However, I have to transform the data before applying the Gaussian distribution. I used this syntax and used log(weight) as: posJy.model<-glm(log(weight) ~ factor(pos), family=gaussian(link='identity'), subset=Soil=="Jy"). This syntax COULD NOT transform the data. But if I transform the

Hello: I'd like to know how and if the GLM convergence problems are addressed in gamlss. For simplicity, let's focus on Normal and Negative Binomial with log link. The convergence issues of the glm() function were alleviated in 2011 when glm2 package was released. Package gamlss was released in 2012, so it might still use the glm-like solution or call glm() directly. Is that the case or

I am attempting to run a glm with a binomial model to analyze proportion data. I have been following Crawley's book closely and am wondering if there is an accepted standard for how much is too much overdispersion? (e.g. change in AIC has an accepted standard of 2). In the example, he fits several models, binomial and quasibinomial and then accepts the quasibinomial. The output for residual

This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c("ZB:X","ZC:X") # and test.terms = ~ ZB:X + ZC:X in

Hi! I would like to perform an F-Test over more than one variable within a generalized mixed model with Gamma-distribution and log-link function. For this purpose, I use the package mgcv. Similar tests may be done using the function "anova", as for example in the case of a normal distributed response. However, if I do so, the error message "error in eval(expr, envir, enclos) :

Hi all I'm getting a NaN returned on using dffits, as explained below. To me, there seems no obvious (or non-obvious reason for that matter) reason why a NaN appears. Before I start digging further, can anyone see why dffits might be failing? Is there a problem with the data? Consider: # Load data dep <-

Ted Harding said: > I can get the estimated RRs from > RRs <- exp(summary(GLM)$coef[,1]) > but do not see how to implement confidence intervals based > on "robust error variances" using the output in GLM. Thanks for the link to the data. Here's my best guess. If you use the following approach, with the HC0 type of robust standard errors in the

Hi i would like to use some graphs or tables to explore the data and make some sensible guesses of what to expect to see in a glm model to assess if toxin concentration and sex have a relationship with the kill rate of rats. But i cant seem to work it out as i have two predictor variables~help?Thanks.:) Here's my data. >

Hi, Suppose I have a vector: > names.select [1] "Idd13" "Idd14" "Idd8.12" "Idd7" automatically generated by some selection criteria. Now, if I have a data frame with many variables, of which the variables in "names.select" are also variables from the data frame. e.g. > all.df[1:5,] Mouse Idd5 Idd6.19.20 Idd13 Idd14 Idd8.12

Dear R-users When running a glm polynomial model with one explanatory variable (example Y~X+X^2), with a poisson or binomial error distribution, the predicted values obtained from using the predict() function and those obtained from using the coefficients from the summary table "as is" in an equation of the form Y=INTERCEPT+ XCoef x X + XCoef x X^2, differ considerably. The former are

Problems with predict and lines in plotting binomial glm Dear R-helpers I have found quite a lot of tips on how to work with glm through this mailing list, but still have a problem that I can't solve. I have got a data set of which the x-variable is count data and the y-variable is proportional data, and I want to know what the relationship between the variables are. The data was