[問題] plot 3D圖如何標示最大值
[問題類型]:
程式諮詢(我想用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://www.ptt.cc/bbs/R_Language/M.1441159158.A.BB8.html
→
09/02 18:23, , 1F
09/02 18:23, 1F
※ 編輯: aee36900 (114.32.117.19), 09/02/2015 20:12:51
→
09/02 20:13, , 2F
09/02 20:13, 2F
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 4 篇)：
R_Language 近期熱門文章
PTT數位生活區 即時熱門文章
