6# basswawa
我都有些佩服你了。
你是想证明sas例题错了?
在proc nlp里,programming statement和其他的statements有timing的问题(谁前谁后)。但其运行机制和在data step中并无大异。有关这些问题好像在proc optmodel(?)的帮助文件有较为翔实的说明。
每一个观察根据是否censored而产生一个log likelihood的描述(s),最后在max 里循环出其和的极值及各个参数的值。
我没有很理解你到底纠结在哪里?
另外那个logf没有实际意义。去掉logf,然后max s;结果应该是一样的,正如我在中文论坛上给出的code。
如果你不相信,你在proc NLP statement 里加上out=选项,比较一下logf 和s,他们应该是一模一样的。