Linear models: Metal ~ Material
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 121.8 126.2 -56.9
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.0001886 2.871
##
## Fixed effects: X27Al ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 4.333 0.8288 17 5.229 1e-04
## materialplasma -6.000 1.1721 17 -5.119 1e-04
## Correlation:
## (Intr)
## materialplasma -0.707
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.1611 -0.4644 -0.2903 0.2322 2.6704
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## -0.3792 3.985 4.19
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.3501 0.119
##
## Fixed effects: log(X52Cr) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 0.3241 0.1470 17 2.205 0.0415
## materialplasma -0.3621 0.0486 17 -7.450 0.0000
## Correlation:
## (Intr)
## materialplasma -0.165
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.61486 -0.42512 0.09768 0.41627 2.12683
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 3.788 8.153 2.106
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.1439 0.1681
##
## Fixed effects: log(X55Mn) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 2.161 0.07620 17 28.36 0
## materialplasma -1.654 0.06863 17 -24.10 0
## Correlation:
## (Intr)
## materialplasma -0.45
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.26974 -0.61864 -0.07546 0.66801 1.57988
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 28.56 32.92 -10.28
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.601 0.2375
##
## Fixed effects: log(X59Co) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) -0.8574 0.25474 17 -3.366 0.0037
## materialplasma -0.4740 0.09697 17 -4.888 0.0001
## Correlation:
## (Intr)
## materialplasma -0.19
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -2.12132 -0.53606 0.09482 0.69190 1.21065
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 39.96 44.32 -15.98
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 1.402 0.26
##
## Fixed effects: log(X75As) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 0.3494 0.5772 17 0.605 0.5530
## materialplasma -0.4236 0.1062 17 -3.990 0.0009
## Correlation:
## (Intr)
## materialplasma -0.092
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.0555 -0.8593 0.1639 0.8013 0.9448
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## -55.24 -50.87 31.62
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.07716 0.03685
##
## Fixed effects: log(X78Se) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 4.647 0.03325 17 139.77 0
## materialplasma -0.095 0.01504 17 -6.34 0
## Correlation:
## (Intr)
## materialplasma -0.226
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.52265 -0.46025 -0.08334 0.33768 1.71503
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 12.03 16.4 -2.017
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.1948 0.1978
##
## Fixed effects: log(X95Mo) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) -0.4611 0.09790 17 -4.710 0.0002
## materialplasma 0.2647 0.08077 17 3.277 0.0044
## Correlation:
## (Intr)
## materialplasma -0.412
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.0809 -0.7181 -0.1076 0.3556 2.4225
##
## Number of Observations: 24
## Number of Groups: 6
## Error: NA/NaN/Inf in foreign function call (arg 1)
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 4.605 8.969 1.698
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.22 0.156
##
## Fixed effects: log(X133Cs) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 0.8588 0.10045 17 8.549 0
## materialplasma -0.6686 0.06367 17 -10.502 0
## Correlation:
## (Intr)
## materialplasma -0.317
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.6784 -0.3913 -0.1980 0.2845 2.0345
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 15.57 19.93 -3.785
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.1454 0.2303
##
## Fixed effects: log(X138Ba) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 0.3137 0.08915 17 3.519 0.0026
## materialplasma 0.3342 0.09404 17 3.554 0.0024
## Correlation:
## (Intr)
## materialplasma -0.527
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -1.44094 -0.71018 0.04561 0.36166 2.43940
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 47.83 52.19 -19.91
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.5085 0.4315
##
## Fixed effects: log(X202Hg) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) -0.341 0.2421 17 -1.408 0.177
## materialplasma -1.721 0.1762 17 -9.767 0.000
## Correlation:
## (Intr)
## materialplasma -0.364
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -2.7234 -0.4373 0.2155 0.5161 1.2640
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## -12.77 -8.406 10.39
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.1751 0.1006
##
## Fixed effects: log(X205Tl) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) -3.676 0.07715 17 -47.64 0
## materialplasma 0.483 0.04108 17 11.76 0
## Correlation:
## (Intr)
## materialplasma -0.266
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -2.07413 -0.59428 0.01372 0.27656 1.95539
##
## Number of Observations: 24
## Number of Groups: 6
## Linear mixed-effects model fit by REML
## Data: mydata
## AIC BIC logLik
## 47.53 51.89 -19.76
##
## Random effects:
## Formula: ~1 | ID
## (Intercept) Residual
## StdDev: 0.553 0.4192
##
## Fixed effects: log(X208Pb) ~ material
## Value Std.Error DF t-value p-value
## (Intercept) 2.556 0.2562 17 9.98 0
## materialplasma -2.383 0.1712 17 -13.92 0
## Correlation:
## (Intr)
## materialplasma -0.334
##
## Standardized Within-Group Residuals:
## Min Q1 Med Q3 Max
## -2.17132 -0.47422 -0.03951 0.42784 1.87022
##
## Number of Observations: 24
## Number of Groups: 6
Boxplots
T.Tests
##
## Welch Two Sample t-test
##
## data: X27Al by material
## t = 5.119, df = 17.81, p-value = 7.418e-05
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 3.536 8.464
## sample estimates:
## mean in group blood mean in group plasma
## 4.333 -1.667
##
## Welch Two Sample t-test
##
## data: X52Cr by material
## t = 1.586, df = 18.7, p-value = 0.1294
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.1123 0.8123
## sample estimates:
## mean in group blood mean in group plasma
## 1.425 1.075
##
## Welch Two Sample t-test
##
## data: X55Mn by material
## t = 9.585, df = 11.07, p-value = 1.073e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 5.676 9.057
## sample estimates:
## mean in group blood mean in group plasma
## 9.033 1.667
##
## Welch Two Sample t-test
##
## data: X59Co by material
## t = 1.262, df = 21.8, p-value = 0.2204
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.08059 0.33059
## sample estimates:
## mean in group blood mean in group plasma
## 0.4667 0.3417
##
## Welch Two Sample t-test
##
## data: X75As by material
## t = 0.3514, df = 21.9, p-value = 0.7287
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -1.798 2.531
## sample estimates:
## mean in group blood mean in group plasma
## 2.567 2.200
##
## Welch Two Sample t-test
##
## data: X78Se by material
## t = 2.861, df = 21, p-value = 0.009366
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 2.616 16.551
## sample estimates:
## mean in group blood mean in group plasma
## 104.67 95.08
##
## Welch Two Sample t-test
##
## data: X95Mo by material
## t = -2.273, df = 21.95, p-value = 0.03313
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.35059 -0.01608
## sample estimates:
## mean in group blood mean in group plasma
## 0.6583 0.8417
##
## Welch Two Sample t-test
##
## data: X111Cd by material
## t = 9.027, df = 11, p-value = 2.035e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.208 0.342
## sample estimates:
## mean in group blood mean in group plasma
## 0.275 0.000
##
## Welch Two Sample t-test
##
## data: X133Cs by material
## t = 5.012, df = 14.68, p-value = 0.0001653
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.6982 1.7351
## sample estimates:
## mean in group blood mean in group plasma
## 2.458 1.242
##
## Welch Two Sample t-test
##
## data: X138Ba by material
## t = -2.872, df = 16.12, p-value = 0.011
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -1.0426 -0.1574
## sample estimates:
## mean in group blood mean in group plasma
## 1.4 2.0
##
## Welch Two Sample t-test
##
## data: X202Hg by material
## t = 6.114, df = 11.85, p-value = 5.51e-05
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.4663 0.9837
## sample estimates:
## mean in group blood mean in group plasma
## 0.8667 0.1417
##
## Welch Two Sample t-test
##
## data: X205Tl by material
## t = -6.209, df = 19.95, p-value = 4.639e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.02115 -0.01051
## sample estimates:
## mean in group blood mean in group plasma
## 0.02583 0.04167
##
## Welch Two Sample t-test
##
## data: X208Pb by material
## t = 3.158, df = 11.01, p-value = 0.009095
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 5.508 30.825
## sample estimates:
## mean in group blood mean in group plasma
## 19.42 1.25
Paired T.Tests
##
## Paired t-test
##
## data: X27Al by material
## t = 4.648, df = 11, p-value = 0.0007079
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 3.159 8.841
## sample estimates:
## mean of the differences
## 6
##
## Paired t-test
##
## data: X52Cr by material
## t = 5.066, df = 11, p-value = 0.0003627
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.1979 0.5021
## sample estimates:
## mean of the differences
## 0.35
##
## Paired t-test
##
## data: X55Mn by material
## t = 9.644, df = 11, p-value = 1.061e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 5.685 9.048
## sample estimates:
## mean of the differences
## 7.367
##
## Paired t-test
##
## data: X59Co by material
## t = 9.574, df = 11, p-value = 1.14e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.09626 0.15374
## sample estimates:
## mean of the differences
## 0.125
##
## Paired t-test
##
## data: X75As by material
## t = 2.438, df = 11, p-value = 0.03296
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.03559 0.69774
## sample estimates:
## mean of the differences
## 0.3667
##
## Paired t-test
##
## data: X78Se by material
## t = 4.881, df = 11, p-value = 0.0004864
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 5.262 13.905
## sample estimates:
## mean of the differences
## 9.583
##
## Paired t-test
##
## data: X95Mo by material
## t = -2.682, df = 11, p-value = 0.02132
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.3338 -0.0329
## sample estimates:
## mean of the differences
## -0.1833
##
## Paired t-test
##
## data: X111Cd by material
## t = 9.027, df = 11, p-value = 2.035e-06
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.208 0.342
## sample estimates:
## mean of the differences
## 0.275
##
## Paired t-test
##
## data: X133Cs by material
## t = 5.832, df = 11, p-value = 0.0001138
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.7575 1.6758
## sample estimates:
## mean of the differences
## 1.217
##
## Paired t-test
##
## data: X138Ba by material
## t = -2.978, df = 11, p-value = 0.01257
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -1.0435 -0.1565
## sample estimates:
## mean of the differences
## -0.6
##
## Paired t-test
##
## data: X202Hg by material
## t = 6.984, df = 11, p-value = 2.318e-05
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 0.4965 0.9535
## sample estimates:
## mean of the differences
## 0.725
##
## Paired t-test
##
## data: X205Tl by material
## t = -10.65, df = 11, p-value = 3.921e-07
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.01911 -0.01256
## sample estimates:
## mean of the differences
## -0.01583
##
## Paired t-test
##
## data: X208Pb by material
## t = 3.2, df = 11, p-value = 0.008463
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 5.669 30.664
## sample estimates:
## mean of the differences
## 18.17