Cálculo del tamaño de la muestra para TOST¡Hola a todos!
Contexto: Tengo un proceso A que es mi referencia, y desarrollé un proceso B que quiero evaluar su comparabilidad con el proceso A. Para hacerlo, planeo usar un TOST (prueba de equivalencia) usando el 3xSTD (desviación estándar) del proceso A como umbral. Ahora, quiero saber cuántas veces debo repetir el proceso B para estar seguro de poder detectar dicha variación (sabiendo que una ejecución es una cantidad significativa de dinero).
Opción 1 : Quería usar la herramienta Tamaño de muestra y poder, con comparación de medias de 2 muestras, con un alfa de 0.1 (2x0.05 para cada proceso), un STD = 1, sin parámetro adicional, Diferencia para detectar = 3 (3x STD de 1), potencia = 0,95 y obtengo 7, por lo que 4 ejecuciones del proceso B para comparar con 4 ejecuciones del proceso A para garantizar la comparabilidad.
opcion 2 : Para usar la distribución t (como se muestra a continuación) y, por lo tanto, sugerir 6 ejecuciones:
¿Qué opinas de las dos opciones, cuál es mejor para este caso? ¿Tienes alguna otra herramienta que no se me haya ocurrido?
Muchas gracias ! 
@martindemel tal vez ?
Modelado Estadístico AvanzadoAnálisis y modelado de datos básicosExploración y visualización de datosLa personalización en masaModelado predictivo y aprendizaje automáticoIngeniería de Calidad y ProcesosRe: Cálculo del tamaño de la muestra para TOSTSí, está en JMP 16. Utilice el menú DOE --> Exploradores de tamaño de muestra --> Potencia para encontrar pruebas de equivalencia para casos de 1 y 2 muestras.
Re: Cálculo del tamaño de la muestra para TOSTEstimado equipo JMP,
¿Se ha agregado esta función en JMP16? ¿Si es así, donde? No he sido capaz de encontrarlo.
gracias de antemano
Ana
Re: Cálculo del tamaño de la muestra para TOSTHola, Juan,
Póngase en contacto con su representante de ventas de JMP SE/JMP para que le proporcione una invitación para el programa de adopción temprana (EA) para JMP 16 (siempre que sea un usuario y esté dispuesto a proporcionar comentarios). En el EA, nuestro desarrollo comparte nuevas funciones potenciales y solicita comentarios antes de que se envíe en la nueva versión. Esta es una gran oportunidad para proporcionar comentarios sobre la "nueva" forma de hacer cálculos de potencia y pruebas de equivalencia, así como si todavía hay algo diferente de lo que esperaría. Y para mediados o finales de mayo lo más probable es que se presente la nueva plataforma de potencia y tamaño de muestra en el EA, incluidas las pruebas de equivalencia hasta donde yo sé. Su retroalimentación es muy apreciada.
Mejor,
Martín
Re: Cálculo del tamaño de la muestra para TOST@calking - ¿Hay alguna posibilidad de que se agreguen pruebas de equivalencia?
Actualmente se puede hacer en un software similar, como Minitab. https://support.minitab.com/en-us/minitab/18/help-and-how-to/statistics/equivalence-tests/supporting-topics/power-for-equivalence-tests/
*se agacha*
Re: Cálculo del tamaño de la muestra para TOST¡Ciertamente!
Para responder a sus preguntas de seguimiento, los valores que elegí para el alfa y la potencia para usar en la Calculadora de potencia en realidad se basaron en los valores de doble alfa y potencia que estaba buscando. Para aclarar, dado que JMP actualmente no tiene una calculadora para una prueba de equivalencia, estamos usando un truco (ver artículo) para aproximar el cálculo de potencia usando la calculadora de potencia de prueba t de dos muestras. El truco consiste en el hecho de que la prueba de equivalencia ha cambiado las hipótesis que normalmente usa en una prueba t de dos muestras, por lo que necesitaremos cambiar algunos valores. Así que para tu caso, procedimos de la siguiente manera:
1. Quieres un alfa de 0.1 (2x0.05). Para el truco, esto significa que establecerá la potencia en la calculadora en 0,9; es decir, 1- (alfa destinada a la prueba de equivalencia).
2. Quieres una potencia de 0,95. Para el truco, esto significa que establecerá el nivel alfa en la calculadora en 0,05; es decir, 1-(potencia destinada a prueba de equivalencia).
Tenga en cuenta que tuvimos que cambiar los valores que usamos para potencia y alfa. Espero que eso ayude a aclarar de dónde vienen esos valores que usé.
Me aseguré de usar el valor alfa total en lugar de los valores individuales para coincidir específicamente con sus valores particulares. Sin embargo, si 0,05 es de hecho el nivel de confianza global que desea lograr, entonces ese es probablemente el valor que debe usar. Un valor de 0,1 dará como resultado un tamaño de muestra total más pequeño, pero solo porque está aflojando la restricción sobre la confianza global (90 % de confianza frente a 95 % de confianza). Además, la prueba t de dos muestras ya tiene en cuenta que hay dos poblaciones y, por lo tanto, ya se ajusta para garantizar que se mantenga el nivel de confianza global.
Ejecuté algunos números rápidos y, de cualquier manera (alfa = 0,1 frente a alfa = 0,05), la diferencia general en el tamaño total de la muestra es solo 8 frente a 9. Así que, en última instancia, todo se reducirá a cómo su gerencia decida sopesar los riesgos y los rendimientos (8 muestras totales para un 90 % de confianza global frente a 9 muestras totales para un 95 % de confianza global).
Re: Cálculo del tamaño de la muestra para TOSTBien genial muchas gracias!
¡No tengo R, así que intentaré esto con el script!Re: Cálculo del tamaño de la muestra para TOSTHola @calking ,
¡Eso suena bien, gracias! Par de preguntas:
1. Sobre el alfa, nos dijeron que como tenemos 2 poblaciones, tenemos que aplicar el alfa de 1-0.95 a ambos, entonces para usar 0.05x2=0.1 en el alfa global de la herramienta, ¿qué piensas al respecto?
2. Sobre la potencia, ¿es normal usar 0.9?
¡Muchas gracias por tu ayuda!
Re: Cálculo del tamaño de la muestra para TOSTPara el caso de una muestra, especifique los valores de estas variables:
n = 15; //total sample size
conf = 0.95; //confidence
mu = 505; //hypothesized mean
muu = 510; //upper equivalence limit
mul = 490; //lower equivalence limit
sigma = 4; //estimated standard deviation
Para el caso de dos muestras, especifique los valores de estas variables:
N = 36; //total sample size
W = 1; //weighting for sample size (e.g. use 2 for 1 of the groups to have 2x as many samples as the other
conf = 0.95; //confidence level
mudiff_l = -2; //lower equivalence limit on mean difference
mudiff_u = 2; //upper equivalence limit on mean difference
mudiff = 0; //estimated mean difference
sigma = 2; //estimated standard deviation
Para encontrar un tamaño de muestra, solo necesitaría implementar un bucle para aumentar N hasta lograr la potencia deseada.
Entiendo su preocupación sobre la necesidad de enviar esto a una autoridad para su revisión.Yo mismo trabajo en una empresa regulada y tenemos un protocolo de validación para el complemento JMP con las calculadoras que construimos en base a este script. La documentación de SAS para PROC POWER muestra las fórmulas para el cálculo de potencia para 1 y 2 muestras TOST.Puede hacer referencia a ese documento para mostrar que los cálculos son precisos, reproducirlo en otro software (por ejemplo, R, tengo un script para eso si lo desea) o confirmar las estimaciones de potencia a través de la simulación.Si tiene SAS en cualquier parte de su organización, puede obtener una copia impresa de la energía del proceso y listo.
Re: Cálculo del tamaño de la muestra para TOSTGracias por esa respuesta, no estoy seguro de haber entendido todo el guión, ¿qué se supone que debo cambiar?
Mi preocupación es que este experimento es "oficial" (se enviará a las autoridades), por lo tanto, no estoy seguro de usar dicho script...Re: Cálculo del tamaño de la muestra para TOSTDesarrollamos una aplicación que realiza cálculos de tamaño de muestra TOST.No puedo enviar eso, pero puedo compartir el script que usé para replicar los cálculos de potencia para TOST en SAS usando Owen's Q.Tiene el código para escenarios de 1 y 2 lados con un ejemplo para cada uno.Tenga en cuenta que esto no manejará bien los escenarios pequeños de señal a ruido.En mi experiencia con esto hasta ahora, solo funciona hasta que nu (parámetro de tamaño de muestra) llega a estar alrededor de 200-250 porque Gamma( nu/2 ) se encuentra con un problema de desbordamiento numérico y devuelve un valor faltante.
//Owen's Q
OwensQ = Function( {t, delta, a, b, nu},
c = Sqrt( 2 * Pi() ) / (Gamma( nu / 2 ) * Power( 2, (nu - 2) / 2 ));
integrand = Expr(
Normal Distribution( ((t * x / Sqrt( nu )) - delta) ) * Power( x, nu - 1 ) * Normal Density( x )
);
c * Integrate( integrand, x, a, b, <<STARTING value=""></STARTING>
Re: Cálculo del tamaño de la muestra para TOSTHola Elofar,
encontré un artículo lo que parecía indicar que podría aproximarse al cálculo de la potencia mediante el uso de una prueba t donde alpha_ttest = 1- power_equiv, power_ttest = 1-alpha_equiv, y la diferencia para detectar coincide con el valor umbral para la prueba de equivalencia. Usando la calculadora Potencia y tamaño de muestra para dos medias de muestra, especificaría alfa=0,05 (1-0,95), StdDev=1, Diferencia para detectar = 3 y Potencia=0,9 (1-0,1). Eso le da un tamaño de muestra total de 8, que parece estar de acuerdo con su Opción 1.
¡Espero que ayude!
