Tenho uma função:
my_gamma = function(x)
{
f = function(t){t^(x-1) * exp(-t)}
integrate(f, 0, Inf)
}
Estou tentando desenhar o seu gráfico com ggplot2
:
plot = ggplot(data.frame(x = seq(0.01, 10)), aes(x=x))
(plot + stat_function(fun = my_gamma))
Mas está me dando o erro:
Warning messages:
1: In t^(x - 1) :
comprimento do objeto maior não é múltiplo do comprimento do objeto menor
2: In t^(x - 1) * exp(-t) :
comprimento do objeto maior não é múltiplo do comprimento do objeto menor
3: Computation failed in `stat_function()`:
evaluation of function gave a result of wrong length
O que estou a fazer de errado?