Re: [問題] 如何用R generate random sample?

看板Statistics作者 (嘿喲~加油!!!)時間18年前 (2007/12/24 19:22), 編輯推噓4(404)
留言8則, 2人參與, 最新討論串2/2 (看更多)
感謝推文的L大和C大 我現在解決第一步了 再請問一下 如果 我現在想造一個像panel data 的 樣本群 \ j=1 2 .... 1000 t=1 z_11 2 . . . . . . . . . . . 500 z_500,1000 一樣是AR(1): Z[t,j] =β* Z[t-1,j] + u[t] u~iid gamma(500,a,b)  (其實就是第一步做一千遍 只是我希望排成array的型式 之後好方便從裡面抽數據) 小弟本來想用兩個loop 可是做不出來 希望懂得板友可以不吝指教 感謝~ -- -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.200.200

12/25 03:10, , 1F
類似我之前的寫法,只要在set.seed上面動點手腳即可
12/25 03:10, 1F

12/25 03:13, , 2F
把原來放在第一行的set.seed(0)搬到原本的for迴圈裡
12/25 03:13, 2F

12/25 03:13, , 3F
並改成set.seed(j)
12/25 03:13, 3F

12/25 03:13, , 4F
然後再從外面加一個迴圈 for(j in 1:1000){ ...... }
12/25 03:13, 4F

12/25 04:12, , 5F
忘了講,你要另外設一個 matrix object 來放模擬的資料
12/25 04:12, 5F

12/25 13:28, , 6F
感謝C大 ^ ^ 我把問題解決了 (C大都這麼晚睡阿XD)
12/25 13:28, 6F

12/25 15:59, , 7F
1.我在米國 2.最好把所有要生成的object都放在for j的迴
12/25 15:59, 7F

12/25 16:01, , 8F
圈裡面,因為error term也需要用不同的seed
12/25 16:01, 8F
文章代碼(AID): #17RvQ3OY (Statistics)
文章代碼(AID): #17RvQ3OY (Statistics)