Puede crear un gráfico utilizando los siguientes scripts que colocan en un nuevo gráfico la línea de referencia y dos líneas de curva ROC.
Los valores de las curvas ROC se obtienen haciendo clic derecho en el gráfico de la curva ROC que obtuvo en la ventana de análisis y eligiendo "personalizar". Luego, en la ventana de script, obtiene el script de la línea. Puede copiarlos y pegarlos en el script a continuación, reemplazando las líneas que están allí actualmente. Repita esto para cada curva ROC que desee agregar.
Espero que esto ayude.
New Window( "ROC Curve",
Graph Box(
FrameSize( 600, 600 ),
X Scale( 0, 1 ),
Y Scale( 0, 1 ),
//add axis titles
xName( "1-Specificity" ),
yName( "Sensitivity" ),
//add reference line
Line Style( 0 );
Y Function( XTemp , XTemp );
//add first line
Line Style( 1 );
Text Size( 12 );
Text({0.45, 0.2}, "......... Group, BMI; AUC = 0.88");
pen size(2);
Line(
[0, 0, 0, 0, 0, 0, 0, 0.153846153846154, 0.230769230769231, 0.230769230769231,
0.307692307692308, 0.384615384615385, 0.384615384615385, 0.461538461538462,
0.461538461538462, 0.538461538461539, 0.615384615384615, 0.692307692307692,
0.769230769230769, 0.846153846153846, 0.923076923076923, 1, 1],
[0, 0.111111111111111, 0.222222222222222, 0.333333333333333, 0.444444444444444,
0.555555555555556, 0.666666666666667, 0.666666666666667, 0.666666666666667,
0.777777777777778, 0.777777777777778, 0.777777777777778, 0.888888888888889,
0.888888888888889, 1, 1, 1, 1, 1, 1, 1, 1, 1]
);
//add second line
Line Style( 4 );
Text Size( 12 );
Text( {0.45, 0.1}, "_.._.._.._ Group, BMI, Endo; AUC = 0.96" );
Line(
[0, 0, 0, 0, 0, 0, 0, 0.0833333333333333, 0.0833333333333333, 0.166666666666667,
0.25, 0.25, 0.333333333333333, 0.416666666666667, 0.5, 0.583333333333333,
0.666666666666667, 0.75, 0.833333333333333, 1, 1],
[0, 0.222222222222222, 0.333333333333333, 0.444444444444444, 0.555555555555556,
0.666666666666667, 0.777777777777778, 0.777777777777778, 0.888888888888889,
0.888888888888889, 0.888888888888889, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
);
);
);Una forma de hacerlo es utilizar la tabla ROC que existe después de elegir crear el gráfico ROC. Básicamente, para cada prueba, convierta la tabla ROC en una tabla de datos, concatene las tablas utilizando la opción para crear una nueva columna que identifique las tablas de origen y, a continuación, utilice el gráfico de superposición con los gráficos agrupados por tabla de origen y active la superposición por grupos. Puede hacerlo de forma interactiva o con algo como el siguiente ejemplo que armé guardando los scripts después de hacerlo de forma interactiva:
dt = open("$DATOS_DE_MUESTRA/Fitness.jmp");
fm = dt << Ajustar modelo(
Y( :Sexo ),
Efectos( :Oxy ),
Personalidad (logística nominal),
Ejecutar modelo(
Positivo
Nivel("M"),
Pruebas de razón de verosimilitud ( 1 ),
Bosque
Pruebas( 0 ),
República de China
Curva( 1 ),
Logístico
Trama( 1 ),
EnviarAlInforme(
Despacho( {"Receptor operativo
Característica"}, "Tabla ROC", OutlineBox, Close( 0 ) )
)
)
);
fmReport = Informe( fm );
fmReport["Tabla ROC"][Cuadro de tabla( 1 )] << convertir en datos
mesa;
tabla de datos(1) << nombre del conjunto("Oxy");
fm = dt << Ajustar modelo(
Y( :Sexo ),
Efectos( :Tiempo de ejecución ),
Personalidad (logística nominal),
Ejecutar modelo(
Positivo
Nivel("M"),
Pruebas de razón de verosimilitud ( 1 ),
Bosque
Pruebas( 0 ),
República de China
Curva( 1 ),
Logístico
Trama( 1 ),
EnviarAlInforme(
Despacho( {"Receptor operativo
Característica"}, "Tabla ROC", OutlineBox, Close( 0 ) )
)
)
);
fmReport = Informe( fm );
fmReport["Tabla ROC"][Cuadro de tabla( 1 )] << convertir en datos
mesa;
tabla de datos(1) << conjunto nombre("Tiempo de ejecución");
Concatenar(
Tabla de datos ( "Tiempo de ejecución")
Tabla de datos ("Oxy"),
Tabla de salida ( "Ambos")
Crear columna de origen
);
Trazado de superposición(
X( :Nombre( "1-Especificidad" ) ),
Y( :Sensibilidad ),
Agrupación( :Tabla de origen ),
Grupos de superposición,
Paso (1),
Mostrar puntos( 0 ),
EnviarAlInforme(
Despacho(
{},
"106",
Caja de escala( 2 ),
{Formato( "Mejor", 10 ), Mostrar cuadrícula principal(
1 )}
),
Despacho(
{}, "101", ScaleBox( 2 ), {Mostrar cuadrícula principal( 1 )} )
)
);
estilo="margen inferior: 0.0001pt; altura de línea: normal;">
estilo="font-size: 10pt; font-family: "Courier New"; color: azul marino;">
estilo="font-size: 10pt; font-family: "Courier New"; color: negro;">
