Uncertain fractional order neutral differential equation is an important model to describe the evolution process of uncertain dynamical system. This paper devotes to studying linear uncertain fractional order neutral differential equations. After providing the analytic solutions for linear uncertain fractional order neutral differential equations by the Mittag-Leffler function, the author investigates the inverse uncertainty distribution of the solution to linear uncertain fractional order neutral differential equation by the α-path and has a discussion of the dependence of the solution on initial function based on the generalized Gronwall inequality.