This function generate an incomplete split-plot design for given number of main plot treatments, number of subplot treatments, number of blocks and number of subplot treatments in each main plot
ispd.cmis(v1,v2,b,k2)
number of main plot treatments
number of subplot treatments
number of block
number of subplot treatments in each main plot
A list containing parameters, design layout and column layout of design