Esta imagen debería volverse a crear como imágenes vectoriales SVG. Esto proporciona muchas ventajas, véase Commons:Media for cleanup (en inglés) para más información. Si ya hay una versión SVG de esta imagen disponible, por favor súbala a Commons. Tras subirla, reemplace esta plantilla con la plantilla {{vector version available|nuevo nombre de imagen.svg}} en esta imagen.
The Matlab script for this example uses
kde.m and is given below.
% Datarandn('seed',1)% Used for reproducibilitydata=[randn(100,1)-10;randn(100,1)+10];% Two Normals mixed% Truephi=@(x)exp(-.5*x.^2)/sqrt(2*pi);% Normal Densitytpdf=@(x)phi(x+10)/2+phi(x-10)/2;% True Density% Kernelh=std(data)*(4/3/numel(data))^(1/5);% Bandwidth estimated by Silverman's Rule of Thumbkernel=@(x)mean(phi((x-data)/h)/h);% Kernel Densitykpdf=@(x)arrayfun(kernel,x);% Elementwise application% Plotfigure(2),clf,holdonx=linspace(-25,+25,1000);% Linear Spaceplot(x,tpdf(x))% Plot True Densityplot(x,kpdf(x))% Plot Kernel Density with Silverman's Rule of Thumbkde(data)% Plot Kernel Density with Solve-the-Equation Bandwidth
#The same code with R language#` Dataset.seed(1)#Used for reproducibilitydata=c(rnorm(100,-10,1),rnorm(100,10,1))#Two Normals mixed#` Truephi=function(x)exp(-.5*x^2)/sqrt(2*pi)#Normal Densitytpdf=function(x)phi(x+10)/2+phi(x-10)/2#True Density#` Kernelh=sd(data)*(4/3/length(data))^(1/5)#Bandwidth estimated by Silverman's Rule of ThumbKernel2=function(x)mean(phi((x-data)/h)/h)#Kernel Densitykpdf=function(x)sapply(x,Kernel2)#Elementwise application#` Plotx=seq(-25,25,length=1000)#Linear Spaceplot(x,tpdf(x),type="l",ylim=c(0,0.23),col="red")#Plot True Densitypar(new=T)plot(x,kpdf(x),type="l",ylim=c(0,0.23),xlab="",ylab="",axes=F)#Plot Kernel Density with Silverman's Rule of Thumb
Licencia
User:Kernel estimator, titular de los derechos de autor de esta obra, la publica en los términos de la siguiente licencia:
La persona que ha asociado una obra a este documento lo dedica al dominio público mediante la cesión mundial de sus derechos bajo la ley de derechos de autor y todos los derechos legales adyacentes propios de dicha, en el ámbito permitido por ley. Puedes copiar, modificar, distribuir y reproducir el trabajo, incluso con objetivos comerciales, sin pedir aprobación del autor.
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse