r画图后如何确定最值点

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问各位大牛，当我画了一个图形之后，我想找到这个图形上的极值点或者最值点，请问是否有函数实现？这个图形本身没有函数存在，我想只能通过自己去找，但是有没有简单一点的程序呢？比如这个程序产生的图形，如何确定最低点是哪个呢？

1. n <- 1000
   
2. set.seed(731)
   
3. age <- 50 + 12*rnorm(n)
   
4. label(age) <- "Age"
   
5. sex <- factor(sample(c('Male','Female'), n,
   
6.                      rep=TRUE, prob=c(.6, .4)))
   
7. cens <- 15*runif(n)
   
8. h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
   
9. dt <- -log(runif(n))/h
   
10. label(dt) <- 'Follow-up Time'
    
11. e <- ifelse(dt <= cens,1,0)
    
12. dt <- pmin(dt, cens)
    
13. units(dt) <- "Year"
    
14. dd <- datadist(age, sex)
    
15. options(datadist='dd')
    
16. S <- Surv(dt,e)
    
17. 
    
18. f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE)
    
19. cox.zph(f, "rank")             # tests of PH
    
20. anova(f)
    
21. plot(Predict(f, age, sex)) # plot age effect, 2 curves for 2 sexes

复制代码

谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：有没有 最低点 如何

或者只看这个图：

1. n <- 1000
   
2. set.seed(731)
   
3. age <- 50 + 12*rnorm(n)
   
4. label(age) <- "Age"
   
5. sex <- factor(sample(c('Male','Female'), n,
   
6.                      rep=TRUE, prob=c(.6, .4)))
   
7. cens <- 15*runif(n)
   
8. h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
   
9. dt <- -log(runif(n))/h
   
10. label(dt) <- 'Follow-up Time'
    
11. e <- ifelse(dt <= cens,1,0)
    
12. dt <- pmin(dt, cens)
    
13. units(dt) <- "Year"
    
14. dd <- datadist(age, sex)
    
15. options(datadist='dd')
    
16. S <- Surv(dt,e)
    
17. 
    
18. f <- cph(S ~ rcs(age,4) , x=TRUE, y=TRUE)
    
19. cox.zph(f, "rank")             # tests of PH
    
20. anova(f)
    
21. plot(Predict(f, age)) # plot age effect, 2 curves for 2 sexes

复制代码

> Predict(f, age)
         age         yhat        lower        upper
1   19.71985 -0.580383381 -1.532437029  0.371670268
2   20.00869 -0.582811675 -1.519170266  0.353546915
3   20.29754 -0.585239970 -1.505919870  0.335439930
4   20.58638 -0.587668265 -1.492686694  0.317350164
5   20.87523 -0.590096559 -1.479471646  0.299278527
6   21.16408 -0.592524854 -1.466275702  0.281225993
.....

mengqinqing 发表于 2015-5-4 13:02
> Predict(f, age)
         age         yhat        lower        upper
1   19.71985 -0.580383381 -1 ...

因为你图中的曲线本身就是R基于Predict(f, age)这组数据进行某种光滑化得到的，所以并不存在数学意义上的解析最小值点。一种最小值点估计方法就是简单地取对应于yhat值最小的age值(

1. xmin <- age[which.min(yhat)]

复制代码

)。从你提供的数据来看， age变量中两点的间隔很小(我推测R其实就是根据这组数据绘出的曲线)，最后得到的估计应该会不错。可以画出来看看。

hugebear 发表于 2015-5-4 13:21
因为你图中的曲线本身就是R基于Predict(f, age)这组数据进行某种光滑化得到的，所以并不存在数学意义上的 ...

谢谢！
可以得到想要的age最小值。
我想把这个点加到图里面去，用points加不上去，不知道为什么？
但是如果用plot(x,y)去画图再加points则可以，这是为什么？

1. n <- 1000
   
2. set.seed(731)
   
3. age <- 50 + 12*rnorm(n)
   
4. label(age) <- "Age"
   
5. sex <- factor(sample(c('Male','Female'), n,
   
6.                      rep=TRUE, prob=c(.6, .4)))
   
7. cens <- 15*runif(n)
   
8. h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
   
9. dt <- -log(runif(n))/h
   
10. label(dt) <- 'Follow-up Time'
    
11. e <- ifelse(dt <= cens,1,0)
    
12. dt <- pmin(dt, cens)
    
13. units(dt) <- "Year"
    
14. dd <- datadist(age, sex)
    
15. options(datadist='dd')
    
16. S <- Surv(dt,e)
    
17. 
    
18. f <- cph(S ~ rcs(age,4) , x=TRUE, y=TRUE)
    
19. cox.zph(f, "rank")             # tests of PH
    
20. anova(f)
    
21. p<- Predict(f, age)
    
22. 
    
23. #用plot(p)画图
    
24. plot(p)
    
25. minpointyhat= min(p$yhat)
    
26. minpointage <- p$age[which.min(p$yhat)]#minpointage=p$age[which(p$yhat == min(p$yhat))]
    
27. points(x=minpointage,y=minpointyhat,col="red",pc=16)
    
28. 
    
29. #用plot自己画图
    
30. plot(p$yhat~p$age,type="l",ylim=c(-2,1.5))
    
31. lines(p$age,p$lower,lty=2)
    
32. lines(p$age,p$upper,lty=2)
    
33. points(x=minpointage,y=minpointyhat,col="red",pch=16)
    
34. text(x=minpointage,y=-0.8, "age=34.7",cex = .8)

复制代码

mengqinqing 发表于 2015-5-4 13:37
谢谢！
可以得到想要的age最小值。
我想把这个点加到图里面去，用points加不上去，不知道为什么？

你的原图看起来是用ggplot2画的，points()不是ggplot2里的绘图命令。加点要用符合ggplot2语法的命令才行。

hugebear 发表于 2015-5-4 21:43
你的原图看起来是用ggplot2画的，points()不是ggplot2里的绘图命令。加点要用符合ggplot2语法的命令才行。 ...

那我怎么样才能把这个图加上一个点呢？或者怎么样才能把后面那个图的置信区间填充上颜色呢？

hugebear 发表于 2015-5-4 21:43
你的原图看起来是用ggplot2画的，points()不是ggplot2里的绘图命令。加点要用符合ggplot2语法的命令才行。 ...

已解决，谢谢！

请问函数是什么，怎么解决的呢，拜托&#129402;