1
A
回答
5
因此,這是你的函數現在(希望你知道怎麼寫的R功能;如果沒有,檢查writing your own function):
f <- function (x) (pi/2) * (1/(1 + 0.25 * x^2))
f
上(-Inf, Inf)
定義,因此在這個範圍內的整合給出了一個不定積分。幸運的是,接近Inf
在x^(-2)
速度,所以整體是明確界定,並可以計算:
C <- integrate(f, -Inf, Inf)
# 9.869604 with absolute error < 1e-09
C <- C$value ## extract integral value
# [1] 9.869604
那麼你一定要規範化f
,因爲我們知道一個概率密度應該集成到1:
f <- function (x) (pi/2) * (1/(1 + 0.25 * x^2))/C
您可以通過繪製其密度:
curve(f, from = -10, to = 10)
1
既然我有可能的分佈函數,我想知道如何使用這個新的分佈函數創建一個隨機樣本
n = 1000
?
一個離題的問題,但沒有你做一個新的線程就可以回答。它很有用,因爲它變得微妙。
比較
set.seed(0); range(simf(1000, 1e-2))
#[1] -56.37246 63.21080
set.seed(0); range(simf(1000, 1e-3))
#[1] -275.3465 595.3771
set.seed(0); range(simf(1000, 1e-4))
#[1] -450.0979 3758.2528
set.seed(0); range(simf(1000, 1e-5))
#[1] -480.5991 8017.3802
所以我覺得e = 1e-2
是合理的。我們可以借鑑的樣本,做一個(縮放)直方圖和覆蓋密度曲線:
set.seed(0); x <- simf(1000)
hist(x, prob = TRUE, breaks = 50, ylim = c(0, 0.16))
curve(f, add = TRUE, col = 2, lwd = 2, n = 201)
相關問題
- 1. 在r中創建函數
- 2. 如何在R中創建一個概率分佈的矩陣
- 3. 函數創建R中
- 4. 在R中創建純函數
- 5. 在r中創建迴文函數
- 6. 如何使用分配gem在Ruby中創建指數分佈?
- 7. 如何使用ggplot和aes函數在R中創建分組條形圖?
- 8. 你如何創建參數分佈?
- 9. 如何在R中創建類似的Python函數?
- 10. 如何在R中創建股票波動函數?
- 11. 如何在R中創建線性變換函數?
- 12. 學習如何在R中創建函數 - 距離公式
- 13. 如何創建一個自定義的隨機分佈函數?
- 14. 如何在R中創建分數階乘設計?
- 15. 如何創建for循環來計算r中分箱數據的gini函數?
- 16. 在數學中創建分佈
- 17. 創建方程函數中的R
- 18. 創建循環函數中的R
- 19. 在R中構建函數
- 20. 如何在R中繪製度分佈
- 21. 如何在C++ 11中創建正常的整數分佈?
- 22. 在SQL中創建分區函數
- 23. 在R Markdown中爲HTML發佈創建循環分頁符
- 24. R中的佈局函數
- 25. Multinomial分佈的dmultinom函數R
- 26. R:分佈函數返回「隨機」值
- 27. 累積分佈函數,閾值,R
- 28. 如何在R中創建聚集布爾變量的圖形?
- 29. 如何在jsp中創建函數
- 30. 如何在Python 2.6中創建函數
像往常一樣,好的選擇。 – akrun
謝謝,再次精彩的迴應!現在我有了可能的分佈函數,我想知道如何使用這個新的分佈函數創建一個說n = 1000的隨機樣本? – Chris95
@ZheyuanLi完成。 – Chris95