Description Details Author(s) References Examples
For the implementation of the step-down or step-up Dunnett test procedures, the package includes R functions to calculate critical constants and R functions to calculate adjusted P-values of test statistics. In addition, the package also contains functions to evaluate testing powers and hence the necessary sample sizes for the classic statistical problem of comparing multiple treatments with a control.
Package: | DunnettTests |
Type: | Package |
Version: | 1.0 |
Date: | 2013-11-11 |
License: | GPL-2 |
FAN XIA Maintainer: FAN XIA <phoebexia@yahoo.com>
Charles W Dunnett and Ajit C Tamhane. Step-down multiple tests for comparing treatments with a control in unbalanced one-way layouts. Statistics in Medicine, 10(6):939-947, 1991.
Charles W. Dunnett and Ajit C. Tamhane. A step-up multiple test procedure. Journal of the American Statistical Association, 87(417):162-170, 1992.
Charles W Dunnett and Ajit C Tamhane. Step-up multiple testing of parameters with unequally correlated estimates. Biometrics, pages 217-227, 1995.
Ajit C Tamhane and Charles W Dunnett. Stepwise multiple test procedures with biometric applications. Journal of statistical planning and inference, 82(1):55-68, 1999.
Koon Shing Kwong and Wei Liu. Calculation of critical values for dunnett and tamhane???s step-up multiple test procedure. Statistics & probability letters, 49(4):411-416, 2000.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | #critical constants
cvSDDT(k=4,alpha=0.05,alternative="U",corr=0.5,df=30)
cvSUDT(k=4,alpha=0.05,alternative="U",corr=0.5,df=30)
#adjusted P-values
qvSDDT(teststats=c(2.20,1.47,2.00),alternative="B",corr=0.5,df=30)
## Not run:
qvSUDT(teststats=c(2.20,1.47,2.00),alternative="B",corr=0.5,df=30)
## End(Not run)
#testing power of SU Dunnett to reject at least one false null
powDT(1,4,mu=22,mu0=20,n=100,n0=80,"means",sigma=5,alpha=0.025,testcall="SU")
#sample size to achieve at least 90% power of rejecting all the
#false nulls for step-down Dunnett test procedure
## Not run:
nvDT(1,0.90,r=4,k=4,mu=0.7,mu0=0.5,"props",dist="zdist",testcall="SD")
## End(Not run)
|
Loading required package: mvtnorm
c1 c2 c3 c4
1.697 1.989 2.148 2.254
c1 c2 c3 c4
1.697 2.008 2.157 2.260
$`ordered test statistics`
H1 H3 H2
2.20 2.00 1.47
$`Adjusted P-values of ordered test statistics`
[1] 0.090 0.098 0.152
$`ordered test statistics`
H2 H3 H1
1.47 2.00 2.20
$`Adjusted P-values of ordered test statistics`
[1] 0.152 0.102 0.092
[1] 0.8574494
$`least sample size required in each treatment groups`
[1] 140
$`least sample size required in the control group`
[1] 140
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.