¿Cómo calcular las estimaciones REML de los componentes de varianza según la validación del ensayo biológico USP <1033>?Hola,
He encontrado un problema al calcular los resultados del análisis REML en base a los datos y resultados presentados en USP <1033>. Tengo problemas para calcular los datos de la Tabla 4 y obtener los mismos resultados que los de la Tabla 9. ¿Hay alguna persona que pueda manejar estos datos y explicar cómo proceder para realizar el análisis REML utilizando este ejemplo?
¡Gracias!
Modelado estadístico avanzadoDiseño de experimentosIngeniería de Calidad y ProcesosRe: ¿Cómo calcular las estimaciones REML de los componentes de varianza según la validación del ensayo biológico USP <1033>?Para mayor información proporciono las fuentes: 1033 VALIDACIÓN DE ENSAYO BIOLÓGICO
Re: ¿Cómo calcular las estimaciones REML de los componentes de varianza según la validación del ensayo biológico USP <1033>?Hola @Luk_s , bienvenido. Puedo reproducir las tablas 5 y 6, pero no la 9. Parece que los resultados de la tabla 9 son un ejemplo que utiliza datos diferentes. Esto se puede ver por la diferencia en Var(error). En la tabla 5 es 0,000766. En la tabla 9 es 0,0022. Deberían ser iguales si ambas utilizan los mismos datos. Por lo tanto, la tabla 9 no se basa en los datos que se muestran en la tabla 4. Es un ejemplo de otra cosa, como lo implica la declaración "Las estimaciones de los componentes de la varianza obtenidas de dicho análisis se presentan en la tabla 9" (la cursiva es mía); observe que no se refiere a los datos de la tabla 4. Y, como señalan correctamente en la tabla 5, "Los analistas y los medios de comunicación no se han incluido debido a la pequeña cantidad de niveles (2 niveles) para cada factor". , por lo que no sería aconsejable realizar el tipo de análisis que se muestra en la tabla 9.
Para que conste, y por más desaconsejable que parezca, obtengo lo siguiente (Tabla 4, Nivel 0.5):
Re: ¿Cómo calcular las estimaciones REML de los componentes de varianza según la validación del ensayo biológico USP <1033>?Hola @MRB3855, muchas gracias por tu atención al problema que te he presentado. Estaba convencido de que estaba haciendo algo mal al desarrollar los datos para el modelo. También hice una comprobación y obtuvimos los mismos resultados, lo que me complace mucho. También traté con las tablas de datos anteriores, sin embargo, me sorprendieron los resultados obtenidos en este ejemplo.
¡Gracias!
