数据分析实验报告

数据分析》实验报告班级：07信计0班学号：姓名：实验日期2010-3-11实验地点：_实验楼505实验名称：样本数据的特征分析使用软件名称：MATLAB实

验

目实

验

目

的2.熟练掌握jbtest与lillietest关于一元数据的正态性检验；3•掌握统计作图方法；4.掌握多元数据的数字特征与相关矩阵的处理方法;安徽省1990-2004年万元工业GDP废气排放量、废水排放量、固体废物排放量以及用于污染治理的投入经费比重见表6.1.1,解决以下问题：表6.1.1废气'废水'固体废物排放量及污染治理的投入经费占GDP比重年份万元工业GDP废气排放量万元工业GDP固体物排放量万元工业GDP废水排放量环境污染治理投资占GDP比重（立方米）（千克）（吨）(%)1990104254.40519.48441.65000476.97398.190.26实199289317.41119.45332.140.23验199363012.4267.93203.910.20内容199445435.047.86128.2004212.45113.390.22199639874.1913.2487.1208537.9776.980.21199835270.7945.3659.680.11199935200.7634.9360.820971.8257.350.19200140348.431.1753.060.11200240392.960.1650.960.12200337237.130.0543.940270.0636.900.131•计算各指标的均值、方差、标准差、变异系数以及相关系数矩阵；2•计算各指标的偏度、峰度、三均值以及极差；3•做出各指标数据直方图并检验该数据是否服从正态分布？若不服从正态分布，利用boxcox变换以后给出该数据的密度函数；4.上网查找1990-2004江苏省万元工业GDP废气排放量，安徽省与江苏省是否服从同样的分布？程序如下：clear;clcformatbank%保留两位小数%%%%%%%%%%%%%%%%安徽省数据%%%%%%%%%%%%%%%%%%A=[104254.40519.48441.650.1894415.00476.97398.190.2689317.41119.45332.140.2363012.4267.93203.910.2045435.047.86128.200.1746383.4212.45113.390.2239874.1913.2487.120.1538412.8537.9776.980.2135270.7945.3659.680.1135200.7634.9360.820.1535848.971.8257.350.1940348.431.1753.060.1140392.960.1650.960.1237237.130.0543.940.1534176.270.0636.900.13]%计算各指标的均值、方差、标准差、变异系数、偏度、峰度以及极差A1=[mean(A);var(A);std(A);std(A)./mean(A);skewness(A,O);kurtosis(A,0)-3;range(A)]%三均值A2=[1/41/21/4]*prctile(A,[255075])%计算各指标的相关系数矩阵A3=corrcoef(A)%做出各指标数据直方图subplot(221),histfit(A(:,1),8)subplot(222),histfit(A(:,2),8)subplot(223),histfit(A(:,3),8)subplot(224),histfit(A(:,4),7)%检验该数据是否服从正态分布fori=1:4[h(i),p(i),lstat(i),cv(i)]=lillietest(A(:,i),0.05);endh,p%计算出前三列不服从正态分布，利用boxcox变换以后给出该数据的密度函数[t1,I1]=boxcox(A(:,1))[t2,I2]=boxcox(A(:,2))[t3,I3]=boxcox(A(:,3))%%%%%%%%%%%%%%%%江苏省数据%%%%%%%%%%%%%%%%%J=[5047.0024.232234.005730.0024.002077.005352.0022.492450.005202.0021.162323.005772.1421.162725.007872.1122.022883.007450.8121.972891.008164.2422.012857.007487.8520.292944.008353.9320.112896.009078.2020.193038.0013343.9927.103553.0014286.4626.303796.0014632.6924.753893.7017818.4226.354672.5320196.5829.635757.37]%检验该数据是否服从正态分布fori=1:3[h1(i),p1(i),lstat1(i),cv1(i)]=lillietest(J(:,i),0.05);endh1,p1%计算出第一、三列不服从正态分布，利用boxcox变换以后给出该数据的密度函数[t11,I11]=boxcox(J(:,1))[t33,I33]=boxcox(J(:,3)1•各指标的均值、方差、标准差、变异系数、偏度、峰度、极差和三均值如下表：指标数字特征万元工业GDP废气排放量(m3)万元工业GDP固体物排放(kg)万元工业GDP废水排放量(t)环境污染治理投资占GDP比重(％)均值51972.0089.26142.950.17方差577053049.9028703.4718645.500.00标准差24021.93169.42136.550.05变异系数0.461.900.960.27偏度1.442.251.440.29

4•经检验可知：江苏省万元工业GDP废气排放量、固体废物排放量数据不服从正态分布，而废水排放量服从正态分布。附：表1990-2005年万元工业GDP废气排放量、废水排放量、固体废物排放量时间工业废气排放量（亿立方米）工业废水排放量（亿吨）工业固体废物排放量（万吨）19905047.0024.232234.0019915730.0024.002077.0019925352.0022.492450.0019935202.0021.162323.0019945772.1421.162725.0019957872.1122.022883.0019967450.8121.972891.0019978164.2422.012857.0019987487.8520.292944.0019998353.9320.112896.00200090