RMark Jolly Seber POPAN Открытое население фиксированных значений

Question

Я хочу запустить несколько моделей POPAN, используя временные, одиночные и фиксированные параметры для моего параметра вероятности захвата (p).У меня есть 6 лет (2012-2017) и у меня есть уникальные фиксированные оценки p для каждого годаКак мне создать уникальные значения р для каждого года?Я могу сделать это в программе MARK, но не в пакете R.

run.Cap1 <- function() {
 Phi.time = list(formula =  ~ time)
Phi.dot = list(formula = ~ 1)
 p.time = list(formula =  ~ time)
p.fixed=list(formula=~time,fixed=list(time=c(1,2,3,4,5,6),value=c(0.0478,0.0478,0.0478,0.0478,0.0556,0.04)))
 pent.time = list(formula =  ~ time)
  N.dot = list(formula =  ~ 1)

  Phi.time.p.time.pent.time.N = mark(
    Cap1.pr,
    Cap1.ddl,
    model.parameters = list(
      Phi = Phi.time,
      p = p.time,
      pent = pent.time,
      N = N.dot))

  Phi.time.p.fixed.pent.time.N = mark(
    Cap1.pr,
    Cap1.ddl,
    model.parameters = list(
      Phi = Phi.time,
      p = p.fixed,
      pent = pent.time,
      N = N.dot))

  Phi.dot.p.time.pent.time.N = mark(
    Cap1.pr,
    Cap1.ddl,
    model.parameters = list(
      Phi = Phi.dot,
      p = p.time,
      pent = pent.time,
      N = N.dot))

 Phi.time.p.time.pent.time.N =
    adjust.parameter.count(Phi.time.p.time.pent.time.N, 16)
  return(collect.models())
}


Cap1.results <- run.Cap1()

options(width = 150)
Cap1.results