[R] Draw a random number about the fitting of a glm model using family=quasi

Hi, I desire to draw a random number about the fitting of a glm model using a quasi linear error distribution. For example, if my family is poisson, then this is what I want: model1 = glm(amount~days,family=poisson(link="log"),data=dataset1) mean1 = predict(model1,newdata =data.frame(days=time),

Re: [R] extracting standard errors in glm.nb

Mea culpa. Thanks! jpl David Winsemius wrote: > > Look more thoroughly? > > You should see this about 2/3 the way down the output: > > $ coefficients : num [1:14, 1:4] 3.019 -0.475 -0.709 -0.724 -0.615 ... >..- attr(*, "dimnames")=List of 2 >

Re: [R] extracting standard errors in glm.nb

Hi Dieter, I don't see anything in the str(model) output that matches the standard errors of the coefficients given in summary(model). Any other suggestion? Thanks. jpl Dieter Menne wrote: > > jpl math.unl.edu> writes: > >> >> How do I extract the standa

[R] extracting standard errors in glm.nb

Hi, How do I extract the standard errors of the coefficients in a glm.nb model? Thanks. -- View this message in context: http://www.nabble.com/extracting-standard-errors-in-glm.nb-tp22019202p22019202.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] Extracting standard error from survreg?

Thanks! This is exactly what I needed. Peter Dalgaard wrote: > > jpl wrote: >> Hello, >> >> I would like to extract the standard error for the coefficients returned >> when using survreg. Does anyone know how to do this? >> >> Thank you. >>

[R] Extracting standard error from survreg?

Hello, I would like to extract the standard error for the coefficients returned when using survreg. Does anyone know how to do this? Thank you. Joan -- View this message in context: http://www.nabble.com/Extracting-standard-error-from-survreg--tp19994742p19994742.html Sent from the R help ma

Re: [R] using nls to fit a curve to data

Thanks Ben! Switching over to the gamma pdf and using the algorithm="plinear" did the trick. jpl -- View this message in context: http://www.nabble.com/using-nls-to-fit-a-curve-to-data-tp19332210p19360761.html Sent from the R help mailing list archive at

[R] using nls to fit a curve to data

Hi, I am trying to fit a curve to data. My command line is: model10=nls(offspring~((A*c^k)/gamma(k))*((degdays-alpha)^(k-1))*exp(-c*(degdays-alpha)), start=list(A=30,k=2,c=.018,alpha=131)) I get the error message: Error in numericDeriv(form[[3]], names(ind), env) : Missing val

[R] question on dispersion parameter

Hi, I'm programming in R and below is a summary of a generalized linear model: ** *** Call: glm(formula = offspring ~ degdays, family = quasi(link = "log", variance = "mu"), data = fecundity) Deviance Residuals: Min 1Q Median 3Q Max -0.76674 -0.291