Re: [問題] plot 3D圖如何標示最大值

看板R_Language作者celestialgod (天)時間10年前 (2015/09/03 01:52)推噓0(0推 0噓 1→)

留言1則, 1人參與討論串3/4 (看更多)

如果要找最大值可以直接使用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://www.ptt.cc/bbs/R_Language/M.1441216376.A.CD2.html