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

With proc glm for Comparison

The MEANS Procedure

The Means Procedure

Summary statistics

Analysis Variable : infection Ave percent surface covered
condition N Obs N Mean Std Dev Minimum Maximum
Control 8 8 22.6250000 8.3637722 15.5600000 35.2600000
Spring300 4 4 16.7475000 10.7824900 5.1300000 30.1200000
Spring600 4 4 18.2500000 4.9273996 11.4000000 22.7500000
Spring1200 4 4 14.2525000 4.8578896 8.1100000 19.8400000
Fall300 4 4 9.4975000 4.9333854 5.0800000 14.8000000
Fall600 4 4 15.5025000 3.7887586 11.1000000 20.2300000
Fall1200 4 4 5.7500000 2.8700639 2.0100000 8.4800000

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 DF Pr > ChiSq
15.8552 6 0.0146

Monte Carlo Estimates for the Exact Test

Monte Carlo Estimates for the Exact Test
Probability Estimate 99% Confidence Limits Samples Seed
Pr >= ChiSq 0.0043 0.0026 0.0060 10000 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 DF Pr > ChiSq
14.3869 6 0.0256

Monte Carlo Estimates for the Exact Test

Monte Carlo Estimates for the Exact Test
Probability Estimate 99% Confidence Limits Samples Seed
Pr >= ChiSq 0.0099 0.0073 0.0125 10000 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 1.2328828

Randomization Tests on the Scab Disease Data

Bigger MC sample size: See Control vs. Fall300

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 200000
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.7871
infection Control vs. Spring600 0.2966 0.9355
infection Control vs. Spring120 0.0520 0.4173
infection Control vs. Fall300 0.0037 0.0515
infection Control vs. Fall600 0.0950 0.6055
infection Control vs. Fall1200 0.0004 0.0074
infection Spring300 vs. Spring6 0.7538 0.9999
infection Spring300 vs. Spring1 0.6032 0.9983
infection Spring300 vs. Fall300 0.1386 0.7319
infection Spring300 vs. Fall600 0.7949 1.0000
infection Spring300 vs. Fall120 0.0288 0.2741
infection Spring600 vs. Spring1 0.4069 0.9791
infection Spring600 vs. Fall300 0.0766 0.5350
infection Spring600 vs. Fall600 0.5673 0.9973
infection Spring600 vs. Fall120 0.0142 0.1576
infection Spring1200 vs. Fall30 0.3253 0.9512
infection Spring1200 vs. Fall60 0.7941 1.0000
infection Spring1200 vs. Fall12 0.0849 0.5679
infection Fall300 vs. Fall600 0.2168 0.8658
infection Fall300 vs. Fall1200 0.4365 0.9852
infection Fall600 vs. Fall1200 0.0502 0.4073
infection Control vs. Spring 0.0532 0.4231
infection Control vs. Fall 0.0004 0.0087
infection Spring vs. Fall 0.0332 0.3047

Randomization Tests on the Scab Disease Data

Another large-sample Z-test of H0: p = 0.05

The IML Procedure

p_Z

p Z
0.0515 3.0779351