热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

R中mgcv包画图——单位x对应y改变量图

我是看了最近发表在《EnvironmentalHealth》上的文章,采用的是R中的mgcv包,目的是探求血压和温度的非线性关系,我利用找到数据在R中学习操作有些问题。ctgam(Blood
我是看了最近发表在《Environmental Health》上的文章,采用的是R中的mgcv包,目的是探求血压和温度的非线性关系,我利用找到数据在R中学习操作有些问题。
ct=gam(Blood pressure~s(Temp,k=5)+s(BMI,k=4)+s(Age,k=6)+s(RH,k=4)+s(Pa,k=5)+Gender+Season,family=gaussian,data=mydata)                                                                                                                                                                                           (公式1)

1. 如何控制个体特征变量:性别(0,1)是分类变量,是否直接将其当成线性关系放进等式里,如公式1所示。年龄是连续性变量,用s平滑。是这样吗?我是把季节(1,2,3,4)当作分类变量,以线性关系放入公式控制的,这样作对吗?

2. 我在R中采用公式1,然后用plot(ct)画图,得到图1,这个纵坐标s(Tmean,3.87)是应变量血压吗?如果图1纵坐标为血压的话,我看文章中画出了降低单位温度血压改变值的图(图2)。想得到纵坐标为血压改变值的图是否就是求图1的导数,但是求导的程序代码一直尝试未果。这个猜想不知道对不对,希望各位大侠给出得到图2的具体代码,不甚感激。

1 个解决方案

#1


请大家帮忙解答,困扰小女好久了。

推荐阅读
author-avatar
1234
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有