R

一元随机数

runif(5,1,3) #生成五个（1，3）之间的均匀分布随机数
runif(3) #默认生成3个（0，1）之间的
rnorm(3,1,5) #产生3个均值=1，标准差为5的正态分布随机数
x=rnorm(100)
hist(x,prob=T,main="normal mu=0,sigma=1") #hist means histogram 直方图
#prob=T??
curve(dnorm(x),add=T)  #在直方图上添加标准正态分布密度函数线

• rnorm 生成正态分布随机数
• dnorm 生成正态分布密度函数
• pnorm 生成正态分布的累计概率密度函数
• qnorm 生成正态分布的分位数函数
• 如： pnorm(2) = P(X<2), P(X<a)=0.95— a=qnorm(0.95）
  

多元随机数

mvtnorm安装包中的多元正态

matrix(c(10,3,3,2),2,2)
x=rmvnorm(1000, mean=rep(0,5),sigma=diag(5)) #产生10个均值为0，方差为1且分量独立的5维正态分布随机向量 
colname(x)<-c("X","Y")  #坐标轴命名
head(x) #显示前面几个
var(x)
colMeans(x) #求均值

MASS程序包的多元正态

#二元正态随机数
mu<-c(0,1)
Sigma<-matrix(c(1,0.5,0.5,1),2,2)
n<-1000
library(MASS)
biv<-mvrnorm(n,mu,Sigma)
colnames(biv)<-c("X","Y")

#参数估计
mu.hat<-colMeans(biv)  #均值X_
Sigma.hat<-cov(biv) #方差S

#常数密度轮廓线
install.packages("mixtools")
library(mixtools)
plot(biv)
ellipse(mu<-colMeans(biv),sigma<-cov(biv),alpha=.05,col='red')  #ellipse椭圆
#alpha:Probability to be excluded from the ellipse.
#The default value is alpha = .05, which results in a 95% ellipse.
points(t(mu),col='red',pch=19) #绘制点t(mu),  pch为点的类型

R语言低级绘图points函数的用法
R语言低级绘图text函数的用法

程序包的安装与调用

install.packages("mvtnorm")  #安装程序包 
library(mvtnorm)  #加载程序包 【每一次使用前都需要加载】
?rmvnorm # 查看程序包的内部具体函数
help(package="程序包名")	#查看程序包的所有内容

*the density function密度函数