看板 Statistics 關於我們 聯絡資訊
感謝推文的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
chien533:類似我之前的寫法,只要在set.seed上面動點手腳即可 12/25 03:10
chien533:把原來放在第一行的set.seed(0)搬到原本的for迴圈裡 12/25 03:13
chien533:並改成set.seed(j) 12/25 03:13
chien533:然後再從外面加一個迴圈 for(j in 1:1000){ ...... } 12/25 03:13
chien533:忘了講,你要另外設一個 matrix object 來放模擬的資料 12/25 04:12
laba1014:感謝C大 ^ ^ 我把問題解決了 (C大都這麼晚睡阿XD) 12/25 13:28
chien533:1.我在米國 2.最好把所有要生成的object都放在for j的迴 12/25 15:59
chien533:圈裡面,因為error term也需要用不同的seed 12/25 16:01