Я хочу реализовать следующую функцию в R.
, где m и o указывают модель и наблюдение соответственно. Вот пример кода с данными и моим кодом:
xm <- runif(100,1,1000)+50
xo <- runif(100,1,1000)
y <- matrix(NA, nrow = length(xo), ncol=1)
for(r in 1:length(xm)){
y[r,1] <- xm[r] + quantile(ecdf(xm)(xo[r]), 0.5) - quantile(ecdf(xm)(xm[r]), 0.5)
}
Я ожидаю, что моя новая "y" имеет меньшую ошибку и значительно отличающиеся значения по сравнению с "xm" в моем наборе данных, но это не так так как значения почти одинаковые.