knit
編譯,產生.html
檔,需印出R
程式碼及執行結果。r1092
。密碼: xxx
。學號-姓名-R-exam2.Rmd
」、「學號-姓名-R-exam2.html
」。 (學號及姓名,改成自己)-2
」、「-3
」,例如: 「學號-姓名-R-exam2-2.Rmd
」等等。You can modify the html file, but please keep the link www.wftpserver.com at least.
」, 請將滑鼠移至「網址列」後,按「Enter」即可。若再不行,請換其它瀏覽器(IE/Edge/Firefox/Chrome)。knit
無法編譯產生文件, 則可以「不執行有錯的程式碼」,但必需列印此段程式碼。助教會依照狀況部份給分。Student’s \(t\)-distribution has the probability density function (pdf) given by \[
p(x)={\frac {\Gamma ({\frac {\nu +1}{2}})}{{\sqrt {\nu \pi }}\,\Gamma ({\frac {\nu }{2}})}}\left(1+{\frac {x^{2}}{\nu }}\right)^{\!-{\frac {\nu +1}{2}}},
\] where \(\nu\) is the number of degrees of freedom and \(\Gamma\) is the gamma function. 試寫一R
函數(命名為t_pdf
,輸入為\(x\)、\(\nu\),輸出為\(p(x)\)), 計算\(t\)分配之機率密度函數值,並利用此函數畫出下圖(注意: \(x-\), \(y-\)軸標號、圖例說明(legend)、線條顏色)。 (註: \(\Gamma\) 函數在R
是什麼指令,請自己查。\(\nu = + \infty\)以\(\nu = 10\)代入。畫圖以R base graphics
套件或ggplot2
皆可。) (圖片來源: https://en.wikipedia.org/wiki/Student%27s_t-distribution)
# your source code here
若隨機變數\(Y\)服從卜瓦松分配\((X\sim Poisson(\lambda))\),其累積機率分佈函數為: \[
F_X(x)=P(X \leq x) =\sum_{k=0}^x \frac{\lambda^ke^{-\lambda}}{k!}。
\] 試寫一R
函數(命名為poisson_cdf
,輸入為\(x\)、\(\lambda\),輸出為\(F_X(x)\)), 計算卜瓦松分配之累積機率分佈函數值, 並利用此函數印出下表(印出紅色框部份的表格即可)。
# your source code here
給定常數\(\lambda\),且 \(X_n \sim Binomial\left(n, \frac{\lambda}{n}\right)\), 則 \[ P(X_n=x) = {C^n_x} \left(\frac{\lambda}{n}\right)^x \left(1-\frac{\lambda}{n}\right)^{n-x} \longrightarrow \frac{e^{-\lambda}\lambda^x}{x!}\quad as \ n \rightarrow \infty \] 亦即,Poisson分配是二項式分配\(Binomial \left(n, \frac{\lambda}{n}\right)\)的一個極限分配。
若\(n=30, \lambda=2.4\),利用R
的內鍵指令(二項式分佈機率質量函數及卜瓦松 分佈機率質量函數), 列出下表(data.frame
)來驗証Poisson極限定理, 其中最一欄位(diff
)為兩者差之絕對值。
x binomial.pmf poisson.pmf diff
1 0 0.0820 0.0907 0.0087517497
2 1 0.2138 0.2177 0.0038982090
3 2 0.2696 0.2613 0.0083375766
4 3 0.2188 0.2090 0.0097959196
5 4 0.1284 0.1254 0.0030235072
6 5 0.0581 0.0602 0.0021224767
7 6 0.0210 0.0241 0.0030372714
8 7 0.0063 0.0083 0.0019823213
9 8 0.0016 0.0025 0.0009083534
10 9 0.0003 0.0007 0.0003270817
11 10 0.0001 0.0002 0.0000976313
# your source code here
有成功編譯出正確的 「學號-姓名-R-exam1.pdf
」及「學號-姓名-R-exam1.doc
」,並上傳。 以下數學式是測試MikTeX/LaTeX,請勿刪。這是常態分佈的機率密度函數: \[
f(x)=\frac{1}{\sqrt{2\pi}\sigma}e^{\frac{(x-\mu)^2}{2\sigma^2}}
\]