В настоящее время я изучаю SDM (моделирование распространения видов), и я знаком с некоторыми из используемых там кодировок. Я сам изучал R, есть некоторые коды, которые я просто не знаю, что они имеют в виду. Может кто-нибудь сказать мне, что означает этот код.
Здесь они иллюстрируют вычисление коэффициента корреляции и AU C (площадь под кривой оператора приемника) с двумя случайными переменными. p (присутствие) имеет более высокие значения и представляет собой прогнозируемое значение для 50 известных случаев (мест), где вид присутствует, а (отсутствие) имеет более низкие значения и представляет прогнозируемое значение для 50 известных случаев (мест), где вид отсутствует. .
>p <- rnorm(50, mean=0.7, sd=0.3)
>a <- rnorm(50, mean=0.4, sd=0.4)
>par(mfrow=c(1, 2))
>plot(sort(p), col='red', pch=21)
>points(sort(a), col='blue', pch=24)
>legend(1, 0.95 * max(a,p), c('presence', 'absence'),pch=c(21,24), col=c('red', 'blue'))
>comb <- c(p,a)
>group <- c(rep('presence', length(p)), rep('absence', length(a)))
>boxplot(comb~group, col=c('blue', 'red'))
Часть, которую я не понимаю, является следующей:
>e <- evaluate(p=p, a=a)
>class(e)
>par(mfrow=c(1, 2))
>density(e)
>boxplot(e, col=c('blue', 'red'))