<title>Background</title> <p>For several decades Risk Assessment has been conducted in the field of human and animal health [e.g. [##UREF##0##1##, ####UREF##1##2##, ##UREF##2##3##, ##UREF##3##4##, ##UREF##4##5####4##5##]]. The need to develop a formal means for Risk Analysis of animal welfare has been recognized at the European level [##UREF##5##6##,##UREF##6##7##]. A recently published report on qualitative Risk Assessment on intensively farmed calves [##UREF##7##8##,##UREF##8##9##] was an important step toward transparent decision making on animal welfare. The methodology, however, was also recognized to be in need of further modification [[##UREF##1##2##,##UREF##7##8##], p.8].</p> <p>In a separate paper [##UREF##9##10##] we reported on a critical analysis from a semantic-modelling perspective, and formulated recommendations for improvement of Risk Assessment applied to animal welfare, as presented in the EFSA report [##UREF##7##8##,##UREF##8##9##]. Semantic modelling is a kind of risk-benefit assessment, i.e. welfare assessment based on a structured analysis of available scientific information [##UREF##10##11##, ####UREF##11##12##, ##REF##12162649##13##, ##UREF##12##14####12##14##]. Several semantic models have successfully been 'validated' against expert opinion [##REF##12162650##15##, ####UREF##13##16##, ##UREF##14##17####14##17##]. In these studies, typically two sets of scores have been requested from experts: welfare scores for housing systems and scores for the importance of welfare-relevant system attributes, which we suggested being the equivalents of the 'hazards' in Risk Assessment [##UREF##9##10##]. Given their value in relation to semantic modelling, these expert-opinion scores probably provide a good starting point for representing expert reasoning about animal welfare. Conceptually, these sets of scores provide the first two steps in backward expert-reasoning from overall scores to the underlying scientific information: welfare scores for housing systems can, in principle, be explained by the attribute (i.e. hazard or risk) scores, which can be explained by the underlying science specifying relationships between two types of attributes, namely design criteria and welfare performance criteria [##UREF##10##11##].</p> <p>In order to check a number of critical points raised in the main study such as the need to specify definitions, the need to include positive (behavioural) aspects of welfare and to complement an assessment of risk components with a perception of overall welfare [##UREF##9##10##], this paper reports on a study comparing the scores for Hazard Characterization (HC), Exposure Assessment (EA) and Risk Characterization (RC) as presented in the EFSA report with semantic-modelling type scores elicited from experts about a selected number of welfare hazards and housing systems for calves. This paper also addresses several additionally suggested points for improvement of Risk Analysis [##UREF##9##10##], including the linking of information (such as between hazards and underlying scientific information, and between HC scores and overall welfare scores), reporting of uncertainty measures, verification of items possibly lacking from EFSA [##UREF##7##8##,##UREF##8##9##] and transparency about underlying values. Finally, the welfare scores given by the experts, which had all been involved or cited in the EFSA report [##UREF##7##8##,##UREF##8##9##], provide complementary information to decision makers on the welfare of calves, and also provides unique information on how groups of experts may differ in assessing animal welfare.</p> <p>The objectives of this paper, therefore, were to elicit expert opinion about calf welfare as part of a semantic modelling-type 'validation' study addressing the above-mentioned aspects of the Risk Assessment (RA) approach developed in the calf EFSA report [##UREF##7##8##,##UREF##8##9##].</p>
<title>Methods</title> <p>A survey was conducted in November-December 2006 by sending an email message to the authors of the EFSA report, to the veterinary experts who had given advice on Exposure Assessment (EA) and to a selected number of applied ethologists, who were the authors of papers cited in EFSA [##UREF##8##9##] (together representing three different roles in the EFSA report). In total 38 experts from 10 different (European and North-American) countries were contacted with the request to assess overall animal welfare of 11 housing systems (on a scale from 0, worst to 10, best) and 18 hazards (also on a scale from 0 to 10, i.e. least to most important for welfare). In the questionnaire, it was emphasized that only welfare was to be assessed, and that welfare could be defined as what matters to the animals from their point of view. The items were presented in a table-format (comparable to Tables ##TAB##1##2## and ##TAB##2##3## below, but providing the full description given in the EFSA report) in a randomized order. In addition, experts were asked to state their professional background and an opinion on the EFSA report [##UREF##7##8##,##UREF##8##9##]. Experts were then classified into those with a background in veterinary science, ethology, or of mixed background, i.e. with a background in both veterinary medicine and ethology. Item descriptions were identical to the ones used in the EFSA report [##UREF##7##8##,##UREF##8##9##], except for two newly added items in each list. White veal in baby boxes and suckler calves at pasture were added as 'controls' to the list of housing systems, and insufficient roughage and insufficient play were added to the list of hazards in order to examine the hypothesis that these are important systems and hazards not adequately addressed in the EFSA report [##UREF##7##8##,##UREF##8##9##] as indicated in Bracke et al. [##UREF##9##10##].</p> <p>Kendall's coefficient of concordance was calculated to determine agreement among experts, and Spearman's correlation coefficients (Rho) were used to determine relationships between median expert scores and Hazard Characterization (HC) scores for hazards, and between expert scores and overall Risk Characterization (RC) scores for housing systems. Hazard scores from the survey were compared with HC scores, because these are indicators of the potential importance of a hazard. In the EFSA report, HC scores were constant across housing systems. RC scores were calculated in the EFSA report by multiplying HC and EA (Exposures Assessment) scores (see Table ##TAB##0##1## for an overview of abbreviations used). RC scores indicate various levels of hazard exposure and risk related to hazards in different housing systems. Overall risk per housing system was calculated from the median and total (i.e. the sum of components) RC scores reported for each housing system in EFSA [##UREF##7##8##,##UREF##8##9##]. These sets of scores differed because not all hazards were scored for all systems. Both scores only give a rough idea of the overall risk, as the underlying scales were not cardinal (i.e. the interval between successive points of the scale may not have been constant).</p> <p>The statistical analyses were done in SPSS 13.0 [##UREF##15##18##]. Correlations for housing systems were expected to be negative, because higher expert scores implied higher welfare, whereas higher RC scores implied more risk for welfare, i.e. lower welfare.</p> <p>To determine main factor effects on the scores given by the experts, a components of variance model was used [##UREF##16##19##], initially ignoring the fact that scores ranged from 0 to 10. The model comprised random effects for experts and fixed effects for Hazard/Housing system, Role and Background as main effects. The additional factor Gender (of the expert) and two-factor interactions were systematically tested, dropping additional factor combinations when not significant. The most relevant models were subsequently analyzed with a threshold model comprising the aforementioned fixed and random effects. The estimation procedure is discussed in Keen and Engel [##UREF##17##20##] where it is shown that this model is appropriate for analyzing ordered scores. In the analyses, the following factors were considered: Housing system (n = 11) or Hazard (n = 18); Background (veterinarian, n = 8; applied ethologist, which often combined the study of animal behaviour and animal science, n = 11; and mixed background, which were mostly veterinarians working as applied ethologist, n = 5); Gender (male, n = 16; female, n = 8); Role (i.e. role of involvement in the writing of the EFSA report [##UREF##8##9##]; these included Working Group member, i.e. authors of the report, n = 3; veterinarian contributing to Exposure Assessment, n = 5; other contacted expert, e.g. by being acknowledged in the report, n = 4; and author of a reference quoted in EFSA [##UREF##8##9##], n = 12 for housing systems and n = 11 for hazards). The interaction between Role and Background could not be examined because they were confounded, e.g. reference authors were all ethologists and exposure assessors were all veterinarians (see Table ##TAB##1##2##).</p> <p>Significance levels were determined with Wald tests employing a chi-square approximation [##UREF##18##21##]. Calculations were performed with GenStat [##UREF##19##22##].</p>
<title>Results</title> <p>The response rate of the questionnaire was 63% (n = 24 respondents for housing systems and n = 23 for hazards), comprising 3 Working Group members, 5 exposure assessors and 16 other scientific experts. In total, ten experts were positive about the Risk Assessment approach in the EFSA report. The other experts either did not respond to this question or stated that they were not familiar with the report. Working Group members (i.e. authors of the report) and exposure assessors generally responded positively, whilst 70% of respondents not personally involved (i.e. only through having a reference cited in the EFSA report) indicated that they were not familiar with the report (that had only recently been published at the time the survey was conducted). Several experts expressed doubt about the scientific value of the questionnaire (e.g. for requesting an instantaneous response without prolonged contemplation). Several experts complained about the vague descriptions of the housing systems, and some experts perceived hazards to be non-uniform (e.g. castration versus floors).</p> <p>Figures ##FIG##0##1## and ##FIG##1##2## give boxplots of the housing and hazard scores given by the experts, grouped by their professional background. Figure ##FIG##0##1##, for example, shows that median welfare scores for the housing system 'baby boxes' were 0.0, 6.0 and 0.0 for ethologists, veterinarians and experts with a mixed (veterinary and ethological) background respectively.</p> <p>Table ##TAB##2##3## shows Kendall's coefficient of concordance (W) for the scores given to housing systems and hazards by background. When considered together, there was low (W = 0.29; W = 0.18), but highly significant (P &lt; 0.001) agreement among the whole group of experts. Agreement was generally less significant when examined within the smaller subgroups of experts with different backgrounds. Table ##TAB##2##3## also shows Spearman's correlation coefficients (Rho) between (group and subgroup) expert opinion scores and EFSA scores (i.e. HC scores for hazards and RC scores for housing systems respectively). Significant correlations were found only for HC scores (reported in EFSA) and (the hazard scores given by) veterinarians (Rho = 0.57), for HC scores and mixed-background experts (Rho = 0.66), and for median RC scores and (the housing system scores given by) veterinarians (Rho = -0.68).</p> <p>Figure ##FIG##2##3## illustrates two relationships found for HC and hazard scores for experts with different backgrounds, namely ethologists (where the relationship was not significant) and veterinarians (where Rho was significant, namely 0.57, see Table ##TAB##2##3##). Figure ##FIG##2##3## shows hazards that received a high HC score in EFSA (2006a, b), but received relatively low expert scores (for both veterinarians and ethologists), such as light (Li) and mixing of calves (Mi). It also illustrates the reverse, especially for access to a natural teat (Te) (particularly for ethologists) and for education (Ed), bedding (Be) and floor (Fl) (both types of expert).</p> <p>In the components of variance models, effects of Gender (main effects and interactions) were neither significant for housing-system scores nor for hazard scores. For hazard scores, no significant interactions were found, resulting in a model with main effects for Hazard (P = 0.00), Role (P = 0.33) and Background (P = 0.08). For housing-system scores, the final model comprised Housing system (P &lt; 0.05), Role (P = 0.01), Background (P &lt; 0.05) and the interaction between Housing system and Background (P &lt; 0.01).</p> <p>In the final threshold model for hazard scores, only Hazard was significant (P &lt; 0.001; see Table ##TAB##4##5##). Role was not significant, and a trend was found for Background (P = 0.06). Respondents with a mixed background tended to give higher hazard scores than veterinarians, and ethologists gave intermediate scores that were closer to the mixed-background group.</p> <p>According to the experts, the least important hazards were insufficient human contact, separation from the dam, overfeeding and lack of maternal care (Table ##TAB##4##5##). These did not significantly differ from each other, and scored significantly lower than all other hazards, except for light which was intermediate. A whole range of hazards with somewhat higher scores did not significantly differ from each other. The 6 most important hazards were underfeeding, inadequate colostrum intake, poor education, insufficient space, inadequate roughage and iron deficiency (in that order, see Table ##TAB##4##5##). In this list, underfeeding scored significantly higher than iron deficiency. The two hazards added for 'validation' purposes, insufficient space to play and inadequate roughage, ended up in the middle and middle-upper range respectively.</p> <p>In the final threshold model for Housing-system scores, the interaction between Housing system and Background failed to reach significance. This left a model with only main effects for Housing system (P &lt; 0.001; see Table ##TAB##3##4##), Role (P = 0.03) and Background (P = 0.03).</p> <p>The various veal and feedlot systems received the lowest scores, with (white veal calves in) Baby boxes (the system added for 'validation' purposes as a negative control) scoring significantly lower than the other systems (see Table ##TAB##3##4##). Pink veal and white veal suckling from a dam scored significantly higher than similar bucket-fed groups of white veal calves. Suckler beef calves kept on pasture (the system added for 'validation purposes as a positive control) scored significantly higher than all other systems.</p> <p>Veterinarians gave significantly higher overall welfare scores for housing systems compared with mixed-background experts and ethologists, but the latter did not differ from each other.</p> <p>Working Group members (i.e. authors of the EFSA report) did not significantly differ from reference authors, but Working Group members did give significantly higher overall welfare scores than veterinary exposure assessors and contacted experts.</p>
<title>Discussion</title> <p>The objectives of this paper were to elicit expert opinion about calf welfare and to verify conclusions from our previous analysis [##UREF##9##10##] of the new Risk Assessment (RA) approach developed in the calf-welfare report of the European Food Safety Authority [##UREF##7##8##,##UREF##8##9##]. This paper reports a first validation-type study of Risk Assessment applied to animal welfare, which is a methodology in need of further refinement [[##UREF##7##8##], p.8; [##UREF##1##2##]] to which end recommendations from a semantic-modelling perspective have been formulated [##UREF##9##10##].</p> <title>Methods</title> <p>A semantic-modelling type questionnaire [##REF##12162650##15##, ####UREF##13##16##, ##UREF##14##17####14##17##] was sent to experts, requesting 'intuitive' welfare scores for housing systems and hazards on scales from 0 to 10. The total number of experts was limited. The experts in this study were all applied ethologists or veterinary scientists that had been involved or cited in the EFSA [##UREF##8##9##] report on the welfare of intensively-reared calves. These scientists had been identified in the EFSA report as <italic>the </italic>experts on this subject in Europe. From a semantic modelling (SM) perspective, however, the term 'experts' must be qualified, because the respondents all had a particular area of expertise (rather than being complete and fully impartial generalists), and few experts had experience with (the technical details of) (semi-)quantified overall welfare assessment as developed, for example, in SM. This may limit the extent to which the survey can be regarded as a 'gold standard'. In the section 'hazards' below this point will be further illustrated with the example of 'underfeeding'.</p> <p>In response to the questionnaire, several experts questioned its scientific value. Perhaps these respondents had not fully realized that in this study <italic>they </italic>were the experimental subjects. By virtue of being knowledgeable experts, their opinion, even when elicited in this stimulus-response like fashion, was inherently valid (by being their expert opinion), especially also because uncertainties about the scores were to become part of the (biological) variation around the group's opinions. These respondents' complaint, however, indicates a legitimate concern about a risk of misinterpretation of the outcomes of this study. When experts believe that one item, housing system or hazard, is better or more important than another, it does not logically follow that it actually is better or more important. For the latter conclusion, further scientific studies are needed, esp. including measurements of animal-based attributes. SM subscribes to that view, but also recognizes that an assessment of animal welfare is always an assessment from a human's point of view [##UREF##20##23##]. It is rarely possible to assess overall welfare within a single scientific study, and it always requires taking a range of (animal- and environment-based) measures that must be selected and interpreted within the context of decades of scientific research. As far as we know, the most structured way available at present to move towards that objective is semantic modelling.</p> <p>Furthermore, the respondents validly complained about the vague and general descriptions provided for hazards and housing systems. Unfortunately, these were unavoidable in this study because they had to be adopted from the EFSA report [##UREF##7##8##,##UREF##8##9##] that was under scrutiny. As indicated in the underlying study [##UREF##9##10##] also from a SM perspective, more detailed descriptions would be required: hazards should be specified in relation to the underlying scientific information and housing systems should be specified using a matrix of welfare-relevant attributes covering the range of conditions prevailing in the housing systems in the assessment domain, including both environment-based inputs and animal-based outcomes covering all welfare-relevant needs [##REF##12162649##13##,##REF##12162650##15##].</p> <p>A further methodological issue concerns the concept of Risk. Risk may differ from welfare assessment in that a risk to welfare may or may not actually compromise welfare, depending on the (negative welfare) effects actually occurring. However, because both survey and EFSA report [##UREF##7##8##,##UREF##8##9##] concerned the European scale, the population of farms was sufficiently large to assume that risks and their associated effects on welfare were (more or less) referring to the same properties of the system. Exactly which properties the respondents considered cannot be determined from this survey. As the respondents were familiar with the housing systems and hazards (as they were experts who had been asked to abstain from scoring when they were not familiar with it) and as they were asked to respond without much contemplation, it is reasonable to assume that in most cases the scores were given for typical, representative examples of systems and hazards.</p> <p>An important caveat with respect to the interpretation of housing-system (and hazard) scores, however, is that the scores were given for the experts' personal interpretation of welfare. Even though welfare was defined in the survey as what matters to the animals from their point of view, differences in interpretation may have contributed to variation in the scores. By contrast, whereas the scores were probably given for 'average' systems, this survey did not address the potentially much larger range of variation existing between individual farms within type of system. In relation to this variation, one expert commented that 'a good farmer can produce good welfare in a poor system'. Though this statement can be challenged, the reverse is certainly true: a bad farmer will cause poor welfare in what is otherwise a good system. Therefore, the scores reported here for the different types of housing systems and hazards cannot be taken to represent welfare scores for all individual cases, and further work is needed to address this point.</p> <p>Finally, welfare scores for housing systems and hazards were expressed on a scale from 0 to 10. A median score of 5 was previously found to be the cut-off point for acceptability proposed by experts who had given welfare scores for enrichment materials for pigs [##UREF##13##16##]. This supports a tentative suggestion to use some score in the middle of the scale (somewhere around 5) as the (implicit) cut-off point for what the experts in the present survey may have considered acceptable/important, also because this would be in accordance with its familiar use as a grading scale, e.g. in schools and psychological tests.</p> <title>General 'validation'</title> <p>Kendall's coefficients of concordance (W) were highest (0.34, P &lt; 0.001) for hazard scores given by mixed-background experts, which is explained by the fact that these were experts that had been involved as authors of the EFSA report [##UREF##8##9##]. Otherwise, W values were low for both housing system and hazard scores (range 0.09 to 0.29, Table ##TAB##2##3##) compared to similar welfare scores for pregnant sows (W = 0.73 and 0.43 for housing systems and attributes respectively, [##REF##12162650##15##]). Nevertheless, W values were highly significant for the whole group of experts, probably due to the larger number of individuals in the dataset. Vague item descriptions and the request to provide intuitive scores may have contributed to this finding. More contemplation about better specified items, e.g. in working group discussions, may improve the level of concordance (but see [##UREF##9##10##]), and this could be monitored with a semantic-modelling type questionnaire. As long as the objective of complete consensus has not been reached, the level of concordance among the experts and the degree of variation in the scores given may provide an entry for specifying the level of uncertainty for scores given in Risk Assessment.</p> <p>Compared to previous studies validating semantic models against expert opinion [##REF##12162650##15##, ####UREF##13##16##, ##UREF##14##17####14##17##], this study yielded moderate correlations for hazards (range 0.28–0.66) and relatively poor and many not significant correlations for housing systems. The correlation for HC scores was highest for experts with a mixed background, followed by veterinarians (Rho: 0.66 and 0.57, both P &lt; 0.05). Surprisingly, the median scores for hazard importance provided by ethologists did not correlate significantly with the HC scores reported in EFSA [##UREF##7##8##,##UREF##8##9##]. This may be explained by the confounding relationship with Role: many ethologists had not directly been involved in the writing of the report (they had only been cited), whereas veterinarians and mixed-background experts in this study had been actively involved as exposure assessors and as authors of the report respectively.</p> <p>Poor correlations between expert opinion and Risk Assessment may reflect the latter's focus on negative hazards, rather than on both negative and positive welfare aspects, and it may reflect the RA's focus on component hazards rather than on overall (risk for poor) welfare (both as indicated in [##UREF##9##10##]). With respect to the representation of overall risk, it may be noted that all reported correlations of median expert scores with total RC scores were lower than the corresponding correlations with median RC scores (see Table ##TAB##2##3##). This may indicate that the procedure followed in the EFSA report [##UREF##7##8##,##UREF##8##9##] of leaving out some hazards for some housing systems reduced its suitability to derive overall welfare, as indicated in this study by expert opinion (but note that this was not an objective of the EFSA report, while it has been proposed from a semantic modelling perspective, [##UREF##9##10##]).</p> <p>Veterinarians were the only group that showed a significant correlation with overall risk related to housing systems, namely -0.68 for median RC scores. Although this may suggest added value of consulting veterinarians in Risk Assessment as described in EFSA [##UREF##7##8##,##UREF##8##9##], it may also simply reflect their involvement in the report or a health-related underlying value in the EFSA [##UREF##7##8##,##UREF##8##9##] report (see [##UREF##9##10##]).</p> <title>Hazards</title> <p>In this study, experts with a mixed background tended to give higher hazard scores than veterinarians, and ethologists gave intermediate scores close to the mixed-background experts. This could well indicate that welfare scientists may attach more importance to animal welfare than veterinarians do.</p> <p>According to the experts, the least important hazards were insufficient human contact, separation from the dam, overfeeding and lack of maternal care. All other hazards had median scores of at least 6.0. The most important hazards (median scores &gt; 6.5) were underfeeding, inadequate colostrum intake, poor education, insufficient space, insufficient roughage and iron deficiency, inadequate/inappropriate ventilation, poor floor conditions and no bedding (in that order). This list only partially confirms the analysis in EFSA ([##UREF##7##8##,##UREF##8##9##]; see also Figure ##FIG##2##3##), especially with respect to the importance of colostrum intake and inadequate ventilation. Insufficient light and mixing of calves were found to be much less important in the survey compared with the HC scores reported in EFSA [##UREF##7##8##,##UREF##8##9##]; e.g. mixing of calves was identified there as a main risk for calf welfare). As can be noted from the Boxplot shown in Figure ##FIG##1##2##, experts with a mixed background gave relatively high scores for these two hazards. Subsequent data exploration (not shown) indicated that Working Group members might have accounted for this difference, indicating that the discrepancy for these two hazards would be even larger if Working Group members who had written the report had been excluded from the analysis. This is in accordance with a previous suggestion [##UREF##9##10##], that the EFSA results may be diverging from current expert opinion. This is also true for several other hazards such as stockman education and to some extent (particularly for ethologists) access to a natural teat, which, conversely, seem to have been considered more important by the consulted experts than indicated by their HC scores reported in EFSA [##UREF##7##8##,##UREF##8##9##]. Furthermore, the median score of 7.0 for 'poor floor conditions' supports its ranking as 4<sup>th </sup>most important hazards-class in Anonymous [##UREF##10##11##] and suggests a higher importance compared to the scores given in the EFSA report, where this hazard had been divided into 5 component hazards (see [##UREF##9##10##]). In addition, roughage was identified by the experts as an important hazard (median: 7.0), especially by experts with a mixed background (see Figure ##FIG##1##2##). This item had been added to the list, because it was considered to be either lacking from the EFSA report, or inadequately referred to by the hazard 'insufficiently balanced solid food' (HC = 3), again confirming our analysis in Bracke et al. [##UREF##9##10##]. The present study, however, did not confirm a similar hypothesis for the added hazard 'space to play' (median score: 6.0), which was rated as of average importance only (though it was still scored as more, but not significantly more, important than insufficient light and mixing of calves). A specific explanation for this finding cannot be provided, because experts did not specify the reasons for their scores (for feasibility reasons).</p> <p>In this survey, underfeeding was the most important hazard. This may not appear to be surprising, because food has been identified as the 'gold standard' resource in consumer demand studies [##UREF##21##24##]; feed refusal is often a first and important sign of illness; and food is a necessary requirement for survival, growth, health and (re-)production. Given these scientific arguments it is surprising that, previously, underfeeding had not been identified as a main hazard by a group of 22 experts [##UREF##10##11##], and that it had received a Hazard Characterization (HC) score of only 4 on the 5 point scale in the EFSA (2006a, b) report. In the report, underfeeding was given the same HC score as, for example, high humidity, poor air quality (ammonia, dust) and continuous restocking (no all-in, all-out), but it received a lower HC score than, for example, inadequate ventilation, poor air quality (H<sub>2</sub>S), insufficient space, insufficient light, social isolation and mixing of calves from different sources. It would seem difficult, if not impossible, to justify these differences based on available scientific evidence. A possible explanation for the absence of underfeeding in Anonymous [##UREF##10##11##] can be found in the EFSA [##UREF##7##8##,##UREF##8##9##] report, where, in accordance with expectation, overall risk associated with underfeeding was low (even classified as 'negligible risk'), because whereas the effect (HC) was reasonably high, the occurrence probability, i.e. Exposure Assessment (EA) scores, were low (1 or 2 on a 5 point scale). In other words, in intensive calf rearing systems, aimed at maximized growth, underfeeding is rare. This corresponds to the procedure described in semantic modelling (SM) to exclude the above mentioned scientific evidence in the formal calculation of the weighting factor for underfeeding (which is equivalent to the HC score in Risk Assessment, see [##UREF##9##10##]), because it does not apply to the assessment domain where good farming practices were assumed. Technically, when the assessment domain only contains housing systems where animals are provided with sufficient food (as is normally the case in modern production systems), underfeeding is in fact much less important than at first suggested.</p> <title>Housing systems</title> <p>The two 'control' systems added to the list of systems taken from EFSA [##UREF##7##8##,##UREF##8##9##] as part of our 'validation' effort [##UREF##9##10##], i.e. baby boxes and suckler beef calves at pasture, indeed obtained the lowest and highest predicted mean overall-welfare score, and they also differed significantly from all other systems. However, whereas the latter could be regarded as a true positive control, defining the upper range of calf welfare, the former system, baby boxes, cannot be regarded as a true negative control, as will be explained below.</p> <p>Our finding that Working Group members did not significantly differ from reference authors indicates that the authors of the EFSA [##UREF##7##8##,##UREF##8##9##] report were 'in line' with the authors of their sources. Significantly higher scores given by Working Group members compared with veterinary exposure assessors and contacted experts were mainly due to higher scores for the 4 high-welfare systems (groups of dairy and sucker beef calves).</p> <p>Veterinarians gave significantly higher scores for housing systems than either mixed-background experts or ethologists. As this was especially the case for the low-welfare (veal and feedlot) systems, the finding may correspond with their lower scores for hazard importance, confirming that veterinarians may have been less concerned by welfare problems in intensive systems for rearing calves than applied ethologists, whether or not they had a veterinary background. Other explanations, however, are also possible, e.g. that different definitions of animal welfare were used (despite the fact that welfare had been defined in the survey's instructions), perhaps involving a different weighting of welfare aspects (e.g. physical versus mental health; physiological versus behavioural needs). Such differences would be expected between experts with different backgrounds, given, for example, the fact that many years of dispute has not yet resulted in a commonly accepted definition of animal welfare among ethologists [##UREF##10##11##,##UREF##22##25##,##UREF##9##10##].</p> <p>Veterinarians gave higher median scores to each of the 5 veal systems, especially for baby boxes, compared with ethologists and experts with a mixed background (see Figure ##FIG##0##1##). Median veterinary scores did not drop below 5.0 for any housing system. Their lowest median scores, given for feedlots, white veal in small groups, and baby boxes, were 5.0, 5.3 and 6.0 respectively. This implied that baby boxes were not a negative control system for veterinarians, because they gave lower (though not significantly lower) scores to the two other systems. By contrast, animal welfare experts, i.e. ethologists and experts with a veterinary background working in applied ethology, were much more negative about these three systems (medians between 0.0 and 4.0), and ethologists gave scores below 6 also to other veal systems, to hutches and to small groups of dairy calves (see Figure ##FIG##0##1##). This apparent difference may be related to differences in professional experience and affinity to health and production in the sector (despite what was claimed about the veterinarians' independence in the EFSA report, see [##UREF##9##10##]). This hypothesis could, for example, explain why veterinarians gave relatively high scores for calves kept in baby boxes (and hutches), as individual housing promotes hygiene. It could, perhaps, also explain why they identified feedlots, an American system which is not prevalent in Europe, as the worst system. Finally, it could explain that veterinarians showed lower median scores for access to a natural teat and for castration/dehorning (see Figures ##FIG##1##2## and ##FIG##2##3##), because natural teat sucking is a typical behavioural requirement and castration/dehorning is very much part of routine veterinary practice.</p>
<title>Conclusion</title> <p>This paper reported a 'validation' study of the EFSA report, comparing its scores for hazard characterization (HC) and Risk Characterization (RC) with semantic-modelling type scores elicited from a limited number of experts about a selected number of welfare hazards and housing systems for calves according to recommendations formulated in the underlying paper [##UREF##9##10##]. Experts included ethologists and veterinarians involved in the publication of the EFSA report, as well as authors of publications cited in that report.</p> <p>Differences in welfare score were found between housing systems, between hazards, between experts with different types of involvement in the EFSA report [##UREF##7##8##,##UREF##8##9##], and between experts with different backgrounds.</p> <p>The poorest welfare systems (with median scores &lt; 5.5) were: white veal housed individually in baby boxes (first 6–8 weeks), bucket fed (i.e. not suckling); white veal in small groups, bucket fed (i.e. not suckling); feed lots (high density groups within outside pen); and white veal in larger groups, with automatic feeding (i.e. not suckling). The best system was: suckler beef calves kept with cows in a herd at pasture.</p> <p>The most important hazards (with median scores &gt; 6.5) were: underfeeding; inadequate colostrum intake; poorly educated stockperson; insufficient floor space allowance; insufficient/inadequate roughage; iron deficiency; inadequate/inappropriate ventilation; poor floor conditions; and no bedding.</p> <p>This study provided only limited support of the scores reported in EFSA [##UREF##7##8##,##UREF##8##9##] and emphasized a number of suggestions identified in Bracke et al. [##UREF##9##10##] to improve the EFSA [##UREF##7##8##,##UREF##8##9##] Risk Assessment from the perspective of semantic modelling. Hazards and housing systems should be specified in more detail in relation to the available scientific information. Risk Characterization (RC) scores could be linked to overall welfare scores for housing systems, which are inherently lacking in Risk Assessment. Measures of uncertainty should be included, such as those based on the level of concordance among the experts or based on the variation in scores. Rules could be specified for including and leaving out items (hazards or housing systems) from a list, e.g. by including positive 'hazards', i.e. benefits for welfare such as play and roughage for calves, and by including items that define the end points of the scale, such as illustrated by the positive and negative 'controls' in the list of housing systems. Finally, transparency about underlying values is called for, as this study confirmed some widely perceived differences in professional backgrounds of experts and the impact of these differences on their perception of animal welfare.</p> <p>These points illustrate that, for further development, Risk Assessment applied to animal welfare can benefit from semantic modelling principles, and that both can benefit from improved definitions of important concepts such as welfare, its components, hazards and risks.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Recently, a Risk Assessment methodology was applied to animal welfare issues in a report of the European Food Safety Authority (EFSA) on intensively housed calves.</p> <title>Methods</title> <p>Because this is a new and potentially influential approach to derive conclusions on animal welfare issues, a so-called semantic-modelling type 'validation' study was conducted by asking expert scientists, who had been involved or quoted in the report, to give welfare scores for housing systems and for welfare hazards.</p> <title>Results</title> <p>Kendall's coefficient of concordance among experts (n = 24) was highly significant (P &lt; 0.001), but low (0.29 and 0.18 for housing systems and hazards respectively). Overall correlations with EFSA scores were significant only for experts with a veterinary or mixed (veterinary and applied ethological) background. Significant differences in welfare scores were found between housing systems, between hazards, and between experts with different backgrounds. For example, veterinarians gave higher overall welfare scores for housing systems than ethologists did, probably reflecting a difference in their perception of animal welfare.</p> <p>Systems with the lowest scores were veal calves kept individually in so-called "baby boxes" (veal crates) or in small groups, and feedlots. A suckler herd on pasture was rated as the best for calf welfare. The main hazards were related to underfeeding, inadequate colostrum intake, poor stockperson education, insufficient space, inadequate roughage, iron deficiency, inadequate ventilation, poor floor conditions and no bedding. Points for improvement of the Risk Assessment applied to animal welfare include linking information, reporting uncertainty and transparency about underlying values.</p> <title>Conclusion</title> <p>The study provides novel information on expert opinion in relation to calf welfare and shows that Risk Assessment applied to animal welfare can benefit from a semantic modelling approach.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>MBMB did most of the data collection and writing. BE en WGB were involved as statisticians. SEA and BA provided most helpful comments and guidance in setting up the work and in getting it published. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>Many thanks to the participating respondents, and to Jos Metz, Jos Noordhuizen, David Morton and Jeff Rushen for their suggestions.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Boxplot of welfare scores for housing systems by background (see also Table 2, n = 24 experts)</bold>. Asterisks and circles indicate two types of outliers identified as standard practice in SPSS. Outliers are scores with values between 1.5 and 3 box lengths from the upper or lower edge of the box. The box length is the interquartile range (i.e. median 25% to 75% of values), while the horizontal line in the box indicates the median value. The two curved lines are connecting median values of ethologists (solid line) and veterinarians (dashed line) respectively.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Boxplot of scores for hazard importance by background (see also Table 3; n = 23 experts)</bold>. Asterisks and circles indicate two types of outliers as standard practice in SPSS. Outliers are scores with values between 1.5 and 3 box lengths from the upper or lower edge of the box. The box length is the interquartile range (i.e. median 25% to 75% of values), while the horizontal line in the box indicates the median value.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Scatter plot of HC scores (horizontal axis) and median hazard scores (y-axis) given by veterinarians (triangles) and ethologists (stars)</bold>. Hazard codes: Be: Bedding; Ca: Castration; Co: Colostrum; Da: Dam; Ed: Education; Fl: Floor; Hb: Haemoglobin; Hu: Human contact; Li: Light; Ma: Maternal care; Mi: Mixing; Of: Overfeeding; Pl: Play; Ro: Roughage; Sp: Space; Te: Teat; Uf: Underfeeding; Ve: Ventilation (see also Table 5).</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Overview of abbreviations used</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Abbreviation</td><td align="left">Meaning</td></tr></thead><tbody><tr><td align="left">EA</td><td align="left">Exposure assessment</td></tr><tr><td align="left">EFSA</td><td align="left">European Food Safety Authority</td></tr><tr><td align="left">HC</td><td align="left">Hazard characterisation</td></tr><tr><td align="left">RC</td><td align="left">Risk characterisation</td></tr><tr><td align="left">SM</td><td align="left">Semantic modelling</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Specification of numbers of respondents according to their background and their role in the writing of the EFSA (2006b) report.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="3">Background</td><td/></tr></thead><tbody><tr><td align="left">Role in EFSA report</td><td align="center">Vet</td><td align="center">Ethol.</td><td align="center">Mixed</td><td align="center"><italic>Total</italic></td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left">Author of EFSA report</td><td align="center">0</td><td align="center">0</td><td align="center">3</td><td align="center">3</td></tr><tr><td align="left">Veterinarian involved in EA</td><td align="center">5</td><td align="center">0</td><td align="center">0</td><td align="center">5</td></tr><tr><td align="left">Contacted expert</td><td align="center">3</td><td align="center">0</td><td align="center">1</td><td align="center">4</td></tr><tr><td align="left">Cited reference author</td><td align="center">0</td><td align="center">11</td><td align="center">1</td><td align="center">12</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="right"><italic>Total</italic></td><td align="center">8</td><td align="center">11</td><td align="center">5</td><td align="center">24</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Agreement among experts (expressed as W, Kendall's coefficients of concordance, for welfare scores given to the 11 housing systems and to the 18 hazards in the questionnaire), and agreement between experts and EFSA report (expressed as Rho, Spearman's rank correlation coefficients, between median expert scores and hazard/risk characterisation)</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left" colspan="3">Housing systems</td><td align="left" colspan="3">Hazards</td><td align="left" colspan="2">Hazards</td><td align="left" colspan="4">Housing systems</td></tr><tr><td/><td colspan="12"><hr/></td></tr><tr><td align="left">Type of expert</td><td/><td/><td/><td/><td/><td/><td align="left" colspan="2">HC</td><td align="left" colspan="2">Median RC</td><td align="left" colspan="2">Total RC</td></tr></thead><tbody><tr><td/><td align="left">W</td><td align="left">P</td><td align="left">n</td><td align="left">W</td><td align="left">P</td><td align="left">n</td><td align="left">Rho</td><td align="left">P</td><td align="left">Rho</td><td align="left">P</td><td align="left">Rho</td><td align="left">P</td></tr><tr><td/><td colspan="12"><hr/></td></tr><tr><td align="left">All experts</td><td align="left">0.29</td><td align="left">0.00</td><td align="left">18</td><td align="left">0.18</td><td align="left">0.00</td><td align="left">21</td><td align="left">0.47</td><td align="left">0.06</td><td align="left">-0.54</td><td align="left">ns</td><td align="left">-0.45</td><td align="left">ns</td></tr><tr><td align="left">Ethologists</td><td align="left">0.23</td><td align="left">0.01</td><td align="left">9</td><td align="left">0.11</td><td align="left">0.04</td><td align="left">9</td><td align="left">0.28</td><td align="left">ns</td><td align="left">-0.03</td><td align="left">ns</td><td align="left">-0.22</td><td align="left">ns</td></tr><tr><td align="left">Veterinarians</td><td align="left">0.10</td><td align="left">ns</td><td align="left">5</td><td align="left">0.09</td><td align="left">ns</td><td align="left">7</td><td align="left">0.57</td><td align="left">0.02</td><td align="left">-0.68</td><td align="left">0.05</td><td align="left">-0.36</td><td align="left">ns</td></tr><tr><td align="left">Mixed background</td><td align="left">0.15</td><td align="left">ns</td><td align="left">4</td><td align="left">0.34</td><td align="left">0.00</td><td align="left">5</td><td align="left">0.66</td><td align="left">0.01</td><td align="left">0.13</td><td align="left">ns</td><td align="left">-0.08</td><td align="left">ns</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Descriptions of housing systems, their median scores and significance levels (Sig.) according to the final threshold model (see text).</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Description of housing systems</td><td align="center">Median</td><td align="left">Sig.</td></tr></thead><tbody><tr><td align="left">White veal housed individually in baby boxes (first 6–8 weeks), bucket fed (ie not suckling)</td><td align="center">3.00</td><td align="left">a</td></tr><tr><td align="left">White veal in small groups, bucket fed (ie not suckling)</td><td align="center">4.00</td><td align="left">bc</td></tr><tr><td align="left">Feed lots (high density groups within outside pen)</td><td align="center">4.00</td><td align="left">c</td></tr><tr><td align="left">White veal in larger groups, with automatic feeding (ie not suckling)</td><td align="center">5.00</td><td align="left">ce</td></tr><tr><td align="left">White veal in small groups, suckling</td><td align="center">5.50</td><td align="left">deg</td></tr><tr><td align="left">Pink veal in small groups, bucket fed + (some) solid fods, not suckling</td><td align="center">6.00</td><td align="left">fg</td></tr><tr><td align="left">Hutches outside with replacement dairy calves, bucket fed (not suckling) + solid foods, weaned at 2–3 months</td><td align="center">6.00</td><td align="left">gi</td></tr><tr><td align="left">Small groups of replacement dairy calves, bucket fed (not suckling) + solid foods, weaned at 2–3 months</td><td align="center">7.00</td><td align="left">hik</td></tr><tr><td align="left">Groups of dairy calves with an automatic feeding system (not suckling) + solid foods, weaned at 2–3 months</td><td align="center">7.00</td><td align="left">ik</td></tr><tr><td align="left">Suckler beef calves in groups kept inside, led twice a day to the dam for suckling up to 6–9 months</td><td align="center">7.00</td><td align="left">jk</td></tr><tr><td align="left">Suckler beef calves kept with cows in a herd at pasture</td><td align="center">9.50</td><td align="left">l</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>Descriptions of hazards, their median scores and significance levels (Sig.) according to the final threshold model (see text).</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Description of hazards</td><td align="left">Median</td><td align="left">Sig.</td></tr></thead><tbody><tr><td align="left">Insufficient contact with humans</td><td align="left">4.00</td><td align="left">a</td></tr><tr><td align="left">Separation from the dam</td><td align="left">4.00</td><td align="left">ac</td></tr><tr><td align="left">Too rich diet (overfeeding)</td><td align="left">5.00</td><td align="left">ac</td></tr><tr><td align="left">Lack of maternal care</td><td align="left">5.00</td><td align="left">ac</td></tr><tr><td align="left">Insufficient light</td><td align="left">6.00</td><td align="left">bce</td></tr><tr><td align="left">Mixing calves from different sources</td><td align="left">6.00</td><td align="left">deg</td></tr><tr><td align="left">No access to natural or artificial teat</td><td align="left">6.50</td><td align="left">eg</td></tr><tr><td align="left">Insufficient space for natural play (eg running and gamboling)</td><td align="left">6.00</td><td align="left">eg</td></tr><tr><td align="left">No bedding</td><td align="left">7.00</td><td align="left">dg</td></tr><tr><td align="left">Poor floor condition</td><td align="left">7.00</td><td align="left">fg</td></tr><tr><td align="left">Castration/dehorning, no anaesthetics</td><td align="left">6.00</td><td align="left">fi</td></tr><tr><td align="left">Inadequate/inappropriate ventilation</td><td align="left">7.00</td><td align="left">fi</td></tr><tr><td align="left">Iron deficiency</td><td align="left">7.00</td><td align="left">fi</td></tr><tr><td align="left">Insufficient/inadequate roughage</td><td align="left">7.00</td><td align="left">gik</td></tr><tr><td align="left">Insufficient floor space allowance</td><td align="left">7.50</td><td align="left">gik</td></tr><tr><td align="left">Poorly educated stockperson</td><td align="left">7.50</td><td align="left">gik</td></tr><tr><td align="left">Inadequate colostrum intake</td><td align="left">8.25</td><td align="left">hik</td></tr><tr><td align="left">Underfeeding</td><td align="left">9.00</td><td align="left">jk</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>Vet: veterinarian; Ethol.: Applied ethologist; Mixed: background both as Vet and as applied ethologist.</p></table-wrap-foot> <table-wrap-foot><p>P: significance level; ns: not significant; n: number of experts without missing values; HC: Hazard Characterisation; RC: Risk Characterization.</p></table-wrap-foot> <table-wrap-foot><p>For significance levels (Sig.), systems without a common letter differ significantly (P &lt; 0.05).</p></table-wrap-foot> <table-wrap-foot><p>For significance levels, systems without a common letter differ significantly (P &lt; 0.05).</p></table-wrap-foot>
<graphic xlink:href="1751-0147-50-29-1"/> <graphic xlink:href="1751-0147-50-29-2"/> <graphic xlink:href="1751-0147-50-29-3"/>
[{"collab": ["CAC (Codex Alimentarius Commission)"], "source": ["Principles and Guidelines for the Conduct of Microbiological Risk Assessment CAC/GL-30"], "year": ["1999"]}, {"collab": ["CAC (Codex Alimentarius Commission)"], "source": ["Principles and guidelines for the conduct of microbiological risk assessment Document CAC/GL 30"], "year": ["2002"]}, {"collab": ["NRC (National Research Council)"], "source": ["Risk assessment in the federal government: Managing the process"], "year": ["1983"], "publisher-name": ["Washington DC, US National Academies Press"]}, {"collab": ["OIE (Office International des Epizooties)"], "source": ["Handbook on Import Risk Analysis for Animals and Animal Products Introduction and qualitative risk analysis"], "year": ["2004"], "volume": ["1"], "publisher-name": ["Paris, OIE"]}, {"collab": ["OIE (Office International des Epizooties)"], "source": ["Handbook on Import Risk Analysis for Animals and Animal Products Quantitative risk assessment"], "year": ["2004"], "volume": ["2"], "publisher-name": ["Paris, OIE"]}, {"collab": ["EC (European Commission)"], "source": ["First report on the harmonisation of risk assessment procedures Part 1: The report of the Scientific Steering Committee's Working Group on Harmonisation of Risk Assessment Procedures in the Scientific Committees advising the European Commission in the area of human and environmental health"], "year": ["2000"], "publisher-name": ["Brussels European Commission"]}, {"collab": ["EFSA (European Food Safety Authority)"], "source": ["Principles of risk assessment of food producing animals: Current and future approaches"], "year": ["2005"], "publisher-name": ["Parma"]}, {"collab": ["EFSA (European Food Safety Authority)"], "article-title": ["Scientific opinion on the risks of poor welfare in intensive calf farming systems. An update of the Scientific Veterinary Committee Report on the Welfare of Calves"], "source": ["EFSA-Q-2005-014 The EFSA Journal"], "year": ["2006"], "volume": ["366"], "fpage": ["1"], "lpage": ["36"]}, {"collab": ["EFSA (European Food Safety Authority)"], "source": ["Scientific report on the risks of poor welfare in intensive calf farming systems An update of the Scientific Veterinary Committee Report on the Welfare of Calves EFSA-Q-2005-014 (144 pp), Annex to EFSA (2006a)"], "year": ["2006"]}, {"surname": ["Bracke", "Edwards", "Metz", "Noordhuizen", "Algers"], "given-names": ["MBM", "SA", "JHM", "JPTM", "B"], "article-title": ["Synthesis of semantic modeling and risk analysis methodology applied to animal welfare"], "source": ["Animal"], "year": ["2008"], "volume": ["2"], "fpage": ["1061"], "lpage": ["1072"], "pub-id": ["10.1017/S1751731108002139"]}, {"collab": ["Anonymous"], "article-title": ["Scientists' assessment of the impact of housing and management on animal welfare"], "source": ["J Appl Animal Welfare Sci"], "year": ["2001"], "volume": ["4"], "fpage": ["3"], "lpage": ["52"], "pub-id": ["10.1207/S15327604JAWS0401_2"]}, {"surname": ["Bracke"], "given-names": ["MBM"], "article-title": ["RICHPIG: a semantic model to assess enrichment materials for pigs"], "source": ["Animal Welfare"], "year": ["2008"], "volume": ["17"], "fpage": ["289"], "lpage": ["304"]}, {"surname": ["Bracke", "Hulsegge", "Keeling", "Blokhuis"], "given-names": ["MBM", "B", "L", "HJ"], "article-title": ["Decision support system with semantic model to assess the risk of tail biting in pigs: 1. Modelling"], "source": ["Appl Anim Behav Sci"], "year": ["2004"], "volume": ["87"], "fpage": ["31"], "lpage": ["44"], "pub-id": ["10.1016/j.applanim.2003.12.005"]}, {"surname": ["Bracke", "Zonderland", "Bleumer"], "given-names": ["MBM", "JJ", "EJB"], "article-title": ["Expert judgement on enrichment materials for pigs validates preliminary RICHPIG Model"], "source": ["Appl Anim Behav Sci"], "year": ["2007"], "volume": ["104"], "fpage": ["1"], "lpage": ["13"], "pub-id": ["10.1016/j.applanim.2006.05.005"]}, {"surname": ["Bracke", "Zonderland", "Bleumer"], "given-names": ["MBM", "JJ", "EJB"], "article-title": ["Expert consultation on weighting factors of criteria for assessing environmental enrichment materials for pigs"], "source": ["Appl Anim Behav Sci"], "year": ["2007"], "volume": ["104"], "fpage": ["14"], "lpage": ["23"], "pub-id": ["10.1016/j.applanim.2006.05.006"]}, {"collab": ["Anonymous"], "source": ["SPSS\u00ae Base 110 User's Guide"], "year": ["2001"], "publisher-name": ["SPSS, Chicago, Illinois, USA, SPSS"]}, {"surname": ["Searle", "Casella", "McCulloch"], "given-names": ["SR", "G", "EC"], "source": ["Variance components"], "year": ["1992"], "publisher-name": ["New Cork, UK, Wiley"]}, {"surname": ["Keen", "Engel"], "given-names": ["A", "B"], "article-title": ["Analysis of a mixed model for ordinal data by iterative re-weighted REML"], "source": ["Statistica Neerlandica"], "year": ["1997"], "volume": ["51"], "fpage": ["129"], "lpage": ["144"], "pub-id": ["10.1111/1467-9574.00044"]}, {"surname": ["Engel", "Buist"], "given-names": ["B", "WG"], "article-title": ["Analysis of a generalized linear model: a case study and simulation results"], "source": ["Biometrics J"], "year": ["1996"], "volume": ["38"], "fpage": ["61"], "lpage": ["80"], "pub-id": ["10.1002/bimj.4710380107"]}, {"collab": ["GenStat Committee"], "source": ["The guide to GenStat"], "year": ["2000"], "publisher-name": ["Oxford, UK, VSN Int"]}, {"surname": ["Bekoff", "Gruen", "Townsend", "Rollin"], "given-names": ["M", "L", "SE", "EB"], "article-title": ["Animals in science: Some areas revisited"], "source": ["Anim Behav"], "year": ["1992"], "volume": ["44"], "fpage": ["473"], "lpage": ["484"], "pub-id": ["10.1016/0003-3472(92)90057-G"]}, {"surname": ["Dawkins"], "given-names": ["MS"], "article-title": ["Battery hens name their price: Consumer demand theory and the measurement of ethological 'needs'"], "source": ["Anim Behav"], "year": ["1983"], "volume": ["31"], "fpage": ["1195"], "lpage": ["1205"], "pub-id": ["10.1016/S0003-3472(83)80026-8"]}, {"surname": ["Bracke", "Spruijt", "Metz"], "given-names": ["MBM", "BM", "JHM"], "article-title": ["Overall welfare assessment reviewed. Part 1: Is it possible?"], "source": ["Netherl J Agric Sci"], "year": ["1999"], "volume": ["47"], "fpage": ["279"], "lpage": ["291"]}]
Acta Vet Scand. 2008 Jul 14; 50(1):29