Равная длина шага в оценке плотности в R - PullRequest
Новая
       21

Равная длина шага в оценке плотности в R

0 голосов
/ 03 апреля 2019

Различают два распределения, но хотят контролировать «длину шага» по значению x графика плотности.В настоящее время R дает то, что кажется оптимизированной или случайной длине шага

Пробовал возиться с пропускной способностью и просматривает сообщения разной плотности, но ни одна из них не решает эту проблему.Справочный документ R также расплывчатый.

Вот код, который я пробовал.

data2 <- rlnormAlt(n = 100000, mean = 5.0, cv = 3.5/5.0)
data1 <- rlnormAlt(n = 100000, mean = 12.3, cv = 19.6/12.3)

x = density(data1, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$x
y = density(data1, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$y -
density(data2, bw = 0.01, kernel = "gaussian", from= range(c(data1, data2))[1], to=range(c(data1, data2))[2], n = 10000)$y

z <- cbind(x,y) %>% as.data.frame() %>% dplyr::filter(., y > 0, x > 4)
sapply(seq(1,10,1), function(i) z$x[i+1] - z$x[i])

Вы можете видеть, что длина шага "иногда" согласована, но я не могу установить это.

Например, я хочу, чтобы длина шага была каждые 1 или 0,5 и т. Д. (Т. Е. Seq (0,50,1), где я могу выбрать это).

Любая помощь будетоценили.

...