作者greatjus (jus)
看板R_Language
标题[问题] 不断将矩阵扩充的办法
时间Mon Jan 9 05:03:46 2017
又来麻烦这个版的高手了,深感抱歉
目前的问题:
我有一个矩阵,始终都只会有12列
一开始只有第一行的值作为起始值,尔後
我每次会将前一行的值经过一个apply的方程式
再度生成新的12个值。
想请问:
有没有能够避免使用rbind,却又能够不停的把新的12个值放入矩阵的下一行方法。
请问高手们有什麽推荐的方程式吗?
(不使用rbind,是因为同样的事情
我可能会重复做20万次甚至到80万次以上,使用bind 笔电好像会爆炸)
PS:这个步骤其实就是在做 MCMC:metropolis
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 223.136.209.189
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1483909428.A.B77.html
1F:→ Edster: 台湾的列是row,行是col, 现在还是这样教吧。第二,先创有 01/09 09:07
2F:→ Edster: 80万列的matrix, 然後把apply改用for就好. 01/09 09:07
3F:→ Edster: 你要做mcmc就直接找package mcmc,这样比较没什麽问题, 01/09 09:12
4F:→ Edster: 不然所有笔的参数初始值都是乱数也好,看你用有方向性的 01/09 09:12
5F:→ Edster: 初始值应该没比较好。 01/09 09:12
6F:→ greatjus: 初始值是刻意要自订的,不过mcmc package 确实有, 01/09 14:37
7F:→ greatjus: 谢谢你的建议,我使用看看 PACKAGE 好了 01/09 14:38
8F:→ greatjus: 可能太累了,我竟然没养到事先创好矩阵就好,谢谢你!! 01/09 14:48