similar to: Extract values from a predict() result... how?

... 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

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 =

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?

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,

Hello, I am trying to figure out how to plot the profile likelihood curve of a GLM parameter with 95% pCI's on the same plot. The example I have been trying with is below. The plots I am getting are not the likelihood curves that I was expecting. The y-axis of the plots is tau and I would like that axis to be the likelihood so that I have a curve that maxes at the parameter estimate. I am

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

> 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

Dear expRts, when I try to build a package by myself, the process fails writing the following messages to 00install.out: * Installing *source* package 'dynalc' ... ** libs WARNING: no source files found chmod: Zugriff auf ?/home/hannes/texte/forschung/clot/programme/dynalc.Rcheck/dynalc/libs/*? nicht m?glich: No such file or directory ** R ** data ** inst ** help Note: removing empty

On 10 September 2010 14:37, Marcus Daniels <mdaniels at lanl.gov> wrote: >  Hello, > > With GCC, it is possible to compile GIMPLE from an object file back to > assembly.  With Dragonegg, it seems to be a NOP.  (See below for a > comparison.) > > This appears it could be a nice way to mix the strengths of GCC and > LLVM.  If it worked..  Should it? I can't

Hello, With GCC, it is possible to compile GIMPLE from an object file back to assembly. With Dragonegg, it seems to be a NOP. (See below for a comparison.) This appears it could be a nice way to mix the strengths of GCC and LLVM. If it worked.. Should it? Thanks, Marcus [mdaniels at dn002 dragonegg]$ cat hello.c #include <stdio.h> int hello () { printf ("Hello

Full_Name: José Otero Version: Version 1.5.0 (2002-04-29) OS: Redhat 7.3 Submission from: (NULL) (192.187.16.164) Hi: Installation of package faraway as root, from tarbal: R CMD INSTALL ./faraway.tar.gz ERROR: cannot extract package from './faraway.tar.gz' idem, from zipped package: R CMD INSTALL faraway.zip gzip: faraway.zip has more than one entry--rest ignored ERROR: cannot

I am trying to plot the results of a random effects model where each subject is allowed to have his, or her, own intercept. I use xyplot to plot the data, lme to compute the regression and then try to put a summary regression line on the xyplot. As can be seen by the output pasted below, I am getting an error message, Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : (list)

Hi list members I have a vector of strings x=c("w","ex","ee") And I want to get a logical vector showing the positions where my search string "e" matches the elements partially, i.e. is at least the left-hand part of the target strings, i.e. I want to get a vector FALSE TRUE TRUE. Any ideas? Thanks Steve Powell proMENTE social research research |

Anyone conceive of any reason why I should have a stream that RealPlayer and iTunes have no problem with, but MP9 refuses to play: just sits there dumbly having loaded the .m3u that points to it, ignoring the play button?

I am experiencing a curious event, and wondering if others have seen this... As well, I have a question related to it. Today, I noticed my puppet summary report from Foreman this morning, that 60 of my 160 hosts all stopped reporting at nearly the exact same time, and have not since restarted. Investigating, it appears that my puppetmaster temporarily ran out of disk space on the /var