作者celestialgod (天)
看板R_Language
标题Re: [问题] plot 3D图如何标示最大值
时间Thu Sep 3 01:52:52 2015
如果要找最大值可以直接使用optim来做
我直接base在andrew大大的程式上写的....
findTheMax = optim(c(1,1),
function(x) -predict(mod, data.frame(x=x[1],y=x[2]), type = "response"),
lower = c(1,1), upper = c(9,9))
findTheMax$par # location of x,y
findTheMax$value # maximum value of z
※ 引述《andrew43 (讨厌有好心推文後删文者)》之铭言：
: ※ 引述《aee36900 (持久战!!)》之铭言：
: : [问题类型]:
: : 程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
: : [软体熟悉度]:
: : 新手(没写过程式，R 是我的第一次)
: : [问题叙述]:
: : 练习dlnm模型
: : 想要在3D图形上将RR最大值的数值及对应的天数找出来
: : 不知道有无function可以直接使用计算出来
: : 图形如下
: : http://imgur.com/6IcvqCW
: : [程式范例]:
: : install.packages("dlnm")
: : library(dlnm)
: : cb3.pm <- crossbasis(chicagoNMMAPS$pm10, lag=1, argvar=list(fun="lin",cen=0),
: : arglag=list(fun="strata"))
: : varknots <- equalknots(chicagoNMMAPS$temp,fun="bs",df=5,degree=2)
: : lagknots <- logknots(30, 3)
: : cb3.temp <- crossbasis(chicagoNMMAPS$temp, lag=30, argvar=list(fun="bs",knots=varknots,cen=21), arglag=list(knots=lagknots))
: : model3 <- glm(death ~ cb3.pm + cb3.temp + ns(time, 7*14) + dow,family=quasipoisson(), chicagoNMMAPS)
: : pred3.temp <- crosspred(cb3.temp, model3, by=1)
: : plot(pred3.temp, xlab="Temperature", zlab="RR", theta=200, phi=40, lphi=30,
: : main="3D graph of temperature effect")
: : [环境叙述]:
: : R version 3.1.3 (2015-03-09)
: : Platform: x86_64-redhat-linux-gnu (64-bit)
: : Running under: CentOS release 6.5 (Final)
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.160.170.238
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1441216376.A.CD2.html
1F：→ andrew43: 多谢补充了。 09/04 04:52