pnorm是dnorm的累计分布函数

Posted 2021-11-07 基督徒Isaac

library(tidyverse)
x <- seq(0,12,0.01)
p <- pnorm(x,4,2)
x %>% plot(.,
           p,
           type="l",
           col=1,
           xlab="工作时间",
           ylab="累计产量",
           xaxt="n",
           yaxt="n"
           )
axis(1,at=seq(0,12,1),
     labels=paste(seq(0,12,1),"h"),las=1)
axis(2,at=seq(0,1,0.2),
     labels=c(seq(0,40,10),"max"),las=2)
d <- dnorm(x,4,2)+1-dnorm(x,4,2) %>% max()
x %>% lines(.,
            d,
            lty=3,
            col=2
            )
legend("bottomright",
       legend=c("工作效率","累计产量"),
       col=c(2,1),
       lty=c(3,1),
       lwd=1
       )

结果：