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
