Distribution Fitting

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

Distribution Fitting

L.Isella
Dear All,
I am experiencing some problems in fitting a trimodal distribution (which should be thought as a sum of three Gaussian distributions) using the nls function for nonlinear fittings.
As a test, just consider the very simple code:

rm(list=ls());
mydata<-rnorm(10000,0,4);
mydens<-density(mydata,kernel="gaussian");
y1<-mydens$y;
x1<-mydens$x;
myfit<-nls(y1~A*exp(-x1^2/sig));


which I use to get the empirical density (as I would from real experimental data) and test it against a Gaussian ansatz.
Well, either R always crashes for a segmentation fault and I have to restart it manually or I get this output:

Error in match.call(definition, call, expand.dots) :
        '.Primitiv�i�d�������������...' is not a function

Am I missing the obvious or is there some bug in my R build?
Many thanks

Lorenzo

_______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-finance