R help - Oct 2007 - estfun & df

Hello EVERYONE, 
 
I need an URGENT help from you please! 

How can I see the "estfun" (empirical estimating function) and
"df" (degree of freedom) from the following mixed-model please?
 
(fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) 
 
Many thanks in advance for your kind help. 
 
Sattar


       
____________________________________________________________________________________

ews, photos & more. 

	[[alternative HTML version deleted]]

On 10/8/07, Abdus Sattar <gsrwork at yahoo.com>
wrote:
> Hello EVERYONE,
> I need an URGENT help from you please!
> How can I see the "estfun" (empirical estimating function) and
"df" (degree of freedom) from the following mixed-model please?
> (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy))
You have asked this question several times on this list and on the
R-SIG-Mixed-Models list.  Repeating the question ad nauseum is
unlikely to generate a helpful answer.

To the best of my knowledge there is no lmer or lmer2 method for the
estfun generic.  Others may be willing to give you some references on
the difficulty in determining degrees of freedom issue - I don't want
to rehash it myself.

So, despite your repeated requests, there are no good answers to your
questions.  You cannot "see" such results because there are, as far as
I know, no functions to provide such results.  You are welcome to
write your own code to provide them.  If that is not suitable you may
need to take advantage of our money-back guarantee and remove R from
your computer to qualify for a full refund of the purchase price.

On Mon, 8 Oct 2007, Abdus Sattar wrote:
> Hello EVERYONE,
>
> I need an URGENT help from you please!
This type of requests is not considered to be very polite, please have a 
look at the posting guide.
> How can I see the "estfun" (empirical estimating function)
I guess (because you are not telling us) that you would like to have an 
estfun() method for "lmer" objects. I don't provide one in my
sandwich
package (where the estfun generic is taken from) and AFAIK there is 
nothing analagous readily available in "lme4"... But I guess that you 
should be able to extract/compute the empirical estimating functions from 
the fitted "lmer" object.
> and "df" (degree of freedom) from the following mixed-model
please?
>
> (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy))
When you compute
   logLik(fm1)
it returns a "logLik" object that has a "df" attribute:
   attr(logLik(fm1), "df")

hth,
Z