“有条件的”地块

发布于 2024-12-03 16:44:52 字数 274 浏览 1 评论 0原文

我很好奇如何绘制这样定义的函数：

 if(x < 1)
   f(x) = x/10 * 1.2
 if(x < 3)
   f(x) = x/12 * 1.7
 ...
 else
   f(x) = x/15 * 2

如果函数很简单，例如 f(x) = x/10 * x/5 ，那么就不会有问题，可以使用 curve() 方法。但是我不确定处理更复杂的函数（如上面的函数）的最佳方法是什么。有什么想法吗？如果可以使用 ggplot() 的话，加分:)

原文

I am curious how to plot function that is defined something like this:

 if(x < 1)
   f(x) = x/10 * 1.2
 if(x < 3)
   f(x) = x/12 * 1.7
 ...
 else
   f(x) = x/15 * 2

If the function was simple, say f(x) = x/10 * x/5 , then there would be no problem, and one could use curve() method. However I am not sure what is the best way to deal with more complex functions, like the one above. Any ideas? Bonus points, if ggplot() could be used :)

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

荒岛晴空 2024-12-10 16:44:52

曲线还是有可能的。（当您阅读统计文献时，这个公式显示为 I[x]，“I”代表“指标”。）

curve( (x <1)*( (x/10)*1.2 ) +       # one line for each case
       (!(x <1)&(x<3) )*(x/12)*1.7 + # logical times (local) function
        (x >=3)*(x/15)*2 ,
        0,4)                         # limits

“在此输入图像描述”

Curve is still a possibility. (And as you read the statistical literature, this formulation shows up as I[x], "I" being for "indicator".)

curve( (x <1)*( (x/10)*1.2 ) +       # one line for each case
       (!(x <1)&(x<3) )*(x/12)*1.7 + # logical times (local) function
        (x >=3)*(x/15)*2 ,
        0,4)                         # limits

enter image description here

回复收藏 0 原文

老子叫无熙 2024-12-10 16:44:52

您在寻找类似 stepfun 的东西吗？

fn <- stepfun(c(1,2,3,4,5),c(0,1,2,3,4,5))
plot(fn,verticals = FALSE)

在此处输入图像描述

指定函数的方式会有点不同，但一旦掌握了就很容易掌握阅读 ?stepfun 并自己绘制一些示例。

Are you looking for something like stepfun?

fn <- stepfun(c(1,2,3,4,5),c(0,1,2,3,4,5))
plot(fn,verticals = FALSE)

enter image description here

The way you specify the function will be a bit different, but it's fairly easy to grasp once you've read ?stepfun and plotted some examples yourself.

回复收藏 0 原文

宁愿没拥抱 2024-12-10 16:44:52

口味问题，但我更喜欢 ifelse 而不是 Dwins 指标（例如mbq 评论）。用于比较：

curve(
    (x <1)           * ( (x/10)*1.2 ) +
    (!(x <1)&(x<3) ) * ( (x/12)*1.7 ) +
    (x >=3)          * ( (x/15)*2   ) ,
    0,4)

# versus

curve(
    ifelse(x < 1, (x/10)*1.2,
    ifelse(x < 3, (x/12)*1.7,
                  (x/15)*2    )),
    0,4)

Matter of taste but I prefer ifelse over Dwins indicators (like in mbq comment). For compare:

curve(
    (x <1)           * ( (x/10)*1.2 ) +
    (!(x <1)&(x<3) ) * ( (x/12)*1.7 ) +
    (x >=3)          * ( (x/15)*2   ) ,
    0,4)

# versus

curve(
    ifelse(x < 1, (x/10)*1.2,
    ifelse(x < 3, (x/12)*1.7,
                  (x/15)*2    )),
    0,4)

回复收藏 0 原文

~没有更多了~