This function generates p-rep designs with parameters:
v = 2ms(s-1) : number of treatments,
b_1 = 2(s-1) : first set of blocks of size k_1 = ms,
b_2 = 2s : second set of blocks of size k_2 = 2m(s-1) and
r = 3 : number of replications.
e = 2 : number of environments,
This function also generates incidence matrix, information matrix, canonical efficiency factor and average variance factor of the generated p-rep design.