thenceforth Fitting distribution in R
you can try here Use fitdistr in MASS package
Example;
fitdistr(as.vector(res.score.m), densefun = 'normal') mean sd -8.0647926561 2.9550064789 ( 0.0005642286) ( 0.0003989698)
http://www.statmethods.net/advgraphs/probability.html