2、抽样与分布

#抽样与分布#
#常用的统计函数#
library(stats)
length(x)#长度
sum(x)#求和
prod(x)#给出 x 中所有元素的乘积
min(x)
max(x)#最大值
which.min(x)
which.max(x)#x中元素最大值所在位置
mean(x)#均值
median(x)#中位数
mad(x)#中位数绝对偏差
var(x)#方差
sd(x)#标准差
range(x) #极差
quantile(x)#分位数
quantile(x,prob) #求任意概率的分位数 自己设定的0~1上的数
fivenum(x)#五分位数
cumsum(x)#累积和
cumprod(x)
cummax(x)
cummin(x)
rev(x)   #取逆序
sort(x,decreasing=F)  
order(x, decreasing=F) 
rank(x)
var(x,y)
cov(x,y)
cor(x,y)
outer(x,y)
#对于矩阵x
var(x)  #协方差阵
cor(x)  #相关系数阵
###随机抽样！！！！！！！！！！！！！！！！！！！！！！！
#等可能的不放回的随机抽样
sample(x, n)
sample(1:52, 4)
#等可能的有放回的随机抽样
sample(x, n, replace=TRUE)
sample(c("H","T"), 10, replace=TRUE)
#不等可能的随机抽样
sample(x, n, replace=TRUE, prob=y)
sample(c("H","T"), 10, replace=T, prob=c(0.8, 0.2))#prob 各个元素被抽到的概率

#stats包中有最常用的18个分布函数名
#离散分布(6个)
binom             二项分布
multinorm       多项分布
nbinorm          负二项分布
geom             几何分布
hyper             超几何分布
pois              泊松分布
#连续分布(12个)
beta               贝塔分布
cauchy           柯西分布
chisq             卡方分布
exp                指数分布
f                     F分布
gamma          伽玛分布
logis               logistic分布
lnorm             对数正态分布
norm              正态分布
t                      t分布
unif                  均匀分布
weib               威布尔分布
#MCMCpack包中的4个分布
dirichlet            Dirichlet分布
invgamma        逆伽玛分布
wish                 Wishart分布   
iwish                逆Wishart分布
#mvtnorm包中的2个分布
mvnorm            多元正态分布
mvt                   多元t分布   
#evir包中的2个分布
gev                   广义极值分布
gpd                   广义Pareto分布     
#################
#若R中分布的函数名为func, 则四类函数的调用格式为
#概率密度函数：dfunc(x, p1, p2,…), x为数值向量，p1为参数
#(累积)分布函数：pfunc(q, p1, p2,…), q为数值向量
#分位数函数：qfunc(p, p1, p2,…), p为由概率构成的向量
#随机数函数：rfunc(n, p1, p2,…), n为生产数据的个数
#其中p1, p2,…是分布的参数值。
#例子：
x <- rnorm(1000)
#x的直方图
hist(x, col="light blue", probability=T, ylim=c(0, max(density(x)$y, 0.4)))
#x的核密度估计值，并连接成线
lines(density(x), col="red", lwd=3)
#标准正态分布的密度函数值
curve(dnorm(x), col="blue", lwd=3, lty=3, add=T)