Analysis of the Fibra data, using Fortaleza as the response variable and Algodon as the single grouping factor.

Prepared by BrailleR

Group summaries

The following objects are masked _by_ .GlobalEnv:

    algodon, fortaleza
Level Mean Standard deviation n Standard error
15% 9.8 3.346640 5 1.4966630
20% 15.4 3.130495 5 1.4000000
25% 17.6 2.073644 5 0.9273618
30% 21.6 2.607681 5 1.1661904
35% 10.8 2.863564 5 1.2806248

The ratio of the largest group standard deviation to the smallest is 1.61

Comparative boxplots

Comparative boxplots

 This graph has 5 boxplots printed vertically
 with the title: 
 "" appears on the x-axis.
 "" appears on the y-axis.
 Tick marks for the y-axis are at: 10, 15, 20, and 25 
 Group 15% has 5 values.
 There are no outliers marked for this group 
 The whiskers extend to 7 and 15 from the ends of the box, 
 which are at 7 and 11 
 The median, 9 is 50 % from the lower end of the box to the upper end.
 The upper whisker is Inf times the length of the lower whisker.
 Group 20% has 5 values.
 There are no outliers marked for this group 
 The whiskers extend to 12 and 18 from the ends of the box, 
 which are at 12 and 18 
 The median, 17 is 83 % from the lower end of the box to the upper end.
 The upper whisker is NaN times the length of the lower whisker.
 Group 25% has 5 values.
 An outlier is marked at: 14 
 The whiskers extend to 18 and 19 from the ends of the box, 
 which are at 18 and 19 
 The median, 18 is 0 % from the lower end of the box to the upper end.
 The upper whisker is NaN times the length of the lower whisker.
 Group 30% has 5 values.
 There are no outliers marked for this group 
 The whiskers extend to 19 and 25 from the ends of the box, 
 which are at 19 and 23 
 The median, 22 is 75 % from the lower end of the box to the upper end.
 The upper whisker is Inf times the length of the lower whisker.
 Group 35% has 5 values.
 An outlier is marked at: 7 15 
 The whiskers extend to 10 and 11 from the ends of the box, 
 which are at 10 and 11 
 The median, 11 is 100 % from the lower end of the box to the upper end.
 The upper whisker is NaN times the length of the lower whisker.

Comparative dotplots

Warning in dotplot.formula(fortaleza ~ algodon, xlab = "Fortaleza", ylab =
"Algodon"): The dotplot command is a wrapper for stripchart, and is not a base R
function.

Comparative dotplots

This graph has 25 dotplots printed horizontally
with the title: 
"" appears on the x-axis.
"" appears on the y-axis.
Tick marks for the x-axis are at: 10, 15, 20, and 25 
the data that range from 7 to 25 has been broken into 6 bins.
The counts are:
15%: 1 0 0 0 0 0 
20%: 0 1 0 0 0 0 
25%: 0 0 1 0 0 0 
30%: 0 0 0 1 0 0 
35%: 1 0 0 0 0 0 
NA: 1 0 0 0 0 0 
NA: 0 1 0 0 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 1 0 0 
NA: 1 0 0 0 0 0 
NA: 1 0 0 0 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 0 1 0 
NA: 0 1 0 0 0 0 
NA: 0 1 0 0 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 0 0 1 
NA: 0 1 0 0 0 0 
NA: 0 0 1 0 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 1 0 0 
NA: 0 0 0 0 0 1 
NA: 0 0 1 0 0 0 

One-way Analysis of Variance


The p value for algodon is 0

The ratios of the group standard deviations to the overall standard deviation 
 (groups ordered by increasing mean) are:
 1.18 1.1 0.73 0.92 1.01

            Df Sum Sq Mean Sq F value   Pr(>F)    
algodon      4  475.8  118.94   14.76 9.13e-06 ***
Residuals   20  161.2    8.06                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual analysis

Tests for homogeneity of Variance


    Bartlett test of homogeneity of variances

data:  fortaleza by algodon
Bartlett's K-squared = 0.93309, df = 4, p-value = 0.9198

    Fligner-Killeen test of homogeneity of variances

data:  fortaleza by algodon
Fligner-Killeen:med chi-squared = 1.7265, df = 4, p-value = 0.7859

Tukey Honestly Significant Difference test

For term algodon the comparisons which are significant at 5% are:
20% and 15% with a difference of 5.6 and P-value of 0.0385 
25% and 15% with a difference of 7.8 and P-value of 0.0026 
30% and 15% with a difference of 11.8 and P-value of 0 
25% and 35% with a difference of 6.8 and P-value of 0.0091 
30% and 35% with a difference of 10.8 and P-value of 1e-04 
30% and 20% with a difference of 6.2 and P-value of 0.0189 
  Tukey multiple comparisons of means
    95% family-wise confidence level
    factor levels have been ordered

Fit: aov(formula = fortaleza ~ algodon, data = fibra)

$algodon
        diff        lwr       upr     p adj
35%-15%  1.0 -4.3729583  6.372958 0.9797709
20%-15%  5.6  0.2270417 10.972958 0.0385024
25%-15%  7.8  2.4270417 13.172958 0.0025948
30%-15% 11.8  6.4270417 17.172958 0.0000190
20%-35%  4.6 -0.7729583  9.972958 0.1162970
25%-35%  6.8  1.4270417 12.172958 0.0090646
30%-35% 10.8  5.4270417 16.172958 0.0000624
25%-20%  2.2 -3.1729583  7.572958 0.7372438
30%-20%  6.2  0.8270417 11.572958 0.0188936
30%-25%  4.0 -1.3729583  9.372958 0.2101089

Plot of Tukey HSD