data("aristolochia")
attach(aristolochia)
#=============================
# Use the DIC function
#=============================
DIC(trat, resp)
#=============================
# Use the DIC function noparametric
#=============================
DIC(trat, resp, test="noparametric")
#=============================
# Use the DIC.glm function
#=============================
resp=resp/4 # total germinated seeds
# the value 25 is the total of seeds in the repetition
DIC.glm(trat, cbind(resp,25-resp), glm.family="binomial")
Run the code above in your browser using DataLab