Salut,
Je cherche à comprendre comment effectuer des tests d'hypothèses avec des données censurées et comment déterminer la taille de l'échantillon pour les études futures.
Mes données sont censurées à gauche et par intervalle. Je souhaite démontrer qu'un « événement » s'est produit pour 95 % des échantillons avant un certain point temporel (limite supérieure).
J'ai utilisé la fonction de distribution de durée de vie de JMP pour créer un graphique d'événements, déterminer le modèle le plus adapté à mes données et identifier le moment où 95 % de mes échantillons présentent un événement. J'ai ensuite comparé ce moment au quantile à 95 % à ma limite supérieure. Cette approche vous semble-t-elle correcte ?
Je souhaite ensuite utiliser les données du modèle issues de cette étude pour déterminer la taille d'échantillon nécessaire à ma prochaine étude. Auriez-vous des suggestions sur la manière de procéder ?
Cordialement
Thor
Vous êtes sur la bonne voie en utilisant l'estimation du quantile à 95 %. De plus, l'estimation du quantile à 95 % n'est qu'une estimation ponctuelle, elle comporte des incertitudes. La valeur réelle avec laquelle vous pouvez comparer n'est généralement pas cette valeur particulière, mais les limites de confiance qui lui sont associées. Mais cela dépend de votre objectif. Par exemple, si votre objectif est de trouver la quantité T indiquant qu'AU MOINS 95 % des événements se produisent avant T, avec un niveau de confiance souhaité L. Ensuite, vous devez rechercher la limite supérieure de la limite de confiance unilatérale au niveau de confiance L. Supposons que le niveau de confiance souhaité est de 90 %. Cette quantité T est la borne supérieure de la borne de confiance unilatérale au niveau de confiance de 90 %. Cette quantité T est supérieure au quantile à 95 %, elle tient donc compte des incertitudes d'estimation. Et si votre point de temps particulier est PLUS GRAND que cette quantité T, vous pouvez fermement affirmer que 95 % des événements se produisent avant ce point de temps. Il y a ici pas mal d'inégalités. Le tout est d'être du bon côté.
Pour trouver le nombre, vous devez savoir où spécifier le niveau de confiance souhaité ; voir cette capture d'écran :
Pour trouver la limite de confiance unilatérale, vous devez trouver "Estimation personnalisée":
Dans Custom Estimation, recherchez la partie Quantile, entrez 0,95 pour la probabilité, choisissez "Upper" pour Side. Et utilisez la "probabilité", au lieu de "Wald" pour votre rapport, au cas où la taille de votre échantillon n'est pas grande.
Pour la question sur la taille de l'échantillon, vous pouvez consulter deux endroits : Tests de durée de vie parmi les explorateurs de taille d'échantillon dans JMP 17 ou Plan de test de fiabilité parmi la taille et la puissance de l'échantillon dans JMP 17 et versions antérieures. Ce dernier est sur le point de prendre sa retraite. Ils sont tous deux sous le menu DOE. Le point principal est : identifiez votre objectif, qui devrait être l'un des deux, et dans votre cas, je suppose que l'objectif est d'estimer le quantile (l'autre est d'estimer la probabilité). À côté de cela, assurez-vous que l'objectif est de trouver la taille de l'échantillon. Spécifiez ensuite les informations restantes.
