作者aee36900 (持久战!!)
看板R_Language
标题[问题] plot 3D图如何标示最大值
时间Wed Sep 2 09:59:15 2015
[问题类型]:
程式谘询(我想用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), 来自: 211.76.246.58
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1441159158.A.BB8.html
1F：→ andrew43: 我会先布下网格直接找，至於精确做法我不知道了。 09/02 18:23
※ 编辑: aee36900 (114.32.117.19), 09/02/2015 20:12:51
2F：→ aee36900: 好建议,但还是不知道怎麽做... 09/02 20:13