Ejemplo de PROC NLIN con Perfilado de Parámetros y Bootstrap

El script ajusta un modelo exponencial a un conjunto de datos. Utiliza el procedimiento NLIN no solo para estimar los parámetros del modelo, sino también para evaluar su estabilidad y sensibilidad. La declaración PROFILE se usa para examinar la verosimilitud del perfil de dos de los parámetros, mientras que la declaración BOOTSTRAP realiza un análisis de remuestreo para proporcionar intervalos de confianza robustos y diagnósticos gráficos sobre los estimadores.

Análisis de datos

Type : CREATION_INTERNE

El conjunto de datos 'clarke1987a' se crea directamente en el script usando un paso DATA y una instrucción 'datalines'.

1 Bloque de código

DATA STEP Data

Explicación :
Este bloque de código crea la tabla SAS 'clarke1987a' en la librería WORK. La tabla contiene dos variables numéricas, 'x' e 'y', cuyos datos se leen directamente desde el flujo de entrada a través de la instrucción 'datalines'.

¡Copiado!

1	DATA clarke1987a;
2	INPUT x y;
3	DATALINES;
4	1 3.183
5	2 3.059
6	3 2.871
7	4 2.622
8	5 2.541
9	6 2.184
10	7 2.110
11	8 2.075
12	9 2.018
13	10 1.903
14	11 1.770
15	12 1.762
16	13 1.550
17	;
18

2 Bloque de código

PROC NLIN

Explicación :
Este bloque activa la generación de gráficos ODS, luego ejecuta el procedimiento de regresión no lineal (NLIN). Define los valores iniciales de los parámetros (parms), especifica un modelo exponencial (model), y solicita dos análisis avanzados: un perfilado de la verosimilitud para 'theta1' y 'theta3' (profile) y un análisis bootstrap con 2000 remuestreos para obtener intervalos de confianza y gráficos de diagnóstico (bootstrap). Los gráficos ODS se desactivan a continuación.

¡Copiado!

1	ods graphics on;
2	PROC NLIN DATA=clarke1987a plots(stats=none)=diagnostics;
3	parms theta1=-0.15
4	theta2=2.0
5	theta3=0.80;
6	profile theta1 theta3 / range = -6 to 2 BY 0.2 all;
7	bootstrap / nsamples = 2000 seed=123 bootplots bootci bootcov;
8	model y = theta3 + theta2exp(theta1x);
9	RUN;
10	ods graphics off;

Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.

Información de copyright : S A S S A M P L E L I B R A R Y

Retour à la liste