I'm running a script to use Rmark package to determine under Pradel models, the survival, recapture and population growth rates. I already use the exact lines in other two dataframes and the script worked correctly, however, with this specific data, its not. I specified all the models I want, the covariates, and groups for the models, however, it appear this message in one of the output files:
`INPUT --- proc stop;
ERROR -- Numerical convergence never reached.
And in the Rstudio console, appear this message:
"MARK did not run properly. If error message was not shown, look at most recent .out file
Model output is not available in R. Check contents of mark035.out
Error in `colnames<-`(`*tmp*`, value = model$design.data[[parameter]]$time[model$pims[[parameter]] [[1]]$pim - :
attempt to set 'colnames' on an object with less than two dimensions"`
Im using this code:
`data <- convert.inp("Mark temporadas.inp", covariates=c("sex", "season"))
data$sex <-factor(data$sex, levels = c(0,1,2), labels = c("males", "females","juven"))
data$season <- factor(data$season, levels = c(0,1,2,3,4), labels = c("1hum21", "1seca","2hum",
"2seca21", "1seca21"))
summary(data)
data.process=process.data(data,model="Pradlambda", groups = c("sex", "season"))
data.dd1=make.design.data(data.process)
run.data <- function()
{
Phi.dot = list(formula = ~ 1)
Phi.season = list(formula = ~ season)
Phi.sex = list(formula = ~ sex)
Phi.time = list(formula = ~ time)
Phi.time.x.season = list(formula = ~ time * season)
Phi.time.x.sex = list(formula = ~ time * sex)
Phi.sex.x.season = list(formula = ~ sex * season)
p.dot = list(formula = ~1)
p.season = list(formula = ~ season)
p.sex = list(formula = ~ sex)
p.time = list(formula = ~ time)
p.time.x.season = list(formula = ~ time * season)
P.time.x.sex = list(formula = ~ time * sex)
P.sex.x.season = list(formula = ~ sex * season)
lambda.dot = list(formula = ~1)
lambda.season = list(formula~ season)
lambda.sex = list(formula ~ sex)
lambda.time = list(formula ~ time)
lambda.time.x.season = list(formula = ~ time * season)
lambda.time.x.sex = list(formula = ~ time * sex)
lambda.sex.x.season = list(formula = ~ sex * season)
model.list = create.model.list("Pradlambda")
model.results = mark.wrapper(model.list, data = data.process, ddl = data.dd1)
return(model.results)}
multigroup.results = run.data()`
Everything goes well until i click on the last line...
I compared the data matrix with the other two matrices I already mentioned which worked well and everything its ok, Idon't understand what is the problem, someone can help me, please??
I really appreciate your help
