Analysis of ordered categorical data with threshold models exemplified by plumage damage scores from laying hens differing in their genotype and rearing environment