Results: scab2.sas

Randomization Tests on the Scab Disease Data

With proc glm for Comparison

The GLM Procedure

 

Dependent Variable: infection Ave percent surface covered

The GLM Procedure

Analysis of Variance

infection

Overall ANOVA

Source DF Sum of Squares Mean Square F Value Pr > F
Model 6 972.413850 162.068975 3.61 0.0103
Error 25 1122.880900 44.915236    
Corrected Total 31 2095.294750      

Randomization Tests on the Scab Disease Data

Classical Kruskal-Wallis Test

The NPAR1WAY Procedure

The Npar1way Procedure

Variable infection

Wilcoxon Analysis

Scores

Wilcoxon Scores (Rank Sums) for Variable infection
Classified by Variable condition
condition N Sum of
Scores
Expected
Under H0
Std Dev
Under H0
Mean
Score
Average scores were used for ties.
Control 8 191.00 132.0 22.976145 23.8750
Spring300 4 69.00 66.0 17.548320 17.2500
Spring600 4 87.00 66.0 17.548320 21.7500
Spring1200 4 60.00 66.0 17.548320 15.0000
Fall300 4 34.50 66.0 17.548320 8.6250
Fall600 4 69.00 66.0 17.548320 17.2500
Fall1200 4 17.50 66.0 17.548320 4.3750

Kruskal-Wallis Test

Kruskal-Wallis Test
Chi-Square 15.8552
DF 6
Pr > Chi-Square 0.0146

Monte Carlo Estimate for the Exact Test

Monte Carlo Estimate for the Exact Test
Pr >= Chi-Square  
Estimate 0.0043
99% Lower Conf Limit 0.0026
99% Upper Conf Limit 0.0060
   
Number of Samples 10000
Initial Seed 77777

Box Plot

Box Plot of Wilcoxon Scores for infection Classified by condition

Randomization Tests on the Scab Disease Data

Randomization Test

The NPAR1WAY Procedure

The Npar1way Procedure

Variable infection

Data Scores Analysis

Scores

Data Scores for Variable infection
Classified by Variable condition
condition N Sum of
Scores
Expected
Under H0
Std Dev
Under H0
Mean
Score
Control 8 181.000 125.2500 20.138047 22.62500
Spring300 4 66.990 62.6250 15.380687 16.74750
Spring600 4 73.000 62.6250 15.380687 18.25000
Spring1200 4 57.010 62.6250 15.380687 14.25250
Fall300 4 37.990 62.6250 15.380687 9.49750
Fall600 4 62.010 62.6250 15.380687 15.50250
Fall1200 4 23.000 62.6250 15.380687 5.75000

One-Way Analysis

Data Scores One-Way Analysis
Chi-Square 14.3869
DF 6
Pr > Chi-Square 0.0256

Monte Carlo Estimate for the Exact Test

Monte Carlo Estimate for the Exact Test
Pr >= Chi-Square  
Estimate 0.0099
99% Lower Conf Limit 0.0073
99% Upper Conf Limit 0.0125
   
Number of Samples 10000
Initial Seed 88888

Box Plot

Box Plot of Data Scores for infection Classified by condition

Randomization Tests on the Scab Disease Data

Randomization multiple comparisons with proc multtest

The Multtest Procedure

The Multtest Procedure

Model Information

Model Information
Test for continuous variables Mean t-test
Degrees of Freedom Method Pooled
Tails for continuous tests Two-tailed
Strata weights None
P-value adjustment Permutation
Center continuous variables No
Number of resamples 20000
Seed 99999

Contrast Coefficients

Contrast Coefficients
Contrast   condition
Control Spring300 Spring600 Spring1200 Fall300 Fall600 Fall1200
Control vs. Spring300 Centered 1 -1 0 0 0 0 0
Control vs. Spring600 Centered 1 0 -1 0 0 0 0
Control vs. Spring120 Centered 1 0 0 -1 0 0 0
Control vs. Fall300 Centered 1 0 0 0 -1 0 0
Control vs. Fall600 Centered 1 0 0 0 0 -1 0
Control vs. Fall1200 Centered 1 0 0 0 0 0 -1
Spring300 vs. Spring6 Centered 0 1 -1 0 0 0 0
Spring300 vs. Spring1 Centered 0 1 0 -1 0 0 0
Spring300 vs. Fall300 Centered 0 1 0 0 -1 0 0
Spring300 vs. Fall600 Centered 0 1 0 0 0 -1 0
Spring300 vs. Fall120 Centered 0 1 0 0 0 0 -1
Spring600 vs. Spring1 Centered 0 0 1 -1 0 0 0
Spring600 vs. Fall300 Centered 0 0 1 0 -1 0 0
Spring600 vs. Fall600 Centered 0 0 1 0 0 -1 0
Spring600 vs. Fall120 Centered 0 0 1 0 0 0 -1
Spring1200 vs. Fall30 Centered 0 0 0 1 -1 0 0
Spring1200 vs. Fall60 Centered 0 0 0 1 0 -1 0
Spring1200 vs. Fall12 Centered 0 0 0 1 0 0 -1
Fall300 vs. Fall600 Centered 0 0 0 0 1 -1 0
Fall300 vs. Fall1200 Centered 0 0 0 0 1 0 -1
Fall600 vs. Fall1200 Centered 0 0 0 0 0 1 -1
Control vs. Spring Centered 3 -1 -1 -1 0 0 0
Control vs. Fall Centered 3 0 0 0 -1 -1 -1
Spring vs. Fall Centered 0 1 1 1 -1 -1 -1

Continuous Variable Tabulations

Continuous Variable Tabulations
Variable condition NumObs Mean Standard Deviation
infection Control 8 22.6250 8.3638
infection Spring300 4 16.7475 10.7825
infection Spring600 4 18.2500 4.9274
infection Spring1200 4 14.2525 4.8579
infection Fall300 4 9.4975 4.9334
infection Fall600 4 15.5025 3.7888
infection Fall1200 4 5.7500 2.8701

P-Values

p-Values
Variable Contrast Raw Permutation
infection Control vs. Spring300 0.1645 0.7877
infection Control vs. Spring600 0.2966 0.9343
infection Control vs. Spring120 0.0520 0.4152
infection Control vs. Fall300 0.0037 0.0519
infection Control vs. Fall600 0.0950 0.6022
infection Control vs. Fall1200 0.0004 0.0068
infection Spring300 vs. Spring6 0.7538 0.9999
infection Spring300 vs. Spring1 0.6032 0.9986
infection Spring300 vs. Fall300 0.1386 0.7307
infection Spring300 vs. Fall600 0.7949 1.0000
infection Spring300 vs. Fall120 0.0288 0.2725
infection Spring600 vs. Spring1 0.4069 0.9791
infection Spring600 vs. Fall300 0.0766 0.5308
infection Spring600 vs. Fall600 0.5673 0.9973
infection Spring600 vs. Fall120 0.0142 0.1578
infection Spring1200 vs. Fall30 0.3253 0.9515
infection Spring1200 vs. Fall60 0.7941 1.0000
infection Spring1200 vs. Fall12 0.0849 0.5645
infection Fall300 vs. Fall600 0.2168 0.8656
infection Fall300 vs. Fall1200 0.4365 0.9855
infection Fall600 vs. Fall1200 0.0502 0.4058
infection Control vs. Spring 0.0532 0.4211
infection Control vs. Fall 0.0004 0.0078
infection Spring vs. Fall 0.0332 0.3031

Randomization Tests on the Scab Disease Data

Testing H0: p = 0.05 with a large-sample Z-test for proportions

The IML Procedure

p_Z

p Z
0.0519 5.6568542