www.phidot.org

Hi everyone,
I'd like to test the influence of the body size on survival of bats.
I wanted to start with two very simple models, one containing body size and one doesn't.

I imported my Dataset, which looks like this ("cov" is the variable for body size and "col" is the information about the colony, the individual bat belongs to, but I would ignore "col" until I can solve my problem with the cov variable)

Code: Select all

> Mnat.size[1:10,]
                                ch     ID  cov col
1  0111111111100000000000000000000 017172 41.4  wt
2  0001000000000000000000000000000 020870 40.4  wt
3  0001000000000000000000000000000 020896 39.3  wt
4  0000110100000000000000000000000 B01602 38.8  wt
5  0000111110100000000000000000000 B01605 39.8  wt
6  0000110000000000000000000000000 B01607 40.2  wt
7  0000111110000000000000000000000 B01611 38.2  wt
8  0000110100000000000000000000000 B01613 40.3  wt
9  0000111110000000000000000000000 B01614 38.4  wt
10 0000111011111110000000000000000 B01615 41.4  wt

processed them, created the design data and two simple models

Code: Select all

Mnat.size$cov<-as.character(Mnat.size$cov)
Mnat.size$cov<-as.numeric(Mnat.size$cov)

Mnat.size.processed=process.data(Mnat.size,model="CJS",begin.time=1989)
Mnat.size.ddl=make.design.data(Mnat.size.processed)

Phi.dot=list(formula=~1)
p.dot=list(formula=~1)
Phi.cov=list(formula=~cov)


Phi.dot.p.dot=make.mark.model(Mnat.size.processed,Mnat.size.ddl,parameters=list(Phi=Phi.dot,p=p.dot))
Phi.cov.p.dot=make.mark.model(Mnat.size.processed,Mnat.size.ddl,parameters=list(Phi=Phi.cov,p=p.dot))

Until now, I got no warning or error messages. In a next step I wanted to compare these models and got this error

Code: Select all

results<-collect.models()
Model 1: Phi.cov.p.dot failed to run.

Model 2: Phi.dot.p.dot failed to run.

Error in result.table[order(result.table$DeltaAICc), ] :
  incorrect number of dimensions
In addition: Warning message:
In min(result.table$AICc) : no non-missing arguments to min; returning Inf

What does this mean and how do I solve the problem?
Thanks a lot!
Bianca

You have only made the models and have not run them. Use mark in place of make.mark.model and change parameters argument to model.parameters.

ok, yes. sometimes it's so easy. I should not work out new things when I'm tired
Thank you!!!!