[R] AICcmodavg

[David Winsemius]

On Sep 13, 2012, at 9:00 AM, David Winsemius wrote:

> 
> On Sep 13, 2012, at 5:52 AM, KerryGK wrote:
> 
>> I am using the AICmodavg package and using R version 2.15.1.  I need help
>> with code that is instead being read as text.  
>> 
>> Below is a subset of code... I actually have 12 models, but I am trying to
>> get this to work for 2 below right now.
>> 
>> Everything 'appears' to work through the line starting with Modnames.  After
>> that the code starting with aictab and beyond is recognized as text and not
>> code.  What am I doing incorrectly with the code?  
>> 
>> buco.models<-list ( )
>> buco.models [[1]] <- glm(P.BUCO~faverage+W15CNT, family=binomial(logit),
>> data=AmphAIC)
>> buco.models [[2]] <- glm(P.BUCO~diam1+W15CNT, family=binomial(logit),
>> data=AmphAIC)
>> Modnames<-paste (“mod” , 1:length(buco.models), sep=")
> ------------------------------------------------------^
> Not at all surprising. The usual response of a programmer faced with code that goes wrong at a particular point should be to _first_ look for missing commas or parentheses. In your case it is a missing comma.

Skitt's Law, again. comma != double-quote.
> 
>> aictab(cand.set = buco.models, modnames = Modnames, sort=TRUE))
>> print(aictab(cand.set = buco.models, modnames = Modnames, sort=TRUE) digits
>> = 4)
>> 
> 
-- 

David Winsemius, MD
Alameda, CA, USA