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?
