R语言提供了summary()函数来计算数据的描述统计,但是其中计算的结果很少。这里老师给大家编写一个自定义函数,可以详细计算数据的描述统计。
- miaoshu<-function(x)
- { #函数体
- n=length(x); #样本量
- avg=mean(x); #样本均值
- med=median(x); #样本中位数
- r=max(x)-min(x); #极差
- s=sd(x);#标准差
- myvar=var(x); #方差
- ups=max(x); #最大值
- lows=min(x); #最小值
- sums=sum(x); #求和
- Q1=quantile(x,0.25);#Q1
- Q3=quantile(x,0.75); #Q3
- s=data.frame("N"=n,"平均数"=avg,"中位数"=med, "方差"=myvar,"标准差"=s,"极差"=r, Q1=Q1,Q3=Q3, "最大值"=ups,"最小值"=lows,"求和"=sums) return(s); #返回值
- }
上面的自定义函数返回值是数据框,下面的小例子,是将该函数与summary()进行对比:
- data=rnorm(100)
- ans1=miaoshu(data)
- ans2=summary(data)
大家自己运行,即可看到结果。