<jupyter_start><jupyter_text>BMI<jupyter_code>lm_INI21002 <- lm (
    INI21002 ~ 0 + age + as.factor(sex) + PC1 + PC2 + PC3 + PC4 + as.factor(PDE3B),
    data_filtered %>% filter(INI21002 != -9 & PDE3B != -9)
)

summary(lm_INI21002)
p_INI21001 <- ggplot(
    data_filtered %>% filter(INI21001 != -9 & PDE3B != -9), 
    aes(x = as.factor(PDE3B), y = INI21001)) +
geom_boxplot(aes(group = factor(PDE3B))) +
labs(y = 'BMI (INI21001)')

p_INI21001
p_INI21001 + facet_grid(sex_label ~ .)
lm_INI21001 <- lm (
    INI21001 ~ 0 + age + as.factor(sex) + PC1 + PC2 + PC3 + PC4 + as.factor(PDE3B),
    data_filtered %>% filter(INI21001 != -9 & PDE3B != -9)
)

summary(lm_INI21001)<jupyter_output><empty_output><jupyter_text>standing height<jupyter_code>p_INI50 <- ggplot(
    data_filtered %>% filter(INI50 != -9 & PDE3B != -9), 
    aes(x = as.factor(PDE3B), y = INI50)) +
geom_boxplot(aes(group = factor(PDE3B))) +
labs(y = 'Standing height (INI50) [cm]')

p_INI50
p_INI50 + facet_grid(sex_label ~ .)
lm_INI50 <- lm (
    INI50 ~ 0 + age + as.factor(sex) + PC1 + PC2 + PC3 + PC4 + as.factor(PDE3B),
    data_filtered %>% filter(INI50 != -9 & PDE3B != -9)
)

summary(lm_INI50)<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.17627401837928153, 'fraction_numerical': 0.1060985797827903, 'mean_word_length': 4.076271186440678, 'pattern_counts': {'":': 0, '<': 14, '<?xml version=': 0, '>': 14, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '2', 'dir': 'jupyter-structured-clean-dedup', 'id': '303999', 'license': 'MIT', 'n_tokens_mistral': 638, 'n_tokens_neox': 550, 'n_words': 91, 'path': 'uk_biobank/DeGAs/notebook/regression_PDE3B.ipynb', 'repo_name': 'bschilder/public-resources', 'size': '1197'}