I suspect I am running into limits of what may be possible in Mark/RMark, but I just wanted to ask your advice in case I am missing something obvious.
I am studying a long-lived species (up to 43 years), and because of historic metal-ringing etc., in many instances we can identify the minimum age of an individual, which ranges from 0 to 39. From a previous post, ages can be automatically converted to an age group and init.age
- Code: Select all
data$MinKnownAge=as.factor(as.integer(data$MinKnownAge))
init.ages=as.integer(levels(data$MinKnownAge))
ms.pr=process.data(data = data, begin.time = c(1983), groups = "MinKnownAge",
initial.age = init.ages, age.var = 1,
model = "MSLiveDead", strata.labels = c("A", "D", "N","W"))
The design data took a long time to create, and no wonder why, the survival design data is 104,000 lines because it includes a line for every initial age, stratum and year (25 years, 1983 - 2018).Psi is 300,000 lines
Is there a way to include a broad array of initial ages without creating a group for every init.age? I have run simpler analyses with just four age groups, but I am investigating senescence and thus it would be great if I could include the true(ish) age of the individual. I later use add.design.data to create agebins, but I would prefer individuals to enter the agebin at the correct time.
Many thanks in advance.
Andy