data("passiflora")
attach(passiflora)
a=with(passiflora, DBC(trat, bloco, MSPA))
b=with(passiflora, DBC(trat, bloco, MSR))
barplot_positive(a, b, var_name = c("DMAP","DRM"), ylab = "Dry root (g)")
a=with(passiflora, DIC(trat, MSPA,test = "noparametric"))
b=with(passiflora, DIC(trat, MSR))
barplot_positive(a, b, var_name = c("DMAP","DRM"), ylab = "Dry root (g)")
Run the code above in your browser using DataLab