1.求平均数。
> x<-c(1,10,20,30,40,50,NA,60)
> xm<-mean(x)
> xm
[1] NA
#na.rm表示允许缺失数据NA
> xm<-mean(x,na.rm=TRUE)
> xm
[1] 30.14286
>> x<-c(1,10,20,30,40,50,60)
> xm<-mean(x)
> xm
[1] 30.14286
#trim表示去两端数据的百分比
> mean(x,trim=0.5)
[1] 30
2.样本方差与标准差
注意,不是总体方差
> x<-c(11,22,33,44,55,66)
> var(x)
[1] 423.5
> sd(x)
[1] 20.57912
>
注意,不是总体方差
> x<-c(11,22,33,44,55,66)
> var(x)
[1] 423.5
> sd(x)
[1] 20.57912
>
3.中位数,分位数,极差
> x<-c(1,10,20,30,40,50,NA,60)
> xmd<-median(x,na.rm=TRUE)
> xmd
[1] 30
probs可以指定分位
> quantile(x,na.rm=TRUE)
0% 25% 50% 75% 100%
1 15 30 45 60
> quantile(x,probs=c(0.1,0.8),na.rm=TRUE)
10% 80%
6.4 48.0
极差
> x=c(123,124,11,90,20,19,26,51)
> max(x)-min(x)
[1] 113
> r=range(x)
> r
[1] 11 124
> r[2]-r[1]
[1] 113
>
半极差
> r=quantile(x,c(1/4,3/4),names=FALSE)
> r
[1] 19.75 98.25
> r[2]-r[1]
[1] 78.5