第八讲统计PPT讲稿.ppt

《第八讲统计PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第八讲统计PPT讲稿.ppt（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第八讲统计第1页，共38页，编辑于2022年，星期三2.参数估计参数估计（分布参数估计和置信区间）（分布参数估计和置信区间）3.描述性统计描述性统计4.线性模型线性模型5.非线性模型非线性模型6.假设检验假设检验7.多元统计多元统计8.统计绘图统计绘图9.统计工序管理统计工序管理10.实验设计实验设计第2页，共38页，编辑于2022年，星期三q 统计工具箱函数统计工具箱函数 数值计算函数数值计算函数 交互式图形工具函数交互式图形工具函数第3页，共38页，编辑于2022年，星期三q 概率密度函数概率密度函数(probability density function)Y=PDF(NAME,X,A,

2、B,C)其中其中NAME为某种概率分布函数的名称，为某种概率分布函数的名称，A,B,C为参数，为参数，X为给定所需计算的自变量的值，为给定所需计算的自变量的值，Y 为返回的函数值为返回的函数值NAME 可为可为 以下典型的概率密度函数标示以下典型的概率密度函数标示第4页，共38页，编辑于2022年，星期三v beta or Beta,v bino or Binomial,v chi2 or Chisquare,v exp or Exponential,v f or F,v gam or Gamma,v geo or Geometric,v hyge or Hypergeometric,v lo

3、gn or Lognormal,v nbin or Negative Binomial,第5页，共38页，编辑于2022年，星期三v ncf or Noncentral F,v nct or Noncentral t,v ncx2 or Noncentral Chi-squarev norm or Normal,v poiss or Poisson,v rayl or Rayleigh,v t or T,v unif or Uniform,v unid or Discrete Uniform,v weib or Weibull.第6页，共38页，编辑于2022年，星期三 同时同时Matlab还

4、提供了其它专用概率还提供了其它专用概率 密度函数密度函数如如：betapdf(a,b,c)等价于等价于 pdf(beta,a,b,c)normpdf(x,a,b)等价于等价于 pdf(norm,x,a,b)第7页，共38页，编辑于2022年，星期三例例：正态分布函数正态分布函数N(0,1)X=-3:0.01:3Y=normpdf(x,0,1)Plot(x,y)第8页，共38页，编辑于2022年，星期三q 累积分布函数和逆累积分布函数累积分布函数和逆累积分布函数(cumulative distribution function)累积分布函数定义累积分布函数定义：即逆累积分布函数：逆累积分布函数：

5、实际上为累积分布函实际上为累积分布函数的逆函数数的逆函数第9页，共38页，编辑于2022年，星期三调用格式调用格式:P=cdf(NAME,X,A,B,C)X=icdf(NAME,P,A,B,C)其中其中NAME的取值与的取值与pdf函数使用方法相同函数使用方法相同另外也拥有专用另外也拥有专用cdf函数，如函数，如normcdf第10页，共38页，编辑于2022年，星期三q 随机数产生函数随机数产生函数产生的基本方法产生的基本方法（直接法，反演法，拒（直接法，反演法，拒绝法）绝法）产生随机数的函数产生随机数的函数Y=random(NAME,A,B,C,m,n)另外也拥有专用随机数产生函数，如另外

6、也拥有专用随机数产生函数，如betarnd(a,b)第11页，共38页，编辑于2022年，星期三例：y=random(norm,0,1,100,80)mesh(y)第12页，共38页，编辑于2022年，星期三q 均值和方差函数（均值和方差函数（20种专用函数）种专用函数）(mean and variance)M,V=NORMSTAT(MU,SIGMA)M,V=WEIBSTAT(A,B)M,V=BETASTAT(A,B)M,V=BIONSTAT(N,P)M,V=EXPSTAT(Mu)M,V=RAYSTAT(B)等等第13页，共38页，编辑于2022年，星期三q 参数估计参数估计 极大似然参数估计

7、方法极大似然参数估计方法（Maximum Likelihood Estimate:MLE)对数似然参数估计方法对数似然参数估计方法（只对部分分布函数（只对部分分布函数提供）提供）第14页，共38页，编辑于2022年，星期三参数估计函数格式参数估计函数格式phat.pci=mle(NAME,Data,alpha,p1)其中其中NAME与前面定义相同（与前面定义相同（如：如：beta)Data为数据样本为数据样本Alpha为用户给定的置信度为用户给定的置信度P1为实验的次数为实验的次数 第15页，共38页，编辑于2022年，星期三参数估计函数格式参数估计函数格式L=normlike(mu,sigm

8、a,Data)Data为数据样本为数据样本Mu,sigma为用户给定的参数为用户给定的参数 Weiblike,gamlike,betalike等等第16页，共38页，编辑于2022年，星期三q 描述性统计描述性统计 目的目的：搜集：搜集 整理整理 加工和分析统计数加工和分析统计数据，使之系统化据，使之系统化 条理化。从中显示数条理化。从中显示数据资料的趋势据资料的趋势 特征和数量关系。特征和数量关系。内容：内容：位置度量位置度量 散布度量散布度量 自助法和自助法和缺失数据情况下的处理办法。缺失数据情况下的处理办法。第17页，共38页，编辑于2022年，星期三样本的几何均值样本的几何均值M=ge

9、omean(X)其中第18页，共38页，编辑于2022年，星期三样本数据的调和均值样本数据的调和均值m=harmmean(X)其中第19页，共38页，编辑于2022年，星期三样本数据的算术平均值样本数据的算术平均值m=mean(X)其中第20页，共38页，编辑于2022年，星期三样本数据的中值（中位数）样本数据的中值（中位数）m=median(X)其中第21页，共38页，编辑于2022年，星期三剔除极端数据的样本中值剔除极端数据的样本中值m=trimmean(X,percent)其中第22页，共38页，编辑于2022年，星期三q 散布度量散布度量 目的目的：描述样本数据中数据偏离数据中描述样本

10、数据中数据偏离数据中心的程度心的程度内容：内容：标准差标准差(std)和方差（和方差（Var)；平均；平均绝对偏差（绝对偏差（mad)；四分位数间距四分位数间距第23页，共38页，编辑于2022年，星期三样本的标准差（样本的标准差（Standard deviation）m=std(X)其中第24页，共38页，编辑于2022年，星期三样本的协方差样本的协方差(Covariance matrix)m=cov(X，Y)其中第25页，共38页，编辑于2022年，星期三样本数据的平均绝对偏差样本数据的平均绝对偏差(mean absolute derive)m=mad(X)其中第26页，共38页，编辑于2

11、022年，星期三样本数据的内四分位数间距样本数据的内四分位数间距m=iqr(X)其中第27页，共38页，编辑于2022年，星期三样本数据的极差样本数据的极差m=range(X)其中m=max(x)-min(x)第28页，共38页，编辑于2022年，星期三q 缺失数据的处理函数缺失数据的处理函数(Standard deviation ignoring NaNs)Nanmax(x)求缺失数据的样本最大值Nanmin(x)求缺失数据的样本最小值Nanmean(x)求缺失数据的样本平均值Nanstd(x)求缺失数据的样本标准差Nansum(x)求缺失数据的样本的和第29页，共38页，编辑于2022年，

12、星期三q 中心矩（中心矩（关于数学期望的矩关于数学期望的矩）定义定义 格式：格式：m=moment(X,order)第30页，共38页，编辑于2022年，星期三q 相关系数相关系数 定义定义：两个随机变量间的线性相依程度两个随机变量间的线性相依程度的度量的度量格式：格式：R=corrcoef(X,Y)第31页，共38页，编辑于2022年，星期三q 样本峰度和样本偏度样本峰度和样本偏度 定义定义：样本峰度样本峰度为单峰分布曲线峰的平坦程度的为单峰分布曲线峰的平坦程度的度量，度量，k=kurtosid(x)；定义定义：样本偏度样本偏度为样本数据围绕其均值的对称情况的为样本数据围绕其均值的对称情况的

13、度量，度量，k=skewness(x),第32页，共38页，编辑于2022年，星期三q 假设检验假设检验 两个总体样本一致性的两个总体样本一致性的Wilcoxon检验检验 p,h=ranksum(x,y,alpha)其中其中h=0 or 1 p=signrank(x,y,alpha)Wilcoxon signed rank 检验检验h,sig,ci=ttest(x,m,alpha)单一样单一样本均值的本均值的t假设检验假设检验第33页，共38页，编辑于2022年，星期三 h,sig,ci=ttest2(x,y,alpha)两个样本均值差异的检验两个样本均值差异的检验h,sig,ci=ztest

14、(x,m,sigma,alpha)已知方差的单样本均值的假设检验已知方差的单样本均值的假设检验第34页，共38页，编辑于2022年，星期三q 线性模型线性模型方差分析方差分析回归分析回归分析第35页，共38页，编辑于2022年，星期三q 多项式回归多项式回归 多项式的创建：多项式的创建：p=x1,x2,x3,xn 高次高次 低次低次 p,s=polyfit(x,y,n).用用n阶多项式拟合阶多项式拟合x和和y polyval(p,x).多项式多项式p在在x上的值上的值Polytool(x,y)拟合多项式及预测的交互式界面拟合多项式及预测的交互式界面 第36页，共38页，编辑于2022年，星期三q 二次响应曲面工具二次响应曲面工具 rstool(x,y,model,alpha,xname,yname)拟合二次曲面拟合二次曲面交互式界面交互式界面 Model 为为interaction,quadratic,purequadratic 第37页，共38页，编辑于2022年，星期三q 非线性回归模型非线性回归模型 Nlinfit非线性最小二乘拟合非线性最小二乘拟合 Nlintool非线性非线性 拟合拟合GUI Nlparci参数置信区间参数置信区间 Nlpredci预测置信区间预测置信区间 Nnls非负最小二乘拟合非负最小二乘拟合第38页，共38页，编辑于2022年，星期三