Я использую RMark для определения численности грызунов по четырем типам растительности. У меня было 26 сайтов с 2 ночами ловушек на каждом сайте. Каждая ночь состояла из 50 ловушек, поэтому 100 ночей ловушки за очко. У меня было всего 3 повторных поимки. Возможно ли провести какой-либо анализ для определения вероятности и численности на разных участках? Когда я попытался запустить модель для «закрытой популяции», на выходе было нулевое значение. Я ценю любые комментарии или рекомендации. Спасибо!
> str(rodentmark2)
'data.frame': 66 obs. of 7 variables:
$ Site : Factor w/ 14 levels "A","B","C","E",..: 1 1 1 2 3 3 3 3 4 4 ...
$ Date : Factor w/ 24 levels "10/12/2019","10/13/2019",..: 12 13 13 14 15 15 16 16 17 17 ...
$ animal_id: int 1 2 3 4 5 6 7 8 9 10 ...
$ ch : chr "10" "01" "01" "10" ...
$ night_1 : int 1 0 0 1 1 1 0 0 1 1 ...
$ night_2 : int 0 1 1 0 0 0 1 1 0 0 ...
$ detect : int 1 1 1 1 1 1 1 1 1 1 ...
> head(rodentmark2)
Site Date animal_id ch night_1 night_2 detect
1 A 5/21/2019 1 10 1 0 1
2 A 5/22/2019 2 01 0 1 1
3 A 5/22/2019 3 01 0 1 1
4 B 5/25/2019 4 10 1 0 1
5 C 5/29/2019 5 10 1 0 1
6 C 5/29/2019 6 10 1 0 1
> dput(rodentmark2)
structure(list(Site = structure(c(1L, 1L, 1L, 2L, 3L, 3L, 3L,
3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 6L, 7L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 11L, 11L, 12L, 12L, 12L, 13L,
13L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L,
14L), .Label = c("A", "B", "C", "E", "F", "K", "O", "Q", "R",
"T", "V", "X", "Y", "Z"), class = "factor"), Date = structure(c(12L,
13L, 13L, 14L, 15L, 15L, 16L, 16L, 17L, 17L, 17L, 17L, 18L, 18L,
18L, 19L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 20L, 20L, 20L,
20L, 20L, 21L, 22L, 23L, 24L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L,
3L, 4L, 4L, 4L, 4L, 5L, 5L, 6L, 7L, 7L, 8L, 9L, 10L, 10L, 10L,
10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L, 11L), .Label = c("10/12/2019",
"10/13/2019", "10/20/2019", "10/21/2019", "11/13/2019", "11/22/2019",
"11/23/2019", "11/25/2019", "11/26/2019", "12/2/2019", "12/3/2019",
"5/21/2019", "5/22/2019", "5/25/2019", "5/29/2019", "5/30/2019",
"6/3/2019", "6/4/2019", "6/6/2019", "6/7/2019", "8/15/2019",
"9/10/2019", "9/24/2019", "9/25/2019"), class = "factor"), animal_id = 1:66,
ch = c("10", "01", "01", "10", "10", "10", "01", "01", "10",
"10", "10", "10", "01", "01", "01", "10", "10", "10", "10",
"10", "01", "01", "01", "01", "01", "01", "01", "01", "01",
"10", "10", "10", "01", "10", "10", "11", "01", "01", "01",
"10", "11", "10", "01", "01", "01", "01", "10", "10", "11",
"01", "01", "10", "01", "10", "10", "10", "10", "10", "10",
"01", "01", "01", "01", "01", "01", "01"), night_1 = c(1L,
0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L,
1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L,
1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L,
1L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L,
0L, 0L, 0L, 0L, 0L), night_2 = c(0L, 1L, 1L, 0L, 0L, 0L,
1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L,
1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L,
0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
), detect = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)), row.names = c(NA,
-66L), class = "data.frame")
rodentmark2$dht$individuals$summary
rodentmark2$dht$Site$summary