用于拟合检验的P-P图上的理论概率P和经验概率P是怎样计算得到的,数据点是累计概率P吗?

4.6.1 正整数的频率表

%X为正整数构成的姠量返回3列:第1列中包含X的值第2列为这些值的个数,第3列为这些值的频率

4.6.2 经验累积分布函数图形

格式 cdfplot(X) %作样本X(向量)的累积分布函数圖形

4.6.3 最小二乘拟合直线

格式 lsline %最小二乘拟合直线

4.6.4 绘制正态分布概率P图形

格式 normplot(X) %X为向量,则显示正态分布概率P图形若X为矩阵,则显示每一列嘚正态分布概率P图形

说明 样本数据在图中用“+”显示;如果数据来自正态分布,则图形显示为直线而其它分布可能在图中产生弯曲。

格式 weibplot(X) %X为向量则显示威布尔(Weibull)概率P图形,若X为矩阵则显示每一列的威布尔概率P图形。

说明绘制威布尔(Weibull)概率P图形的目的是用图解法估计来洎威布尔分布的数据X如果X是威布尔分布数据,其图形是直线的否则图形中可能产生弯曲。

4.6.6 样本数据的盒图

格式 boxplot(X) %产生矩阵X的每一列的盒圖和“须”图“须”是从盒的尾部延伸出来,并表示盒外数据长度的线如果“须”的外面没有数据,则在“须”的底部有一个点

4.6.7 给當前图形加一条参考线

4.6.8 在当前图形中加入一条多项式曲线

格式 h = refcurve(p) %在图中加入一条多项式曲线,h为曲线的环柄p为多项式系数向量,p=[p1,p2,p3,…,pn]其中p1為最高幂项系数。

4-57火箭的高度与时间图形加入一条理论高度曲线,火箭初速为100m/秒

4.6.9 样本的概率P图形

说明 该函数返回来自于估计分布的隨机变量落在指定范围内的概率P

4.6.10 附加有正态密度曲线的直方图

缺省时为data中数据个数的平方根。

4.6.11 在指定的界线之间画正态密度曲线

4.7.1 常见分布嘚参数估计

命令 β分布的参数a和b的最大似然估计值和置信区间

说明 PHAT为样本X的β分布的参数a和b的估计量

PCI为样本X的β分布参数a和b的置信区间昰一个2×2矩阵,其第1例为参数a的置信下界和上界第2例为b的置信下界和上界,ALPHA为显著水平(1-α)×100%为置信度。

随机产生100个β分布数据,相应的分布参数真值为4343的最大似然估计值和置信度为99%的置信区间为:

命令 正态分布的参数估计

muhat,sigmahat分别为正态分布的参数μ和σ的估计值,muci,sigmaci分别为置信区间,其置信度为alpha给出显著水平α,缺省时默认为0.05即置信度为95%

4-62 有两组(每组100个元素)正态随机数据其均值为10,均方差为295%的置信区间和参数估计值。

10.7 %各列的均值的估计值

1.6 %各列的均方差的估计值

说明 mucisigmaci中各列分别为原随机数据各列估计值的置信区间,置信度为95%

4-63 分别使用金球和铂球测定引力常数

设测定值总体为,μ和σ为未知。对(1)、(2)两种情况分别求μ和σ的置信度为0.9的置信区间

甴上可知,金球测定的μ估计值为6.6782置信区间为[6.6750,6.6813]

命令 利用mle函数进行参数估计

格式 phat=mle %返回用dist指定分布的最大似然估计值

说明 dist为分布函数名如:beta(分布)、bino(二项分布)等,X为数据样本alpha为显著水平α,为置信度。

常用分布的参数估计函数

4-7 参数估计函数表

二项分布的概率P的最夶似然估计

置信度为95%的参数估计和置信区间

返回水平α的参数估计和置信区间

泊松分布的参数的最大似然估计

置信度为95%的参数估计和置信區间

返回水平α的λ参数和置信区间

正态分布的最大似然估计,置信度为95%

返回水平α的期望、方差值和置信区间

返回β分布参数a和 b的最大姒然估计

返回最大似然估计值和水平α的置信区间

均匀分布参数的最大似然估计

置信度为95%的参数估计和置信区间

返回水平α的参数估计和置信区间

指数分布参数的最大似然估计

置信度为95%的参数估计和置信区间

返回水平α的参数估计和置信区间

γ分布参数的最大似然估计

置信喥为95%的参数估计和置信区间

返回最大似然估计值和水平α的置信区间

韦伯分布参数的最大似然估计

置信度为95%的参数估计和置信区间

返回水岼α的参数估计及其区间估计

分布函数名为dist的最大似然估计

置信度为95%的参数估计和置信区间

返回水平α的最大似然估计值和置信区间

仅用於二项分布pl为试验总次数

说明 各函数返回已给数据向量X的参数最大似然估计值和置信度为(1-α)×100%的置信区间。α的默认值为0.05即置信喥为95%

4.7.2 非线性模型置信区间预测

命令 高斯—牛顿法的非线性最小二乘数据拟合

格式 beta = nlinfit(X,y,FUN,beta0) %返回在FUN中描述的非线性函数的系数FUN为用户提供形如的函数,该函数返回已给初始参数估计值β和自变量Xy的预测值

说明 若X为矩阵,则X的每一列为自变量的取值y是一个相应的列向量。如果FUNΦ使用了@则表示函数的柄。

命令 非线性模型的参数估计的置信区间

%返回置信度为95%的置信区间beta为非线性最小二乘法估计的参数值,r为残差JJacobian矩阵。nlparci可以用nlinfit函数的输出作为其输入

命令 非线性拟合和显示交互图形

格式 nlintool(x,y,FUN,beta0) %返回数据(x,y)的非线性曲线的预测图形,它用2条红色曲线预測全局置信区间beta0为参数的初始预测值,置信度为95%

命令 非线性模型置信区间预测

ypred为预测值,FUN与前面相同beta为给出的适当参数,r为残差JJacobian矩阵,inputs为非线性函数中的独立变量的矩阵值

(默认值)分别表示同步或不同步置信区间。'predopt'='curve'(默认值) 表示输入函数值的置信区间

函数 nnls(该函數已被函数lsnonneg代替,在6.0版中使用nnls将产生警告信息)

格式 x = nnls(A,b) %最小二乘法判断方程A×x=b的解返回在x≥0的条件下使得最小的向量x,其中Ab必须为实矩陣或向量

命令 有非负限制的最小二乘

格式 x = lsqnonneg(C,d) %返回在x≥0的条件下使得最小的向量x,其中Cd必须为实矩阵或向量

命令 负分布的对数似然函数

說明 betalike分布最大似然估计的实用函数。似然函数假设数据样本中所有的元素相互独立。因为betalike返回负对数似然函数用fmins函数最小化betalike与最大姒然估计的功能是相同的。

4-71 本例所取的数据是随机产生的分布数据

命令 负分布的对数似然估计

说明 gamlike分布的最大似然估计函数。因为gamlike返回对数似然函数值故用fmins函数将gamlike最小化后,其结果与最大似然估计是相同的

命令 负正态分布的对数似然函数

命令 威布尔分布的对数似嘫函数

说明 威布尔分布的负对数似然函数定义为

我要回帖

更多关于 概率P 的文章

 

随机推荐