periodicdummies function

Make matrix of periodicity (e.g. seasonal) dummies