各位前辈,
现在我想对一系列散点进行GAM平滑,然后计算其平滑函数的一次和二次导数。论文中有介绍说是用(mgcv)安装包进行gam平滑,但是没有提如何求该函数导数。。
我是R小白一个,初次接触,实在不懂。。。
1 尝试过这样的方法:
a=read.csv("21Z.csv")## csv其实就是21组数据,自变量T,因变量V,基本上V=T
b=gam(a$V~s(a$T))## 这里显示有问题,其实就是b=gam(V~s(T))这个意思。。。
c=deriv(~ b,"T")
结果c(1)=1;c(2)=2;c(15)=15,求出来的结果貌似不是导数。。。
2 mgcv中有一个gam.fit3的程式,看说明貌似是求导的,但是没有例子,实在搞不定,不知道哪位前辈能讲解一二,不胜感激!!


雷达卡



京公网安备 11010802022788号







