n = 24 full Protestant editions; ten are KJV descendants. Fitted separately by testament, because the two have rates two orders of magnitude apart and pooling them lets the Old Testament dominate. ================ Old Testament ================ --- broad ~ year Call: lm(formula = broad ~ year, data = d) Residuals: Min 1Q Median 3Q Max -0.051957 -0.034216 -0.017802 0.007624 0.230103 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.4761150 0.2467336 1.930 0.0666 . year 0.0001542 0.0001272 1.212 0.2382 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.0616 on 22 degrees of freedom Multiple R-squared: 0.06263, Adjusted R-squared: 0.02002 F-statistic: 1.47 on 1 and 22 DF, p-value: 0.2382 --- broad ~ year + family fixed effects Call: lm(formula = broad ~ year + family, data = d) Residuals: Min 1Q Median 3Q Max -0.057796 -0.004336 0.000025 0.006425 0.073413 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.8090170 0.2079248 3.891 0.0013 ** year -0.0000265 0.0001075 -0.246 0.8085 familyLiteral 0.0364024 0.0172887 2.106 0.0514 . familyModern 0.0363935 0.0224742 1.619 0.1249 familyReformation -0.0959743 0.0487167 -1.970 0.0664 . familySimplified 0.2493197 0.0331458 7.522 1.22e-06 *** familyTyndale-KJV -0.0155079 0.0332544 -0.466 0.6473 familyWEB-family -0.0054189 0.0224198 -0.242 0.8121 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.03155 on 16 degrees of freedom Multiple R-squared: 0.8211, Adjusted R-squared: 0.7429 F-statistic: 10.49 on 7 and 16 DF, p-value: 6.043e-05 --- one observation per family, weighted Call: lm(formula = broad ~ year, data = byfam, weights = n) Weighted Residuals: 1 2 3 4 5 6 7 0.044110 -0.048672 -0.022195 0.228934 -0.003087 -0.016471 -0.074849 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.2731624 0.5653126 0.483 0.649 year 0.0002589 0.0002915 0.888 0.415 Residual standard error: 0.1123 on 5 degrees of freedom Multiple R-squared: 0.1363, Adjusted R-squared: -0.03646 F-statistic: 0.7889 on 1 and 5 DF, p-value: 0.4151 ================ Gospels ================ --- broad ~ year Call: lm(formula = broad ~ year, data = d) Residuals: Min 1Q Median 3Q Max -18.581 -1.140 1.132 4.337 10.001 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 17.57384 28.24639 0.622 0.5402 year 0.02557 0.01456 1.756 0.0929 . --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 7.052 on 22 degrees of freedom Multiple R-squared: 0.123, Adjusted R-squared: 0.08311 F-statistic: 3.085 on 1 and 22 DF, p-value: 0.09294 --- broad ~ year + family fixed effects Call: lm(formula = broad ~ year + family, data = d) Residuals: Min 1Q Median 3Q Max -16.3349 -0.8606 1.3411 2.8204 6.8692 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 82.515953 41.961701 1.966 0.0668 . year -0.007273 0.021700 -0.335 0.7419 familyLiteral -0.460544 3.489058 -0.132 0.8966 familyModern 4.593988 4.535563 1.013 0.3262 familyReformation -23.665422 9.831615 -2.407 0.0285 * familySimplified -0.267320 6.689208 -0.040 0.9686 familyTyndale-KJV -0.271421 6.711135 -0.040 0.9682 familyWEB-family -6.175051 4.524591 -1.365 0.1912 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 6.368 on 16 degrees of freedom Multiple R-squared: 0.4799, Adjusted R-squared: 0.2524 F-statistic: 2.109 on 7 and 16 DF, p-value: 0.1026 --- one observation per family, weighted Call: lm(formula = broad ~ year, data = byfam, weights = n) Weighted Residuals: 1 2 3 4 5 6 7 2.0944 5.1566 2.9415 0.4555 3.6079 -4.7852 -14.9249 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -19.31810 38.88158 -0.497 0.6404 year 0.04461 0.02005 2.225 0.0767 . --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 7.727 on 5 degrees of freedom Multiple R-squared: 0.4975, Adjusted R-squared: 0.397 F-statistic: 4.95 on 1 and 5 DF, p-value: 0.07665 ================ New Testament, rest ================ --- broad ~ year Call: lm(formula = broad ~ year, data = d) Residuals: Min 1Q Median 3Q Max -6.0069 -0.0888 0.5807 0.9716 1.5512 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 26.753539 6.964763 3.841 0.000888 *** year 0.003794 0.003589 1.057 0.301952 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 1.739 on 22 degrees of freedom Multiple R-squared: 0.04834, Adjusted R-squared: 0.005078 F-statistic: 1.117 on 1 and 22 DF, p-value: 0.302 --- broad ~ year + family fixed effects Call: lm(formula = broad ~ year + family, data = d) Residuals: Min 1Q Median 3Q Max -5.752 0.000 0.489 0.857 1.403 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 39.628431 12.166557 3.257 0.00495 ** year -0.002768 0.006292 -0.440 0.66591 familyLiteral -0.353239 1.011633 -0.349 0.73151 familyModern 0.379367 1.315061 0.288 0.77668 familyReformation -4.400626 2.850621 -1.544 0.14220 familySimplified 1.423442 1.939498 0.734 0.47361 familyTyndale-KJV 0.356018 1.945856 0.183 0.85713 familyWEB-family -0.189546 1.311879 -0.144 0.88692 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 1.846 on 16 degrees of freedom Multiple R-squared: 0.2196, Adjusted R-squared: -0.1218 F-statistic: 0.6433 on 7 and 16 DF, p-value: 0.7147 --- one observation per family, weighted Call: lm(formula = broad ~ year, data = byfam, weights = n) Weighted Residuals: 1 2 3 4 5 6 7 -0.3917 0.7107 0.8969 1.4668 -0.3870 -0.7293 -1.3485 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 19.382218 5.564687 3.483 0.0176 * year 0.007598 0.002869 2.648 0.0455 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 1.106 on 5 degrees of freedom Multiple R-squared: 0.5838, Adjusted R-squared: 0.5005 F-statistic: 7.013 on 1 and 5 DF, p-value: 0.04553