美赛中的统计方法简介.pptVIP

余吉昌 yujc@znufe.edu.cn 中南财经政法大学统计与数学学院 数字特征---协方差、相关系数 协方差、相关系数 cov(x) cov(x,y) corrcoef(x,y) corrcoef(X) 常见分布中参数的极大似然估计 函数名以 ‘fit’ 结尾 正态总体的均值和标准差的极大似然估计和置信区间 Normfit(x,alpha) Example: x=[15.14 14.81 15.11 15.26 15.08 15.17 15.12 14.95 15.05 14.97] [muhat, sigmahat,muci,sigmaci]=normfit(x,0.05) 常见分布中参数的极大似然估计 函数名 调 用 形 式 函 数 说 明 binofit [PHAT, PCI] = binofit (X, N, ALPHA) 返回 水平的参数估计和置信区间 poissfit [LAMBDAHAT, LAMBDACI]= poissfit (X,alpha) 返回 水平的 参数和置信区间 unifit [AHAT, BHAT, ACI, BCI] = unifit (X, ALPHA) 返回 水平的参数估计和置信区间 expfit [MUHAT, MUCI] = expfit (X, ALPHA) 返回 水平的参数估计和置信区间 mle [PHAT, PCI] = mle (DIST, DATA, ALPHA) 返回最大似然估计值和 水平的置信区间 对数似然函数值 [logL, info] = betalike (params, data) logL 返回 分布对数似然函数负值。Info为Fisher信息阵，其对角线元素为相应参数的渐进方差 normlike: 正态分布负对数似然函数 gemlike： 分布负对数似然函数 ? 正态总体参数的检验 总体标准差已知时的单个正态总体均值的U检验 某切割机正常工作时，切割的金属棒的长度服从正态分布N(100,4)，从一批金属棒中随机抽取15根，测得它们的长度如下; X=[97 102 105 112 99 103 102 94 100 95 105 98 102 100 103] 检验 ? ? 正态总体参数的检验 [h,p,muci,zval]=ztest(x,100,2,0.05) ：双边检验 [h,p,muci,zval]=ztest(X,100,2,0.05,right) ? 总体标准差未知时单个正态总体的参数的检验 X=[49.4 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9] [h,p,muci,stats]=ttest(X,50,0.05) 矩阵元素赋值 大矩阵可以把小矩阵作为其元素 例： A=[A ; 11 12 13] 在原矩阵的下方加一行 如何在原矩阵的右边添加一列？ Matlab---数据操作 定义矩阵：直接输入法 矩阵用方括号 “[ ]” 括起 例： A = [1 2 3; 4 5 6; 7 8 9] 矩阵同一行中的元素之间用 空格 或 逗号 分隔 矩阵行与行之间用 分号 分开 直接输入法中，分号可以用 回车 代替 例： Matlab 的操作对象是 矩阵 ：命令提示符，不用输入 回车 ：运行所输入的命令 矩阵元素的引用 单个元素的引用 例： A(2,3) 多个元素的引用：冒号的特殊用法 利用小括弧和元素所在的位置(下标) x ( i ) ：向量 x 中的第 i 个元素 A ( i, j ) ：矩阵 A 中的第 i 行，第 j 列元素 a:b:c 产生一个由等差序列组成的向量； a 是首项，b 是公差，c 确定最后一项；若 b=1，则 b 可以省略。 矩阵元素的引用 例： x(1:3) A(3,1:3) A(i:j, m:n) 表示由矩阵 A 的第 i 到第 j 行和第 m 到第 n 列交叉线上的元素组成的子矩阵。 可利用冒号提取矩阵 的整行或整列。 例： A(1, :) A(:, 1:3) A(:, :) 4．建立矩阵的函数 常用函数有： eye(size(A)) 产生与A矩阵同阶的单位矩阵 zeros(m,n) 产生0矩阵 ones(m,n) 产生1矩阵 rand (m,n) 产生随机元素的矩阵 Size(a) 返回包含两个元素的向量。 Length(a) 返回向量的长度。 常见矩阵生成函数 zeros(m,n) 生成一个 m 行 n 列的零矩阵，m=n 时可简写为 zeros(n) ones(m,n)