求教 sas拟合多水平模型初始值设定问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位好：
    我是一名统计新人，最近在学习使用sas中 proc nlmixed拟合多水平logistic模型，该命令中parms语句能够为模型中各参数设定初始值，我的问题如下：

1.模型中设定初始值的意义是什么，有什么具体作用。

2.该设为多少，parms默认是设置为1，但有人说这会影响估计的精度。

3.以前拟合多水平logistic模型使用MLwin，很方便不涉及初始值设定问题，只要设定初次拟合MQL，进而PQL,即可得到结果，但MLwin不能提供-2log 似然值，无法进行嵌套模型比较。

我想能否这样，用MLWIN拟合的参数估计结果作为初始值，在proc mlmixed中  parms语句对各待估参数进行初始值设定。另一方面，将parms参数初始值都设定为1，然后比较哪个结果更接近mlwin。

4.MLWIN中多水平logistic模型算法，与proc nlmixed是否一致。

       不知这方法是否可行，请各位不吝赐教，谢谢。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：多水平模型 多水平 初始值 proc nlmixed Logistic模型 模型 统计 影响

这个影响有多大，正常情况下应该是从0开始，因为你应该假设他们为0，根据数据来推翻这个假设

那么设定初始值，可否利用mlwin拟合的结果呢

freetest 发表于 2010-10-18 10:26
各位好：
    我是一名统计新人，最近在学习使用sas中 proc nlmixed拟合多水平logistic模型，该命令中parms语句能够为模型中各参数设定初始值，我的问题如下：

1.模型中设定初始值的意义是什么，有什么具体作用。

2.该设为多少，parms默认是设置为1，但有人说这会影响估计的精度。

3.以前拟合多水平logistic模型使用MLwin，很方便不涉及初始值设定问题，只要设定初次拟合MQL，进而PQL,即可得到结果，但MLwin不能提供-2log 似然值，无法进行嵌套模型比较。

我想能否这样，用MLWIN拟合的参数估计结果作为初始值，在proc mlmixed中  parms语句对各待估参数进行初始值设定。另一方面，将parms参数初始值都设定为1，然后比较哪个结果更接近mlwin。

4.MLWIN中多水平logistic模型算法，与proc nlmixed是否一致。

       不知这方法是否可行，请各位不吝赐教，谢谢。

It depents on your original problem. If your problem has one and one only maximum value. The initial values only affect the speed of convergence. If your problem size is big, then a proper initial value setting (using your prior information about parameters) will save your big time.

On the other hand, if the problem is highly nonlinear, it may have many LOCAL maximum points. Different initial value may converge to different local maximum instead of a global one. Search for global maximum value is really cumbersome.

感谢楼上指点。但我不明白您所说的many LOCAL maximum points 指什么。
我的目的是分析，地区和个人经济水平对一种疾病患病情况的影响。结局变量为疾病患病率，高水平自变量为地区财政收入，低水平自变量为个人收入、文化程度。使用sas分析主要是为了得到-2log 似然值进行嵌套比较。

freetest 发表于 2010-10-19 13:51
感谢楼上指点。但我不明白您所说的many LOCAL maximum points 指什么。
我的目的是分析，地区和个人经济水平对一种疾病患病情况的影响。结局变量为疾病患病率，高水平自变量为地区财政收入，低水平自变量为个人收入、文化程度。使用sas分析主要是为了得到-2log 似然值进行嵌套比较。

You may think about a highly nonlear problem as you face many mountains, what you need is to find highest  one among them. The numerical algorithm may find a peak(one) and stop there because the first order derivative is zero and the second order derivative is negative, but that is not the highest one.  That is what I mean a LOCAL maximum.

However some problems only have one peak/mountain, for example a simple logistic model. In this case the solution is global.

感谢您的指点，如您所说可能会出现多个LOCAL maximum points ，那么如何确定哪个是最适合的结果呢。

freetest 发表于 2010-10-20 20:02
感谢您的指点，如您所说可能会出现多个LOCAL maximum points ，那么如何确定哪个是最适合的结果呢。

If your problem is really complicated and may have many LOCAL maximum points, then you may try different sets of initial values and pick one with the highest log-likelihood value.

那么我能否利用mlwin拟合的结果作为proc nlmixed  parms语句中的初始值进行拟合呢？
我曾想尝试：分别用0，1，mlwin结果分别作为初始值，进行比较。

freetest 发表于 2010-10-21 11:25
那么我能否利用mlwin拟合的结果作为proc nlmixed  parms语句中的初始值进行拟合呢？
我曾想尝试：分别用0，1，mlwin结果分别作为初始值，进行比较。

Why not.