无法制作 for 循环来模拟我的马尔可夫链 10000 次，并在表格中每次获取 c 的频率以制作直方图

问题描述

这是我的链：@Override

我的链条没有问题，

然而，我真的在为 for 循环而苦苦挣扎。这是我目前所拥有的：

chain[i+1]<-sample(nucleotides,1,p=P[chain[i],])

freqc 定义为：

for (q in 1:nsim)
{
chain <- rep(NA,LengthOfChain)
chain[i+1]<-sample(nucleotides,])
}  
table(freqc)

所以我知道 freqc 可以从我的链的 ONE 模拟中提取 c 的频率，但我不知道如何对所有这些进行此操作，然后将其生成到表格中，以便我以后可以制作一个它的直方图。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）