#R语言# 生成随机数

生成随机数

生成随机数的函数名的格式为r+分布名缩写，例如正态分布(norm)对应于rnorm，均匀分布(unif)对应于runif
若想生成整数，在前面加上 options(digits=1)

1. 正态分布

# rnorm(n,mean = 0,sd = 1) 生成服从正态分布 N(mean,sd^2) 的随机变量
# n为随机数数目
> rnorm(5)
[1] -0.5504527 -0.5495538  0.7608639 -0.6813592 -1.0085563

2. 均匀分布

# runif(n, min = 0, max = 1) 生成服从均匀分布 U(min,max) 随机变量
# n为随机数数目
> runif(5)
[1] 0.8125779 0.2392292 0.9477868 0.4024770 0.6578547

3. 指数分布

# rexp(n,rate = 1) 生成服从指数分布（率参数为rate）的随机变量
# n为随机数数目
> rexp(5)
[1] 1.4806552 1.1144723 0.9901113 0.6998010 0.4210202

指数分布的概率密度函数为
$$f(x)=\{\begin{array}{ll}\lambda e^{-\lambda x}& x>0\\ 0& x\le 0\end{array}$$

4. F分布

# rf(n,df1,df2,ncp) 生成服从F分布(自由度n1 = df1，n2 = df2)的随机变量
# df1和df2必须输入
# n为随机数数目，ncp为非中心参数
> rf(5,3,5,1)
[1] 0.1720921 1.1186105 1.7349433 3.5147749 0.6602628

F分布的概率密度函数为
$$f\left(x,n_1,n_2\right)=\{\begin{array}{l}\frac{{\left(n_1/n_2\right)}^{\frac{n_1}{2}}}{B\left(n_1/2,n_2/2\right)}{x}^{\frac{n_1}{2}-1}{\left(1+\frac{n_1}{n_2}x\right)}^{-\frac{n_1+n_2}{2}},x>0\\ 0,x\le 0\end{array}$$

若总体$X\sim N(\mu ,1)$与总体$Y\sim N(\mu ,1)$独立，$\left(X_1,X_2,\dots ,X_{n_1}\right)$为来自$X$的一个样本，$\left(Y_1,Y_2,\dots ,Y_{n_2}\right)$为来自$Y$的一个样本，则统计量
$$F=\frac{{\sum }_{i=1}^{n_1}{X}_i^2}{n_1}/\frac{{\sum }_{i=1}^{n_2}{Y}_i^2}{n_2}$$
服从自由度为$n_1$和$n_2$，非中心参数$\delta =n{\mu }^2$的非中心F分布，记为$F\sim F\left(n_1,n_2,\delta \right)$。

---

以下函数有空再补充

函数名	分布名
rgamma	Gamma分布
rgeom	几何分布
rhyper	超几何分布
rlnorm	对数正态分布
rlogis	Logistic分布
rmultinom	多项分布
rnbinom	负二项分布
rpois	泊松分布
rt	t分布
rchisq	卡方分布