2011-05-13 30 views
0

我有一個如下的矩陣[500,2],所以我們有500行和2列,左邊給出了X觀察值的索引,右邊給出了這個X實現的概率,所以 - 一個典型的概率密度關係。R關於以正確的方式繪製概率/密度直方圖的問題

所以,我的問題是,如何繪製直方圖的正確方式,以便X軸是X指標,而Y軸是密度(0.01-1.00)。估計器的帶寬是0.33。

在此先感謝!

整個數據的結尾是這樣的:只是一點點方向

[490,] 2.338260830 0.04858685 
[491,] 2.347839477 0.04797310 
[492,] 2.357418125 0.04736149 
[493,] 2.366996772 0.04675206 
[494,] 2.376575419 0.04614482 
[495,] 2.386154067 0.04553980 
[496,] 2.395732714 0.04493702 
[497,] 2.405311361 0.04433653 
[498,] 2.414890008 0.04373835 
[499,] 2.424468656 0.04314252 
[500,] 2.434047303 0.04254907 

@everyone, 是的,我以前的估計做,所以..的帶寬就是我提到的,數據按從低到高的順序排列,因此相應的開始概率爲0,22,峯值約爲0.48,結束時爲0.15。

密度線被繪製爲一個魅力,但我必須做的另外是繪製直方圖!所以,我該如何做到這一點,正確地命令塊(ho數據被分割成盒子等)。

有什麼建議嗎?

這裏是估計後數據的一部分,所有值都是離散的,所以我認爲直方圖可以創建..希望。

[491,] 4.956164 0.2618131 
[492,] 4.963014 0.2608723 
[493,] 4.969863 0.2599309 
[494,] 4.976712 0.2589889 
[495,] 4.983562 0.2580464 
[496,] 4.990411 0.2571034 
[497,] 4.997260 0.2561599 
[498,] 5.004110 0.2552159 
[499,] 5.010959 0.2542716 
[500,] 5.017808 0.2533268 
[501,] 5.024658 0.2523817 

最好的問候, 欣賞的快速響應!(弓)

什麼將做的工作是創建一個直方圖只爲指數,在某種程度上X25/X50每個,例如將它們分組...並計算每個25或50/100/150/200/250等盒子的平均概率..?

回答

1

假設行是爲了從最低到x的最高價值,因爲他們似乎是,你可以使用默認的繪圖命令,你需要的唯一的變化是類型:

plot(your.data, type = 'l') 

編輯:

好吧,我不知道這是不是密度圖更好,但它可以做到:

x = dnorm(seq(-1, 1, length = 500)) 
x.bins = rep(1:50, each = 10) 
bars = aggregate(x, by = list(x.bins), FUN = sum)[,2] 
barplot(bars) 

在你的情況下,從矩陣的第二列的概率替換x。

EDIT2:

關於第二個想法,這纔有意義,如果你的500行表示離散事件。如果它們是連續分佈函數中的點,則將它們加在一起,因爲我已經做了不正確的事情。在數學上,我不認爲你只能使用該範圍內的幾個點來產生範圍的分類概率。

+0

@Tyler,這個工作就像一個魅力,非常感謝,我仍然需要做的一件事是將兩列數據放在一個直方圖中,x軸=觀察值,y軸=概率/密度。 .. – 2011-05-13 21:47:47

+0

好的,我在下半場添加了一些代碼。 – Tyler 2011-05-13 22:00:32

+0

嗨,爲了更加準確,我會在評估之後發佈這些數據,第一篇文章中的數據只是一個正態分佈,所以現在我將發佈估計完全離散的數據。將在2分鐘內編輯我的問題。 – 2011-05-13 22:09:07

0

假設M是矩陣。是不是這樣:

plot(x=M[ , 1], y = M[ , 2]) 

您已經完成密度估計,因爲這不是原始數據。