#=======================================
# Using the simulate1 dataset
#=======================================
data("simulate1")
# Converting to readable format for function
dados=cbind(simulate1[simulate1$tempo==1,3],
simulate1[simulate1$tempo==2,3],
simulate1[simulate1$tempo==3,3],
simulate1[simulate1$tempo==4,3],
simulate1[simulate1$tempo==5,3],
simulate1[simulate1$tempo==6,3])
colnames(dados)=c(1,2,3,4,5,6)
dados
# Creating the treatment vector
resp=aacp(dados)
trat=simulate1$trat[simulate1$tempo==1]
# Analyzing by DIC function
DIC(trat,resp)
Run the code above in your browser using DataLab