Bonjour,
j'essaie de programmer un design d'essai randomisé avec N(groupe expé) = 2*N(groupe réf), et une analyse intermédiaire. Le critère de jugement est binaire, on attend p=0.3 dans le groupe de réf et p=0.55 dans le groupe expérimental. La puissance de l'essai est fixée à 80%, alpha bilatéral 5%. Pour vérifier que je comprends bien la proc seqdesign, je l'ai d'abord progammée pour un essai sans analyse intermédiaire, et j'ai comparé avec le NSN que j'obtiendrais avec la proc power. Et j'ai pas le même : 126 avec la proc seqdesign, 138 avec la proc power. Cette différence me semble trop importante pour ne pas signifier que je rate un truc quelque part ... qu'en pensez-vous ?

voici les codes :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
 
ods graphics on;
proc seqdesign
	/*altref=0.25*/
	plots=boundary(hscale=samplesize)
	;
TwoSidedOBrienFleming: design
	beta=0.2
	nstages=1
	method=obf
	stop=both
	;
samplesize
	model=twosamplefreq(nullprop = 0.3
						prop=0.55/*	 groupe A de la procédure : c'est le groupe expérimental
							pas besoin de le spécifier si altref est précisé */
						test = prop /* précise que nos hypothèses sont données sous forme de proportions et pas d'OR */
						ref = prop /*prop ? ou nullprop ? défaut = prop */
						weight = 2 /* poids du groupe expérimental */ )
						;
ods output Boundary=Bnd ;
run;
ods graphics off ;
 
proc power;
twosamplefreq test=pchi
groupproportions   = (0.3 0.55)
nullproportiondiff = 0
alpha              = 0.05
sides              = 2
power              = 0.8
gweights		   = (1 2)
ntotal             = .;
run;