library("haven")
library("tidyverse")
library("dplyr")
library("foreign")
library("survey")
library("ggplot2")
library("car")
library("rms")
library("SciViews")#list variable
colnames(Fulldat_Pheno)#Perfluorohexane_sulfonic_acid_comment #Perfluorononanoic_acid_comment #perfluorooctanoic_acid_comment #perfluorooctane_sulfonic_acid_comment
#check effect modifiers Subgroup analysis for gender, race, BMI, income, smoking and cancer #for subgroup analysis by gender
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * Gender, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.5952 0.1960 8.141 4.24e-16 ***
Perfluorohexane_sulfonic_acid_comment 6.1603 1.6604 3.710 0.000208 ***
Gender -1.3197 0.1229 -10.737 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment:Gender 1.1009 0.9862 1.116 0.264266
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.18197)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 820059 on 14861 degrees of freedom
AIC: 101809
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desmen, family = "gaussian", data = men)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = men)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.3721 0.1071 -3.475 0.000691 ***
Perfluorohexane_sulfonic_acid_comment 7.4917 1.8308 4.092 7.37e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.27316)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.870465 11.112951
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = deswomen, family = "gaussian", data = women)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = women)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.671 0.112 -14.927 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 9.278 1.655 5.605 1.15e-07 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 50.39522)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
6.003902 12.552164
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * Gender, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.6089 0.1966 8.184 2.96e-16 ***
Perfluorononanoic_acid_comment 4.9167 1.6769 2.932 0.00337 **
Gender -1.3125 0.1233 -10.642 < 2e-16 ***
Perfluorononanoic_acid_comment:Gender 0.8658 0.9886 0.876 0.38113
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.54255)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 825418 on 14861 degrees of freedom
AIC: 101906
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desmen, family = "gaussian", data = men)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = men)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.3858 0.1092 -3.534 0.000564 ***
Perfluorononanoic_acid_comment 6.9288 2.0675 3.351 0.001048 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.22506)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.839324 11.018355
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = deswomen, family = "gaussian", data = women)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = women)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.6684 0.1152 -14.483 < 2e-16 ***
Perfluorononanoic_acid_comment 6.4530 1.7351 3.719 0.000294 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 50.76961)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.020963 9.885070
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * Gender, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.44030 0.24219 1.818 0.0691 .
perfluorooctanoic_acid_comment 3.33485 0.42925 7.769 8.5e-15 ***
Gender -1.40799 0.15232 -9.244 < 2e-16 ***
perfluorooctanoic_acid_comment:Gender 0.08009 0.26876 0.298 0.7657
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.17327)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 665457 on 13004 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88110
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desmen, family = "gaussian", data = men)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = men)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.7337 0.1155 -15.01 <2e-16 ***
perfluorooctanoic_acid_comment 3.3703 0.2809 12.00 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 37.60543)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.814148 3.926515
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = deswomen, family = "gaussian", data = women)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = women)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.1552 0.1243 -25.38 <2e-16 ***
perfluorooctanoic_acid_comment 3.5584 0.2362 15.06 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.44878)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.090583 4.026248
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * Gender, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.4522 0.2039 7.122 1.12e-12 ***
perfluorooctane_sulfonic_acid_comment 5.5355 2.8571 1.937 0.0527 .
Gender -1.3870 0.1280 -10.837 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment:Gender 4.2539 1.7364 2.450 0.0143 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.93476)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 688364 on 13004 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88551
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desmen, family = "gaussian", data = men)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = men)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6023 0.1154 -5.220 7.77e-07 ***
perfluorooctane_sulfonic_acid_comment 9.0445 2.0534 4.405 2.34e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.85434)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.978269 13.110728
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = deswomen, family = "gaussian", data = women)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = women)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.9927 0.1154 -17.268 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 13.5428 3.1089 4.356 2.84e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 47.34834)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
7.38634 19.69926 #for subgroup analysis by Race
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * as.factor(Race),
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7115 0.1515 -4.698 2.65e-06 ***
Perfluorohexane_sulfonic_acid_comment 4.9492 0.9304 5.320 1.06e-07 ***
as.factor(Race)2 0.1864 0.2546 0.732 0.46417
as.factor(Race)3 -0.2032 0.1783 -1.139 0.25455
as.factor(Race)4 1.9672 0.2004 9.816 < 2e-16 ***
as.factor(Race)5 -0.4655 0.2409 -1.932 0.05339 .
Perfluorohexane_sulfonic_acid_comment:as.factor(Race)2 5.4430 1.8136 3.001 0.00269 **
Perfluorohexane_sulfonic_acid_comment:as.factor(Race)3 3.2250 1.3206 2.442 0.01461 *
Perfluorohexane_sulfonic_acid_comment:as.factor(Race)4 3.5303 1.2820 2.754 0.00590 **
Perfluorohexane_sulfonic_acid_comment:as.factor(Race)5 3.1356 1.9077 1.644 0.10028
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 54.77517)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 813685 on 14855 degrees of freedom
AIC: 101705
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_Me, family = "gaussian", data = Mexican)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Mexican)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.5365 0.2011 -2.667 0.0090 **
Perfluorohexane_sulfonic_acid_comment 4.6642 2.5276 1.845 0.0681 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.2044)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-0.3543727 9.6828627
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_Hispanic, family = "gaussian", data = Other_Hispanic)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Hispanic)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6022 0.2706 -2.225 0.02838 *
Perfluorohexane_sulfonic_acid_comment 10.4945 3.6972 2.838 0.00552 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.33981)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.156444 17.832475
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_white, family = "gaussian", data = Non_Hispanic_white)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_white)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.4415 0.1229 -11.724 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 7.5792 1.8626 4.069 8.17e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.34425)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.894114 11.264323
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_Black, family = "gaussian", data = Non_Hispanic_Black)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_Black)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.1240 0.1785 6.297 5.62e-09 ***
Perfluorohexane_sulfonic_acid_comment 10.1317 2.4251 4.178 5.73e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 66.70799)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.328519 14.934890
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_Other, family = "gaussian", data = Other_Race)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Race)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.4384 0.2539 -5.665 1.16e-07 ***
Perfluorohexane_sulfonic_acid_comment 9.9477 3.9640 2.510 0.0135 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.115)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.093614 17.801795
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * as.factor(Race), data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6794 0.1518 -4.477 7.64e-06 ***
Perfluorononanoic_acid_comment 3.8545 0.9470 4.070 4.72e-05 ***
as.factor(Race)2 0.1909 0.2548 0.749 0.453695
as.factor(Race)3 -0.2459 0.1788 -1.375 0.169025
as.factor(Race)4 1.9910 0.2006 9.923 < 2e-16 ***
as.factor(Race)5 -0.4431 0.2419 -1.832 0.066956 .
Perfluorononanoic_acid_comment:as.factor(Race)2 7.3313 2.0437 3.587 0.000335 ***
Perfluorononanoic_acid_comment:as.factor(Race)3 3.1463 1.2526 2.512 0.012022 *
Perfluorononanoic_acid_comment:as.factor(Race)4 3.7969 1.3773 2.757 0.005846 **
Perfluorononanoic_acid_comment:as.factor(Race)5 -1.1542 1.7024 -0.678 0.497803
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.04869)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 817748 on 14855 degrees of freedom
AIC: 101779
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_Me, family = "gaussian", data = Mexican)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Mexican)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.5032 0.2005 -2.510 0.0138 *
Perfluorononanoic_acid_comment 2.3205 1.3354 1.738 0.0855 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.5349)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-0.3309801 4.9719834
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_Hispanic, family = "gaussian", data = Other_Hispanic)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Hispanic)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.5743 0.2737 -2.098 0.0385 *
Perfluorononanoic_acid_comment 9.3529 3.9552 2.365 0.0200 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.82539)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.502917 17.202791
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_white, family = "gaussian", data = Non_Hispanic_white)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_white)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.4775 0.1237 -11.95 < 2e-16 ***
Perfluorononanoic_acid_comment 7.7234 2.2192 3.48 0.000684 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.06343)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.332568 12.114179
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_Black, family = "gaussian", data = Non_Hispanic_Black)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_Black)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.2074 0.1727 6.991 1.85e-10 ***
Perfluorononanoic_acid_comment 7.2312 2.0300 3.562 0.000535 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 67.97717)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.210544 11.251917
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_Other, family = "gaussian", data = Other_Race)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Race)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.3690 0.2628 -5.209 8.74e-07 ***
Perfluorononanoic_acid_comment 2.1736 1.3238 1.642 0.103
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 53.12703)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-0.4493478 4.7965911
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * as.factor(Race), data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.92370 0.17904 -10.744 < 2e-16 ***
perfluorooctanoic_acid_comment 3.56561 0.33891 10.521 < 2e-16 ***
as.factor(Race)2 -0.21326 0.32388 -0.658 0.5103
as.factor(Race)3 -0.09587 0.21103 -0.454 0.6496
as.factor(Race)4 1.73436 0.24456 7.092 1.39e-12 ***
as.factor(Race)5 -0.83532 0.33485 -2.495 0.0126 *
perfluorooctanoic_acid_comment:as.factor(Race)2 0.41861 0.54610 0.767 0.4434
perfluorooctanoic_acid_comment:as.factor(Race)3 -0.20350 0.40480 -0.503 0.6152
perfluorooctanoic_acid_comment:as.factor(Race)4 0.08332 0.43750 0.190 0.8490
perfluorooctanoic_acid_comment:as.factor(Race)5 -0.78879 0.51996 -1.517 0.1293
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 50.84169)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 660840 on 12998 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88032
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_Me, family = "gaussian", data = Mexican)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Mexican)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.1765 0.2389 -9.110 2.94e-14 ***
perfluorooctanoic_acid_comment 3.7655 0.4314 8.729 1.75e-13 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.81887)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.907902 4.623026
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_Hispanic, family = "gaussian", data = Other_Hispanic)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Hispanic)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.2713 0.4004 -5.673 1.94e-07 ***
perfluorooctanoic_acid_comment 3.7574 0.5971 6.293 1.35e-08 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 48.02612)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.570036 4.944711
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_white, family = "gaussian", data = Non_Hispanic_white)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_white)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.7940 0.1080 -25.87 <2e-16 ***
perfluorooctanoic_acid_comment 3.3689 0.2895 11.64 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 35.99008)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.795409 3.942384
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_Black, family = "gaussian", data = Non_Hispanic_Black)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_Black)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.2256 0.1727 -1.306 0.194
perfluorooctanoic_acid_comment 3.4169 0.3153 10.837 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 63.31689)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.791492 4.042238
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_Other, family = "gaussian", data = Other_Race)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Race)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.0761 0.4160 -7.394 5.08e-11 ***
perfluorooctanoic_acid_comment 3.3334 0.5165 6.453 4.31e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 50.67045)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.308229 4.358550
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * as.factor(Race),
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.9815 0.1550 -6.332 2.50e-10 ***
perfluorooctane_sulfonic_acid_comment 10.5764 2.1921 4.825 1.42e-06 ***
as.factor(Race)2 0.3552 0.2638 1.346 0.1783
as.factor(Race)3 -0.2276 0.1835 -1.240 0.2150
as.factor(Race)4 2.0094 0.2062 9.746 < 2e-16 ***
as.factor(Race)5 -0.4196 0.2539 -1.652 0.0985 .
perfluorooctane_sulfonic_acid_comment:as.factor(Race)2 10.2731 4.2426 2.421 0.0155 *
perfluorooctane_sulfonic_acid_comment:as.factor(Race)3 2.5188 2.7538 0.915 0.3604
perfluorooctane_sulfonic_acid_comment:as.factor(Race)4 -0.1072 2.5763 -0.042 0.9668
perfluorooctane_sulfonic_acid_comment:as.factor(Race)5 -0.1326 3.2695 -0.041 0.9676
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.59557)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 683637 on 12998 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88473
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_Me, family = "gaussian", data = Mexican)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Mexican)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.8310 0.2261 -3.675 0.000412 ***
perfluorooctane_sulfonic_acid_comment 5.8035 1.9844 2.925 0.004410 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 47.83446)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.858603 9.748373
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_Hispanic, family = "gaussian", data = Other_Hispanic)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Hispanic)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6994 0.2940 -2.379 0.0196 *
perfluorooctane_sulfonic_acid_comment 20.3798 7.8008 2.613 0.0106 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 49.15044)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.867019 35.892669
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_white, family = "gaussian", data = Non_Hispanic_white)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_white)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.7480 0.1299 -13.457 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 10.6027 2.9420 3.604 0.000466 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.0873)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.774701 16.430738
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_Black, family = "gaussian", data = Non_Hispanic_Black)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Non_Hispanic_Black)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.9148 0.1857 4.925 3.26e-06 ***
perfluorooctane_sulfonic_acid_comment 10.1224 3.8238 2.647 0.0094 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 64.86676)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.537853 17.706861
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_Other, family = "gaussian", data = Other_Race)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Other_Race)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.6442 0.2717 -6.051 2.7e-08 ***
perfluorooctane_sulfonic_acid_comment 18.3124 8.3141 2.203 0.03 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.09565)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.811096 34.813654 ##for subgroup analysis by BMI
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * BMI_cat, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.9372 0.1071 -18.082 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 8.4423 0.7958 10.609 < 2e-16 ***
BMI_catobesity 3.3788 0.1469 23.003 < 2e-16 ***
BMI_catoverweight 0.9736 0.1514 6.429 1.33e-10 ***
Perfluorohexane_sulfonic_acid_comment:BMI_catobesity 0.4993 1.1155 0.448 0.6544
Perfluorohexane_sulfonic_acid_comment:BMI_catoverweight -2.8420 1.1963 -2.376 0.0175 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 53.47542)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 794591 on 14859 degrees of freedom
AIC: 101344
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desBMI_1, family = "gaussian", data = BMI_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.6351 0.1391 -18.946 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 10.5537 2.5787 4.093 7.35e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.75155)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.453193 15.654257
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desBMI_2, family = "gaussian", data = BMI_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.5966 0.1158 -13.786 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 6.2200 1.9369 3.211 0.00166 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.83277)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.388968 10.051035
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desBMI_3, family = "gaussian", data = BMI_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.9077 0.1447 6.274 4.58e-09 ***
Perfluorohexane_sulfonic_acid_comment 8.2549 1.7589 4.693 6.60e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.15996)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.775821 11.734040
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * BMI_cat, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.9302 0.1074 -17.978 < 2e-16 ***
Perfluorononanoic_acid_comment 9.1166 0.8483 10.746 < 2e-16 ***
BMI_catobesity 3.4390 0.1473 23.343 < 2e-16 ***
BMI_catoverweight 0.9673 0.1519 6.368 1.97e-10 ***
Perfluorononanoic_acid_comment:BMI_catobesity -4.6463 1.1334 -4.099 4.17e-05 ***
Perfluorononanoic_acid_comment:BMI_catoverweight -4.0169 1.2077 -3.326 0.000883 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 53.81901)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 799697 on 14859 degrees of freedom
AIC: 101439
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desBMI_1, family = "gaussian", data = BMI_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.680 0.134 -20.006 < 2e-16 ***
Perfluorononanoic_acid_comment 11.039 2.671 4.133 6.29e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.21026)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.756221 16.320828
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desBMI_2, family = "gaussian", data = BMI_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.6002 0.1206 -13.274 <2e-16 ***
Perfluorononanoic_acid_comment 5.3868 2.3739 2.269 0.0249 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.8722)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
0.6914221 10.0822367
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desBMI_3, family = "gaussian", data = BMI_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.9459 0.1476 6.407 2.37e-09 ***
Perfluorononanoic_acid_comment 3.6160 1.5055 2.402 0.0177 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.80188)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
0.6382136 6.5937903
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * BMI_cat, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.9651 0.1307 -22.679 < 2e-16 ***
perfluorooctanoic_acid_comment 3.1587 0.2401 13.156 < 2e-16 ***
BMI_catobesity 2.9621 0.1836 16.131 < 2e-16 ***
BMI_catoverweight 0.8305 0.1850 4.490 7.17e-06 ***
perfluorooctanoic_acid_comment:BMI_catobesity 0.4028 0.3213 1.254 0.210
perfluorooctanoic_acid_comment:BMI_catoverweight -0.2321 0.3395 -0.684 0.494
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 49.8447)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 648081 on 13002 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 87770
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desBMI_1, family = "gaussian", data = BMI_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.8210 0.1631 -23.42 <2e-16 ***
perfluorooctanoic_acid_comment 3.3857 0.3378 10.02 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.66721)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.716865 4.054572
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desBMI_2, family = "gaussian", data = BMI_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.8212 0.1401 -20.14 <2e-16 ***
perfluorooctanoic_acid_comment 2.9324 0.2647 11.08 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 34.08267)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.408208 3.456646
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desBMI_3, family = "gaussian", data = BMI_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6855 0.1522 -4.505 1.57e-05 ***
perfluorooctanoic_acid_comment 3.4187 0.3102 11.021 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 47.44875)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.804389 4.032947
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * BMI_cat, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.1173 0.1114 -19.002 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 16.7276 1.5294 10.937 < 2e-16 ***
BMI_catobesity 3.3194 0.1527 21.731 < 2e-16 ***
BMI_catoverweight 0.8415 0.1576 5.341 9.42e-08 ***
perfluorooctane_sulfonic_acid_comment:BMI_catobesity -2.7488 1.9874 -1.383 0.167
perfluorooctane_sulfonic_acid_comment:BMI_catoverweight -14.3483 2.2793 -6.295 3.17e-10 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.18889)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 665558 on 13002 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88116
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desBMI_1, family = "gaussian", data = BMI_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.8288 0.1498 -18.883 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 17.2257 3.4739 4.959 2.4e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.98791)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
10.34654 24.10491
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desBMI_2, family = "gaussian", data = BMI_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.8699 0.1259 -14.847 <2e-16 ***
perfluorooctane_sulfonic_acid_comment 2.7231 1.4013 1.943 0.0544 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 35.95703)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-0.0518824 5.4980597
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desBMI_3, family = "gaussian", data = BMI_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = BMI_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.5861 0.1534 3.821 0.000213 ***
perfluorooctane_sulfonic_acid_comment 13.2282 4.0348 3.279 0.001371 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 49.19045)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.238259 21.218215 #for subgroup analysis by income
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * as.factor(income_cat),
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.18813 0.08854 -2.125 0.033618 *
Perfluorohexane_sulfonic_acid_comment 7.26575 0.67448 10.772 < 2e-16 ***
as.factor(income_cat)poor 0.53733 0.14424 3.725 0.000196 ***
as.factor(income_cat)rich -1.60891 0.15549 -10.348 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment:as.factor(income_cat)poor -0.07016 1.01161 -0.069 0.944710
Perfluorohexane_sulfonic_acid_comment:as.factor(income_cat)rich 3.16739 1.75781 1.802 0.071581 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 54.99103)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 817112 on 14859 degrees of freedom
AIC: 101760
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desincome_1, family = "gaussian", data = income_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.1731 0.1886 0.918 0.360349
Perfluorohexane_sulfonic_acid_comment 7.6148 2.1202 3.592 0.000461 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 60.64966)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.421199 11.808492
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desincome_2, family = "gaussian", data = income_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7170 0.1034 -6.933 1.62e-10 ***
Perfluorohexane_sulfonic_acid_comment 7.7398 1.5914 4.864 3.20e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 48.43931)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.592036 10.887476
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desincome_3, family = "gaussian", data = income_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.2168 0.1361 -16.288 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 10.0628 2.6317 3.824 0.000201 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 32.43671)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.857289 15.268237
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * as.factor(income_cat), data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.14671 0.08883 -1.652 0.098645 .
Perfluorononanoic_acid_comment 4.53065 0.65323 6.936 4.21e-12 ***
as.factor(income_cat)poor 0.48775 0.14456 3.374 0.000743 ***
as.factor(income_cat)rich -1.62086 0.15600 -10.390 < 2e-16 ***
Perfluorononanoic_acid_comment:as.factor(income_cat)poor 3.77794 1.02664 3.680 0.000234 ***
Perfluorononanoic_acid_comment:as.factor(income_cat)rich 0.24162 1.62927 0.148 0.882107
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.28401)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 821465 on 14859 degrees of freedom
AIC: 101839
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desincome_1, family = "gaussian", data = income_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.0537 0.1732 0.310 0.757025
Perfluorononanoic_acid_comment 10.3622 2.6900 3.852 0.000181 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 59.02526)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.041388 15.683037
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desincome_2, family = "gaussian", data = income_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6808 0.1073 -6.343 3.27e-09 ***
Perfluorononanoic_acid_comment 3.7769 1.2646 2.987 0.00336 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 48.9918)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.275669 6.278137
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desincome_3, family = "gaussian", data = income_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.2010 0.1368 -16.094 <2e-16 ***
Perfluorononanoic_acid_comment 4.2937 1.8683 2.298 0.0231 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 32.76351)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
0.5982343 7.9892044
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * as.factor(income_cat), data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.36874 0.10863 -12.600 <2e-16 ***
perfluorooctanoic_acid_comment 3.21697 0.19472 16.521 <2e-16 ***
as.factor(income_cat)poor 0.25458 0.17963 1.417 0.1564
as.factor(income_cat)rich -1.69968 0.19261 -8.824 <2e-16 ***
perfluorooctanoic_acid_comment:as.factor(income_cat)poor 0.68111 0.31318 2.175 0.0297 *
perfluorooctanoic_acid_comment:as.factor(income_cat)rich 0.02193 0.34220 0.064 0.9489
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 50.97822)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 662819 on 13002 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88063
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desincome_1, family = "gaussian", data = income_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.5909 0.1669 -9.531 2.59e-16 ***
perfluorooctanoic_acid_comment 4.2412 0.3792 11.184 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.0874)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.490244 4.992167
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desincome_2, family = "gaussian", data = income_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.9587 0.1236 -15.85 <2e-16 ***
perfluorooctanoic_acid_comment 3.1547 0.2507 12.59 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 46.06662)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.658291 3.651106
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desincome_3, family = "gaussian", data = income_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.6876 0.1228 -30.02 <2e-16 ***
perfluorooctanoic_acid_comment 3.4211 0.2552 13.41 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 28.47216)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.915777 3.926497
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * as.factor(income_cat),
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.44758 0.09195 -4.868 1.14e-06 ***
perfluorooctane_sulfonic_acid_comment 13.21151 1.18128 11.184 < 2e-16 ***
as.factor(income_cat)poor 0.57565 0.15001 3.838 0.000125 ***
as.factor(income_cat)rich -1.60155 0.16224 -9.872 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment:as.factor(income_cat)poor 0.51567 1.92380 0.268 0.788666
perfluorooctane_sulfonic_acid_comment:as.factor(income_cat)rich -9.03725 2.49091 -3.628 0.000287 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.70492)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 685269 on 13002 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88496
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desincome_1, family = "gaussian", data = income_1)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_1)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.09178 0.19781 -0.464 0.64352
perfluorooctane_sulfonic_acid_comment 14.19665 4.91564 2.888 0.00461 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 54.26434)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.462334 23.930958
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desincome_2, family = "gaussian", data = income_2)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_2)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.0058 0.1139 -8.829 1.16e-14 ***
perfluorooctane_sulfonic_acid_comment 13.1078 3.0187 4.342 3.00e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 47.19879)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
7.129859 19.085693
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desincome_3, family = "gaussian", data = income_3)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = income_3)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.4692 0.1323 -18.671 <2e-16 ***
perfluorooctane_sulfonic_acid_comment 2.3929 1.8054 1.325 0.188
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 31.18655)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-1.182281 5.968159 #for subgroup analysis by cancer
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * had_cancer, data = Fulldat_Pheno)
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.8760 0.2063 4.247 2.18e-05 ***
Perfluorohexane_sulfonic_acid_comment 17.9480 1.9284 9.307 < 2e-16 ***
had_cancer2 -1.3311 0.2168 -6.141 8.42e-10 ***
had_cancer9 6.9978 2.6335 2.657 0.00789 **
had_cancerNone -2.4697 0.3222 -7.665 1.90e-14 ***
Perfluorohexane_sulfonic_acid_comment:had_cancer2 -10.2671 1.9955 -5.145 2.71e-07 ***
Perfluorohexane_sulfonic_acid_comment:had_cancer9 NA NA NA NA
Perfluorohexane_sulfonic_acid_comment:had_cancerNone -17.8615 2.7306 -6.541 6.30e-11 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.14143)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 819291 on 14858 degrees of freedom
AIC: 101801
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_cancer, family = "gaussian", data = cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.05352 0.26966 0.198 0.843
Perfluorohexane_sulfonic_acid_comment 22.74850 3.94539 5.766 6e-08 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.74015)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
14.94007 30.55692
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.14044 0.08711 -13.09 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 7.84032 1.23862 6.33 3.48e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.69392)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
5.390383 10.290249
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * had_cancer, data = Fulldat_Pheno)
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.0110 0.2076 4.870 1.13e-06 ***
Perfluorononanoic_acid_comment 4.6107 1.6808 2.743 0.00609 **
had_cancer2 -1.4492 0.2181 -6.644 3.16e-11 ***
had_cancer9 6.8628 2.6454 2.594 0.00949 **
had_cancerNone -2.6160 0.3237 -8.082 6.89e-16 ***
Perfluorononanoic_acid_comment:had_cancer2 2.0383 1.7577 1.160 0.24621
Perfluorononanoic_acid_comment:had_cancer9 NA NA NA NA
Perfluorononanoic_acid_comment:had_cancerNone -3.7188 2.6771 -1.389 0.16481
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.6396)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 826693 on 14858 degrees of freedom
AIC: 101935
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_cancer, family = "gaussian", data = cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.1764 0.2777 0.635 0.527
Perfluorononanoic_acid_comment 3.2018 2.2714 1.410 0.161
(Dispersion parameter for gaussian family taken to be 59.31664)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-1.293469 7.697137
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.15765 0.09109 -12.709 < 2e-16 ***
Perfluorononanoic_acid_comment 6.98972 1.64974 4.237 4.21e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.66237)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.726592 10.252840
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * had_cancer, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.4511 0.2591 -1.741 0.081669 .
perfluorooctanoic_acid_comment 4.2892 0.4503 9.525 < 2e-16 ***
had_cancer2 -1.2946 0.2721 -4.759 1.97e-06 ***
had_cancer9 -0.4836 3.2128 -0.151 0.880347
had_cancerNone -2.0274 0.3865 -5.245 1.59e-07 ***
perfluorooctanoic_acid_comment:had_cancer2 -0.9392 0.4731 -1.985 0.047139 *
perfluorooctanoic_acid_comment:had_cancer9 19.2001 5.2488 3.658 0.000255 ***
perfluorooctanoic_acid_comment:had_cancerNone -1.6807 0.7498 -2.241 0.025015 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.27506)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 666576 on 13000 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88140
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_cancer, family = "gaussian", data = cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.3175 0.3620 -3.640 0.000416 ***
perfluorooctanoic_acid_comment 3.8986 0.5411 7.204 7.45e-11 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 57.28193)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.826291 4.970889
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.56511 0.09009 -28.47 <2e-16 ***
perfluorooctanoic_acid_comment 3.40019 0.20672 16.45 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.88168)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.990827 3.809553
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * had_cancer, data = Fulldat_Pheno)
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.9100 0.2160 4.212 2.55e-05 ***
perfluorooctane_sulfonic_acid_comment 13.3994 3.2649 4.104 4.08e-05 ***
had_cancer2 -1.6323 0.2269 -7.192 6.71e-13 ***
had_cancer9 6.9638 2.5845 2.694 0.00706 **
had_cancerNone -2.7921 0.3353 -8.327 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment:had_cancer2 -1.2122 3.3865 -0.358 0.72039
perfluorooctane_sulfonic_acid_comment:had_cancer9 NA NA NA NA
perfluorooctane_sulfonic_acid_comment:had_cancerNone -12.6722 7.9869 -1.587 0.11262
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 53.06482)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 689896 on 13001 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88586
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_cancer, family = "gaussian", data = cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.1067 0.3097 0.344 0.731
perfluorooctane_sulfonic_acid_comment 9.8745 8.6906 1.136 0.258
(Dispersion parameter for gaussian family taken to be 60.36244)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-7.346495 27.095469
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.44500 0.09115 -15.853 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 11.99581 2.33803 5.131 1.15e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.82679)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
7.365872 16.625754 #for subgroup analysis by smoking
Call:
glm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Perfluorohexane_sulfonic_acid_comment * now_smoke, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.82409 0.07648 -10.775 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 6.61939 0.59571 11.112 < 2e-16 ***
now_smoke 0.46001 0.04970 9.256 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment:now_smoke 1.36515 0.39941 3.418 0.000633 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.21313)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 820522 on 14861 degrees of freedom
AIC: 101818
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = descurrent_smokers, family = "gaussian", data = current_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = current_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.4881 0.1640 2.977 0.00347 **
Perfluorohexane_sulfonic_acid_comment 9.0906 4.0359 2.252 0.02597 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.86239)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.106084 17.075116
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desformer_smokers, family = "gaussian", data = former_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = former_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7992 0.1597 -5.005 1.74e-06 ***
Perfluorohexane_sulfonic_acid_comment 9.4329 2.5838 3.651 0.000375 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.0045)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.322236 14.543536
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = desnon_smokers, family = "gaussian", data = non_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.5451 0.1086 -14.228 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 7.8065 1.4971 5.214 6.88e-07 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.58836)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
4.845293 10.767756
Call:
glm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Perfluorononanoic_acid_comment * now_smoke, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.79949 0.07674 -10.418 < 2e-16 ***
Perfluorononanoic_acid_comment 5.10597 0.59562 8.573 < 2e-16 ***
now_smoke 0.46052 0.04987 9.234 < 2e-16 ***
Perfluorononanoic_acid_comment:now_smoke 1.20718 0.39648 3.045 0.00233 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 55.58049)
Null deviance: 840986 on 14864 degrees of freedom
Residual deviance: 825982 on 14861 degrees of freedom
AIC: 101916
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = descurrent_smokers, family = "gaussian", data = current_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = current_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.5144 0.1617 3.181 0.00183 **
Perfluorononanoic_acid_comment 5.4740 3.4318 1.595 0.11313
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.42867)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-1.315472 12.263425
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desformer_smokers, family = "gaussian", data = former_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = former_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7605 0.1604 -4.740 5.41e-06 ***
Perfluorononanoic_acid_comment 5.4194 2.1733 2.494 0.0139 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 52.66954)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
1.120809 9.718025
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = desnon_smokers, family = "gaussian", data = non_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.574 0.111 -14.173 < 2e-16 ***
Perfluorononanoic_acid_comment 7.267 2.206 3.295 0.00126 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.43543)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.904712 11.629794
Call:
glm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
perfluorooctanoic_acid_comment * now_smoke, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.11342 0.09537 -22.161 < 2e-16 ***
perfluorooctanoic_acid_comment 3.27963 0.16651 19.697 < 2e-16 ***
now_smoke 0.45774 0.06129 7.468 8.64e-14 ***
perfluorooctanoic_acid_comment:now_smoke 0.20797 0.10903 1.907 0.0565 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.21673)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 666022 on 13004 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88121
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = descurrent_smokers, family = "gaussian", data = current_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = current_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7812 0.1588 -4.920 2.91e-06 ***
perfluorooctanoic_acid_comment 3.8372 0.3970 9.666 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.06004)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.050853 4.623523
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desformer_smokers, family = "gaussian", data = former_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = former_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.1898 0.1768 -12.39 <2e-16 ***
perfluorooctanoic_acid_comment 3.6643 0.3610 10.15 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 49.5252)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.949335 4.379225
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = desnon_smokers, family = "gaussian", data = non_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.1600 0.1205 -26.22 <2e-16 ***
perfluorooctanoic_acid_comment 3.5151 0.2387 14.73 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.43287)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
3.042374 3.987740
Call:
glm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
perfluorooctane_sulfonic_acid_comment * now_smoke, data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.09641 0.07976 -13.746 <2e-16 ***
perfluorooctane_sulfonic_acid_comment 11.55955 1.13208 10.211 <2e-16 ***
now_smoke 0.49016 0.05172 9.476 <2e-16 ***
perfluorooctane_sulfonic_acid_comment:now_smoke 0.47998 0.69708 0.689 0.491
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 53.04973)
Null deviance: 705266 on 13007 degrees of freedom
Residual deviance: 689859 on 13004 degrees of freedom
(因為不存在，1857 個觀察量被刪除了)
AIC: 88579
Number of Fisher Scoring iterations: 2
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = descurrent_smokers, family = "gaussian", data = current_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = current_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.2582 0.1687 1.531 0.1286
perfluorooctane_sulfonic_acid_comment 6.4072 3.3876 1.891 0.0611 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.77507)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
-0.3029717 13.1173688
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desformer_smokers, family = "gaussian", data = former_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = former_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.9907 0.1669 -5.935 3.01e-08 ***
perfluorooctane_sulfonic_acid_comment 10.7524 4.1875 2.568 0.0115 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 51.92854)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
2.460115 19.044762
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = desnon_smokers, family = "gaussian", data = non_smokers)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_smokers)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.8882 0.1182 -15.97 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 14.8581 2.9598 5.02 1.85e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.48088)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
8.996779 20.719365 #sensitivity analysis without cancer patient only for pfas
#subset cancer
cancer <- Fulldat_Pheno[Fulldat_Pheno$had_cancer == 1, ]
cancer$had_cancer[is.na(cancer$had_cancer)] <- 0
cancer <- cancer[cancer$had_cancer != 0, , drop = FALSE]
non_cancer <- Fulldat_Pheno[Fulldat_Pheno$had_cancer == 2, ]
non_cancer$had_cancer[is.na(non_cancer$had_cancer)] <- 0
non_cancer <- non_cancer[non_cancer$had_cancer != 0, , drop = FALSE]
des_cancer <- svydesign(id =~ psu, strata =~ Strata, weights =~ weight_2, nest = TRUE, data = cancer)
des_non_cancer <- svydesign(id =~ psu, strata =~ Strata, weights =~ weight_2, nest = TRUE, data = non_cancer) #sensitivity Perfluorohexane_sulfonic_acid
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.14044 0.08711 -13.09 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 7.84032 1.23862 6.33 3.48e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.69392)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.312744 -0.9681331
Perfluorohexane_sulfonic_acid_comment 5.390383 10.2902490
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.05895 0.35397 5.817 4.93e-08 ***
Perfluorohexane_sulfonic_acid_comment 7.50347 1.21185 6.192 8.30e-09 ***
Gender -1.53289 0.14751 -10.392 < 2e-16 ***
Race 0.15978 0.06895 2.317 0.02216 *
Marital_Status2 2.24997 0.30033 7.492 1.19e-11 ***
Marital_Status3 0.97063 0.29537 3.286 0.00133 **
Marital_Status4 0.97625 0.40354 2.419 0.01703 *
Marital_Status5 -0.51639 0.20075 -2.572 0.01130 *
Marital_Status6 -0.81617 0.29542 -2.763 0.00662 **
Marital_Status77 -2.17975 2.55510 -0.853 0.39528
Marital_Status99 3.10636 6.56941 0.473 0.63716
Marital_StatusNone 2.22194 1.41199 1.574 0.11816
Ratio_income_poverty -0.55301 0.04307 -12.839 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.72106)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.3582391 2.7596610
Perfluorohexane_sulfonic_acid_comment 5.1044851 9.9024556
Gender -1.8249023 -1.2408806
Race 0.0232826 0.2962780
Marital_Status2 1.6554473 2.8444966
Marital_Status3 0.3859060 1.5553520
Marital_Status4 0.1774070 1.7750883
Marital_Status5 -0.9137866 -0.1189905
Marital_Status6 -1.4009732 -0.2313620
Marital_Status77 -7.2378283 2.8783324
Marital_Status99 -9.8984366 16.1111536
Marital_StatusNone -0.5732363 5.0171182
Ratio_income_poverty -0.6382807 -0.4677452
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8.412e-01 7.566e-01 1.112 0.268612
Perfluorohexane_sulfonic_acid_comment 7.892e+00 1.232e+00 6.405 3.74e-09 ***
Gender -1.579e+00 1.439e-01 -10.971 < 2e-16 ***
Race 1.472e-01 6.548e-02 2.248 0.026559 *
Marital_Status2 2.205e+00 3.103e-01 7.105 1.22e-10 ***
Marital_Status3 8.878e-01 2.821e-01 3.147 0.002120 **
Marital_Status4 7.953e-01 3.327e-01 2.390 0.018528 *
Marital_Status5 -2.866e-01 1.831e-01 -1.565 0.120349
Marital_Status6 -6.250e-01 2.759e-01 -2.265 0.025425 *
Marital_Status77 -4.807e+00 1.975e+00 -2.434 0.016516 *
Marital_Status99 3.861e+00 5.477e+00 0.705 0.482352
Marital_StatusNone 5.480e+00 4.727e+00 1.159 0.248810
Ratio_income_poverty -4.491e-01 4.406e-02 -10.195 < 2e-16 ***
BMI 2.045e-01 1.187e-02 17.226 < 2e-16 ***
sleep_disorders2 -7.747e-01 2.145e-01 -3.611 0.000459 ***
sleep_disorders7 2.130e+01 2.579e+00 8.259 3.44e-13 ***
sleep_disorders9 -5.840e-02 5.155e+00 -0.011 0.990982
sleep_disordersNone -2.021e+00 2.227e-01 -9.074 4.81e-15 ***
Smoked_days -2.217e+00 2.989e-01 -7.418 2.55e-11 ***
now_smoke -6.723e-01 1.179e-01 -5.702 9.93e-08 ***
quit_smoking 1.406e-05 1.470e-05 0.956 0.341040
Avg_alcohol_drinks2 5.152e-01 2.316e-01 2.224 0.028161 *
Avg_alcohol_drinks9 2.858e+00 2.263e+00 1.263 0.209150
Avg_alcohol_drinksNone 2.707e-01 2.746e-01 0.986 0.326324
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.53696)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -6.580239e-01 2.340411e+00
Perfluorohexane_sulfonic_acid_comment 5.450433e+00 1.033394e+01
Gender -1.864222e+00 -1.293820e+00
Race 1.743972e-02 2.769609e-01
Marital_Status2 1.589807e+00 2.819510e+00
Marital_Status3 3.287237e-01 1.446859e+00
Marital_Status4 1.359356e-01 1.454611e+00
Marital_Status5 -6.493425e-01 7.619911e-02
Marital_Status6 -1.171596e+00 -7.831478e-02
Marital_Status77 -8.720918e+00 -8.940723e-01
Marital_Status99 -6.992130e+00 1.471333e+01
Marital_StatusNone -3.886746e+00 1.484711e+01
Ratio_income_poverty -5.364431e-01 -3.618469e-01
BMI 1.809979e-01 2.280521e-01
sleep_disorders2 -1.199847e+00 -3.496214e-01
sleep_disorders7 1.618831e+01 2.640913e+01
sleep_disorders9 -1.027398e+01 1.015718e+01
sleep_disordersNone -2.462213e+00 -1.579607e+00
Smoked_days -2.809605e+00 -1.625008e+00
now_smoke -9.059769e-01 -4.386635e-01
quit_smoking -1.507381e-05 4.318887e-05
Avg_alcohol_drinks2 5.619396e-02 9.741917e-01
Avg_alcohol_drinks9 -1.625433e+00 7.342356e+00
Avg_alcohol_drinksNone -2.733872e-01 8.148028e-01
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorohexane_sulfonic_acid),
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.6624 0.0998 -6.637 7.41e-10 ***
ln(Perfluorohexane_sulfonic_acid) -1.1008 0.1097 -10.033 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.40329)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -0.8597935 -0.4649820
ln(Perfluorohexane_sulfonic_acid) -1.3178516 -0.8838188
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorohexane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.52671 0.34921 10.099 < 2e-16 ***
ln(Perfluorohexane_sulfonic_acid) -1.37166 0.11238 -12.205 < 2e-16 ***
Gender -2.30387 0.15904 -14.486 < 2e-16 ***
Race 0.18730 0.06976 2.685 0.00826 **
Marital_Status2 2.55124 0.31221 8.172 3.25e-13 ***
Marital_Status3 0.93353 0.29982 3.114 0.00230 **
Marital_Status4 1.07402 0.41284 2.602 0.01043 *
Marital_Status5 -0.42486 0.20113 -2.112 0.03670 *
Marital_Status6 -0.79174 0.28042 -2.823 0.00555 **
Marital_Status77 -1.60660 2.21289 -0.726 0.46922
Marital_Status99 4.46715 5.21019 0.857 0.39291
Marital_StatusNone 2.46760 1.52312 1.620 0.10779
Ratio_income_poverty -0.49943 0.04100 -12.181 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.96024)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 2.83541650 4.21800449
ln(Perfluorohexane_sulfonic_acid) -1.59413114 -1.14918255
Gender -2.61871387 -1.98903445
Race 0.04920198 0.32539784
Marital_Status2 1.93319697 3.16928038
Marital_Status3 0.34000299 1.52706345
Marital_Status4 0.25675917 1.89128150
Marital_Status5 -0.82302940 -0.02669751
Marital_Status6 -1.34686089 -0.23662512
Marital_Status77 -5.98722766 2.77403220
Marital_Status99 -5.84693429 14.78124048
Marital_StatusNone -0.54757388 5.48277342
Ratio_income_poverty -0.58059068 -0.41826537
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorohexane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.368e+00 7.906e-01 2.996 0.003379 **
ln(Perfluorohexane_sulfonic_acid) -1.279e+00 1.100e-01 -11.632 < 2e-16 ***
Gender -2.281e+00 1.577e-01 -14.465 < 2e-16 ***
Race 1.746e-01 6.627e-02 2.635 0.009618 **
Marital_Status2 2.474e+00 3.227e-01 7.668 7.18e-12 ***
Marital_Status3 8.562e-01 2.848e-01 3.007 0.003268 **
Marital_Status4 8.542e-01 3.462e-01 2.467 0.015141 *
Marital_Status5 -2.175e-01 1.869e-01 -1.164 0.246989
Marital_Status6 -6.566e-01 2.660e-01 -2.469 0.015090 *
Marital_Status77 -3.970e+00 1.953e+00 -2.033 0.044426 *
Marital_Status99 5.216e+00 4.256e+00 1.226 0.222954
Marital_StatusNone 5.356e+00 4.861e+00 1.102 0.272960
Ratio_income_poverty -4.073e-01 4.232e-02 -9.624 2.64e-16 ***
BMI 1.983e-01 1.189e-02 16.672 < 2e-16 ***
sleep_disorders2 -7.604e-01 2.100e-01 -3.621 0.000443 ***
sleep_disorders7 2.067e+01 2.558e+00 8.079 8.71e-13 ***
sleep_disorders9 1.721e+00 4.446e+00 0.387 0.699444
sleep_disordersNone -1.581e+00 2.276e-01 -6.947 2.68e-10 ***
Smoked_days -2.275e+00 2.988e-01 -7.614 9.44e-12 ***
now_smoke -6.535e-01 1.186e-01 -5.508 2.37e-07 ***
quit_smoking 1.717e-05 1.484e-05 1.156 0.249977
Avg_alcohol_drinks2 5.442e-01 2.396e-01 2.271 0.025047 *
Avg_alcohol_drinks9 3.220e+00 2.453e+00 1.312 0.192069
Avg_alcohol_drinksNone 1.777e-01 2.664e-01 0.667 0.506154
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.06252)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 8.016655e-01 3.934812e+00
ln(Perfluorohexane_sulfonic_acid) -1.496841e+00 -1.061086e+00
Gender -2.593080e+00 -1.968225e+00
Race 4.330006e-02 3.059400e-01
Marital_Status2 1.834754e+00 3.113526e+00
Marital_Status3 2.918953e-01 1.420514e+00
Marital_Status4 1.681631e-01 1.540192e+00
Marital_Status5 -5.877204e-01 1.527955e-01
Marital_Status6 -1.183631e+00 -1.295357e-01
Marital_Status77 -7.839664e+00 -1.006741e-01
Marital_Status99 -3.217602e+00 1.364992e+01
Marital_StatusNone -4.277192e+00 1.498911e+01
Ratio_income_poverty -4.911761e-01 -3.234411e-01
BMI 1.747134e-01 2.218459e-01
sleep_disorders2 -1.176533e+00 -3.443411e-01
sleep_disorders7 1.559809e+01 2.573608e+01
sleep_disorders9 -7.089016e+00 1.053081e+01
sleep_disordersNone -2.032185e+00 -1.130111e+00
Smoked_days -2.867441e+00 -1.683152e+00
now_smoke -8.885998e-01 -4.183860e-01
quit_smoking -1.224764e-05 4.657981e-05
Avg_alcohol_drinks2 6.945557e-02 1.018980e+00
Avg_alcohol_drinks9 -1.641599e+00 8.081783e+00
Avg_alcohol_drinksNone -3.502257e-01 7.056312e-01#sensitivity “Perfluorononanoic_acid” “Perfluorononanoic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.15765 0.09109 -12.709 < 2e-16 ***
Perfluorononanoic_acid_comment 6.98972 1.64974 4.237 4.21e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.66237)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.337828 -0.9774805
Perfluorononanoic_acid_comment 3.726592 10.2528395
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.02575 0.34127 5.936 2.82e-08 ***
Perfluorononanoic_acid_comment 6.47270 1.58464 4.085 7.93e-05 ***
Gender -1.52942 0.14723 -10.388 < 2e-16 ***
Race 0.16090 0.06909 2.329 0.02151 *
Marital_Status2 2.18296 0.30124 7.247 4.24e-11 ***
Marital_Status3 0.91661 0.28292 3.240 0.00154 **
Marital_Status4 1.05297 0.40934 2.572 0.01130 *
Marital_Status5 -0.43412 0.20413 -2.127 0.03546 *
Marital_Status6 -0.78872 0.29774 -2.649 0.00914 **
Marital_Status77 -2.16140 2.55508 -0.846 0.39925
Marital_Status99 3.12148 6.56507 0.475 0.63530
Marital_StatusNone 2.23611 1.41214 1.583 0.11590
Ratio_income_poverty -0.55026 0.04117 -13.364 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.74044)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.35016765 2.70133201
Perfluorononanoic_acid_comment 3.33575118 9.60965632
Gender -1.82087235 -1.23796947
Race 0.02413385 0.29766241
Marital_Status2 1.58663050 2.77929519
Marital_Status3 0.35655486 1.47667366
Marital_Status4 0.24264574 1.86329266
Marital_Status5 -0.83822015 -0.03001848
Marital_Status6 -1.37812148 -0.19931678
Marital_Status77 -7.21942559 2.89663428
Marital_Status99 -9.87472998 16.11769583
Marital_StatusNone -0.55935940 5.03157318
Ratio_income_poverty -0.63177289 -0.46875538
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 7.737e-01 7.445e-01 1.039 0.300961
Perfluorononanoic_acid_comment 6.303e+00 1.709e+00 3.689 0.000351 ***
Gender -1.581e+00 1.461e-01 -10.821 < 2e-16 ***
Race 1.478e-01 6.480e-02 2.280 0.024489 *
Marital_Status2 2.139e+00 3.128e-01 6.836 4.61e-10 ***
Marital_Status3 8.432e-01 2.730e-01 3.088 0.002544 **
Marital_Status4 8.607e-01 3.419e-01 2.517 0.013250 *
Marital_Status5 -2.100e-01 1.868e-01 -1.124 0.263369
Marital_Status6 -6.115e-01 2.790e-01 -2.192 0.030482 *
Marital_Status77 -4.760e+00 1.969e+00 -2.418 0.017234 *
Marital_Status99 3.911e+00 5.479e+00 0.714 0.476779
Marital_StatusNone 5.338e+00 4.682e+00 1.140 0.256743
Ratio_income_poverty -4.500e-01 4.197e-02 -10.722 < 2e-16 ***
BMI 2.034e-01 1.190e-02 17.086 < 2e-16 ***
sleep_disorders2 -7.872e-01 2.056e-01 -3.829 0.000213 ***
sleep_disorders7 2.109e+01 2.544e+00 8.289 2.94e-13 ***
sleep_disorders9 -7.027e-02 5.186e+00 -0.014 0.989213
sleep_disordersNone -1.853e+00 2.187e-01 -8.473 1.13e-13 ***
Smoked_days -2.189e+00 2.890e-01 -7.574 1.16e-11 ***
now_smoke -6.599e-01 1.168e-01 -5.648 1.26e-07 ***
quit_smoking 1.185e-05 1.495e-05 0.793 0.429478
Avg_alcohol_drinks2 5.542e-01 2.399e-01 2.310 0.022732 *
Avg_alcohol_drinks9 2.828e+00 2.212e+00 1.279 0.203679
Avg_alcohol_drinksNone 3.292e-01 2.704e-01 1.217 0.226108
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.66472)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -7.015651e-01 2.248933e+00
Perfluorononanoic_acid_comment 2.916816e+00 9.688869e+00
Gender -1.871024e+00 -1.291831e+00
Race 1.936924e-02 2.761732e-01
Marital_Status2 1.518691e+00 2.758422e+00
Marital_Status3 3.021650e-01 1.384246e+00
Marital_Status4 1.832079e-01 1.538251e+00
Marital_Status5 -5.802606e-01 1.601931e-01
Marital_Status6 -1.164272e+00 -5.865080e-02
Marital_Status77 -8.661671e+00 -8.592209e-01
Marital_Status99 -6.944966e+00 1.476747e+01
Marital_StatusNone -3.940565e+00 1.461635e+01
Ratio_income_poverty -5.331216e-01 -3.668056e-01
BMI 1.797817e-01 2.269525e-01
sleep_disorders2 -1.194623e+00 -3.798747e-01
sleep_disorders7 1.604680e+01 2.612959e+01
sleep_disorders9 -1.034690e+01 1.020635e+01
sleep_disordersNone -2.286821e+00 -1.419964e+00
Smoked_days -2.761329e+00 -1.616095e+00
now_smoke -8.914125e-01 -4.283907e-01
quit_smoking -1.776448e-05 4.146831e-05
Avg_alcohol_drinks2 7.881362e-02 1.029582e+00
Avg_alcohol_drinks9 -1.554839e+00 7.211835e+00
Avg_alcohol_drinksNone -2.067012e-01 8.650220e-01
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorononanoic_acid),
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.41625 0.09072 -15.61 <2e-16 ***
ln(Perfluorononanoic_acid) -1.76367 0.14027 -12.57 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.4502)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.595689 -1.236820
ln(Perfluorononanoic_acid) -2.041110 -1.486222
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorononanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.7172 0.3368 5.098 1.27e-06 ***
ln(Perfluorononanoic_acid) -1.7542 0.1398 -12.548 < 2e-16 ***
Gender -1.8261 0.1488 -12.270 < 2e-16 ***
Race 0.2766 0.0657 4.210 4.91e-05 ***
Marital_Status2 2.1462 0.3084 6.959 1.85e-10 ***
Marital_Status3 0.7238 0.2884 2.510 0.01340 *
Marital_Status4 1.1204 0.3790 2.956 0.00374 **
Marital_Status5 -0.3204 0.1980 -1.619 0.10813
Marital_Status6 -0.6157 0.2860 -2.153 0.03330 *
Marital_Status77 -1.4656 2.4903 -0.589 0.55727
Marital_Status99 2.8258 5.7961 0.488 0.62676
Marital_StatusNone 3.4975 1.2537 2.790 0.00612 **
Ratio_income_poverty -0.5006 0.0401 -12.484 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.51127)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.0503295 2.38398653
ln(Perfluorononanoic_acid) -2.0308975 -1.47743685
Gender -2.1207235 -1.53150738
Race 0.1465516 0.40668803
Marital_Status2 1.5356595 2.75672157
Marital_Status3 0.1528728 1.29475871
Marital_Status4 0.3701282 1.87061925
Marital_Status5 -0.7123563 0.07148234
Marital_Status6 -1.1817923 -0.04954469
Marital_Status77 -6.3952752 3.46414780
Marital_Status99 -8.6481484 14.29964751
Marital_StatusNone 1.0157338 5.97921893
Ratio_income_poverty -0.5800033 -0.42123710
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorononanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.965e-01 7.271e-01 0.545 0.586665
ln(Perfluorononanoic_acid) -1.633e+00 1.458e-01 -11.197 < 2e-16 ***
Gender -1.860e+00 1.451e-01 -12.820 < 2e-16 ***
Race 2.565e-01 6.129e-02 4.184 5.75e-05 ***
Marital_Status2 2.107e+00 3.179e-01 6.627 1.28e-09 ***
Marital_Status3 6.843e-01 2.743e-01 2.495 0.014061 *
Marital_Status4 9.118e-01 3.263e-01 2.794 0.006134 **
Marital_Status5 -1.187e-01 1.829e-01 -0.649 0.517561
Marital_Status6 -4.860e-01 2.676e-01 -1.816 0.072060 .
Marital_Status77 -4.223e+00 1.717e+00 -2.460 0.015455 *
Marital_Status99 3.753e+00 4.781e+00 0.785 0.434181
Marital_StatusNone 6.118e+00 4.449e+00 1.375 0.171816
Ratio_income_poverty -4.083e-01 4.119e-02 -9.914 < 2e-16 ***
BMI 1.986e-01 1.197e-02 16.591 < 2e-16 ***
sleep_disorders2 -7.961e-01 2.051e-01 -3.882 0.000177 ***
sleep_disorders7 2.167e+01 2.325e+00 9.319 1.33e-15 ***
sleep_disorders9 -1.137e-01 5.390e+00 -0.021 0.983210
sleep_disordersNone -1.406e+00 2.293e-01 -6.132 1.36e-08 ***
Smoked_days -2.133e+00 2.872e-01 -7.427 2.44e-11 ***
now_smoke -6.105e-01 1.141e-01 -5.351 4.75e-07 ***
quit_smoking 1.407e-05 1.513e-05 0.930 0.354467
Avg_alcohol_drinks2 6.224e-01 2.340e-01 2.659 0.008985 **
Avg_alcohol_drinks9 2.250e+00 1.606e+00 1.401 0.164081
Avg_alcohol_drinksNone 4.569e-01 2.675e-01 1.708 0.090440 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.70616)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.044382e+00 1.837343e+00
ln(Perfluorononanoic_acid) -1.921558e+00 -1.343691e+00
Gender -2.147662e+00 -1.572639e+00
Race 1.349976e-01 3.779082e-01
Marital_Status2 1.477051e+00 2.737057e+00
Marital_Status3 1.408590e-01 1.227749e+00
Marital_Status4 2.651336e-01 1.558454e+00
Marital_Status5 -4.811173e-01 2.436693e-01
Marital_Status6 -1.016257e+00 4.428887e-02
Marital_Status77 -7.624653e+00 -8.205670e-01
Marital_Status99 -5.721790e+00 1.322779e+01
Marital_StatusNone -2.697222e+00 1.493362e+01
Ratio_income_poverty -4.899157e-01 -3.266903e-01
BMI 1.748445e-01 2.222737e-01
sleep_disorders2 -1.202494e+00 -3.896809e-01
sleep_disorders7 1.705893e+01 2.627312e+01
sleep_disorders9 -1.079477e+01 1.056738e+01
sleep_disordersNone -1.860534e+00 -9.517407e-01
Smoked_days -2.701900e+00 -1.563733e+00
now_smoke -8.365672e-01 -3.844202e-01
quit_smoking -1.591723e-05 4.406283e-05
Avg_alcohol_drinks2 1.586491e-01 1.086179e+00
Avg_alcohol_drinks9 -9.330698e-01 5.433523e+00
Avg_alcohol_drinksNone -7.319217e-02 9.869533e-01#sensitivity “perfluorooctanoic_acid” “perfluorooctanoic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.56511 0.09009 -28.47 <2e-16 ***
perfluorooctanoic_acid_comment 3.40019 0.20672 16.45 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.88168)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -2.743518 -2.386709
perfluorooctanoic_acid_comment 2.990827 3.809553
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.75689 0.33538 2.257 0.02605 *
perfluorooctanoic_acid_comment 3.22550 0.20582 15.671 < 2e-16 ***
Gender -1.53814 0.14501 -10.607 < 2e-16 ***
Race 0.11634 0.07174 1.622 0.10784
Marital_Status2 1.81547 0.31135 5.831 5.92e-08 ***
Marital_Status3 0.76375 0.27971 2.731 0.00740 **
Marital_Status4 1.05077 0.34826 3.017 0.00319 **
Marital_Status5 0.10850 0.18382 0.590 0.55625
Marital_Status6 -0.24953 0.29982 -0.832 0.40711
Marital_Status77 -1.01787 3.37154 -0.302 0.76331
Marital_Status99 1.36573 6.57411 0.208 0.83582
Marital_StatusNone 3.66221 1.41276 2.592 0.01087 *
Ratio_income_poverty -0.55241 0.03997 -13.819 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.08648)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 0.09203898 1.4217502
perfluorooctanoic_acid_comment 2.81747539 3.6335169
Gender -1.82559856 -1.2506784
Race -0.02588656 0.2585636
Marital_Status2 1.19824992 2.4326950
Marital_Status3 0.20926233 1.3182362
Marital_Status4 0.36039545 1.7411498
Marital_Status5 -0.25589901 0.4729084
Marital_Status6 -0.84389672 0.3448306
Marital_Status77 -7.70154757 5.6658077
Marital_Status99 -11.66668693 14.3981379
Marital_StatusNone 0.86158416 6.4628409
Ratio_income_poverty -0.63165624 -0.4731704
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -4.928e-01 6.884e-01 -0.716 0.47586
perfluorooctanoic_acid_comment 3.374e+00 2.098e-01 16.077 < 2e-16 ***
Gender -1.549e+00 1.434e-01 -10.803 < 2e-16 ***
Race 1.033e-01 6.546e-02 1.578 0.11775
Marital_Status2 1.791e+00 3.181e-01 5.630 1.79e-07 ***
Marital_Status3 7.398e-01 2.690e-01 2.750 0.00712 **
Marital_Status4 8.314e-01 3.067e-01 2.711 0.00795 **
Marital_Status5 2.499e-01 1.751e-01 1.427 0.15686
Marital_Status6 -2.938e-01 2.891e-01 -1.016 0.31203
Marital_Status77 -4.037e+00 2.644e+00 -1.527 0.13009
Marital_Status99 2.558e+00 5.607e+00 0.456 0.64931
Marital_StatusNone 5.517e+00 4.968e+00 1.110 0.26957
Ratio_income_poverty -4.671e-01 4.109e-02 -11.367 < 2e-16 ***
BMI 1.968e-01 1.171e-02 16.807 < 2e-16 ***
sleep_disorders2 -5.646e-01 2.156e-01 -2.619 0.01025 *
sleep_disorders7 2.359e+01 3.086e+00 7.643 1.60e-11 ***
sleep_disorders9 -3.913e+00 2.273e+00 -1.721 0.08845 .
sleep_disordersNone 1.234e-01 2.339e-01 0.528 0.59893
Smoked_days -2.462e+00 2.863e-01 -8.600 1.50e-13 ***
now_smoke -6.878e-01 1.132e-01 -6.078 2.45e-08 ***
quit_smoking -3.628e-06 1.265e-05 -0.287 0.77492
Avg_alcohol_drinks2 5.654e-01 2.161e-01 2.616 0.01033 *
Avg_alcohol_drinks9 3.829e+00 2.673e+00 1.432 0.15528
Avg_alcohol_drinksNone 6.014e-01 2.762e-01 2.177 0.03191 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 36.37583)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.859225e+00 8.737181e-01
perfluorooctanoic_acid_comment 2.957074e+00 3.790118e+00
Gender -1.833425e+00 -1.264257e+00
Race -2.661341e-02 2.332746e-01
Marital_Status2 1.159612e+00 2.422554e+00
Marital_Status3 2.058595e-01 1.273685e+00
Marital_Status4 2.226623e-01 1.440126e+00
Marital_Status5 -9.773791e-02 5.975097e-01
Marital_Status6 -8.675834e-01 2.800060e-01
Marital_Status77 -9.285568e+00 1.211421e+00
Marital_Status99 -8.572890e+00 1.368866e+01
Marital_StatusNone -4.344557e+00 1.537834e+01
Ratio_income_poverty -5.486638e-01 -3.855294e-01
BMI 1.735470e-01 2.200292e-01
sleep_disorders2 -9.925864e-01 -1.366969e-01
sleep_disorders7 1.746207e+01 2.971359e+01
sleep_disorders9 -8.425427e+00 5.999165e-01
sleep_disordersNone -3.409201e-01 5.878162e-01
Smoked_days -3.030787e+00 -1.894104e+00
now_smoke -9.124406e-01 -4.631984e-01
quit_smoking -2.874502e-05 2.148834e-05
Avg_alcohol_drinks2 1.364257e-01 9.943880e-01
Avg_alcohol_drinks9 -1.477165e+00 9.135325e+00
Avg_alcohol_drinksNone 5.311342e-02 1.149603e+00
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctanoic_acid),
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.4894 0.1560 3.137 0.00216 **
ln(perfluorooctanoic_acid) -2.0684 0.1223 -16.910 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.78158)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 0.1804373 0.7983283
ln(perfluorooctanoic_acid) -2.3106715 -1.8262166
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.22545 0.37542 11.255 < 2e-16 ***
ln(perfluorooctanoic_acid) -2.14248 0.13357 -16.040 < 2e-16 ***
Gender -2.21051 0.15228 -14.516 < 2e-16 ***
Race 0.19358 0.07466 2.593 0.010856 *
Marital_Status2 2.45472 0.32271 7.607 1.16e-11 ***
Marital_Status3 0.95025 0.27599 3.443 0.000822 ***
Marital_Status4 1.00373 0.38448 2.611 0.010333 *
Marital_Status5 -0.05240 0.19392 -0.270 0.787535
Marital_Status6 -0.43735 0.29498 -1.483 0.141113
Marital_Status77 -0.73238 2.87674 -0.255 0.799532
Marital_Status99 3.02410 5.78699 0.523 0.602355
Marital_StatusNone 3.77963 1.10343 3.425 0.000872 ***
Ratio_income_poverty -0.43051 0.04171 -10.321 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.6031)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 3.4812174 4.9696866
ln(perfluorooctanoic_acid) -2.4072634 -1.8776913
Gender -2.5123853 -1.9086407
Race 0.0455633 0.3415886
Marital_Status2 1.8149875 3.0944533
Marital_Status3 0.4031339 1.4973645
Marital_Status4 0.2415524 1.7659124
Marital_Status5 -0.4368298 0.3320368
Marital_Status6 -1.0221251 0.1474211
Marital_Status77 -6.4351852 4.9704305
Marital_Status99 -8.4479268 14.4961217
Marital_StatusNone 1.5922186 5.9670460
Ratio_income_poverty -0.5132010 -0.3478220
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.991e+00 6.901e-01 4.334 3.61e-05 ***
ln(perfluorooctanoic_acid) -2.168e+00 1.381e-01 -15.702 < 2e-16 ***
Gender -2.187e+00 1.479e-01 -14.788 < 2e-16 ***
Race 1.949e-01 6.843e-02 2.848 0.005376 **
Marital_Status2 2.425e+00 3.294e-01 7.363 6.18e-11 ***
Marital_Status3 9.041e-01 2.619e-01 3.452 0.000829 ***
Marital_Status4 7.780e-01 3.361e-01 2.314 0.022776 *
Marital_Status5 7.733e-02 1.827e-01 0.423 0.673107
Marital_Status6 -5.046e-01 2.859e-01 -1.765 0.080742 .
Marital_Status77 -3.719e+00 2.353e+00 -1.581 0.117226
Marital_Status99 4.109e+00 4.796e+00 0.857 0.393639
Marital_StatusNone 5.685e+00 4.578e+00 1.242 0.217388
Ratio_income_poverty -3.568e-01 4.021e-02 -8.873 3.91e-14 ***
BMI 1.980e-01 1.190e-02 16.646 < 2e-16 ***
sleep_disorders2 -6.768e-01 2.071e-01 -3.269 0.001500 **
sleep_disorders7 2.328e+01 2.791e+00 8.342 5.33e-13 ***
sleep_disorders9 -4.221e+00 3.123e+00 -1.352 0.179676
sleep_disordersNone -1.191e-01 2.241e-01 -0.531 0.596466
Smoked_days -2.371e+00 2.836e-01 -8.361 4.87e-13 ***
now_smoke -6.575e-01 1.174e-01 -5.598 2.05e-07 ***
quit_smoking 5.903e-06 1.569e-05 0.376 0.707534
Avg_alcohol_drinks2 2.326e-01 2.134e-01 1.090 0.278380
Avg_alcohol_drinks9 3.035e+00 2.290e+00 1.325 0.188212
Avg_alcohol_drinksNone 2.816e-01 2.727e-01 1.033 0.304303
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 35.90674)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.620984e+00 4.360464e+00
ln(perfluorooctanoic_acid) -2.442505e+00 -1.894260e+00
Gender -2.480425e+00 -1.893343e+00
Race 5.908115e-02 3.307428e-01
Marital_Status2 1.771414e+00 3.079137e+00
Marital_Status3 3.842142e-01 1.424022e+00
Marital_Status4 1.107383e-01 1.445199e+00
Marital_Status5 -2.854041e-01 4.400700e-01
Marital_Status6 -1.072183e+00 6.289301e-02
Marital_Status77 -8.388563e+00 9.509954e-01
Marital_Status99 -5.409970e+00 1.362871e+01
Marital_StatusNone -3.403302e+00 1.477310e+01
Ratio_income_poverty -4.366430e-01 -2.770010e-01
BMI 1.744032e-01 2.216300e-01
sleep_disorders2 -1.087823e+00 -2.658065e-01
sleep_disorders7 1.774318e+01 2.882312e+01
sleep_disorders9 -1.041893e+01 1.977791e+00
sleep_disordersNone -5.639924e-01 3.258387e-01
Smoked_days -2.933915e+00 -1.808096e+00
now_smoke -8.906308e-01 -4.243756e-01
quit_smoking -2.523550e-05 3.704108e-05
Avg_alcohol_drinks2 -1.909356e-01 6.561489e-01
Avg_alcohol_drinks9 -1.510534e+00 7.580051e+00
Avg_alcohol_drinksNone -2.596379e-01 8.228519e-01#sensitivity “perfluorooctane_sulfonic_acid” “perfluorooctane_sulfonic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.44500 0.09115 -15.853 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 11.99581 2.33803 5.131 1.15e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.82679)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.625501 -1.264492
perfluorooctane_sulfonic_acid_comment 7.365872 16.625754
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.94063 0.37759 5.139 1.25e-06 ***
perfluorooctane_sulfonic_acid_comment 11.65931 2.34984 4.962 2.65e-06 ***
Gender -1.61964 0.14268 -11.352 < 2e-16 ***
Race 0.14688 0.07699 1.908 0.05909 .
Marital_Status2 2.33293 0.31736 7.351 4.14e-11 ***
Marital_Status3 1.07554 0.28213 3.812 0.00023 ***
Marital_Status4 0.67182 0.40820 1.646 0.10274
Marital_Status5 -0.49588 0.20130 -2.463 0.01536 *
Marital_Status6 -0.93247 0.30717 -3.036 0.00301 **
Marital_Status77 -1.33411 3.01920 -0.442 0.65947
Marital_Status99 3.47466 6.58801 0.527 0.59899
Marital_StatusNone 2.51239 1.40115 1.793 0.07578 .
Ratio_income_poverty -0.56328 0.04277 -13.171 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.68774)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.192091751 2.68916615
perfluorooctane_sulfonic_acid_comment 7.001034132 16.31759525
Gender -1.902482501 -1.33679575
Race -0.005735359 0.29950196
Marital_Status2 1.703802079 2.96205646
Marital_Status3 0.516254576 1.63481740
Marital_Status4 -0.137399701 1.48103061
Marital_Status5 -0.894932879 -0.09682369
Marital_Status6 -1.541401869 -0.32353481
Marital_Status77 -7.319332496 4.65110353
Marital_Status99 -9.585299654 16.53462520
Marital_StatusNone -0.265229376 5.29000456
Ratio_income_poverty -0.648063586 -0.47850617
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.231e-01 7.274e-01 0.857 0.393804
perfluorooctane_sulfonic_acid_comment 1.147e+01 2.576e+00 4.453 2.29e-05 ***
Gender -1.592e+00 1.395e-01 -11.412 < 2e-16 ***
Race 1.398e-01 7.200e-02 1.942 0.055016 .
Marital_Status2 2.267e+00 3.233e-01 7.012 3.28e-10 ***
Marital_Status3 9.923e-01 2.702e-01 3.673 0.000395 ***
Marital_Status4 5.168e-01 3.390e-01 1.524 0.130677
Marital_Status5 -3.226e-01 1.890e-01 -1.707 0.091099 .
Marital_Status6 -8.672e-01 2.977e-01 -2.913 0.004458 **
Marital_Status77 -4.231e+00 2.574e+00 -1.644 0.103509
Marital_Status99 4.235e+00 5.571e+00 0.760 0.448986
Marital_StatusNone 5.317e+00 4.836e+00 1.099 0.274384
Ratio_income_poverty -4.822e-01 4.293e-02 -11.232 < 2e-16 ***
BMI 2.066e-01 1.208e-02 17.111 < 2e-16 ***
sleep_disorders2 -6.633e-01 2.212e-01 -2.998 0.003460 **
sleep_disorders7 2.142e+01 3.028e+00 7.075 2.44e-10 ***
sleep_disorders9 -3.571e+00 3.302e+00 -1.081 0.282214
sleep_disordersNone -1.312e+00 2.286e-01 -5.741 1.10e-07 ***
Smoked_days -2.283e+00 2.966e-01 -7.697 1.24e-11 ***
now_smoke -6.911e-01 1.183e-01 -5.840 7.09e-08 ***
quit_smoking 9.211e-06 1.631e-05 0.565 0.573606
Avg_alcohol_drinks2 1.385e-01 2.247e-01 0.617 0.538970
Avg_alcohol_drinks9 3.565e+00 2.849e+00 1.251 0.213838
Avg_alcohol_drinksNone 2.905e-01 2.802e-01 1.037 0.302535
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 37.76279)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -8.207566e-01 2.066869e+00
perfluorooctane_sulfonic_acid_comment 6.358214e+00 1.658637e+01
Gender -1.868832e+00 -1.315055e+00
Race -3.063835e-03 2.827601e-01
Marital_Status2 1.625246e+00 2.908774e+00
Marital_Status3 4.559677e-01 1.528565e+00
Marital_Status4 -1.561039e-01 1.189653e+00
Marital_Status5 -6.978752e-01 5.259367e-02
Marital_Status6 -1.458184e+00 -2.761872e-01
Marital_Status77 -9.340268e+00 8.784676e-01
Marital_Status99 -6.823003e+00 1.529318e+01
Marital_StatusNone -4.283423e+00 1.491659e+01
Ratio_income_poverty -5.673826e-01 -3.969615e-01
BMI 1.826603e-01 2.306002e-01
sleep_disorders2 -1.102472e+00 -2.241199e-01
sleep_disorders7 1.541361e+01 2.743600e+01
sleep_disorders9 -1.012558e+01 2.983606e+00
sleep_disordersNone -1.766213e+00 -8.586239e-01
Smoked_days -2.872039e+00 -1.694424e+00
now_smoke -9.259536e-01 -4.562004e-01
quit_smoking -2.316702e-05 4.158849e-05
Avg_alcohol_drinks2 -3.074798e-01 5.845714e-01
Avg_alcohol_drinks9 -2.089847e+00 9.219482e+00
Avg_alcohol_drinksNone -2.657794e-01 8.467610e-01
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctane_sulfonic_acid),
design = des_non_cancer, family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.39439 0.24297 5.739 7.52e-08 ***
ln(perfluorooctane_sulfonic_acid) -1.29961 0.09325 -13.937 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.86485)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 0.9132372 1.875540
ln(perfluorooctane_sulfonic_acid) -1.4842737 -1.114953
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des_non_cancer,
family = "gaussian", data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.51886 0.45001 12.264 < 2e-16 ***
ln(perfluorooctane_sulfonic_acid) -1.45093 0.10267 -14.132 < 2e-16 ***
Gender -2.25366 0.15991 -14.093 < 2e-16 ***
Race 0.27189 0.07808 3.482 0.000722 ***
Marital_Status2 2.47441 0.33663 7.351 4.15e-11 ***
Marital_Status3 0.77935 0.28205 2.763 0.006741 **
Marital_Status4 0.74595 0.37966 1.965 0.052030 .
Marital_Status5 -0.45480 0.19690 -2.310 0.022816 *
Marital_Status6 -0.70561 0.29950 -2.356 0.020295 *
Marital_Status77 -0.69088 2.82973 -0.244 0.807581
Marital_Status99 3.66399 5.72973 0.639 0.523884
Marital_StatusNone 3.74761 1.51976 2.466 0.015254 *
Ratio_income_poverty -0.50339 0.04096 -12.291 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.42912)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 4.626763065 6.4109626
ln(perfluorooctane_sulfonic_acid) -1.654455621 -1.2473974
Gender -2.570654656 -1.9366576
Race 0.117093709 0.4266769
Marital_Status2 1.807089518 3.1417381
Marital_Status3 0.220218531 1.3384734
Marital_Status4 -0.006675801 1.4985841
Marital_Status5 -0.845129094 -0.0644763
Marital_Status6 -1.299332064 -0.1118939
Marital_Status77 -6.300490124 4.9187252
Marital_Status99 -7.694532745 15.0225148
Marital_StatusNone 0.734866913 6.7603456
Ratio_income_poverty -0.584580716 -0.4222025
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks, design = des_non_cancer, family = "gaussian",
data = non_cancer)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = non_cancer)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.051e+00 7.338e-01 5.521 2.87e-07 ***
ln(perfluorooctane_sulfonic_acid) -1.547e+00 1.170e-01 -13.229 < 2e-16 ***
Gender -2.300e+00 1.576e-01 -14.588 < 2e-16 ***
Race 2.824e-01 7.246e-02 3.898 0.00018 ***
Marital_Status2 2.476e+00 3.441e-01 7.196 1.37e-10 ***
Marital_Status3 7.493e-01 2.670e-01 2.806 0.00607 **
Marital_Status4 5.177e-01 3.301e-01 1.568 0.12009
Marital_Status5 -3.335e-01 1.849e-01 -1.804 0.07432 .
Marital_Status6 -7.636e-01 2.914e-01 -2.620 0.01022 *
Marital_Status77 -3.407e+00 2.300e+00 -1.482 0.14172
Marital_Status99 4.859e+00 4.712e+00 1.031 0.30503
Marital_StatusNone 5.403e+00 4.773e+00 1.132 0.26046
Ratio_income_poverty -4.297e-01 4.019e-02 -10.692 < 2e-16 ***
BMI 1.993e-01 1.194e-02 16.686 < 2e-16 ***
sleep_disorders2 -6.340e-01 2.193e-01 -2.891 0.00474 **
sleep_disorders7 2.164e+01 2.751e+00 7.865 5.47e-12 ***
sleep_disorders9 -2.900e+00 2.277e+00 -1.274 0.20583
sleep_disordersNone 2.508e-01 2.530e-01 0.991 0.32414
Smoked_days -2.188e+00 2.864e-01 -7.638 1.65e-11 ***
now_smoke -6.290e-01 1.157e-01 -5.435 4.15e-07 ***
quit_smoking 1.310e-05 1.647e-05 0.795 0.42846
Avg_alcohol_drinks2 3.639e-01 2.146e-01 1.696 0.09321 .
Avg_alcohol_drinks9 4.258e+00 2.353e+00 1.810 0.07350 .
Avg_alcohol_drinksNone 2.769e-01 2.779e-01 0.996 0.32153
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 36.7221)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 2.594792e+00 5.507911e+00
ln(perfluorooctane_sulfonic_acid) -1.779507e+00 -1.315146e+00
Gender -2.612403e+00 -1.986621e+00
Race 1.386023e-01 4.262525e-01
Marital_Status2 1.793395e+00 3.159573e+00
Marital_Status3 2.192931e-01 1.279262e+00
Marital_Status4 -1.375426e-01 1.173002e+00
Marital_Status5 -7.004616e-01 3.339734e-02
Marital_Status6 -1.342084e+00 -1.850802e-01
Marital_Status77 -7.971845e+00 1.157617e+00
Marital_Status99 -4.494152e+00 1.421227e+01
Marital_StatusNone -4.071613e+00 1.487818e+01
Ratio_income_poverty -5.094488e-01 -3.499093e-01
BMI 1.755518e-01 2.229596e-01
sleep_disorders2 -1.069233e+00 -1.987391e-01
sleep_disorders7 1.617764e+01 2.709947e+01
sleep_disorders9 -7.419072e+00 1.619238e+00
sleep_disordersNone -2.514945e-01 7.530530e-01
Smoked_days -2.756338e+00 -1.619197e+00
now_smoke -8.587464e-01 -3.992709e-01
quit_smoking -1.959961e-05 4.579843e-05
Avg_alcohol_drinks2 -6.212325e-02 7.899698e-01
Avg_alcohol_drinks9 -4.129324e-01 8.929024e+00
Avg_alcohol_drinksNone -2.746704e-01 8.284324e-01