Problem in fitting of erlang distribution and estimation of initial or starting values

Question

I have fitted many distributions to my data but I am facing difficulty in fitting Erlang distribution to the data in R. I also used plotdist function to find starting or ini