Hi Kamil,

(I have replied to the forum but my posts as replies never seem to get
accepted, so here is the text in an email also. Apologies if it comes
through twice.)

Thanks so much for the prompt reply. You've solved it! Actually I had the
latest version of the package, but an older version of R that wasn't
supporting the package evidently... when I updated it everything ran with
the "AICc" rank. However, when I run it with "QAICc" (which i think Iwant
to use), R returns the error "Error in chat < 1 : 'chat' is missing" and
does not conduct the dredge. My chat (variance inflation factor) is 0.597.
As this is <1, should I just rank by AICc? I still have a small sample size
necessitating AICc.

For interest, I can make the dredge run with rank="QAICc", but I need to
include the chat, and a error message then tells me "In rank(x, chat =
0.597604265592854) : 'chat' given is < 1, increased to 1". The script for
this is below:

##Information Theoretic Approach
globalmodel<- glm(TB~lat+protocol+tested+streams+goats+hay+cattle+deer,
family="binomial")
chat<- deviance(globalmodel)/59 #There we 59 residual degrees of freedom in
this model.
models<- dredge(globalmodel, beta=FALSE, evaluate=TRUE, rank="QAICc",
chat=chat, fixed=NULL, trace=FALSE)
models
warnings()

I think I will use just AICc, but do let me know if you think this is
wrong! Many thanks again for your help,
Cat


On 1 April 2013 22:19, Kamil Barton [via R] <
ml-node+s789695n4662958h84@n4.nabble.com> wrote:

> Hi Cat,
> are you using some very old version of MuMIn? That would explain the
> missing 'QAICc'.
> As for the error message about 'logLik', it usually occurs when there
> are some misspelled arguments (that go into "..." and are passed to the
> rank function, 'AICc' in your case). Check if there is some argument (of
> type "logical") in your call to 'dredge' that is not its formal argument.
>
> kamil
>
> On 2013-03-30 03:13, [hidden email]<http://user/SendEmail.jtp?type=node&node=4662958&i=0>wrote:
>
> > Hi Kamil,
> >
> > Thanks for your help. I do want to use rank="QAICc", but I when I try it
> I get the following error message:
> >
> > Error in get(as.character(FUN), mode = "function", envir = envir) :
> >    object 'QAICc' of mode 'function' was not found
> >
> > Do I need to install another package that would recognise this?
> >
> > When I leave the script with rank="AICc" (to see if it will run), and
> having removed "chat" (which I put in after previously recieving an error
> message about it being missing from the formula), I still get my original
> error message:
> >
> > Error in UseMethod("logLik") :
> >    no applicable method for 'logLik' applied to an object of class
> "logical"
> >
> > Any further thoughts?
> > Thanks, Cat
> >
> >
> > <quote author='Kamil Barton'>
> > 'rank' should be "QAICc".
> > AICc does not have argument 'chat', hence the error.
> >
> > kamil
> >
> >
> >
> > Hi all,
> >
> > I'm having trouble with the model generating 'dredge' function in the
> MuMIn
> > 'Multi-model Inference' package.
> >
> > Here's the script:
> >
> > globalmodel<- glm(TB~lat+protocol+tested+
> > streams+goats+hay+cattle+deer,
> > family="binomial")
> > chat<- deviance(globalmodel)/59 #There we 59 residual degrees of freedom
> in
> > this global model.
> > models<- dredge(globalmodel, beta=FALSE, evaluate=TRUE, rank="AICc",
> > chat=chat, fixed=NULL, trace=FALSE)
> >
> > And the error message is:
> >
> > Error in UseMethod("logLik") :
> >    no applicable method for 'logLik' applied to an object of class
> "logical"
> >
> > I have trawled the literature and it seems to be ok to use a binary GLM
> as
> > the global model - could this be the problem? The variables are a mix of
> > binary and continuous data.
> >
> > Any thoughts?
> > Thanks,
> > Cat
> >
> >
> > </quote>
> > Quoted from:
> >
> http://r.789695.n4.nabble.com/Error-message-in-dredge-function-MuMIn-package-used-with-binary-GLM-tp4662842p4662881.html
> >
>
>
>
> The University of Aberdeen is a charity registered in Scotland, No
> SC013683.
>
> ______________________________________________
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4662958&i=1>mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
>
> http://r.789695.n4.nabble.com/Error-message-in-dredge-function-MuMIn-package-used-with-binary-GLM-tp4662842p4662958.html
>  To unsubscribe from Error message in dredge function (MuMIn package) used
> with binary GLM, click here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4662842&code=Y2F0LmUuY293aWVAZ21haWwuY29tfDQ2NjI4NDJ8MTMyMDc2OTMxNw==>
> .
> NAML<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>




--
View this message in context: http://r.789695.n4.nabble.com/Error-message-in-dredge-function-MuMIn-package-used-with-binary-GLM-tp4662842p4663033.html
Sent from the R help mailing list archive at Nabble.com.
	[[alternative HTML version deleted]]

