Смотрите участок рядом с перекрестком - PullRequest
0 голосов
/ 13 июня 2018
## Define the moment generating function of Weibull distribution
scale <- 1/2
shape <- 1 
lambda <- 2 
beta <- 0.1

## I have specified nmax=160 since I cant perform the sum until infinity
mgfw <- function(x){
  nmax <- 160
  scale <- scale
  shape <- shape
  suma <- 0
  for(n in 0:nmax){
    suma <- suma + ((x^n) * ((scale)^n)) * gamma(1 + (n/shape)) / factorial(n)
    }   
  return(suma)
  }  

curve(mgfw, from=0.1, 0.25, ylim=c(1, 1.2))
mu <- (scale) * gamma(1 + (1 / shape))

fun2 <- function(x) 1 + x * (1 + beta) * mu

x <- seq(0, 10, length.out=100)  
y <- fun2(x)

curve(fun2, from=0, 10, add=TRUE)
grid()

Решая предыдущие уравнения, я получил следующие результаты:

library(rootSolve)
r <- uniroot.all(function(x) mgfw(x) - fun2(x), c(0.1, 0.185))
r
abline(v=r)

Я получил график, подобный этому

enter image description here

Пересечение обеих линий дается вертикальной линией.Но я хотел бы получить сюжет, где пересечение может быть четким в сюжете.Как изменить размер сюжета?Или видите в области с разным масштабом?

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...