Я хочу выполнить загрузку моих данных с помощью функции popbio boot.transitions.Тем не менее, это как-то выдает странную ошибку со столбцом фертильности:
stagefate = read.table("file.csv", header = TRUE, sep=",")
stagefatenames<- c("individual", "stage", "fate","fertility")
trans2.3 <- stagefate[stagefatenames]
boot2.3 <- boot.transitions(trans2.3, 100)
ci2.3 <- quantile(boot2.3$lambda, c(0.025, 0.975))
hist(boot2.3$lambda, col = "grey80", xlab = "Lambda", main = "")
abline(v = ci2.3, lty = 3)
dput(trans2.3)
50 + ошибки:
Missing a fertility column with individual fertility rates
Что я оба на самом деле не понимаю.Я добавил столбец «Плодородие», не знаю, как его не увидеть и как это исправить.
Кстати, моя таблица выглядит так:
individual stage fate fertility
1 1 infant infant 0
2 2 infant infant 0
3 3 infant juvenile 0
4 4 infant juvenile 0
5 5 infant juvenile 0
6 6 juvenile adult 0
7 7 adult adult 0
8 8 adult adult 1
9 9 adult adult 1
dput (trans2.3)
structure(list(individual = 1:9, stage = structure(c(2L, 2L,
2L, 2L, 2L, 3L, 1L, 1L, 1L), .Label = c("adult", "infant", "juvenile"
), class = "factor"), fate = structure(c(2L, 2L, 3L, 3L, 3L,
1L, 1L, 1L, 1L), .Label = c("adult", "infant", "juvenile"), class =
"factor"),
fertility = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L)), row.names = c(NA,
-9L), class = "data.frame")
edit: данные были не в правильном формате, и это исправило одну предыдущую ошибку (матрица сортировки ошибок).
Однако ошибка рождаемости все еще сохраняется.В любом случае, это дает некоторый результат, который я нахожу действительно странным и, вероятно, будет вызван ошибкой рождаемости.Доверительные интервалы после начальной загрузки для любых данных, которые я вставляю, всегда меньше 1, что странно, поскольку наши лямбды, которые мы вычисляем, никогда не бывают такими низкими.