Я хотел бы разделить два распределения плотности (или две гистограммы) в R. Я бы даже взял «вычитание одного из другого с помощью оператора» - но я не вижу очевидного способа делайте это, кроме выборки и вычитайте/разделите длинный путь.Разделите распределения плотности в R
Есть ли пакет функций/R, который позволяет управлять плотностями R без отбора проб?
Я новичок в R, я - трансплантат корня CERN, если это актуально.
Трюк с операциями с двумя (или более) плотностями заключается в том, чтобы их диапазоны были идентичными. Тогда векторы «x» будут одинаковыми, а векторы «y» будут хорошо вести себя с помощью обычных арифметических операций.
> xx <- rnorm(100)
> yy <- rnorm(100)
> yd <- density(yy, from=-3, to=3)
> xd <- density(xx, from=-3, to=3)
> plot(xd, col="red")
> str(xd)
List of 7
$ x : num [1:512] -3 -2.99 -2.98 -2.96 -2.95 ...
$ y : num [1:512] 0.0122 0.0126 0.0131 0.0136 0.0141 ...
$ bw : num 0.377
$ n : int 100
$ call : language density.default(x = xx, from = -3, to = 3)
$ data.name: chr "xx"
$ has.na : logi FALSE
- attr(*, "class")= chr "density"
> lines(yd$x, yd$y, col="red")
# difference of the densities
> lines(yd$x, xd$y-yd$y, col="green")
Я, кажется, повторяюсь: http://stackoverflow.com/questions/8808751/difference-between-two-density-plots –
Я не знаю, если это ответ. Я думаю, что они хотят построить плотность частного (или различия) двух случайных величин, учитывая плотность каждого из них. –
ОП запросил соотношение плотностей, ... не плотность отношения. –
Только для того, чтобы быть четким: у вас есть распределение плотности по двум случайным величинам x, y (по наиболее вероятному образцу) и хотите построить плотность отношения x/y? –
Если вы хотите работать с теоретическими дистрибутивами, для этого есть семейство пакетов, начиная с букв «до». –
Рикардо, да, это именно то, что я хотел бы. @DWin, спасибо, я рассмотрю пакеты (однако, дистрибутивы - экспериментальные выборки, так что это может не сработать). – user2533409