[问题类型]:
程式咨询(我想用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 这些资料 又觉得并没有我想的复杂
但还是目前无法动弹
不知道有没有大大愿意给小鲁一点指点
感谢不尽～
也希望大家台风天都一切平安
[关键字]:
选择性，也许未来有用

sample(3)[1]==sample(3)[1] 是{1,2,3}抽出且放回二次的sample(3)[1]==sample(3)[1] 的意义和三选一独立进行二次，看二次的结果是否相同。我看不懂这和问题的关系是什么...

啊是这样吗？谢谢您～我再来改改因为我是找老师的code来改 我也不甚理解其中的意义..

用排列组合方式算完再丢不就没training 的意义吗……原本只是要你们制造一个机率函数 然后用seed去喂n次再来filtering 相同出现的机率吧

他应该是要我们先算theoretical再算empirical但我真的很难推敲出R code要怎写 无奈

你下的nrorm是常态mu=0,s=1乱数。uniform用runifrnorm

那方便问一下a大 如何判断这函数uniform distribution

不懂你的问题。

就是为何要用uniform的function做这题呢？不好意思>"<

喔，我以为这题就是要用uniform来做建议你一次提一个问题，这样别人要回答的话较方便不然你觉得什么要用normal?

我后来看了一下资料 这是一个multivariate 用MC的uniform去approxiate 这是我的理解～不好意思思绪有点乱