you can use the fitdistr() function in MASS package.

An example is :

library(MASS)

sample<-rexp(n=50,rate=3)

result<-fitdistr(sample,"exponential")

result$estimate
    rate
2.669378

result$loglik
[1] -0.9077293