[R] pglm package: fitted values and residuals
pauljohn32 at gmail.com
Thu Apr 25 21:14:25 CEST 2013
On Wed, Apr 24, 2013 at 4:37 PM, Achim Zeileis <Achim.Zeileis at uibk.ac.at> wrote:
> On Wed, 24 Apr 2013, Paul Johnson wrote:
>> On Wed, Apr 24, 2013 at 3:11 AM, <alfonso.carfora at uniparthenope.it>
>>> I'm using the package pglm and I'have estimated a "random probit model".
>>> I need to save in a vector the fitted values and the residuals of the model
>>> but I can not do it.
>>> I tried with the command fitted.values using the following procedure
>>> without results:
>> This is one of those "ask the pglm authors" questions. You should take it
>> up with the authors of the package. There is a specialized email list
>> R-sig-mixed where you will find more people working on this exact same
>> pglm looks like fun to me, but it is not quite done, so far as I can tell.
> I'm sure that there are many. One of my attempts to write up a list is in
> Table 1 of vignette("betareg", package = "betareg").
Yes! That's exactly the list I was thinking of. It was driving me
crazy I could not find it.
Thanks for the explanation. I don't think I should have implied that
the pglm author must actually implement all the methods, it is
certainly acceptable to leverage the methods that exist. It just
happened that the ones I tested were not implemented by any of the
But this thread leads me to one question I've wondered about recently.
Suppose I run somebody's regression function and out comes an object.
Do we have a way to ask that object "what are all of the methods that
might apply to you?" Here's why I wondered. You've noticed that
predict.lm has the interval="confidence" argument, but predict.glm
does not. So if I receive a regression model, I'd like to say to it
"do you have a predict method" and if I could get that predict method,
I could check to see if there is a formal argument interval. If it
does not, maybe I'd craft one for them.
> Personally, I don't write anova() methods for my model objects because I can
> leverage lrtest() and waldtest() from "lmtest" and linearHypothesis() and
> deltaMethod() from "car" as long as certain standard methods are available,
> including coef(), vcov(), logLik(), etc.
> Similarly, an AIC() method is typically not needed as long as logLik() is
> available. And BIC() works if nobs() is available in addition.
>>> m1_S<-pglm(Feed ~ Cons_PC_1 + imp_gen_1 + LGDP_PC_1 + lnEI_1 +
>>> Can someone help me about it?
>> R-help at r-project.org mailing list
>> PLEASE do read the posting guide
>> and provide commented, minimal, self-contained, reproducible code.
Paul E. Johnson
Professor, Political Science Assoc. Director
1541 Lilac Lane, Room 504 Center for Research Methods
University of Kansas University of Kansas
More information about the R-help