1. sd()
求一组数据的标准差
> x = rep(1,15)
> x
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
> sd(x)
[1] 0
2.var()
求一组数据的方差
> x = runif(5) #生成均匀分布的随机数 runif(num,min = 0,max = 1)
> sd(x)
[1] 0.1953973
> var(x)
[1] 0.03818011
3.随机数与抽样模拟(正态分布、二项分布、指数分布。) https://blog.****.net/lilanfeng1991/article/details/18505723
rbinom(n,size,p) size表示每次进行伯努利试验的次数,p表示每次实验的成功概率,n表示共进行了n次此二项分布试验,返回的是每次二项分布试验中成功的次数。
par(mfrow=c(1,3)) #图形以一行三列显示
p=0.25
for( n in c(10,20,50))
{ x=rbinom(100,n,p)
hist(x,prob=T,main=paste("n =",n))
#hist函数绘制出柱状图
xvals=0:n
points(xvals,dbinom(xvals,n,p),type="h",lwd=3) #lwd是线条的宽度,默认为1
#绘制柱状图的竖线
}
par(mfrow=c(1,1))
4.par() https://blog.****.net/heidlyn/article/details/54581605
5.predict 函数
predict(fm,data.frame(xx = 58),interval = "confidence",level = 0.95)
#求置信水平为0.95的置信区间