R mgcv GAM画图 计算RR

请问各位大神，我现在想用R mgcv包做GAM，广义加性模型。我的数据资料是每个病例信息，以及所各自对应的pm2.5暴露，要研究这种疾病与空气暴露的关系。请问，mgcv如何做这种图，并且输出的是RR值（相对危险度）！

纵坐标是RR值，而不是s（x）！

您好  我想和您交流怎么做GAM  请您加我QQ 2027660071 谢谢

请问怎么解决？急求~

陈超123 发表于 2017-10-31 19:32
纵坐标是RR值，而不是s（x）！

你好，请问您的问题解决了吗？我也想知道如何求RR

1. library(epiR) # Used for Risk ratio
   
2. library(MASS) # Used for data
   
3. 
   
4. dat1 <- birthwt; head(dat1)
   
5. 
   
6. ## Generate a table of cell frequencies. First set the levels of the outcome
   
7. ## and the exposure so the frequencies in the 2 by 2 table come out in the
   
8. ## conventional format:
   
9. dat1$low <- factor(dat1$low, levels = c(1,0))
   
10. dat1$smoke <- factor(dat1$smoke, levels = c(1,0))
    
11. dat1$race <- factor(dat1$race, levels = c(1,2,3))
    
12. ## Generate the 2 by 2 table. Exposure (rows) = smoke. Outcome (columns) = low.
    
13. tab1 <- table(dat1$smoke, dat1$low, dnn = c("Smoke", "Low BW"))
    
14. print(tab1)
    
15. ## Compute the incidence risk ratio and other measures of association:
    
16. epi.2by2(dat = tab1, method = "cohort.count",
    
17. conf.level = 0.95, units = 100, outcome = "as.columns")

复制代码

源答案见这里：
https://stats.stackexchange.com/questions/180135/calculating-a-risk-ratio-for-specific-x-values-from-a-gam-model-using-the-mgcv-p