¿Código SAS en JMP?Me gustaría saber si hay alguna manera de ejecutar un código SAS (abajo) en JMP. O cómo puedo emular el código SAS a continuación en JMP.
proc mixed
Class year week treatment rep;
model stdmisummax=year|week|treatment/ddfm=kr residual;
random rep(year) treatment*rep(year);
repeated week/type=ar(1) subject=treatment*rep(year);
run;
Mis factores son: Tratamiento, semana (medida repetida), año, rep (repetición).
Gracias de antemano.
PD: He intentado ejecutar algo similar en JMP usando el complemento de medidas repetidas, pero no funcionará, porque no puedo cruzar tres factores para hacer que la identificación del sujeto.
Creo que mi versión JMP no es pro o es una especie de versión pro para la universidad.
No tengo el software SAS.
Re: Codigo SAS en JMP?Sí, el efecto aleatorio, Year(Subject), o cualquiera que sea el primero, es la variación estimada entre sujetos .
Re: Codigo SAS en JMP?@jiancao Gracias por su explicación.
Cada vez que incluyo más de un efecto aleatorio, obtengo "0" para Estimación, Error estándar (ver captura de pantalla) para uno de los efectos aleatorios. ¿Por qué sería eso? ¿Es por la anidación?
Gracias por su ayuda,
isadora
Re: Codigo SAS en JMP?Puede agregar más de un efecto aleatorio en el análisis de medidas repetidas. Sin embargo, recibirá una nota que dice " La columna de tema repetido se ha anidado en todos los efectos aleatorios .” Vea la captura de pantalla a continuación. Creo que es por eso que Adam sugirió dejar Rep(Año) porque anidará en Tratamiento.
¿Has probado a incluir dos efectos aleatorios? ¿Qué han visto?
Re: Codigo SAS en JMP?Hola @AdamMorris y @DonMcCormack
Me gustaría saber si me podrías ayudar a entender un poco mejor por qué usaste solo un término como efecto aleatorio. ¿Es porque soy ¿Se me permite usar solo un término aleatorio en Efectos aleatorios cuando uso AR1?
Gracias,
isadora
Re: Codigo SAS en JMP?@DonMcCormack ¡Eso es muy útil!
Gracias por compartirlo conmigo y por darme su opinión sobre qué solución es la mejor. Realmente aprecio eso.
isadora
Re: Codigo SAS en JMP?Mi solución está a continuación:
Fit Model(
Y( :Response ),
Effects(
:Treatment,
:Week,
:Treatment * :Week,
:Year,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Rep[:Treatment, :Year] ),
NoBounds( 1 ),
Personality( "Mixed Model" ),
Subject( :Treatment, :Year, :Rep ),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Run( Variogram( AR1( 1 ) ) )
);
Solo hay dos diferencias. Para el efecto aleatorio adicional, anidé Rep en Year. Usé Tratamiento, Año y Rep para definir implícitamente el Sujeto. Sin embargo, encuentro que el enfoque de Adam de definir explícitamente el Sujeto es más atractivo intuitivamente.
Re: Codigo SAS en JMP?@AdamMorris , ¡Has sido de mucha ayuda!
¡Muchas gracias por compartir tus conocimientos y toda esta buena información! Soy consciente de que.
Re: Codigo SAS en JMP?@DonMcCormack Muchas gracias por su respuesta.
¿Te importaría compartir tu solución también? Me gustaría compararlos porque me interesa aprender más sobre la codificación en JMP y otros enfoques para mis datos.
¡Gracias!
isadora
Re: Codigo SAS en JMP?Hola Isadora,
Para compartir con otros que lean este hilo, la solución propuesta (del soporte técnico de JMP) es crear una nueva columna Asunto que es la concatenación de las 3 variables que identifican a los sujetos (Tratamiento, Año y Rep).
También es necesario crear una copia continua de la columna Semana para utilizarla como variable Repetida.
JMP Pro no puede ajustarse exactamente al mismo modelo que especifica el código SAS (para los modelos AR(1), JMP anida los efectos aleatorios en la variable sujeto).Sin embargo, JMP puede ajustarse al modelo omitiendo el término Rep(Year).La siguiente secuencia de comandos muestra cómo se puede configurar esto.
Gracias,
Adán
New Column("Subject", Character, Formula( :Treatment || ", " || Char( :Year ) || ", " || Char( :Rep ) ));
New Column("Week2", numeric, continuous, formula(:Week));
Fit Model(
Y( :response ),
Effects(
:Treatment,
:Week,
:Year,
:Treatment * :Week,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Treatment*:Year*:Rep ),
Subject(:Subject),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Personality( "Mixed Model" )
);Re: Codigo SAS en JMP?isadora,
Hablé con Adam Morris sobre su pregunta. Independientemente, obtuve los mismos resultados que él, aunque configuré el modelo de una manera ligeramente diferente. Me gustó un poco más su configuración, así que usaría esa. Por favor, hágamelo saber si tiene preguntas adicionales.
Don
Re: Codigo SAS en JMP?
Hola @DonMcCormack , he adjuntado una captura de pantalla de cómo se ve la tabla de datos para cada sujeto.
Tengo 4 repeticiones, 6 tratamientos y 2 años. Las muestras se recolectaron en 3 días diferentes después de que se aplicaron los tratamientos. Llamo a estos días de muestreo de No. de muestreo en la tabla de datos y se denominan "semana" en el script SAS.
Tengo JMP Pro 13.
Gracias de antemano,
isadora
Re: Codigo SAS en JMP?Sería útil saber un poco más sobre la estructura de la tabla de datos. ¿Cómo sería para un sujeto? ¿Se producen repeticiones dentro de la semana? Necesitará JMP Pro para usar la estructura de covarianza AR(1).
Re: Codigo SAS en JMP?No soy un experto en modelos de efectos mixtos en general ni en esta plataforma JMP Pro en particular. Tal vez alguien que sea tan experto se una a la discusión. De lo contrario, puede enviarlo al soporte técnico de JMP para obtener más sugerencias.
Re: Codigo SAS en JMP?@Mark_Bailey Gracias por compartir esta nota. Lo he leído antes.
Tengo la Plataforma Mixta en mi versión JMP. El problema es que no puedo estipular correctamente el ID del sujeto en esta plataforma, porque mi sujeto es una interacción de tratamiento*rep(año). Esa plataforma no me permite cruzar y anidar para agregar el asunto "tratamiento*rep(año)" en el cuadro ID del sujeto. Le permite agregar solo un factor.
Gracias
Re: Codigo SAS en JMP?Perdón por la confusion. No, no es posible ejecutar programas SAS en JMP.
JMP es capaz de realizar algunos análisis simples de medidas repetidas. Mira esto nota .Situaciones más complejas (es decir, estructuras de covarianza no simples) requieren la plataforma de modelo mixto en JMP Pro.
Re: Codigo SAS en JMP?@Mark_Bailey Gracias por su respuesta. El problema es que no quiero ejecutar los análisis en SAS. Quiero saber si es posible ejecutar este código en JMP de alguna manera o cómo puedo emular ese código en un script JMP o en Fit Model.
Gracias.
Re: Codigo SAS en JMP?En caso de que esto aplique a su situación, debe saber que JMP es un cliente completo de SAS. Entonces, así como SAS se ejecuta como un cliente, también lo puede hacer JMP. Requiere una conexión de servidor SAS. Ese servidor puede ejecutarse en la misma computadora, como es el caso de mi computadora portátil de trabajo, o puede ser remoto a través de una conexión de red. ¿Tiene acceso a un servidor SAS remoto (empresa, universidad)? Si es así, pegue su código SAS en un editor de secuencias de comandos JMP. Haga clic con el botón derecho en el editor y seleccione Enviar SAS.
Seleccione Archivo > Preferencias y configure primero la conexión al servidor SAS.
Re: Codigo SAS en JMP?@uday_guntupalli Gracias por su respuesta. Sí, lo sé, pero no quiero entrar en SAS. Obtuve este código de un colega, porque no pude o no pude averiguar cómo ejecutar los mismos análisis con los factores estipulados de esa manera en JMP.
Gracias 
Re: Codigo SAS en JMP?@bordini ,
Sé que dices que no tienes SAS, pero ¿sabías que puedes descargar la edición SAS University de forma gratuita? Si realmente necesita probar el código SAS, puede descargar la edición universitaria de SAS y probarlo.
https://www.sas.com/en_us/software/university-edition/download-software.html
No sé cómo realizar el paso equivalente en JMP, por lo que también buscaré otras respuestas
