代码:
DATA:
FILE IS C:\Users\hp\Desktop\LPA.dat;
VARIABLE:
MISSING ARE ALL (-99);
NAMES ARE NUM ql bql bh pg;
USEVARIABLES ARE QL BQL BH PG;
CLASSES = c (5) ; ! 设定潜类别个数,从1个类别开始,依次增加;
ANALYSIS : TYPE = MIXTURE ;
MITERATIONS=5000
STARTS = 1000 250; ! 避免局部最大化解,增加随机起始值数
LRTSTARTS= 0 0 500 200;
OUTPUT: TECH11 TECH14;
SAVEDATA: FILE = zc5.txt;
Save = cprob ; ! 保存后验分类概率
Plot: ! 通过绘图命令 可以获得描述性统计图和条件概 示意图
type is plot3 ;
series = QL BQL BH PG ( * );
比书里多的两行代码是因为下面的报错,修改后报错就没有了。ps.不管有没有这行代码tech11结果基本不变
WARNING: OF THE 5 BOOTSTRAP DRAWS, 4 DRAWS HAD BOTH A SMALLER LRT VALUE THAN THE
OBSERVED LRT VALUE AND NOT A REPLICATED BEST LOGLIKELIHOOD VALUE FOR THE 5-CLASS MODEL.
THIS MEANS THAT THE P-VALUE MAY NOT BE TRUSTWORTHY DUE TO LOCAL MAXIMA.
INCREASE THE NUMBER OF RANDOM STARTS USING THE LRTSTARTS OPTION.
tech11结果:TECHNICAL 11 OUTPUT
TECHNICAL 11 OUTPUT
Random Starts Specifications for the k-1 Class Analysis Model
Number of initial stage random starts 20
Number of final stage optimizations 4
VUONG-LO-MENDELL-RUBIN LIKELIHOOD RATIO TEST FOR 4 (H0) VERSUS 5 CLASSES
H0 Loglikelihood Value -9690.979
2 Times the Loglikelihood Difference 11603.859
Difference in the Number of Parameters 6
Mean 156127.952
Standard Deviation 221672.534
P-Value 1.0000
LO-MENDELL-RUBIN ADJUSTED LRT TEST
Value 11393.281
P-Value 1.0000
tech14结果
TECHNICAL 14 OUTPUT
Random Starts Specifications for the k-1 Class Analysis Model
Number of initial stage random starts 20
Number of final stage optimizations 4
Random Starts Specification for the k-1 Class Model for Generated Data
Number of initial stage random starts 0
Number of final stage optimizations for the
initial stage random starts 0
Random Starts Specification for the k Class Model for Generated Data
Number of initial stage random starts 1000
Number of final stage optimizations 500
Number of bootstrap draws requested Varies
PARAMETRIC BOOTSTRAPPED LIKELIHOOD RATIO TEST FOR 4 (H0) VERSUS 5 CLASSES
H0 Loglikelihood Value -9690.979
2 Times the Loglikelihood Difference 11603.859
Difference in the Number of Parameters 6
Approximate P-Value 0.0000
Successful Bootstrap Draws 5
请问为什么p值会等于1,5分类模型除了tech11结果之外其他指标和理论解释性都很好,是否可以接受其作为最优模型,如果可以在论文中应该如何解释p=1以及可以选定其作为最优模型?
感谢各位大佬!
代码:
DATA:
FILE IS C:\Users\hp\Desktop\LPA.dat;
VARIABLE:
MISSING ARE ALL (-99);
NAMES ARE NUM ql bql bh pg;
USEVARIABLES ARE QL BQL BH PG;
CLASSES = c (5) ; ! 设定潜类别个数,从1个类别开始,依次增加;
ANALYSIS : TYPE = MIXTURE ;
MITERATIONS=5000
STARTS = 1000 250; ! 避免局部最大化解,增加随机起始值数
LRTSTARTS= 0 0 500 200;
OUTPUT: TECH11 TECH14;
SAVEDATA: FILE = zc5.txt;
Save = cprob ; ! 保存后验分类概率
Plot: ! 通过绘图命令 可以获得描述性统计图和条件概 示意图
type is plot3 ;
series = QL BQL BH PG ( * );
比书里多的两行代码是因为下面的报错,修改后报错就没有了。ps.不管有没有这行代码tech11结果基本不变
WARNING: OF THE 5 BOOTSTRAP DRAWS, 4 DRAWS HAD BOTH A SMALLER LRT VALUE THAN THE
OBSERVED LRT VALUE AND NOT A REPLICATED BEST LOGLIKELIHOOD VALUE FOR THE 5-CLASS MODEL.
THIS MEANS THAT THE P-VALUE MAY NOT BE TRUSTWORTHY DUE TO LOCAL MAXIMA.
INCREASE THE NUMBER OF RANDOM STARTS USING THE LRTSTARTS OPTION.
tech11结果:TECHNICAL 11 OUTPUT
TECHNICAL 11 OUTPUT
Random Starts Specifications for the k-1 Class Analysis Model
Number of initial stage random starts 20
Number of final stage optimizations 4
VUONG-LO-MENDELL-RUBIN LIKELIHOOD RATIO TEST FOR 4 (H0) VERSUS 5 CLASSES
H0 Loglikelihood Value -9690.979
2 Times the Loglikelihood Difference 11603.859
Difference in the Number of Parameters 6
Mean 156127.952
Standard Deviation 221672.534
P-Value 1.0000
LO-MENDELL-RUBIN ADJUSTED LRT TEST
Value 11393.281
P-Value 1.0000
tech14结果
TECHNICAL 14 OUTPUT
Random Starts Specifications for the k-1 Class Analysis Model
Number of initial stage random starts 20
Number of final stage optimizations 4
Random Starts Specification for the k-1 Class Model for Generated Data
Number of initial stage random starts 0
Number of final stage optimizations for the
initial stage random starts 0
Random Starts Specification for the k Class Model for Generated Data
Number of initial stage random starts 1000
Number of final stage optimizations 500
Number of bootstrap draws requested Varies
PARAMETRIC BOOTSTRAPPED LIKELIHOOD RATIO TEST FOR 4 (H0) VERSUS 5 CLASSES
H0 Loglikelihood Value -9690.979
2 Times the Loglikelihood Difference 11603.859
Difference in the Number of Parameters 6
Approximate P-Value 0.0000
Successful Bootstrap Draws 5
请问为什么p值会等于1,5分类模型除了tech11结果之外其他指标和理论解释性都很好,是否可以接受其作为最优模型,如果可以在论文中应该如何解释p=1以及可以选定其作为最优模型?
感谢各位大佬!