最大似然威布尔模型（nlpse06）
/***************************************************************//* *//*S A S S A M P L E L I B R A R Y*//* *//*名称：nlpse06*//*标题：最大似然威布尔模型（nlpse06）*//*产品：或*//*系统：所有*//*关键字：或*//*程序：OPTMODEL*//*数据：*//* *//*支持：更新：*//*参考：*//*MISC：非线性规划求解器中的示例6*//*数学规划章节*//* *//***************************************************************/数据枪；输入天数cens@@；数据线；143 0 164 0 188 0 188 0190 0 192 0 206 0 209 0213 0 216 0 220 0 227 0230 0 234 0 246 0 265 0304 0 216 1 244 1;过程操作模型；设置OBS；num天{OBS}；num cens{OBS}；将数据插入OBS=[_N_]天cens；var sig>=1.0e-6初始值10；var c>=1.0e-6初始值10；var theta>=0<=min{i在OBS:cens[i]=0}天[i]init 10；OBS中的impvar fi｛i｝=（如果cens[i]=0，则log（c）-c*log（sig）+（c-1）*log（天[i]-θ）)-（（天[i]-θ）/sig）^c；max logf=OBS中的总和{i}fi[i]；设置VARS=1.._NVAR_；VARS中的num mycov{i，1..i}中的j；用NLP/covest=（cov=2 covot=mycov）求解；打印符号cθ；打印支原体；从VARS中的[ij]={i、1..i}中的j创建数据covdatavar_i=_var_[i].name var_j=_var _[j].name mycov；退出；