作者AmigoSafin ()
看板R_Language
标题[问题] 用R plot all the possibilities
时间Sun Sep 16 07:59:07 2018
[问题类型]:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
入门(写过其他程式，只是对语法不熟悉)
[问题叙述]:
大家好
我现在遇到一个问题 有４个蓝色玩具 B1 B2 B3 B4
另有6个红色玩具 R1 R2 ... R6
现在要看随机抽出三个玩具 颜色相同的机率
我的做法是 (2C1*4C3 + 2C1*6C3)/10C3 不知道这样对吗？
然後现在我遇到一个问题是
write a piece of R code to empirically compute the prob. of the three toys
in the same color, let n be the number of random samples in the bag.
我用的R code如下 当n=50,000
n <- 50000
f <- 0
for(i in 1:n)
f <- f + ifelse((sample(3)[1]==sample(3)[1]),1,0)
print(f/n)
但现在问题要我 Plot your estimated probabilities as a function of n (log scale again here,
please) 根据 n=50,100,200,.....,50000 用R绘制出机率的图表
我在想这是不是要先把取出三种同颜色的玩具的C.D.F.先求出来
再用plot C.D.F的方式去绘制呢？
想了两天都还是没头绪
老师有给蒙地卡罗模拟的参考资料
但我觉得这好像不是
另外还有buffon needle probability的参考资料
但我都无法refer 这些资料 又觉得并没有我想的复杂
但还是目前无法动弹
不知道有没有大大愿意给小鲁一点指点
感谢不尽～
也希望大家台风天都一切平安
[关键字]:
选择性，也许未来有用
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 8.41.66.215
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1537055950.A.AD8.html
1F：→ andrew43: sample(3)[1]==sample(3)[1] 是{1,2,3}抽出且放回二次的 09/16 08:33
2F：→ andrew43: sample(3)[1]==sample(3)[1] 的意义和三选一独立 09/16 08:35
3F：→ andrew43: 进行二次，看二次的结果是否相同。我看不懂这和问题的关 09/16 08:35
4F：→ andrew43: 系是什麽... 09/16 08:35
5F：→ AmigoSafin: 啊是这样吗？谢谢您～我再来改改 09/16 10:19
6F：→ AmigoSafin: 因为我是找老师的code来改 我也不甚理解其中的意义.. 09/16 10:33
7F：推 july81212: 用排列组合方式算完再丢不就没training 的意义吗…… 09/19 22:30
8F：推 july81212: 原本只是要你们制造一个机率函数 然後用seed去喂n次 09/19 22:33
9F：→ july81212: 再来filtering 相同出现的机率吧 09/19 22:33
大大真是高手
我回去找老师给的资料
有看到以下说明：
# The randomness generated a computer is not entirely random indeed
# It is referred to as pseudo-randomness
set.seed(8262009)
这就是您说的吗？
其实我还是不甚明白啊 因为我没有学过这个/_\
那就藉此再问一下
Compute Pr[1 X 3] where X FX(x) with fX(x) = 3x2 in [0,1].
我用的code:
n <- 10000
z <- rnorm(n)
proba <- length(which(abs(z)<1/4))/n
print(proba)
根据的是老师给的资料
Probability as expectation
（因为Monte Carlo Simulation有两种 Probability as expectation/
Probability using the uniform density as the weighting/importance function
我用的是前者
不知道这样对吗？
谢谢大家～～
※ 编辑: AmigoSafin (8.41.66.215), 09/21/2018 10:27:39
10F：→ AmigoSafin: 他应该是要我们先算theoretical再算empirical但我真的 09/21 10:28
11F：→ AmigoSafin: 很难推敲出R code要怎写 无奈 09/21 10:29
12F：→ andrew43: 你下的nrorm是常态mu=0,s=1乱数。uniform用runif 09/21 14:56
13F：→ andrew43: rnorm 09/21 14:56
14F：→ AmigoSafin: 那方便问一下a大 如何判断这函数uniform distribution 09/21 19:41
15F：→ andrew43: 不懂你的问题。 09/21 21:22
16F：→ AmigoSafin: 就是为何要用uniform的function做这题呢？不好意思>"< 09/21 22:24
17F：→ andrew43: 喔，我以为这题就是要用uniform来做 09/21 23:52
18F：→ andrew43: 建议你一次提一个问题，这样别人要回答的话较方便 09/21 23:54
19F：→ andrew43: 不然你觉得什麽要用normal? 09/22 00:04
20F：→ AmigoSafin: 我後来看了一下资料 这是一个multivariate 用MC的uni 09/22 08:34
21F：→ AmigoSafin: form去approxiate 这是我的理解～不好意思思绪有点乱 09/22 08:35