Model fixed effects | Plant | Zone and department | UJC | β Coefficient | SE |

Intercept | 2.636 | 0.632 | |||

Follow-up year (Year-1935) | −0.058 | 0.012 | |||

Spline 1 (Year-46, If >0)* | −0.016 | 0.016 | |||

Spline 2 (Year-64, If >0) | 0.029 | 0.012 | |||

Spline 3 (Year-71, If >0) | −0.210 | 0.025 | |||

Spline 4 (Year-76, If >0) | 0.092 | 0.028 | |||

UJCs nested within zones and plants | 1 | 1-Preparation | A | 1.281 | 0.675 |

1 | 1-Preparation | B | 1.659 | 0.650 | |

1 | 1-Preparation | D | 2.113 | 0.735 | |

1 | 2-Carding | A | 1.658 | 0.639 | |

1 | 2-Carding | B | 1.613 | 0.670 | |

1 | 2-Carding | D | 2.198 | 0.767 | |

1 | 3-Spinning | A | 0.603 | 0.640 | |

1 | 3-Spinning | B | 1.200 | 0.648 | |

1 | 4-Twisting | A | 0.758 | 0.646 | |

1 | 4-Twisting | B | 1.344 | 0.653 | |

1 | 5-Weaving | A | −0.329 | 0.651 | |

1 | 5-Weaving | B | −0.123 | 0.644 | |

3 | 1-Preparation | A | 1.241 | 0.634 | |

3 | 1-Preparation | B | 1.621 | 0.634 | |

3 | 1-Preparation | D | 1.371 | 0.650 | |

3 | 2-Carding | A | 1.595 | 0.632 | |

3 | 2-Carding | B | 1.916 | 0.631 | |

3 | 2-Carding | D | 1.610 | 0.659 | |

3 | 3-Spinning | A | 0.292 | 0.637 | |

3 | 3-Spinning | B | 0.920 | 0.631 | |

3 | 4-Twisting | A | 0.307 | 0.641 | |

3 | 4-Twisting | B | 1.553 | 0.632 | |

3 | 5-Weaving | A | 1.604 | 0.640 | |

3 | 5-Weaving | B | 1.033 | 0.629 | |

4 | 1-Preparation | A | 1.608 | 0.670 | |

4 | 1-Preparation | B | 1.683 | 0.639 | |

4 | 1-Preparation | D | 2.231 | 0.723 | |

4 | 2-Carding | A | 1.610 | 0.635 | |

4 | 2-Carding | B | 1.696 | 0.636 | |

4 | 2-Carding | D | 2.908 | 0.705 | |

4 | 3-Spinning | A | 1.913 | 0.649 | |

4 | 3-Spinning | B | 1.938 | 0.640 | |

4 | 4-Twisting | A | 2.279 | 0.661 | |

4 | 4-Twisting | B | 1.926 | 0.643 | |

4 | 5-Weaving | A | 1.279 | 0.643 | |

4 | 5-Weaving | B | 1.573 | 0.631 |

Example arithmetic mean PCM concentration calculation for a weaver (UJC B) in plant 1 in 1955. Model residual variance = 1.169.

AM = exp[(intercept+(β

_{year})×(1955–1935)+(β_{spline_1})×(1955–1946)+β_{UJC zone plant}+0.5×(residual model variance)]AM = exp[2.636+(−0.058×20)+(−0.016×9)+(−0.123)+(0.5×1.169)]

AM = 6.01 fibres>5 μm/ml

*In linear spline time trends the coefficients for the second term through remaining terms are constrained to apply only when the variable (year-boundary) is a positive value.

A, arithmetic mean; PCM, phase contract microscopy; SE, standard error; UJC, uniform job category.