统计数据基本特征的把握.pptVIP

* * * * * * * * * * * * * * * data=read.csv(某公司05至07年两项数据月度表.csv,head=TRUE);summary(data) 均值=mean(data[,2:3]);方差=var(data[,2:3])#此命令给出的是样本修正方差 标准差=sd(data[,2:3]) 全距=c(max(data[,2])-min(data[,2]),max(data[,3]-min(data[,3]))) 中位数=c(median(data[,2]),median(data[,3])) 离散系数=标准差/均值 第五节 统计软件使用 data=read.csv(某矿井2006年混煤煤质每日数据.csv,head=TRUE) data=data[,-1] #将回归中用不到的第一列（日期序号）删除 attributes(data) attr(data,names)=c(水分,灰分,硫分,发热量) #由于原数据变量名中含有%致使导入的数据变量名后均含有圆点（看来R不能识别带%的变量名）所以重新为变量命名以便于使用 attach(data) plot.new() figureparameter=par(mfrow=c(2,2),bg=light green,cex.lab=1,font=3) #设定图形参数：画2行2列的散点图（背景为light green，主标题大小1.5磅，坐标轴标签大小1磅，字体为粗体） hist(水分);hist(灰分);hist(硫分);hist(发热量) plot.new() figureparameter=par(mfrow=c(2,2),bg=light green,cex.lab=1,font=3) #设定图形参数：画2行2列的散点图（背景为light green，主标题大小1.5磅，坐标轴标签大小1磅，字体为粗体） plot(硫分,pch=4,col=red) #画硫分的散点图（图标为4，颜色为red） plot(水分,pch=8,col=red) #画水分的散点图（图标为8） plot(灰分,pch=23,col=red,bg=yellow) #画灰分的散点图（图标为23，颜色为red，图标的填充色（背景色）为yellow） plot(发热量,pch=23,col=red,bg=yellow) #画发热量的散点图（图标为23，颜色为red，图标的填充色（背景色）为yellow） 第五节 统计软件使用 x=as.vector(data$发热量);summary(x) 单价=numeric(length(x)) #建立单价作为相应发热量混煤的单价并设定单价的属性(应注意numeric与as.numeric不同) for(i in 1:length(x)){ if(x[i]=5900) 单价[i]=105 #发热量5900以上的单价为105 if(5900x[i]x[i]=5800) 单价[i]=102 #发热量5800—5900的单价为102 if(5800x[i]x[i]=5700) 单价[i]=99 #发热量5700—5800的单价为99 if(5700x[i]x[i]=5600) 单价[i]=96 #发热量5600—5700的单价为96 if(5600x[i]x[i]=5500) 单价[i]=93 #发热量5500—5600的单价为93 if(x[i]5500) 单价[i]=88 #发热量5500的单价为88 } 等级=numeric(length(x)) #建立向量以便把不同发热量的煤分成不同等级 for(i in 1:length(x)){ if(x[i]=5900) 等级[i]=1 #发热量5900以上的等级为1 if(5900x[i]x[i]=5800) 等级[i]=2 #发热量5800—5900的等级为2 if(5800x[i]x[i]=5700) 等级[i]=3 #发热量5700—5800的等级为3 if(5700x[i]x[i]=5600) 等级[i]=4 #发热量5600—5700的等级为4 if(5600x[i]x[i]=5500) 等级[i]=5 #发热量5500—5600的等级为5 if(x[i]5500) 等级[i]=6 #发热量5500的等级为6 } data=cbind(data,单价,等级);data 第五节 统计软件使用 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 未分组数据四分位数位置的确定： 下四分位数(QL)位置： N+1 4 上四分