Analysis: Dunnett test
dunnett.RdThe function performs the Dunnett test
dunnett(
trat,
resp,
control,
model = "DIC",
block = NA,
column = NA,
line = NA,
alpha.t = 0.05,
pointsize = 5,
pointshape = 21,
linesize = 1,
labelsize = 4,
textsize = 12,
errorsize = 1,
widthsize = 0.2,
label = "Response",
fontfamily = "sans"
)Arguments
- trat
Numerical or complex vector with treatments
- resp
Numerical vector containing the response of the experiment.
- control
Treatment considered control (write identical to the name in the vector)
- model
Experimental design (DIC, DBC or DQL)
- block
Numerical or complex vector with blocks
- column
Numerical or complex vector with columns
- line
Numerical or complex vector with lines
- alpha.t
Significance level (default is 0.05)
- pointsize
Point size
- pointshape
Shape
- linesize
Line size
- labelsize
Label size
- textsize
Font size
- errorsize
Errorbar size
- widthsize
Width errorbar
- label
Variable label
- fontfamily
font family
Value
I return the Dunnett test for experiments in a completely randomized design, randomized blocks or Latin square.
Note
Do not use the "-" symbol or space in treatment names
Examples
#====================================================
# complete randomized design
#====================================================
data("pomegranate")
with(pomegranate,dunnett(trat=trat,resp=WL,control="T1"))
#> Estimate IC-lwr IC-upr t value p-value sig
#> T1 - T2 -0.1650 -0.6433104 0.3133104 -0.9528 0.8035 ns
#> T1 - T3 0.7750 0.2966896 1.2533104 4.4753 0.0013 *
#> T1 - T4 0.7775 0.2991896 1.2558104 4.4898 0.0013 *
#> T1 - T5 0.7950 0.3166896 1.2733104 4.5908 0.0011 *
#> T1 - T6 0.3200 -0.1583104 0.7983104 1.8479 0.2667 ns
#====================================================
# randomized block design in factorial double
#====================================================
library(AgroR)
data(cloro)
attach(cloro)
#> The following object is masked from simulate3:
#>
#> resp
#> The following object is masked from simulate1:
#>
#> resp
#> The following object is masked from aristolochia (pos = 5):
#>
#> resp
#> The following objects are masked from simulate2:
#>
#> bloco, resp
#> The following objects are masked from laranja:
#>
#> bloco, resp
#> The following object is masked from aristolochia (pos = 8):
#>
#> resp
#> The following objects are masked from cloro (pos = 9):
#>
#> bloco, f1, f2, resp
#> The following object is masked from passiflora:
#>
#> bloco
respAd=c(268, 322, 275, 350, 320)
a=FAT2DBC.ad(f1, f2, bloco, resp, respAd,
ylab="Number of nodules",
legend = "Stages",mcomp="sk")
#>
#> -----------------------------------------------------------------
#> Normality of errors
#> -----------------------------------------------------------------
#> Method Statistic p.value
#> Shapiro-Wilk normality test(W) 0.9548911 0.1117923
#>
#> As the calculated p-value is greater than the 5% significance level, hypothesis H0 is not rejected. Therefore, errors can be considered normal
#>
#> -----------------------------------------------------------------
#> Homogeneity of Variances
#> -----------------------------------------------------------------
#> Method Statistic p.value
#> Bartlett test(Bartlett's K-squared) 16.11086 0.02412261
#>
#> As the calculated p-value is less than the 5% significance level, H0 is rejected. Therefore, the variances are not homogeneous
#>
#> -----------------------------------------------------------------
#> Independence from errors
#> -----------------------------------------------------------------
#> Method Statistic p.value
#> Durbin-Watson test(DW) 2.047899 0.1769663
#>
#> As the calculated p-value is greater than the 5% significance level, hypothesis H0 is not rejected. Therefore, errors can be considered independent
#>
#> -----------------------------------------------------------------
#> Additional Information
#> -----------------------------------------------------------------
#>
#> CV (%) = 28.29
#> Mean Factorial = 218.35
#> Median Factorial = 185
#> Mean Aditional = 307
#> Median Aditional = 320
#> Possible outliers = No discrepant point
#>
#> -----------------------------------------------------------------
#> Analysis of Variance
#> -----------------------------------------------------------------
#> Df Sum Sq Mean.Sq F value Pr(F)
#> Fator1 1 16160.400 16160.400 3.8772062 5.765683e-02
#> Fator2 3 116554.500 38851.500 9.3212591 1.409219e-04
#> block 4 7122.311 1780.578 0.4271966 7.878578e-01
#> Fator1:Fator2 3 452096.200 150698.733 36.1556682 2.154175e-10
#> Ad x Factorial 1 34928.100 34928.100 8.3799563 6.781309e-03
#> Residuals 32 133377.689 4168.053
#>
#>
#> Your analysis is not valid, suggests using a non-parametric test and try to transform the data
#>
#> -----------------------------------------------------------------
#> Significant interaction: analyzing the interaction
#> -----------------------------------------------------------------
#> Df Sum Sq Mean Sq F value Pr(>F)
#> Fator2 3 116554 38851 9.3213 0.0001409 ***
#> block 4 11614 2903 0.6966 0.5999233
#> Fator2:Fator1 4 468257 117064 28.0861 4.580e-10 ***
#> Fator2:Fator1: Plantio 1 26112 26112 6.2648 0.0176146 *
#> Fator2:Fator1: R1+15 1 70896 70896 17.0095 0.0002467 ***
#> Fator2:Fator1: V1+15 1 258888 258888 62.1125 5.419e-09 ***
#> Fator2:Fator1: V3+15 1 112360 112360 26.9574 1.137e-05 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> -----------------------------------------------------------------
#> Analyzing F2 inside of the level of F1
#> -----------------------------------------------------------------
#>
#> Df Sum Sq Mean Sq F value Pr(>F)
#> Fator1 1 16160 16160 3.8772 0.057657 .
#> block 4 11614 2903 0.6966 0.599923
#> Fator1:Fator2 6 568651 94775 22.7385 2.974e-10 ***
#> Fator1:Fator2: IN 3 75470 25157 6.0356 0.002232 **
#> Fator1:Fator2: NI 3 493181 164394 39.4414 7.343e-11 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> -----------------------------------------------------------------
#> Final table
#> -----------------------------------------------------------------
#> Plantio R1+15 V1+15 V3+15
#> IN 272.8 aA 236.6 aA 140.4 bB 304.0 aA
#> NI 170.6 bB 68.2 bC 462.2 aA 92.0 bC
#>
#>
#> Averages followed by the same lowercase letter in the column and
#> uppercase in the row do not differ by the sk (p< 0.05 )
data=rbind(data.frame(trat=paste(f1,f2,sep = ""),bloco=bloco,resp=resp),
data.frame(trat=c("Test","Test","Test","Test","Test"),
bloco=unique(bloco),resp=respAd))
with(data,dunnett(trat = trat,
resp = resp,
control = "Test",
block=bloco,model = "DBC"))
#> Estimate IC-lwr IC-upr t value p-value sig
#> Test - INPlantio -34.2 -148.72047 80.32047 -0.8376 0.9478 ns
#> Test - INV1+15 -166.6 -281.12047 -52.07953 -4.0802 0.0019 *
#> Test - INV3+15 -3.0 -117.52047 111.52047 -0.0735 1.0000 ns
#> Test - INR1+15 -70.4 -184.92047 44.12047 -1.7242 0.4035 ns
#> Test - NIPlantio -136.4 -250.92047 -21.87953 -3.3405 0.0138 *
#> Test - NIV1+15 155.2 40.67953 269.72047 3.8010 0.0043 *
#> Test - NIV3+15 -215.0 -329.52047 -100.47953 -5.2655 0.0001 *
#> Test - NIR1+15 -238.8 -353.32047 -124.27953 -5.8484 0.0000 *
