核密度中的多个列并将其与另一个图组合绘制

Question

我有一个由模拟值组成的数据。数据集是一个95 * 525矩阵，下面我打印了一小部分数据样本。

0.01076     0.01076     0.01076     0.01076     0.01076
0.0105259   0.010586823 0.010968957 0.010919492 0.010685372
0.010461666 0.010271274 0.01124356  0.010362122 0.010699493
0.010621132 0.009699747 0.01143447  0.010169863 0.010920117
0.01090401  0.009537089 0.011461832 0.010613986 0.010964597
0.010876345 0.008634657 0.011353565 0.0111837   0.010456199
0.010245947 0.009195114 0.012275926 0.011230311 0.010093507
0.010056854 0.009557125 0.01199076  0.010158304 0.009759821
0.008903207 0.009537452 0.011057318 0.010154896 0.00985161
0.008249492 0.010232035 0.01133204  0.009719779 0.00995871
0.006726547 0.010441961 0.010392739 0.009166036 0.01043827
0.007159706 0.010254763 0.009552638 0.009371703 0.010261791
0.007692135 0.01099915  0.009913062 0.008190215 0.009553324

为了绘制内核密度，我已经下载了包KernSmooth，它运行良好。但是，我想创建包含所有模拟的一个内核密度图。也就是说，我想为 95列绘制一个内核密度，而不是仅包含一列的内核密度。

我用来获得单列内核的代码如下

> kernel=bkde(cir[,1])
> plot(kernel, type="l")

如果我忽略[，1]使其包含所有列，我会收到错误。

此外，我在一张图中为我的模拟值绘制了不同的路径。我接下来要做的是将内核密度图旋转90度，然后将其与模拟路径图的结尾相结合。如果我对此的解释是漫反射的，我可以尝试提供图片作为示例。

Answer 1

使用

kernel <- bkde(as.vector(cir))

获取所有数据的内核密度估计值。

要将图形旋转90度，只需在x命令中交换y和plot坐标，例如在

中

x <- seq(-3,3,.1)
plot(x, dnorm(x))
plot(dnorm(x), x)

希望有所帮助。