Suite aux récentes questions que nous avions ici .
Je sautais pour savoir si quelqu'un était tombé sur ou pouvait partager du code R pour effectuer une analyse de puissance personnalisée basée sur la simulation d'un modèle linéaire?
Plus tard, je voudrais évidemment l'étendre à des modèles plus complexes, mais lm semble être le bon endroit pour commencer. Merci.
r
power
power-analysis
Tal Galili
la source
la source
Voici quelques sources de code de simulation dans R. Je ne sais pas s'il en existe spécifiquement sur les modèles linéaires, mais ils fournissent peut-être suffisamment d'exemples pour obtenir l'essentiel:
Il existe deux autres exemples de simulation sur les sites suivants:
la source
Adapté de Bolker 2009 Ecological Models and Data in R. Vous devez déclarer la force de la tendance (c'est-à-dire la pente) que vous souhaitez tester. Intuitivement, une tendance forte et une faible variabilité nécessiteront une petite taille d'échantillon, une tendance faible et une grande variabilité nécessiteront une grande taille d'échantillon.
Vous pouvez également simuler la tendance minimale que vous pourriez tester pour une taille d'échantillon donnée, comme indiqué dans le livre
la source