У меня есть фрейм данных присутствия / отсутствия, упрощенно здесь:
A <- c(1,1,0,0)
B <- c(1,0,1,0)
C <- c(1,1,0,0)
D <- c(1,1,0,1)
E <- c(1,0,1,0)
df <- cbind(A,B,C,D,E)
A B C D E
[1,] 1 1 1 1 1
[2,] 1 0 1 1 0
[3,] 0 1 0 0 1
[4,] 0 0 0 1 0
Я хочу вычислить изоклину на основе уравнения, которое включает в себя количество столбцов и строк. Уравнение для изоклины основано на опубликованной литературе (Rodríguez-Gironés & Santamaría 2006) и имеет вид:
y <- ( 0.5 / n ) + ( (n-1) / n) * (1 - (1 - ((m*x - 0.5) / (m - 1) ))^p )^(1-p)
, где
m = количество столбцов
n = количество строк
p = диапазон значений, подходящих для изоклины
Я начал с нескольких значений для p . Значение p является диапазоном и должно быть таким, чтобы число отсутствий (0) над изоклиной было равно числу присутствий (1) ниже изоклины.
Как проанализировать диапазон значений p с помощью этого уравнения?
Как мне кодировать так, чтобы я решал для значения p , которое удовлетворяет условиям равных чисел присутствия ниже и отсутствия над изоклиной?
Этот расчет является частью более крупного процесса моделирования, однако я застрял как в математическом процессе, так и в кодировании для производства изоклина.