library("haven")警告： 套件 ‘haven’ 是用 R 版本 4.3.3 來建造的library("tidyverse")警告： 套件 ‘tidyverse’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘ggplot2’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘tidyr’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘readr’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘purrr’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘forcats’ 是用 R 版本 4.3.3 來建造的警告： 套件 ‘lubridate’ 是用 R 版本 4.3.3 來建造的── Attaching core tidyverse packages ──────────────────────────────────────────────────────────────── tidyverse 2.0.0 ──
✔ dplyr 1.1.4 ✔ readr 2.1.5
✔ forcats 1.0.0 ✔ stringr 1.5.1
✔ ggplot2 3.5.0 ✔ tibble 3.2.1
✔ lubridate 1.9.3 ✔ tidyr 1.3.1
✔ purrr 1.0.2 ── Conflicts ────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag() masks stats::lag()
ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errorslibrary("dplyr")
library("foreign")
library("survey")警告： 套件 ‘survey’ 是用 R 版本 4.3.3 來建造的載入需要的套件：grid
載入需要的套件：Matrix
載入套件：‘Matrix’
下列物件被遮斷自 ‘package:tidyr’:
expand, pack, unpack
載入需要的套件：survival
載入套件：‘survival’
下列物件被遮斷由 ‘.GlobalEnv’:
cancer
載入套件：‘survey’
下列物件被遮斷自 ‘package:graphics’:
dotchartlibrary("ggplot2")
library("car")警告： 套件 ‘car’ 是用 R 版本 4.3.3 來建造的載入需要的套件：carData
警告： 套件 ‘carData’ 是用 R 版本 4.3.3 來建造的
載入套件：‘car’
下列物件被遮斷自 ‘package:dplyr’:
recode
下列物件被遮斷自 ‘package:purrr’:
somelibrary("rms")警告： 套件 ‘rms’ 是用 R 版本 4.3.3 來建造的載入需要的套件：Hmisc
警告： 套件 ‘Hmisc’ 是用 R 版本 4.3.3 來建造的Registered S3 method overwritten by 'htmlwidgets':
method from
print.htmlwidget tools:rstudio
載入套件：‘Hmisc’
下列物件被遮斷自 ‘package:survey’:
deff
下列物件被遮斷自 ‘package:dplyr’:
src, summarize
下列物件被遮斷自 ‘package:base’:
format.pval, units
警告： "replValueSp" 類別的子類別 "ndiMatrix" 沒有定義；因此沒有更新
載入套件：‘rms’
下列物件被遮斷自 ‘package:car’:
Predict, vif
下列物件被遮斷自 ‘package:survey’:
calibratelibrary("SciViews")警告： 套件 ‘SciViews’ 是用 R 版本 4.3.3 來建造的#list variable
Fulldat_mediation_pfas <- Fulldat_mediation_pfas %>% mutate(accelerated_age = Phenotypic_Age-chronological_age)
colnames(Fulldat_mediation_pfas) [1] "SEQN" "chronological_age"
[3] "Gender" "Race"
[5] "Pregnancy" "Marital_Status"
[7] "Ratio_income_poverty" "Interview_Weight"
[9] "MEC_Weight" "psu"
[11] "Strata" "BMI"
[13] "Vitamin_A" "Vitamin_C"
[15] "Vitamin_E" "Zinc"
[17] "Selenium" "sleep_disorders"
[19] "Smoked_days" "now_smoke"
[21] "quit_smoking" "Avg_alcohol_drinks"
[23] "equipment_walk" "walk_difficulty"
[25] "had_cancer" "weight_2"
[27] "Perfluorohexane_sulfonic_acid" "Perfluorohexane_sulfonic_acid_comment"
[29] "Perfluorononanoic_acid" "Perfluorononanoic_acid_comment"
[31] "perfluorooctanoic_acid" "perfluorooctanoic_acid_comment"
[33] "perfluorooctane_sulfonic_acid" "perfluorooctane_sulfonic_acid_comment"
[35] "White_blood_cell_count" "Lymphocyte_percent"
[37] "Red_cell_distribution_width" "Mean_cell_volume"
[39] "Albumin" "Creatinine"
[41] "Glucose_serum" "Alkaline_phosphotase"
[43] "xb" "Phenotypic_Age"
[45] "cate_age" "age_binary"
[47] "BMI_cat" "income_cat"
[49] "triglycerides" "fastglucose"
[51] "TriGlu_BMI" "pfas_comment"
[53] "accelerated_age" Fulldat_Pheno <- Fulldat_mediation_pfas#Examine the pfas and Phenotypic_Age
ggplot(Fulldat_Pheno, aes(x = Perfluorohexane_sulfonic_acid)) +
geom_histogram(binwidth = 100, color = "skyblue", fill = "red", alpha = 0.7) +
labs(title = "Distribution of Perfluorohexane_sulfonic_acid",
x = "Perfluorohexane_sulfonic_acid",
y = "Frequency") +
theme_minimal()
ggplot(Fulldat_Pheno, aes(x = Perfluorohexane_sulfonic_acid, y = Phenotypic_Age)) +
geom_point() +
geom_smooth(method = "lm", se = FALSE) +
labs(x = "Perfluorohexane_sulfonic_acid", y = "Phenotypic_Age", title = "Scatter Plot Perfluorohexane_sulfonic_acid vs Phenotypic_Age with Regression Line")
# Define breaks for age groups
breaks <- c(20, 45, 65, Inf) # Breaks represent the age boundaries
# Define labels for the age groups
labels <- c("1", "2", "3")
# Categorize chronological_age into groups and assign custom labels
Fulldat_Pheno$cate_age <- cut(Fulldat_Pheno$chronological_age, breaks = breaks, labels = labels, include.lowest = TRUE)#sample density curves of pfas concentrations among accelerated and delayed age
Fulldat_Pheno <- Fulldat_Pheno %>% mutate(age_binary = case_when(
Phenotypic_Age-chronological_age >= 0 ~"accelerated",
Phenotypic_Age-chronological_age < 0 ~"delayed"
))
library(ggplot2)
cols <- c("#F76D5E", "#72D8FF")
#Perfluorohexane_sulfonic_acid
ggplot(Fulldat_Pheno, aes(x = Perfluorohexane_sulfonic_acid, colour = age_binary)) +
geom_density(lwd = 1.2, linetype = 1) +
scale_color_manual(values = cols)
#Perfluorohexane_sulfonic_acid_comment
ggplot(Fulldat_Pheno, aes(x = Perfluorononanoic_acid, colour = age_binary)) +
geom_density(lwd = 1.2, linetype = 1) +
scale_color_manual(values = cols)
#Perfluorononanoic_acid
ggplot(Fulldat_Pheno, aes(x = perfluorooctanoic_acid, colour = age_binary)) +
geom_density(lwd = 1.2, linetype = 1) +
scale_color_manual(values = cols)
#Perfluorononanoic_acid_comment
ggplot(Fulldat_Pheno, aes(x = perfluorooctane_sulfonic_acid, colour = age_binary)) +
geom_density(lwd = 1.2, linetype = 1) +
scale_color_manual(values = cols)
#Main model of regression for association, and adjust for covariates (Table 2) #Perfluorohexane_sulfonic_acid
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment,
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.04263 0.08603 -12.119 < 2e-16 ***
Perfluorohexane_sulfonic_acid_comment 8.44505 1.23200 6.855 2.42e-10 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 46.40257)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.212798 -0.87247
Perfluorohexane_sulfonic_acid_comment 6.008196 10.88191
Call:
svyglm(formula = accelerated_age ~ Perfluorohexane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.20781 0.36200 6.099 1.30e-08 ***
Perfluorohexane_sulfonic_acid_comment 8.11426 1.20702 6.723 6.09e-10 ***
Gender -1.61239 0.14535 -11.093 < 2e-16 ***
Race 0.19115 0.07252 2.636 0.009480 **
Marital_Status2 2.33524 0.28239 8.270 1.92e-13 ***
Marital_Status3 0.84089 0.27994 3.004 0.003236 **
Marital_Status4 0.95429 0.41698 2.289 0.023826 *
Marital_Status5 -0.73707 0.18788 -3.923 0.000145 ***
Marital_Status6 -0.97283 0.27703 -3.512 0.000625 ***
Marital_Status77 -1.28210 2.49494 -0.514 0.608267
Marital_Status99 2.96413 6.55118 0.452 0.651743
Marital_StatusNone -0.80414 0.34708 -2.317 0.022178 *
Ratio_income_poverty -0.54192 0.04133 -13.111 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.35492)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.49120172 2.9244182
Perfluorohexane_sulfonic_acid_comment 5.72483435 10.5036797
Gender -1.90011544 -1.3246601
Race 0.04759596 0.3346992
Marital_Status2 1.77622177 2.8942561
Marital_Status3 0.28672259 1.3950672
Marital_Status4 0.12883014 1.7797498
Marital_Status5 -1.10900900 -0.3651383
Marital_Status6 -1.52123412 -0.4244286
Marital_Status77 -6.22108130 3.6568891
Marital_Status99 -10.00458185 15.9328486
Marital_StatusNone -1.49121336 -0.1170608
Ratio_income_poverty -0.62374489 -0.4600962
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.888e+00 7.721e-01 2.445 0.016116 *
Perfluorohexane_sulfonic_acid_comment 8.646e+00 1.224e+00 7.062 1.67e-10 ***
Gender -1.660e+00 1.435e-01 -11.567 < 2e-16 ***
Race 1.808e-01 6.906e-02 2.618 0.010126 *
Marital_Status2 2.191e+00 2.813e-01 7.788 4.44e-12 ***
Marital_Status3 7.460e-01 2.665e-01 2.800 0.006064 **
Marital_Status4 7.550e-01 3.613e-01 2.090 0.039009 *
Marital_Status5 -3.671e-01 1.786e-01 -2.056 0.042208 *
Marital_Status6 -7.161e-01 2.620e-01 -2.734 0.007324 **
Marital_Status77 -3.378e+00 2.207e+00 -1.530 0.128894
Marital_Status99 3.848e+00 5.508e+00 0.699 0.486317
Marital_StatusNone 5.388e+00 4.599e+00 1.171 0.244020
Ratio_income_poverty -4.612e-01 4.301e-02 -10.722 < 2e-16 ***
BMI 2.024e-01 1.095e-02 18.486 < 2e-16 ***
sleep_disorders2 -7.687e-01 2.049e-01 -3.751 0.000285 ***
sleep_disorders7 1.927e+01 2.699e+00 7.139 1.14e-10 ***
sleep_disorders9 -6.431e-02 5.138e+00 -0.013 0.990036
sleep_disordersNone -2.027e+00 2.258e-01 -8.977 9.80e-15 ***
Smoked_days -2.113e+00 2.854e-01 -7.403 3.07e-11 ***
now_smoke -6.086e-01 1.136e-01 -5.356 4.85e-07 ***
quit_smoking 2.666e-05 1.263e-05 2.111 0.037109 *
Avg_alcohol_drinks2 5.285e-01 2.110e-01 2.505 0.013733 *
Avg_alcohol_drinks9 5.713e-01 2.813e+00 0.203 0.839455
Avg_alcohol_drinksNone 2.801e-01 2.544e-01 1.101 0.273261
had_cancer2 -1.135e+00 2.444e-01 -4.642 9.77e-06 ***
had_cancer9 2.449e+00 4.059e+00 0.603 0.547481
had_cancerNone -4.985e+00 4.595e+00 -1.085 0.280349
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.26516)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 3.570977e-01 3.418061e+00
Perfluorohexane_sulfonic_acid_comment 6.218952e+00 1.107231e+01
Gender -1.944512e+00 -1.375567e+00
Race 4.387910e-02 3.176402e-01
Marital_Status2 1.633141e+00 2.748356e+00
Marital_Status3 2.178049e-01 1.274216e+00
Marital_Status4 3.878656e-02 1.471150e+00
Marital_Status5 -7.210613e-01 -1.315757e-02
Marital_Status6 -1.235295e+00 -1.968159e-01
Marital_Status77 -7.752860e+00 9.976661e-01
Marital_Status99 -7.070234e+00 1.476598e+01
Marital_StatusNone -3.728993e+00 1.450400e+01
Ratio_income_poverty -5.464661e-01 -3.759453e-01
BMI 1.807031e-01 2.241089e-01
sleep_disorders2 -1.174950e+00 -3.625388e-01
sleep_disorders7 1.391730e+01 2.461639e+01
sleep_disorders9 -1.024806e+01 1.011944e+01
sleep_disordersNone -2.475043e+00 -1.579714e+00
Smoked_days -2.678722e+00 -1.547265e+00
now_smoke -8.338251e-01 -3.833433e-01
quit_smoking 1.623133e-06 5.169953e-05
Avg_alcohol_drinks2 1.103369e-01 9.467409e-01
Avg_alcohol_drinks9 -5.004802e+00 6.147378e+00
Avg_alcohol_drinksNone -2.241250e-01 7.844211e-01
had_cancer2 -1.619078e+00 -6.500888e-01
had_cancer9 -5.596086e+00 1.049465e+01
had_cancerNone -1.409250e+01 4.122322e+00
Call:
svyglm(formula = Phenotypic_Age ~ ln(Perfluorohexane_sulfonic_acid),
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 44.9501 0.2767 162.427 < 2e-16 ***
ln(Perfluorohexane_sulfonic_acid) 1.7040 0.2593 6.571 1.04e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 374.8053)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 44.402761 45.49753
ln(Perfluorohexane_sulfonic_acid) 1.191031 2.21689
Call:
svyglm(formula = Phenotypic_Age ~ ln(Perfluorohexane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 48.0439 0.9771 49.172 < 2e-16 ***
ln(Perfluorohexane_sulfonic_acid) 1.4648 0.2652 5.523 1.91e-07 ***
Gender -1.1838 0.4710 -2.514 0.01325 *
Race 1.0604 0.1538 6.896 2.54e-10 ***
Marital_Status2 18.6290 0.8280 22.499 < 2e-16 ***
Marital_Status3 -3.0927 0.7319 -4.226 4.62e-05 ***
Marital_Status4 -3.5611 1.0720 -3.322 0.00118 **
Marital_Status5 -18.7433 0.5347 -35.052 < 2e-16 ***
Marital_Status6 -15.1011 0.7862 -19.208 < 2e-16 ***
Marital_Status77 3.2268 6.9485 0.464 0.64320
Marital_Status99 23.0357 3.4388 6.699 6.86e-10 ***
Marital_StatusNone -31.6577 0.5538 -57.160 < 2e-16 ***
Ratio_income_poverty -0.4140 0.1353 -3.060 0.00272 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 268.7147)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 46.1096892 49.9780639
ln(Perfluorohexane_sulfonic_acid) 0.9398148 1.9897695
Gender -2.1161471 -0.2515361
Race 0.7560248 1.3648672
Marital_Status2 16.9898922 20.2681162
Marital_Status3 -4.5414770 -1.6438950
Marital_Status4 -5.6833547 -1.4389165
Marital_Status5 -19.8018521 -17.6847432
Marital_Status6 -16.6574133 -13.5448147
Marital_Status77 -10.5285297 16.9821114
Marital_Status99 16.2282848 29.8432140
Marital_StatusNone -32.7541129 -30.5613166
Ratio_income_poverty -0.6818136 -0.1462218
Call:
svyglm(formula = Phenotypic_Age ~ ln(Perfluorohexane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty + BMI +
sleep_disorders + Smoked_days + now_smoke + quit_smoking +
Avg_alcohol_drinks + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.206e+01 1.708e+00 30.470 < 2e-16 ***
ln(Perfluorohexane_sulfonic_acid) 1.451e+00 2.589e-01 5.604 1.62e-07 ***
Gender -1.356e+00 4.396e-01 -3.085 0.002585 **
Race 9.981e-01 1.442e-01 6.920 3.35e-10 ***
Marital_Status2 1.622e+01 7.728e-01 20.989 < 2e-16 ***
Marital_Status3 -2.732e+00 7.164e-01 -3.814 0.000228 ***
Marital_Status4 -2.975e+00 9.932e-01 -2.996 0.003396 **
Marital_Status5 -1.606e+01 4.949e-01 -32.445 < 2e-16 ***
Marital_Status6 -1.326e+01 7.753e-01 -17.101 < 2e-16 ***
Marital_Status77 -3.179e+00 4.984e+00 -0.638 0.524949
Marital_Status99 2.692e+01 1.635e+00 16.459 < 2e-16 ***
Marital_StatusNone 7.971e+00 4.931e+00 1.617 0.108894
Ratio_income_poverty -3.883e-01 1.299e-01 -2.990 0.003457 **
BMI 1.884e-01 2.478e-02 7.601 1.14e-11 ***
sleep_disorders2 -3.521e+00 4.906e-01 -7.178 9.44e-11 ***
sleep_disorders7 2.578e+01 5.924e+00 4.351 3.08e-05 ***
sleep_disorders9 -8.691e+00 1.111e+01 -0.782 0.435890
sleep_disordersNone -4.079e+00 7.252e-01 -5.625 1.47e-07 ***
Smoked_days 1.542e+00 6.087e-01 2.533 0.012748 *
now_smoke 2.403e+00 2.806e-01 8.563 8.39e-14 ***
quit_smoking 2.066e-04 2.881e-05 7.173 9.64e-11 ***
Avg_alcohol_drinks2 4.757e+00 4.732e-01 10.052 < 2e-16 ***
Avg_alcohol_drinks9 4.971e+00 5.943e+00 0.837 0.404681
Avg_alcohol_drinksNone -2.132e-01 5.747e-01 -0.371 0.711337
had_cancer2 -1.340e+01 5.452e-01 -24.570 < 2e-16 ***
had_cancer9 -6.369e+00 4.197e+00 -1.518 0.132030
had_cancerNone -4.885e+01 4.984e+00 -9.801 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 243.4841)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 4.866969e+01 5.544259e+01
ln(Perfluorohexane_sulfonic_acid) 9.377143e-01 1.964246e+00
Gender -2.227671e+00 -4.848634e-01
Race 7.122201e-01 1.283995e+00
Marital_Status2 1.468770e+01 1.775125e+01
Marital_Status3 -4.152113e+00 -1.312013e+00
Marital_Status4 -4.943942e+00 -1.006615e+00
Marital_Status5 -1.703658e+01 -1.507481e+01
Marital_Status6 -1.479508e+01 -1.172161e+01
Marital_Status77 -1.305766e+01 6.700152e+00
Marital_Status99 2.367587e+01 3.015924e+01
Marital_StatusNone -1.802881e+00 1.774486e+01
Ratio_income_poverty -6.457039e-01 -1.308644e-01
BMI 1.392669e-01 2.375231e-01
sleep_disorders2 -4.493391e+00 -2.548635e+00
sleep_disorders7 1.403328e+01 3.751731e+01
sleep_disorders9 -3.071748e+01 1.333621e+01
sleep_disordersNone -5.516646e+00 -2.641730e+00
Smoked_days 3.352220e-01 2.748148e+00
now_smoke 1.846845e+00 2.959428e+00
quit_smoking 1.495331e-04 2.637292e-04
Avg_alcohol_drinks2 3.818748e+00 5.694655e+00
Avg_alcohol_drinks9 -6.807792e+00 1.675051e+01
Avg_alcohol_drinksNone -1.352444e+00 9.259578e-01
had_cancer2 -1.447581e+01 -1.231452e+01
had_cancer9 -1.468743e+01 1.949527e+00
had_cancerNone -5.872693e+01 -3.896839e+01
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorohexane_sulfonic_acid),
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.54233 0.09871 -5.494 1.93e-07 ***
ln(Perfluorohexane_sulfonic_acid) -1.10924 0.09760 -11.365 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 46.18725)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -0.7375652 -0.3470871
ln(Perfluorohexane_sulfonic_acid) -1.3022918 -0.9161789
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorohexane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.67530 0.35524 10.346 < 2e-16 ***
ln(Perfluorohexane_sulfonic_acid) -1.38245 0.09977 -13.856 < 2e-16 ***
Gender -2.35649 0.15783 -14.930 < 2e-16 ***
Race 0.21842 0.07377 2.961 0.003688 **
Marital_Status2 2.62626 0.29456 8.916 5.79e-15 ***
Marital_Status3 0.81096 0.28388 2.857 0.005034 **
Marital_Status4 1.02410 0.42443 2.413 0.017314 *
Marital_Status5 -0.62672 0.18874 -3.321 0.001185 **
Marital_Status6 -0.96358 0.26234 -3.673 0.000357 ***
Marital_Status77 -0.90615 2.14377 -0.423 0.673264
Marital_Status99 4.29297 5.18784 0.828 0.409566
Marital_StatusNone -0.95876 0.35144 -2.728 0.007311 **
Ratio_income_poverty -0.49186 0.03887 -12.654 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.65106)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 2.97207302 4.3785322
ln(Perfluorohexane_sulfonic_acid) -1.57995716 -1.1849405
Gender -2.66893663 -2.0440456
Race 0.07238704 0.3644613
Marital_Status2 2.04314050 3.2093751
Marital_Status3 0.24898353 1.3729318
Marital_Status4 0.18390163 1.8643040
Marital_Status5 -1.00034352 -0.2530978
Marital_Status6 -1.48290907 -0.4442510
Marital_Status77 -5.14995975 3.3376508
Marital_Status99 -5.97688074 14.5628211
Marital_StatusNone -1.65447234 -0.2630462
Ratio_income_poverty -0.56880741 -0.4149184
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.512e+00 7.945e-01 4.421 2.35e-05 ***
ln(Perfluorohexane_sulfonic_acid) -1.315e+00 1.025e-01 -12.837 < 2e-16 ***
Gender -2.352e+00 1.596e-01 -14.732 < 2e-16 ***
Race 2.074e-01 7.036e-02 2.948 0.003922 **
Marital_Status2 2.448e+00 2.944e-01 8.315 3.00e-13 ***
Marital_Status3 7.219e-01 2.681e-01 2.693 0.008209 **
Marital_Status4 7.860e-01 3.765e-01 2.088 0.039146 *
Marital_Status5 -3.055e-01 1.834e-01 -1.666 0.098685 .
Marital_Status6 -7.554e-01 2.513e-01 -3.006 0.003292 **
Marital_Status77 -2.786e+00 2.032e+00 -1.371 0.173238
Marital_Status99 5.213e+00 4.263e+00 1.223 0.223986
Marital_StatusNone 5.275e+00 4.719e+00 1.118 0.266150
Ratio_income_poverty -4.220e-01 4.105e-02 -10.278 < 2e-16 ***
BMI 1.964e-01 1.101e-02 17.837 < 2e-16 ***
sleep_disorders2 -7.379e-01 2.053e-01 -3.595 0.000491 ***
sleep_disorders7 1.878e+01 2.551e+00 7.363 3.75e-11 ***
sleep_disorders9 1.768e+00 4.409e+00 0.401 0.689235
sleep_disordersNone -1.571e+00 2.345e-01 -6.701 9.74e-10 ***
Smoked_days -2.158e+00 2.862e-01 -7.539 1.56e-11 ***
now_smoke -5.910e-01 1.148e-01 -5.149 1.19e-06 ***
quit_smoking 2.927e-05 1.248e-05 2.346 0.020814 *
Avg_alcohol_drinks2 5.552e-01 2.204e-01 2.519 0.013227 *
Avg_alcohol_drinks9 1.380e+00 2.678e+00 0.515 0.607466
Avg_alcohol_drinksNone 1.924e-01 2.482e-01 0.775 0.439930
had_cancer2 -1.254e+00 2.409e-01 -5.207 9.26e-07 ***
had_cancer9 2.322e+00 4.068e+00 0.571 0.569302
had_cancerNone -5.308e+00 4.713e+00 -1.126 0.262595
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.8335)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.937486e+00 5.087259e+00
ln(Perfluorohexane_sulfonic_acid) -1.518618e+00 -1.112376e+00
Gender -2.668007e+00 -2.035212e+00
Race 6.794517e-02 3.468734e-01
Marital_Status2 1.864453e+00 3.031540e+00
Marital_Status3 1.905794e-01 1.253305e+00
Marital_Status4 3.985247e-02 1.532243e+00
Marital_Status5 -6.690642e-01 5.805905e-02
Marital_Status6 -1.253573e+00 -2.572940e-01
Marital_Status77 -6.813853e+00 1.242105e+00
Marital_Status99 -3.235949e+00 1.366228e+01
Marital_StatusNone -4.079502e+00 1.462978e+01
Ratio_income_poverty -5.033413e-01 -3.405911e-01
BMI 1.745965e-01 2.182525e-01
sleep_disorders2 -1.144826e+00 -3.310237e-01
sleep_disorders7 1.372457e+01 2.383564e+01
sleep_disorders9 -6.972207e+00 1.050830e+01
sleep_disordersNone -2.035914e+00 -1.106357e+00
Smoked_days -2.724964e+00 -1.590336e+00
now_smoke -8.184647e-01 -3.634588e-01
quit_smoking 4.536427e-06 5.400049e-05
Avg_alcohol_drinks2 1.183544e-01 9.920712e-01
Avg_alcohol_drinks9 -3.928747e+00 6.688296e+00
Avg_alcohol_drinksNone -2.995973e-01 6.844212e-01
had_cancer2 -1.731435e+00 -7.766140e-01
had_cancer9 -5.740774e+00 1.038456e+01
had_cancerNone -1.465048e+01 4.034654e+00#“Perfluorononanoic_acid” “Perfluorononanoic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment,
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.04695 0.08925 -11.731 < 2e-16 ***
Perfluorononanoic_acid_comment 6.52553 1.50079 4.348 2.71e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 46.57207)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.223484 -0.870423
Perfluorononanoic_acid_comment 3.557023 9.494046
Call:
svyglm(formula = accelerated_age ~ Perfluorononanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.20968 0.34819 6.346 3.92e-09 ***
Perfluorononanoic_acid_comment 6.06875 1.43479 4.230 4.55e-05 ***
Gender -1.61547 0.14754 -10.949 < 2e-16 ***
Race 0.19337 0.07216 2.680 0.008381 **
Marital_Status2 2.27204 0.28400 8.000 8.13e-13 ***
Marital_Status3 0.79886 0.26922 2.967 0.003616 **
Marital_Status4 1.03196 0.42287 2.440 0.016108 *
Marital_Status5 -0.65848 0.19174 -3.434 0.000812 ***
Marital_Status6 -0.95223 0.27868 -3.417 0.000861 ***
Marital_Status77 -1.28245 2.49409 -0.514 0.608045
Marital_Status99 2.96548 6.55601 0.452 0.651835
Marital_StatusNone -0.88175 0.34221 -2.577 0.011167 *
Ratio_income_poverty -0.54438 0.03902 -13.950 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.55117)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.52040934 2.8989519
Perfluorononanoic_acid_comment 3.22843140 8.9090611
Gender -1.90753384 -1.3233963
Race 0.05053368 0.3362113
Marital_Status2 1.70983891 2.8342356
Marital_Status3 0.26591464 1.3318025
Marital_Status4 0.19485087 1.8690673
Marital_Status5 -1.03804065 -0.2789181
Marital_Status6 -1.50389450 -0.4005599
Marital_Status77 -6.21974616 3.6548449
Marital_Status99 -10.01280087 15.9437626
Marital_StatusNone -1.55918627 -0.2043207
Ratio_income_poverty -0.62163212 -0.4671290
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.847e+00 7.592e-01 2.433 0.016620 *
Perfluorononanoic_acid_comment 5.854e+00 1.569e+00 3.732 0.000305 ***
Gender -1.665e+00 1.476e-01 -11.280 < 2e-16 ***
Race 1.824e-01 6.855e-02 2.661 0.008974 **
Marital_Status2 2.129e+00 2.842e-01 7.491 1.98e-11 ***
Marital_Status3 7.123e-01 2.588e-01 2.752 0.006943 **
Marital_Status4 8.234e-01 3.716e-01 2.216 0.028789 *
Marital_Status5 -3.001e-01 1.822e-01 -1.647 0.102371
Marital_Status6 -7.095e-01 2.634e-01 -2.694 0.008198 **
Marital_Status77 -3.357e+00 2.201e+00 -1.525 0.130080
Marital_Status99 3.885e+00 5.522e+00 0.704 0.483257
Marital_StatusNone 5.238e+00 4.557e+00 1.149 0.252897
Ratio_income_poverty -4.663e-01 4.046e-02 -11.525 < 2e-16 ***
BMI 2.010e-01 1.108e-02 18.139 < 2e-16 ***
sleep_disorders2 -7.703e-01 1.982e-01 -3.887 0.000176 ***
sleep_disorders7 1.909e+01 2.675e+00 7.136 1.16e-10 ***
sleep_disorders9 -8.500e-02 5.168e+00 -0.016 0.986907
sleep_disordersNone -1.856e+00 2.230e-01 -8.321 2.91e-13 ***
Smoked_days -2.090e+00 2.776e-01 -7.526 1.66e-11 ***
now_smoke -5.948e-01 1.129e-01 -5.269 7.07e-07 ***
quit_smoking 2.549e-05 1.286e-05 1.982 0.050035 .
Avg_alcohol_drinks2 5.720e-01 2.225e-01 2.571 0.011490 *
Avg_alcohol_drinks9 5.468e-01 2.758e+00 0.198 0.843206
Avg_alcohol_drinksNone 3.382e-01 2.504e-01 1.351 0.179622
had_cancer2 -1.126e+00 2.470e-01 -4.560 1.36e-05 ***
had_cancer9 2.433e+00 4.057e+00 0.600 0.549990
had_cancerNone -4.952e+00 4.556e+00 -1.087 0.279468
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.60545)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 3.421950e-01 3.351937e+00
Perfluorononanoic_acid_comment 2.744589e+00 8.964205e+00
Gender -1.957478e+00 -1.372336e+00
Race 4.654551e-02 3.182948e-01
Marital_Status2 1.565359e+00 2.691929e+00
Marital_Status3 1.992953e-01 1.225234e+00
Marital_Status4 8.688706e-02 1.559970e+00
Marital_Status5 -6.612210e-01 6.097658e-02
Marital_Status6 -1.231610e+00 -1.873871e-01
Marital_Status77 -7.719851e+00 1.005296e+00
Marital_Status99 -7.061334e+00 1.483136e+01
Marital_StatusNone -3.794297e+00 1.426988e+01
Ratio_income_poverty -5.465130e-01 -3.861124e-01
BMI 1.790347e-01 2.229630e-01
sleep_disorders2 -1.163154e+00 -3.774343e-01
sleep_disorders7 1.378899e+01 2.439523e+01
sleep_disorders9 -1.032822e+01 1.015822e+01
sleep_disordersNone -2.298197e+00 -1.413961e+00
Smoked_days -2.639843e+00 -1.539179e+00
now_smoke -8.185551e-01 -3.710377e-01
quit_smoking -3.935925e-09 5.098847e-05
Avg_alcohol_drinks2 1.310772e-01 1.013021e+00
Avg_alcohol_drinks9 -4.920050e+00 6.013723e+00
Avg_alcohol_drinksNone -1.581060e-01 8.344356e-01
had_cancer2 -1.615765e+00 -6.366579e-01
had_cancer9 -5.609262e+00 1.047524e+01
had_cancerNone -1.398370e+01 4.078719e+00
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorononanoic_acid),
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.30480 0.08697 -15.00 <2e-16 ***
ln(Perfluorononanoic_acid) -1.72937 0.12661 -13.66 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 45.33779)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.476821 -1.132771
ln(Perfluorononanoic_acid) -1.979805 -1.478930
Call:
svyglm(formula = accelerated_age ~ ln(Perfluorononanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.89209 0.34008 5.564 1.59e-07 ***
ln(Perfluorononanoic_acid) -1.72159 0.12478 -13.797 < 2e-16 ***
Gender -1.88323 0.14656 -12.849 < 2e-16 ***
Race 0.30261 0.06955 4.351 2.83e-05 ***
Marital_Status2 2.20962 0.29177 7.573 7.77e-12 ***
Marital_Status3 0.60765 0.27099 2.242 0.02674 *
Marital_Status4 1.06306 0.39945 2.661 0.00883 **
Marital_Status5 -0.55146 0.18556 -2.972 0.00357 **
Marital_Status6 -0.81255 0.26782 -3.034 0.00295 **
Marital_Status77 -0.61904 2.43177 -0.255 0.79949
Marital_Status99 2.65720 5.79835 0.458 0.64757
Marital_StatusNone -1.18459 0.35135 -3.372 0.00100 **
Ratio_income_poverty -0.49432 0.03819 -12.942 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.3101)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.21886250 2.5653244
ln(Perfluorononanoic_acid) -1.96860118 -1.4745813
Gender -2.17336738 -1.5930981
Race 0.16492570 0.4403040
Marital_Status2 1.63202381 2.7872158
Marital_Status3 0.07120738 1.1440962
Marital_Status4 0.27231674 1.8538131
Marital_Status5 -0.91879128 -0.1841302
Marital_Status6 -1.34272215 -0.2823706
Marital_Status77 -5.43296834 4.1948883
Marital_Status99 -8.82120453 14.1356061
Marital_StatusNone -1.88011514 -0.4890595
Ratio_income_poverty -0.56992508 -0.4187091
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.544e+00 7.345e-01 2.102 0.037875 *
ln(Perfluorononanoic_acid) -1.595e+00 1.336e-01 -11.935 < 2e-16 ***
Gender -1.924e+00 1.443e-01 -13.332 < 2e-16 ***
Race 2.843e-01 6.588e-02 4.315 3.55e-05 ***
Marital_Status2 2.068e+00 2.901e-01 7.128 1.21e-10 ***
Marital_Status3 5.542e-01 2.580e-01 2.148 0.033953 *
Marital_Status4 8.427e-01 3.601e-01 2.340 0.021108 *
Marital_Status5 -2.133e-01 1.793e-01 -1.189 0.236860
Marital_Status6 -6.011e-01 2.546e-01 -2.361 0.019995 *
Marital_Status77 -2.843e+00 2.025e+00 -1.404 0.163283
Marital_Status99 3.717e+00 4.835e+00 0.769 0.443660
Marital_StatusNone 6.012e+00 4.318e+00 1.392 0.166658
Ratio_income_poverty -4.238e-01 3.984e-02 -10.638 < 2e-16 ***
BMI 1.966e-01 1.110e-02 17.715 < 2e-16 ***
sleep_disorders2 -7.835e-01 2.017e-01 -3.885 0.000177 ***
sleep_disorders7 1.962e+01 2.501e+00 7.844 3.34e-12 ***
sleep_disorders9 -1.456e-01 5.370e+00 -0.027 0.978415
sleep_disordersNone -1.434e+00 2.346e-01 -6.110 1.61e-08 ***
Smoked_days -2.025e+00 2.762e-01 -7.334 4.34e-11 ***
now_smoke -5.549e-01 1.102e-01 -5.035 1.93e-06 ***
quit_smoking 2.666e-05 1.262e-05 2.113 0.036929 *
Avg_alcohol_drinks2 6.530e-01 2.155e-01 3.031 0.003054 **
Avg_alcohol_drinks9 3.867e-01 2.164e+00 0.179 0.858524
Avg_alcohol_drinksNone 4.608e-01 2.497e-01 1.846 0.067698 .
had_cancer2 -1.216e+00 2.396e-01 -5.075 1.62e-06 ***
had_cancer9 2.180e+00 3.991e+00 0.546 0.586009
had_cancerNone -6.514e+00 4.320e+00 -1.508 0.134511
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.60734)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 8.803074e-02 2.999665e+00
ln(Perfluorononanoic_acid) -1.859933e+00 -1.330124e+00
Gender -2.210096e+00 -1.637963e+00
Race 1.536872e-01 4.148445e-01
Marital_Status2 1.492886e+00 2.643069e+00
Marital_Status3 4.276970e-02 1.065546e+00
Marital_Status4 1.289406e-01 1.556461e+00
Marital_Status5 -5.686358e-01 1.421200e-01
Marital_Status6 -1.105731e+00 -9.655666e-02
Marital_Status77 -6.857930e+00 1.171720e+00
Marital_Status99 -5.866531e+00 1.330146e+01
Marital_StatusNone -2.546425e+00 1.457054e+01
Ratio_income_poverty -5.028132e-01 -3.448631e-01
BMI 1.746018e-01 2.185968e-01
sleep_disorders2 -1.183308e+00 -3.837513e-01
sleep_disorders7 1.465862e+01 2.457238e+01
sleep_disorders9 -1.078984e+01 1.049858e+01
sleep_disordersNone -1.898647e+00 -9.685162e-01
Smoked_days -2.572787e+00 -1.477968e+00
now_smoke -7.734086e-01 -3.364522e-01
quit_smoking 1.647164e-06 5.166534e-05
Avg_alcohol_drinks2 2.259159e-01 1.080117e+00
Avg_alcohol_drinks9 -3.902858e+00 4.676202e+00
Avg_alcohol_drinksNone -3.411327e-02 9.556832e-01
had_cancer2 -1.691079e+00 -7.411789e-01
had_cancer9 -5.730335e+00 1.009036e+01
had_cancerNone -1.507704e+01 2.049060e+00#“perfluorooctanoic_acid” “perfluorooctanoic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment,
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.46556 0.09133 -27.00 <2e-16 ***
perfluorooctanoic_acid_comment 3.46625 0.19785 17.52 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.09919)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -2.646417 -2.284706
perfluorooctanoic_acid_comment 3.074454 3.858039
Call:
svyglm(formula = accelerated_age ~ perfluorooctanoic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.93744 0.34435 2.722 0.00757 **
perfluorooctanoic_acid_comment 3.28329 0.19443 16.887 < 2e-16 ***
Gender -1.62598 0.14847 -10.952 < 2e-16 ***
Race 0.15170 0.07711 1.967 0.05173 .
Marital_Status2 1.94547 0.29158 6.672 1.15e-09 ***
Marital_Status3 0.62453 0.26706 2.339 0.02122 *
Marital_Status4 1.13962 0.40884 2.787 0.00629 **
Marital_Status5 -0.07484 0.17676 -0.423 0.67284
Marital_Status6 -0.43036 0.28660 -1.502 0.13615
Marital_Status77 0.04600 3.19895 0.014 0.98855
Marital_Status99 1.15275 6.56907 0.175 0.86103
Marital_StatusNone -1.11365 0.36164 -3.079 0.00264 **
Ratio_income_poverty -0.54875 0.03697 -14.843 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.22419)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 0.254806565 1.6200638
perfluorooctanoic_acid_comment 2.897854995 3.6687201
Gender -1.920296974 -1.3316626
Race -0.001155511 0.3045598
Marital_Status2 1.367437737 2.5234993
Marital_Status3 0.095123030 1.1539361
Marital_Status4 0.329139647 1.9500997
Marital_Status5 -0.425244397 0.2755602
Marital_Status6 -0.998515141 0.1378024
Marital_Status77 -6.295536468 6.3875450
Marital_Status99 -11.869665809 14.1751597
Marital_StatusNone -1.830553052 -0.3967509
Ratio_income_poverty -0.622042978 -0.4754617
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.044e+00 7.532e-01 1.386 0.16900
perfluorooctanoic_acid_comment 3.415e+00 2.063e-01 16.549 < 2e-16 ***
Gender -1.659e+00 1.480e-01 -11.205 < 2e-16 ***
Race 1.361e-01 7.256e-02 1.876 0.06385 .
Marital_Status2 1.811e+00 2.871e-01 6.309 9.39e-09 ***
Marital_Status3 5.770e-01 2.543e-01 2.268 0.02562 *
Marital_Status4 9.094e-01 3.788e-01 2.401 0.01835 *
Marital_Status5 1.649e-01 1.749e-01 0.943 0.34820
Marital_Status6 -4.043e-01 2.784e-01 -1.452 0.14981
Marital_Status77 -2.211e+00 2.905e+00 -0.761 0.44844
Marital_Status99 2.528e+00 5.644e+00 0.448 0.65526
Marital_StatusNone 5.439e+00 4.884e+00 1.114 0.26832
Ratio_income_poverty -4.824e-01 3.933e-02 -12.265 < 2e-16 ***
BMI 1.928e-01 1.103e-02 17.474 < 2e-16 ***
sleep_disorders2 -5.809e-01 2.155e-01 -2.696 0.00832 **
sleep_disorders7 2.143e+01 3.288e+00 6.519 3.62e-09 ***
sleep_disorders9 -3.994e+00 2.251e+00 -1.774 0.07931 .
sleep_disordersNone 8.943e-02 2.457e-01 0.364 0.71670
Smoked_days -2.412e+00 2.818e-01 -8.558 2.29e-13 ***
now_smoke -6.400e-01 1.118e-01 -5.725 1.26e-07 ***
quit_smoking 8.107e-06 1.220e-05 0.665 0.50785
Avg_alcohol_drinks2 6.463e-01 2.071e-01 3.120 0.00241 **
Avg_alcohol_drinks9 9.679e-02 3.806e+00 0.025 0.97976
Avg_alcohol_drinksNone 6.053e-01 2.612e-01 2.317 0.02267 *
had_cancer2 -1.413e+00 2.525e-01 -5.598 2.18e-07 ***
had_cancer9 2.668e+00 3.767e+00 0.708 0.48064
had_cancerNone -6.779e+00 4.879e+00 -1.389 0.16809
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.69223)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -4.516569e-01 2.539940e+00
perfluorooctanoic_acid_comment 3.004927e+00 3.824414e+00
Gender -1.952646e+00 -1.364704e+00
Race -7.994957e-03 2.801766e-01
Marital_Status2 1.241241e+00 2.381464e+00
Marital_Status3 7.187530e-02 1.082027e+00
Marital_Status4 1.572089e-01 1.661497e+00
Marital_Status5 -1.824418e-01 5.122987e-01
Marital_Status6 -9.570645e-01 1.485455e-01
Marital_Status77 -7.980134e+00 3.557329e+00
Marital_Status99 -8.679645e+00 1.373558e+01
Marital_StatusNone -4.259931e+00 1.513771e+01
Ratio_income_poverty -5.604703e-01 -4.042697e-01
BMI 1.708574e-01 2.146699e-01
sleep_disorders2 -1.008790e+00 -1.530943e-01
sleep_disorders7 1.490495e+01 2.796366e+01
sleep_disorders9 -8.464326e+00 4.764921e-01
sleep_disordersNone -3.984812e-01 5.773391e-01
Smoked_days -2.971341e+00 -1.852063e+00
now_smoke -8.619358e-01 -4.179898e-01
quit_smoking -1.611082e-05 3.232492e-05
Avg_alcohol_drinks2 2.349557e-01 1.057660e+00
Avg_alcohol_drinks9 -7.460367e+00 7.653940e+00
Avg_alcohol_drinksNone 8.663384e-02 1.123970e+00
had_cancer2 -1.914780e+00 -9.120281e-01
had_cancer9 -4.813431e+00 1.014890e+01
had_cancerNone -1.646800e+01 2.910993e+00
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctanoic_acid),
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.6548 0.1554 4.213 4.95e-05 ***
ln(perfluorooctanoic_acid) -2.0976 0.1195 -17.547 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 42.06655)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 0.3470523 0.9625496
ln(perfluorooctanoic_acid) -2.3343112 -1.8608543
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctanoic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.42308 0.37398 11.827 < 2e-16 ***
ln(perfluorooctanoic_acid) -2.17327 0.12684 -17.134 < 2e-16 ***
Gender -2.26057 0.14901 -15.170 < 2e-16 ***
Race 0.22911 0.08026 2.855 0.005176 **
Marital_Status2 2.59830 0.30741 8.452 1.57e-13 ***
Marital_Status3 0.77380 0.26592 2.910 0.004398 **
Marital_Status4 1.01117 0.43815 2.308 0.022932 *
Marital_Status5 -0.25463 0.18405 -1.383 0.169395
Marital_Status6 -0.64657 0.27559 -2.346 0.020813 *
Marital_Status77 0.07015 2.69445 0.026 0.979278
Marital_Status99 2.78997 5.77125 0.483 0.629781
Marital_StatusNone -1.42827 0.38964 -3.666 0.000386 ***
Ratio_income_poverty -0.42590 0.03892 -10.943 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.77172)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 3.68171644 5.1644409
ln(perfluorooctanoic_acid) -2.42471566 -1.9218226
Gender -2.55596792 -1.9651682
Race 0.07000069 0.3882121
Marital_Status2 1.98889063 3.2077171
Marital_Status3 0.24664795 1.3009482
Marital_Status4 0.14259707 1.8797487
Marital_Status5 -0.61948011 0.1102255
Marital_Status6 -1.19289872 -0.1002373
Marital_Status77 -5.27128843 5.4115875
Marital_Status99 -8.65085286 14.2307875
Marital_StatusNone -2.20069286 -0.6558561
Ratio_income_poverty -0.50305690 -0.3487423
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.662e+00 7.654e-01 6.091 2.50e-08 ***
ln(perfluorooctanoic_acid) -2.193e+00 1.327e-01 -16.524 < 2e-16 ***
Gender -2.256e+00 1.475e-01 -15.288 < 2e-16 ***
Race 2.275e-01 7.489e-02 3.038 0.00309 **
Marital_Status2 2.441e+00 3.024e-01 8.073 2.38e-12 ***
Marital_Status3 7.127e-01 2.502e-01 2.848 0.00542 **
Marital_Status4 7.810e-01 4.073e-01 1.918 0.05823 .
Marital_Status5 -1.152e-02 1.820e-01 -0.063 0.94965
Marital_Status6 -6.252e-01 2.714e-01 -2.303 0.02350 *
Marital_Status77 -2.211e+00 2.471e+00 -0.895 0.37321
Marital_Status99 4.049e+00 4.830e+00 0.838 0.40399
Marital_StatusNone 5.617e+00 4.445e+00 1.264 0.20954
Ratio_income_poverty -3.718e-01 3.843e-02 -9.675 9.95e-16 ***
BMI 1.932e-01 1.133e-02 17.047 < 2e-16 ***
sleep_disorders2 -6.818e-01 2.048e-01 -3.328 0.00125 **
sleep_disorders7 2.123e+01 2.834e+00 7.493 3.83e-11 ***
sleep_disorders9 -4.311e+00 3.089e+00 -1.396 0.16612
sleep_disordersNone -1.577e-01 2.291e-01 -0.688 0.49290
Smoked_days -2.285e+00 2.767e-01 -8.255 9.90e-13 ***
now_smoke -6.115e-01 1.132e-01 -5.403 5.01e-07 ***
quit_smoking 2.053e-05 1.408e-05 1.459 0.14805
Avg_alcohol_drinks2 3.057e-01 2.023e-01 1.511 0.13416
Avg_alcohol_drinks9 1.945e-01 3.033e+00 0.064 0.94900
Avg_alcohol_drinksNone 2.749e-01 2.572e-01 1.069 0.28801
had_cancer2 -1.602e+00 2.667e-01 -6.008 3.63e-08 ***
had_cancer9 1.481e+00 3.713e+00 0.399 0.69091
had_cancerNone -6.662e+00 4.454e+00 -1.496 0.13811
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 38.20793)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 3.142374e+00 6.182307e+00
ln(perfluorooctanoic_acid) -2.456849e+00 -1.929703e+00
Gender -2.548568e+00 -1.962583e+00
Race 7.883374e-02 3.762655e-01
Marital_Status2 1.840656e+00 3.041524e+00
Marital_Status3 2.157259e-01 1.209609e+00
Marital_Status4 -2.778026e-02 1.589774e+00
Marital_Status5 -3.729075e-01 3.498625e-01
Marital_Status6 -1.164247e+00 -8.615870e-02
Marital_Status77 -7.118510e+00 2.696099e+00
Marital_Status99 -5.541735e+00 1.363941e+01
Marital_StatusNone -3.210232e+00 1.444323e+01
Ratio_income_poverty -4.481541e-01 -2.955178e-01
BMI 1.706837e-01 2.156917e-01
sleep_disorders2 -1.088544e+00 -2.750142e-01
sleep_disorders7 1.560637e+01 2.686175e+01
sleep_disorders9 -1.044391e+01 1.822450e+00
sleep_disordersNone -6.125660e-01 2.971809e-01
Smoked_days -2.834181e+00 -1.735046e+00
now_smoke -8.362689e-01 -3.867898e-01
quit_smoking -7.422054e-06 4.848869e-05
Avg_alcohol_drinks2 -9.603065e-02 7.073813e-01
Avg_alcohol_drinks9 -5.828217e+00 6.217271e+00
Avg_alcohol_drinksNone -2.359365e-01 7.857150e-01
had_cancer2 -2.131670e+00 -1.072537e+00
had_cancer9 -5.892805e+00 8.854957e+00
had_cancerNone -1.550660e+01 2.182708e+00#“perfluorooctane_sulfonic_acid” “perfluorooctane_sulfonic_acid_comment”
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment,
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.31799 0.09158 -14.391 < 2e-16 ***
perfluorooctane_sulfonic_acid_comment 11.58481 2.27727 5.087 1.39e-06 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 44.1984)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) -1.499348 -1.136631
perfluorooctane_sulfonic_acid_comment 7.075198 16.094419
Call:
svyglm(formula = accelerated_age ~ perfluorooctane_sulfonic_acid_comment +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.14660 0.38529 5.571 1.90e-07 ***
perfluorooctane_sulfonic_acid_comment 11.31755 2.31384 4.891 3.55e-06 ***
Gender -1.71136 0.14561 -11.753 < 2e-16 ***
Race 0.18051 0.08064 2.238 0.027264 *
Marital_Status2 2.46913 0.30262 8.159 7.05e-13 ***
Marital_Status3 0.87989 0.27333 3.219 0.001703 **
Marital_Status4 0.72153 0.44874 1.608 0.110806
Marital_Status5 -0.71565 0.19104 -3.746 0.000291 ***
Marital_Status6 -1.10517 0.28924 -3.821 0.000223 ***
Marital_Status77 -0.38439 2.86905 -0.134 0.893672
Marital_Status99 3.29881 6.57580 0.502 0.616938
Marital_StatusNone -1.05698 0.39697 -2.663 0.008950 **
Ratio_income_poverty -0.55551 0.04155 -13.370 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 41.95277)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.38281156 2.9103967
perfluorooctane_sulfonic_acid_comment 6.73062767 15.9044663
Gender -2.00002565 -1.4227022
Race 0.02064643 0.3403654
Marital_Status2 1.86920800 3.0690437
Marital_Status3 0.33803957 1.4217340
Marital_Status4 -0.16805174 1.6111148
Marital_Status5 -1.09436890 -0.3369401
Marital_Status6 -1.67854691 -0.5317949
Marital_Status77 -6.07195319 5.3031785
Marital_Status99 -9.73695632 16.3345719
Marital_StatusNone -1.84392449 -0.2700392
Ratio_income_poverty -0.63788084 -0.4731456
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.092e+00 8.023e-01 2.607 0.010626 *
perfluorooctane_sulfonic_acid_comment 1.102e+01 2.512e+00 4.388 3.02e-05 ***
Gender -1.693e+00 1.426e-01 -11.876 < 2e-16 ***
Race 1.729e-01 7.674e-02 2.254 0.026573 *
Marital_Status2 2.284e+00 2.970e-01 7.691 1.49e-11 ***
Marital_Status3 7.847e-01 2.590e-01 3.030 0.003163 **
Marital_Status4 5.606e-01 3.985e-01 1.407 0.162849
Marital_Status5 -4.073e-01 1.874e-01 -2.173 0.032299 *
Marital_Status6 -9.643e-01 2.830e-01 -3.407 0.000971 ***
Marital_Status77 -2.581e+00 2.730e+00 -0.945 0.346887
Marital_Status99 4.200e+00 5.602e+00 0.750 0.455291
Marital_StatusNone 5.222e+00 4.722e+00 1.106 0.271574
Ratio_income_poverty -4.969e-01 4.266e-02 -11.650 < 2e-16 ***
BMI 2.025e-01 1.148e-02 17.648 < 2e-16 ***
sleep_disorders2 -6.891e-01 2.156e-01 -3.196 0.001905 **
sleep_disorders7 1.927e+01 3.124e+00 6.170 1.76e-08 ***
sleep_disorders9 -3.652e+00 3.283e+00 -1.112 0.268832
sleep_disordersNone -1.376e+00 2.364e-01 -5.822 8.23e-08 ***
Smoked_days -2.202e+00 2.868e-01 -7.678 1.59e-11 ***
now_smoke -6.353e-01 1.149e-01 -5.529 2.94e-07 ***
quit_smoking 2.547e-05 1.477e-05 1.724 0.087950 .
Avg_alcohol_drinks2 1.902e-01 2.108e-01 0.902 0.369375
Avg_alcohol_drinks9 3.701e-01 3.557e+00 0.104 0.917360
Avg_alcohol_drinksNone 2.651e-01 2.625e-01 1.010 0.315255
had_cancer2 -1.373e+00 2.708e-01 -5.070 2.02e-06 ***
had_cancer9 2.341e+00 4.009e+00 0.584 0.560631
had_cancerNone -4.923e+00 4.724e+00 -1.042 0.300003
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.19745)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 4.987494e-01 3.685130e+00
perfluorooctane_sulfonic_acid_comment 6.033605e+00 1.600969e+01
Gender -1.976313e+00 -1.410066e+00
Race 2.054718e-02 3.253173e-01
Marital_Status2 1.694320e+00 2.873844e+00
Marital_Status3 2.704974e-01 1.298969e+00
Marital_Status4 -2.307771e-01 1.351941e+00
Marital_Status5 -7.793978e-01 -3.513537e-02
Marital_Status6 -1.526353e+00 -4.023267e-01
Marital_Status77 -8.001532e+00 2.839942e+00
Marital_Status99 -6.924181e+00 1.532437e+01
Marital_StatusNone -4.154254e+00 1.459913e+01
Ratio_income_poverty -5.816557e-01 -4.122371e-01
BMI 1.797465e-01 2.253265e-01
sleep_disorders2 -1.117305e+00 -2.609026e-01
sleep_disorders7 1.307101e+01 2.547724e+01
sleep_disorders9 -1.017042e+01 2.867086e+00
sleep_disordersNone -1.845662e+00 -9.068490e-01
Smoked_days -2.771831e+00 -1.632634e+00
now_smoke -8.634520e-01 -4.070888e-01
quit_smoking -3.859547e-06 5.478984e-05
Avg_alcohol_drinks2 -2.284839e-01 6.088154e-01
Avg_alcohol_drinks9 -6.693937e+00 7.434154e+00
Avg_alcohol_drinksNone -2.562396e-01 7.863673e-01
had_cancer2 -1.910510e+00 -8.350522e-01
had_cancer9 -5.620018e+00 1.030274e+01
had_cancerNone -1.430305e+01 4.456955e+00
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctane_sulfonic_acid),
design = des, family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.50630 0.24622 6.118 1.27e-08 ***
ln(perfluorooctane_sulfonic_acid) -1.28738 0.09427 -13.656 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 43.22353)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 1.018714 1.993888
ln(perfluorooctane_sulfonic_acid) -1.474062 -1.100699
Call:
svyglm(formula = accelerated_age ~ ln(perfluorooctane_sulfonic_acid) +
Gender + Race + Marital_Status + Ratio_income_poverty, design = des,
family = "gaussian", data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.77446 0.44804 12.888 < 2e-16 ***
ln(perfluorooctane_sulfonic_acid) -1.46710 0.10211 -14.368 < 2e-16 ***
Gender -2.32962 0.15876 -14.674 < 2e-16 ***
Race 0.30484 0.08224 3.707 0.000334 ***
Marital_Status2 2.59025 0.31888 8.123 8.48e-13 ***
Marital_Status3 0.60744 0.27059 2.245 0.026833 *
Marital_Status4 0.77685 0.42804 1.815 0.072342 .
Marital_Status5 -0.63695 0.18719 -3.403 0.000940 ***
Marital_Status6 -0.92620 0.28173 -3.287 0.001368 **
Marital_Status77 0.10791 2.66148 0.041 0.967735
Marital_Status99 3.45240 5.71447 0.604 0.547023
Marital_StatusNone -1.82759 0.40028 -4.566 1.34e-05 ***
Ratio_income_poverty -0.49880 0.03967 -12.575 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 40.62025)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 4.88627186 6.6626561
ln(perfluorooctane_sulfonic_acid) -1.66951751 -1.2646882
Gender -2.64434362 -2.0148883
Race 0.14181745 0.4678717
Marital_Status2 1.95810182 3.2223927
Marital_Status3 0.07102802 1.1438599
Marital_Status4 -0.07169366 1.6253956
Marital_Status5 -1.00802989 -0.2658703
Marital_Status6 -1.48470254 -0.3676907
Marital_Status77 -5.16816805 5.3839830
Marital_Status99 -7.87587330 14.7806653
Marital_StatusNone -2.62108909 -1.0340846
Ratio_income_poverty -0.57743098 -0.4201637
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 + had_cancer, design = des, family = "gaussian",
data = Fulldat_Pheno)
Survey design:
svydesign(id = ~psu, strata = ~Strata, weights = ~weight_2, nest = TRUE,
data = Fulldat_Pheno)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.812e+00 8.193e-01 7.094 2.52e-10 ***
ln(perfluorooctane_sulfonic_acid) -1.571e+00 1.195e-01 -13.149 < 2e-16 ***
Gender -2.398e+00 1.592e-01 -15.061 < 2e-16 ***
Race 3.148e-01 7.794e-02 4.039 0.00011 ***
Marital_Status2 2.453e+00 3.152e-01 7.781 9.68e-12 ***
Marital_Status3 5.583e-01 2.532e-01 2.205 0.02993 *
Marital_Status4 5.467e-01 3.972e-01 1.376 0.17206
Marital_Status5 -4.212e-01 1.842e-01 -2.287 0.02446 *
Marital_Status6 -9.005e-01 2.779e-01 -3.240 0.00166 **
Marital_Status77 -1.964e+00 2.402e+00 -0.818 0.41572
Marital_Status99 4.826e+00 4.742e+00 1.018 0.31146
Marital_StatusNone 5.317e+00 4.643e+00 1.145 0.25508
Ratio_income_poverty -4.461e-01 3.955e-02 -11.278 < 2e-16 ***
BMI 1.946e-01 1.135e-02 17.151 < 2e-16 ***
sleep_disorders2 -6.407e-01 2.184e-01 -2.933 0.00422 **
sleep_disorders7 1.966e+01 2.794e+00 7.036 3.30e-10 ***
sleep_disorders9 -2.962e+00 2.241e+00 -1.322 0.18937
sleep_disordersNone 2.211e-01 2.582e-01 0.856 0.39400
Smoked_days -2.096e+00 2.789e-01 -7.517 3.42e-11 ***
now_smoke -5.742e-01 1.123e-01 -5.115 1.68e-06 ***
quit_smoking 3.080e-05 1.532e-05 2.011 0.04726 *
Avg_alcohol_drinks2 4.240e-01 2.031e-01 2.088 0.03956 *
Avg_alcohol_drinks9 1.180e+00 3.227e+00 0.366 0.71547
Avg_alcohol_drinksNone 2.560e-01 2.609e-01 0.981 0.32896
had_cancer2 -1.630e+00 2.722e-01 -5.988 3.97e-08 ***
had_cancer9 2.192e+00 3.712e+00 0.591 0.55619
had_cancerNone -6.551e+00 4.655e+00 -1.407 0.16268
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 39.04171)
Number of Fisher Scoring iterations: 2
2.5 % 97.5 %
(Intercept) 4.185095e+00 7.439105e+00
ln(perfluorooctane_sulfonic_acid) -1.808281e+00 -1.333762e+00
Gender -2.714681e+00 -2.082212e+00
Race 1.600435e-01 4.696086e-01
Marital_Status2 1.826876e+00 3.078799e+00
Marital_Status3 5.548551e-02 1.061205e+00
Marital_Status4 -2.421512e-01 1.335516e+00
Marital_Status5 -7.869896e-01 -5.547473e-02
Marital_Status6 -1.452394e+00 -3.486225e-01
Marital_Status77 -6.734951e+00 2.806765e+00
Marital_Status99 -4.590525e+00 1.424181e+01
Marital_StatusNone -3.903329e+00 1.453774e+01
Ratio_income_poverty -5.246094e-01 -3.675216e-01
BMI 1.721013e-01 2.171733e-01
sleep_disorders2 -1.074396e+00 -2.069414e-01
sleep_disorders7 1.411189e+01 2.520927e+01
sleep_disorders9 -7.411818e+00 1.487086e+00
sleep_disordersNone -2.916103e-01 7.338243e-01
Smoked_days -2.650225e+00 -1.542578e+00
now_smoke -7.970607e-01 -3.512448e-01
quit_smoking 3.804027e-07 6.121391e-05
Avg_alcohol_drinks2 2.070006e-02 8.273155e-01
Avg_alcohol_drinks9 -5.228885e+00 7.588956e+00
Avg_alcohol_drinksNone -2.620109e-01 7.739953e-01
had_cancer2 -2.170744e+00 -1.089482e+00
had_cancer9 -5.178775e+00 9.563720e+00
had_cancerNone -1.579415e+01 2.693010e+00#Perfluorohexane_sulfonic_acid
library("rcssci")警告： 套件 ‘rcssci’ 是用 R 版本 4.3.3 來建造的Fulldat_Pheno <- Fulldat_Pheno %>% mutate(Ln_Perfluorohexane_sulfonic_acid = ln(Perfluorohexane_sulfonic_acid))
rcssci_linear(data = Fulldat_Pheno, y = "accelerated_age", x = "Ln_Perfluorohexane_sulfonic_acid", covs=c("Gender", "Race", "BMI","had_cancer"), prob = 0.1,ref.zero=FALSE,
filepath = "C:/Users/HKUSCM/Documents")$aics
d.f. d.f. d.f. d.f. d.f.
100770.4 100746.5 100739.1 100739.8 100735.5
$kn
[1] 7
$Q20
0% 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55%
-2.6592600 -1.2039728 -0.9162907 -0.6931472 -0.3566749 -0.2231436 -0.1053605 0.0000000 0.1823216 0.2623643 0.3435897 0.4700036
60% 65% 70% 75% 80% 85% 90% 95% 100%
0.5877867 0.6931472 0.8241754 0.9321641 1.0647107 1.2527630 1.4586150 1.8245493 4.4067192
$lshapcicrossNA#Perfluorononanoic_acid
library("rcssci")
Fulldat_Pheno <- Fulldat_Pheno %>% mutate(Ln_Perfluorononanoic_acid = ln(Perfluorononanoic_acid))
rcssci_linear(data = Fulldat_Pheno, y = "accelerated_age", x = "Perfluorononanoic_acid", covs=c("Gender", "Race", "BMI","had_cancer"), prob = 0.1,ref.zero=FALSE,
filepath = "C:/Users/HKUSCM/Documents")$aics
d.f. d.f. d.f. d.f. d.f.
100818.5 100754.6 100748.2 100745.0 100739.4
$kn
[1] 7
$Q20
0% 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70% 75% 80% 85% 90% 95%
0.058 0.200 0.300 0.400 0.492 0.500 0.600 0.656 0.700 0.800 0.890 0.902 1.000 1.100 1.230 1.394 1.500 1.722 2.100 2.800
100%
80.770
$lshapcicrossNA#perfluorooctanoic_acid
library("rcssci")
Fulldat_Pheno <- Fulldat_Pheno %>% mutate(Ln_perfluorooctanoic_acid = ln(perfluorooctanoic_acid))
rcssci_linear(data = Fulldat_Pheno, y = "accelerated_age", x = "perfluorooctanoic_acid", covs=c("Gender", "Race", "BMI","had_cancer"), prob = 0.1,ref.zero=FALSE,
filepath = "C:/Users/HKUSCM/Documents")$aics
d.f. d.f. d.f. d.f. d.f.
87276.61 87179.36 87119.81 87098.44 87084.35
$kn
[1] 7
$Q20
0% 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70%
0.0700 0.6035 0.8700 1.0700 1.2700 1.4700 1.6700 1.8700 2.0800 2.3000 2.5000 2.8000 3.0700 3.4000 3.7990
75% 80% 85% 90% 95% 100%
4.2000 4.7000 5.3700 6.3000 7.9000 104.0000
$lshapcicrossNA#perfluorooctane_sulfonic_acid
library("rcssci")
Fulldat_Pheno <- Fulldat_Pheno %>% mutate(Ln_perfluorooctane_sulfonic_acid = ln(perfluorooctane_sulfonic_acid))
rcssci_linear(data = Fulldat_Pheno, y = "accelerated_age", x = "perfluorooctane_sulfonic_acid", covs=c("Gender", "Race", "BMI","had_cancer"), prob = 0.1,ref.zero=FALSE,
filepath = "C:/Users/HKUSCM/Documents")$aics
d.f. d.f. d.f. d.f. d.f.
87679.27 87602.10 87564.83 87549.13 87535.16
$kn
[1] 7
$Q20
0% 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70% 75% 80% 85% 90% 95%
0.14 1.60 2.40 3.10 3.86 4.70 5.43 6.30 7.30 8.30 9.40 10.60 12.00 13.60 15.70 17.80 20.70 24.20 29.60 39.70
100%
435.00
$lshapcicrossNA