Publicado el :
Estadística CREATION_INTERNE

Contrastes isotónicos con PROC MIXED

Este código también está disponible en: Deutsch English Français
En espera de validación
El script comienza con un DATA STEP que genera el conjunto de datos 'FerriteCores'. Este conjunto de datos simula mediciones de 'MagneticForce' para cuatro niveles de 'Temp' (temperatura), con cinco repeticiones por nivel, los datos se proporcionan a través de DATALINES. Luego, se invoca PROC MIXED para realizar un análisis de varianza sobre 'MagneticForce' en función de 'Temp'. La declaración CLASS especifica 'Temp' como variable categórica. La declaración LSMESTIMATE se utiliza para definir y probar contrastes isotónicos específicos, lo que permite comparar las medias de los niveles de 'Temp' de manera ordenada. Las opciones '/ adjust=simulate(seed=1) cl upper' se aplican para el ajuste de los valores p, el cálculo de los intervalos de confianza y la definición de un límite superior. Finalmente, se utiliza ODS SELECT LSMestimates para mostrar solo los resultados relevantes de las estimaciones de mínimos cuadrados.
Análisis de datos

Type : CREATION_INTERNE


Los datos se generan directamente en el script a través de un DATA STEP y DATALINES para el conjunto de datos 'FerriteCores'. El campo 'MagneticForce' se lee de los DATALINES.

1 Bloque de código
DATA STEP Data
Explicación :
Este bloque DATA STEP crea el conjunto de datos 'FerriteCores'. Inicia un bucle para generar una variable 'Temp' (temperatura) de 1 a 4, y un bucle interno para 5 repeticiones por 'Temp'. La variable 'MagneticForce' se lee de los DATALINES, simulando mediciones para diferentes condiciones de temperatura.
¡Copiado!
1DATA FerriteCores;
2 DO Temp = 1 to 4;
3 DO rep = 1 to 5; drop rep;
4 INPUT MagneticForce;
5 OUTPUT;
6 END;
7 END;
8 DATALINES;
910.8 9.9 10.7 10.4 9.7
1010.7 10.6 11.0 10.8 10.9
1111.9 11.2 11.0 11.1 11.3
1211.4 10.7 10.9 11.3 11.7
13;
14 
2 Bloque de código
PROC MIXED
Explicación :
Este bloque utiliza PROC MIXED para analizar el conjunto de datos 'FerriteCores'. 'Temp' se define como una variable de clasificación. Se especifica un modelo con 'MagneticForce' como respuesta y 'Temp' como predictor. La declaración LSMESTIMATE se utiliza para probar contrastes isotónicos, evaluando hipótesis ordenadas sobre las medias de los niveles de 'Temp'. Las opciones incluyen el ajuste de los valores p por simulación, el cálculo de los intervalos de confianza y la especificación de un límite superior. ODS SELECT LSMestimates limita la salida a los resultados de las estimaciones de mínimos cuadrados.
¡Copiado!
1PROC MIXED DATA=FerriteCores;
2 class Temp;
3 model MagneticForce = Temp;
4 lsmestimate Temp
5 'avg(1:1)<avg(2:4)' -3 1 1 1 divisor=3,
6 'avg(1:2)<avg(3:4)' -1 -1 1 1 divisor=2,
7 'avg(1:3)<avg(4:4)' -1 -1 -1 3 divisor=3
8 / adjust=simulate(seed=1) cl upper;
9 ods select LSMestimates;
10RUN;
11 
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