一、分析单变量常用的函数
用R 自带的数据BJsales
#基本数据查看
#均值
View(BJsales)
mean(BJsales)
mean(BJsales,trim = 0.1)
#中位数
median(BJsales)
#分位数
quantile(BJsales)
quantile(BJsales,c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9))
#分布图
par(mfrow= c(2,2))
hist(BJsales)
#离散度
#极值
range(BJsales)
#方差(方差越大越离散)平方
var(BJsales)
#标准差(方差的标准化)开平方
sd(BJsales)
#标准化(标准化以后可以比较两个变量的方差和标准差)
scale(BJsales)
#标准误(和标准差相比消除了长度的影响)
error <- sqrt(var(BJsales)/length(BJsales))
error
#偏度峰度函数
mystats <- function(x,na.omit=F){
if (na.omit) x <- x[!is.na(x)]
m <- mean(x)
n <- length(x)
s <- sd(x)
skew <- sum((x-m)^3/s^3)/n
kurt <- sum((x-m)^4/s^4)/n-3
return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))
}
mystats(BJsales)
二、R 提供的统计汇总函数
用R 自带的数据attenu
#查看数据集(地震观测站的信息)
#event:编号mag:级数station:站台号dist: