0

在我开始之前，可以在这里找到这个问题的基本答案： 正确绑定系数到汇总表

这个问题的不同之处在于我需要根据放置结的位置将正确的系数正确连接到汇总表中的正确位置。我使用该I(pmax(0, variable - knot))技术来放置我的样条线。最终结果是每个变量的唯一值表、汇总度量和正确的模型统计信息（请参阅下面示例代码中的最终（尚未完成）表）。 

library(tidyverse)
library(broom)

#pull in and gather data
mtcars1 <- as_tibble(mtcars)
mtcars1$cyl <- as.factor(mtcars$cyl)
#run model and produce model-summary table
model <- glm(mpg ~ cyl + hp + I(pmax(0, hp - 100)), data = mtcars1)

model_summary <- tidy(model)

#produce final summary table
summary_table <- mtcars1 %>%
  select(cyl, hp, wt) %>%
  gather(key = variable, level, - wt) %>%
  group_by(variable, level) %>%
  summarise("sum_wt" = sum(wt)) %>%
  mutate(term = paste0(variable, level)) %>%
  left_join(model_summary, by = c("term" = "term"))

挑战在于采用表中的I(pmax(0, hp -100))术语model_summary并将估计值、std.error、statistic 和 p.value 正确加入到summary_table<= 100 中的每个 hp 观察值中，此外还将其他hp估计统计量加入到 hp 观察值中summary_table即 > 100。

4

0 回答 0