replace for loops with apply in internal functions
update odds ratio calculation

Deconstruct glm and glmer approaches into separate functions;
rename pairwise.odds.ratios() to spaa();
change spaa() behaviour to allow OR specification using ‘method’;
clean object returned by spaa()
change arrangement of functions within .R files
clean help files, generate S3 objects

spaa() given plot, print and summary methods