184128_zh-CN

liliconcom

考克斯-斯内尔 R2

大家好，

我正在估计逻辑回归并且需要得到考克斯-斯内尔右²它隐藏在任何地方还是我必须手动计算它？有人手边有脚本吗？

默认情况下，JMP 仅提供 McFadden 和 Nagelkerke。

这些是我遵循的定义：

https://stats.idre.ucla.edu/other/mult-pkg/faq/general/faq-what-are-pseudo-r-squareds/

 

谢谢你！

罗恩

回复：考克斯-斯内尔 R2谢谢马克。工作完美。

罗恩
回复：考克斯-斯内尔 R2

毕竟剧本还不错。 此示例展示了一种实现方法：

 

Names Default to Here( 1 );

// example case: Big Class
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );

// fit binary logistic model
lr = dt << Fit Model(
 Y( :sex ),
 Effects( :age, :height, :weight ),
 Personality( "Nominal Logistic" ),
 Run( Likelihood Ratio Tests( 1 ), Wald Tests( 0 ) )
);

// access Whole Model Test report
lr rep = (lr << Report)["Whole Model Test"];

// obtain likelihood for reduced and full models
L intercept = Exp( -(lr rep[NumberColBox(1)] << Get( 3 )) );
L full      = Exp( -(lr rep[NumberColBox(1)] << Get( 2 )) );
N           = lr rep[NumberColBox(8)] << Get( 1 );

// compute Cox-Snell R square
r sqr cs    = 1 - (L intercept/L full)^(2/N);

// copy original report
label = lr rep[StringColBox(2)] << Insert Row( 2, { "RSquare (C-S)" } ) << Clone Box;
r sqr u = lr rep[NumberColBox(5)] << Get( 1 );
aicc    = lr rep[NumberColBox(6)] << Get( 1 );
bic     = lr rep[NumberColBox(7)] << Get( 1 );

// replace with new report
lr rep[Tablebox(2)] << Delete;
lr rep << Append(
 TableBox(
  label,
  Number Col Box( "", { r sqr u, r sqr cs, aicc, bic, N } )
 )
);

回复：考克斯-斯内尔 R2

您必须使用您引用的公式根据整个模型测试报告中提供的结果来计算它。 我不知道用于此计算的脚本，但您可以搜索此社区的文件交换区域。 不知道是否需要脚本。 公式非常简单。