Erro tentando desenhar grafico: “comprimento do objeto maior não é múltiplo do comprimento do objeto menor”

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?