作者s3011 (薛丁格算子)
看板D95323XXX
標題Re: [計量] hist指令裡的 freq 參數
時間Thu Mar 8 12:23:53 2007
※ 引述《obscure (馬不停蹄的憂傷)》之銘言：
: ※ 引述《s3011 (薛丁格算子)》之銘言：
: : 今天測試的結果發現 R 2.4版hist指令有bug
: : 有些情況freq不管選TRUE or FALSE都無法正確把圖形畫成density
: : 造成各位花時間try指令非常抱歉
: 助教 請問你是指
: normal_random_variable <- rnorm(n,0,1) 所產生的樣本無法畫成normal density嗎?
不是耶 ^^;
而是hist指令有些情況無法以正確尺度畫出圖形
導致後續無法把normal density畫上去
: 可以載入程式套件中的grid和lattice
: 然後用densityplot(normal_random_variable)
: 就可以畫出樣本的density
: 然後再用這個圖和模擬CLT的圖相互比較
: 請問用此方法比較接不接近normal是否可行?
當然可以啊
不過你要畫的normal density
要用dnorm畫出來的 才是"理論"上的常態分配
用densityplot畫出來的 是"模擬"出來的資料 透過無母數方式畫出來的
所以修正一下你的敘述
CLT的資料 用densityplot畫
接著再用dnorm配合plot指令繪出常態分配的函數圖形
比如說像這樣
> mydata<-seq(-3,3,0.01)
> plot(mydata,dnorm(mydata),type="l")
--
不要花太多時間畫圖啦^^; 以後可以慢慢練
習題6.11頗難才是要花最多時間的地方
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.166.127.117