R: 模拟相关的硬币投掷

huangapple
go评论101阅读模式
英文:

R: Simulating Correlated Coin Flips

问题

我正在使用R编程语言。

我想要模拟硬币翻转，规则如下：

  • 如果是正面（heads），下一个翻转是正面的概率为0.6，是反面的概率为0.4。
  • 如果是反面（tails），下一个翻转是反面的概率为0.6，是正面的概率为0.4。

在R中，使用'markovchain'包，我可以这样做：

  1. library(markovchain)
  3. # 转移矩阵
  4. P <- matrix(c(0.6, 0.4, 0.4, 0.6), byrow = TRUE, nrow = 2)
  5. rownames(P) <- colnames(P) <- c("H", "T")
  7. mc <- new("markovchain", states = c("H", "T"), transitionMatrix = P)
  9. # 生成状态
  10. states <- rmarkovchain(n = 100, object = mc, t0 = "H")
  12. # 打印
  13. table(states)

输出结果类似于：

  1. > states
  2. [1] "H" "T" "T" "H" "T" "T" "T" "H" "H"

我的问题： 有人能展示如何在基本的R中实现这个吗？

我认为我需要：

  • 创建一个大小为'n'的空列表
  • 以0.5的概率为n[1]分配值为H或T
  • 编写一个IFELSE语句，以n[i] = ifelse(n[i-1] == "H", sample(c("H", "T"), prob = c(0.6, 0.4), sample(c("H", "T"), prob = c(0.4, 0.6))为条件

但我不确定如何实现这一点。

有人可以展示给我如何做吗？

谢谢！

英文:

I am working with the R programming language.

I want to simulate coin flips such that:

  • If heads, then next head with p = 0.6 and tail = 0.4
  • if tails, the next tails with p = 0.6 and heads = 0.4

Using the 'markovchain' package in R, I did this as follows:

  1. library(markovchain)
  3. # transition matrix
  4. P &lt;- matrix(c(0.6, 0.4, 0.4, 0.6), byrow = TRUE, nrow = 2)
  5. rownames(P) &lt;- colnames(P) &lt;- c(&quot;H&quot;, &quot;T&quot;)
  8. mc &lt;- new(&quot;markovchain&quot;, states = c(&quot;H&quot;, &quot;T&quot;), transitionMatrix = P)
  10. # Generate states
  11. states &lt;- rmarkovchain(n = 100, object = mc, t0 = &quot;H&quot;)
  13. # Print 
  14. table(states)

The output looks something like this:

  1. &gt; states
  2.   [1] &quot;H&quot; &quot;T&quot; &quot;T&quot; &quot;H&quot; &quot;T&quot; &quot;T&quot; &quot;T&quot; &quot;H&quot; &quot;H&quot;

My Question: Can someone please show me how I can do this in base R?

I think I need to:

  • create an empty list of size "n"
  • assign n[1] = H or T with prob 0.5
  • write an IFELSE statement that says n[i] = ifelse(n[i-1] == &quot;H&quot;, sample(c(&quot;H&quot;, &quot;T&quot;), prob = c(0.6, 0.4), sample(c(&quot;H&quot;, &quot;T&quot;), prob = c(0.4, 0.6))

But I am not sure how to do this.

Can someone please show me how to do this?

Thanks!

答案1

得分: 1

这是代码的翻译部分：

  1. 这是代码的翻译部分：
  3. ```R
  4. 这是代码的翻译部分：
  6. ```R
  7. 这是代码的翻译部分：
  9. ```R
  10. 这是代码的翻译部分：
  12. ```R
  13. 这是代码的翻译部分：
  15. ```R
  16. 这是代码的翻译部分：
  18. ```R
  19. 这是代码的翻译部分：
  21. ```R
  22. 这是代码的翻译部分：
  24. ```R
  25. 这是代码的翻译部分：
  27. ```R
  28. 这是代码的翻译部分：
  30. ```R
  31. 这是代码的翻译部分：
  33. ```R
  34. 这是代码的翻译部分：
  36. ```R
  37. 这是代码的翻译部分：
  39. ```R
  40. 这是代码的翻译部分：
  42. ```R
  43. 这是代码的翻译部分：
  45. ```R
  46. 这是代码的翻译部分：
  48. ```R
  49. 这是代码的翻译部分：
  51. ```R
  52. 这是代码的翻译部分：
  54. ```R
  55. 这是代码的翻译部分：
  57. ```R
  58. 这是代码的翻译部分：
  60. ```R
  61. 这是代码的翻译部分：
  63. ```R
  64. 这是代码的翻译部分：
  66. ```R
  67. 这是代码的翻译部分：
  69. ```R
  70. 这是代码的翻译部分：
  72. ```R
  73. 这是代码的翻译部分：
  75. ```R
  76. 这是代码的翻译部分：
  78. ```R
  79. 这是代码的翻译部分：
  81. ```R
  82. 这是代码的翻译部分：
  84. ```R
  85. 这是代码的翻译部分：
  87. ```R
  88. 这是代码的翻译部分：
  90. ```R
  91. 这是代码的翻译部分：
  93. ```R
  94. 这是代码的翻译部分：
  96. ```R
  97. 这是代码的翻译部分：
  99. ```R
  100. 这是代码的翻译部分：
  102. ```R
  103. 这是代码的翻译部分：
  105. ```R
  106. 这是代码的翻译部分：
  108. ```R
  109. 这是代码的翻译部分：
  111. ```R
  112. 这是代码的翻译部分：
  114. ```R
  115. 这是代码的翻译部分：
  117. ```R
  118. 这是代码的翻译部分：
  120. ```R
  121. 这是代码的翻译部分：
  123. ```R
  124. 这是代码的翻译部分：
  126. ```R
  127. 这是代码的翻译部分：
  129. ```R
  130. 这是代码的翻译部分：
  132. ```R
  133. 这是代码的翻译部分：
  135. ```R
  136. 这是代码的翻译部分：
  138. ```R
  139. 这是代码的翻译部分：
  141. ```R
  142. 这是代码的翻译部分：
  144. ```R
  145. 这是代码的翻译部分：
  147. ```R
  148. 这是代码的翻译部分：
  150. ```R
  151. 这是代码的翻译部分：
  153. ```R
  154. 这是代码的翻译部分：
  156. ```R
  157. 这是代码的翻译部分：
  159. ```R
  160. 这是代码的翻译部分：
  162. ```R
  163. 这是代码的翻译部分：
  165. ```R
  166. 这是代码的翻译部分：
  168. ```R
  169. 这是代码的翻译部分：
  171. ```R
  172. 这是代码的翻译部分：
  174. ```R
  175. 这是代码的翻译部分：
  177. ```R
  178. 这是代码的翻译部分：
  180. ```R
  181. 这是代码的翻译部分：
  183. ```R
  184. 这是代码的翻译部分：
  186. ```R
  187. 这是代码的翻译部分：
  189. ```R
  190. 这是代码的翻译部分：
  192. ```R
  193. 这是代码的翻译部分：
  195. ```R
  196. 这是代码的翻译部分：
  198. ```R
  199. 这是代码的翻译部分：
  201. ```R
  202. 这是代码的翻译部分：
  204. ```R
  205. 这是代码的翻译部分：
  207. ```R
  208. 这是代码的翻译部分：
  210. ```R
  211. 这是代码的翻译部分：
  213. ```R
  214. 这是代码的翻译部分：
  216. ```R
  217. 这是代码的翻译部分：
  219. ```R
  220. 这是代码的翻译部分：
  222. ```R
  223. 这是代码的翻译部分：
  225. ```R
  226. 这是代码的翻译部分：
  228. ```R
  229. 这是代码的翻译部分：
  231. ```R
  232. 这是代码的翻译部分：
  234. ```R
  235. 这是代码的翻译部分：
  237. ```R
  238. 这是代码的翻译部分：
  240. ```R
  241. 这是代码的翻译部分：
  243. ```R
  244. 这是代码的翻译部分：
  246. ```R
  247. 这是代码的翻译部分：
  249. ```R
  250. 这是代码的翻译部分：
  252. ```R
  253. 这是代码的翻译部分：
  255. ```R
  256. 这是代码的翻译部分：
  258. ```R
  259. 这是代码的翻译部分：
  261. ```R
  262. 这是代码的翻译部分：
  264. ```R
  265. 这是代码的翻译部分：
  267. ```R
  268. 这是代码的翻译部分：
  270. ```R
  271. 这是代码的翻译部分：
  273. ```R
  274. 这是代码的翻译部分：
  276. ```R
  277. 这是代码的翻译部分：
  279. ```R
  280. 这是代码的翻译部分：
  282. ```R
  283. 这是代码的翻译部分：
  285. ```R
  286. 这是代码的翻译部分：
  288. ```R
  289. 这是代码的翻译部分：
  291. ```R
  292. 这是代码的翻译部分：
  294. ```R
  295. 这是代码的翻译部分：
  297. ```R
  298. 这是代码的翻译部分：
  300. ```R
  301. 这是代码的翻译部分：
  303. ```R
  304. 这是代码的
  306. <details>
  307. <summary>英文:</summary>
  309. Here it is.
  311.     Fair=function() sample( c(rep(&#39;H&#39;,5), rep(&#39;T&#39;,5)),1)
  312.     A=function() sample( c(rep(&#39;H&#39;,6), rep(&#39;T&#39;,4)),1)
  313.     B=function() sample( c(rep(&#39;H&#39;,4), rep(&#39;T&#39;,6)),1)
  314.     x=Fair()
  315.     for (i in 1:10) x=c(x,ifelse(tail(x,n=1)==&#39;H&#39;, A(), B()))
  316.     print(x)
  320. </details>

huangapple
  • 本文由 发表于 2023年2月19日 06:58:53
  • 转载请务必保留本文链接：https://go.coder-hub.com/75496896.html
  • r
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定