Archivo:Kernel density estimation, comparison between rule of thumb and solve-the-equation bandwidth.png

Contenido de la página no disponible en otros idiomas.
De Wikipedia, la enciclopedia libre

Kernel_density_estimation,_comparison_between_rule_of_thumb_and_solve-the-equation_bandwidth.png(733 × 561 píxeles; tamaño de archivo: 28 kB; tipo MIME: image/png)

Resumen

Descripción
Fecha
Fuente
Matlab Logo 
Este diagrama fue creado con MATLAB.
Autor Kernel estimator
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.

% Data
randn('seed',1)                            % Used for reproducibility
data = [randn(100,1)-10; randn(100,1)+10]; % Two Normals mixed
% True
phi = @(x) exp(-.5*x.^2)/sqrt(2*pi);       % Normal Density
tpdf = @(x) phi(x+10)/2+phi(x-10)/2;       % True Density
% Kernel
h = std(data)*(4/3/numel(data))^(1/5);     % Bandwidth estimated by Silverman's Rule of Thumb
kernel = @(x) mean(phi((x-data)/h)/h);     % Kernel Density
kpdf = @(x) arrayfun(kernel,x);            % Elementwise application
% Plot
figure(2), clf, hold on
x = linspace(-25,+25,1000);                % Linear Space
plot(x,tpdf(x))                            % Plot True Density
plot(x,kpdf(x))                            % Plot Kernel Density with Silverman's Rule of Thumb
kde(data)                                  % Plot Kernel Density with Solve-the-Equation Bandwidth
#The same code with R language
#` Data
set.seed(1)                                                     #Used for reproducibility
data = c(rnorm(100,-10,1),rnorm(100,10,1))                      #Two Normals mixed
#` True
phi = function(x) exp(-.5*x^2)/sqrt(2*pi)                       #Normal Density
tpdf = function(x) phi(x+10)/2+phi(x-10)/2                      #True Density
#` Kernel
h = sd(data)*(4/3/length(data))^(1/5)                           #Bandwidth estimated by Silverman's Rule of Thumb
Kernel2 = function(x) mean(phi((x-data)/h)/h)                   #Kernel Density
kpdf = function(x) sapply(x,Kernel2)                            #Elementwise application
#` Plot
x=seq(-25,25,length=1000)                                       #Linear Space
plot(x,tpdf(x),type="l",ylim=c(0,0.23),col="red")               #Plot True Density
par(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:
Creative Commons CC-Zero Este archivo está disponible bajo la licencia Creative Commons Dedicación de Dominio Público CC0 1.0 Universal.
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.

Registro original de carga

Transferido desde en.wikipedia a Commons usando For the Common Good.

Aquí se muestra la página de descripción original. Los siguientes nombres de usuario se refieren a en.wikipedia.
Fecha y hora Dimensiones Usuario Comentario
15:47, 19 April 2016 733 × 561 (28,809 bytes) w:en:Kernel estimator (discusión | contribs.) author: kernel estimator source: created by kernel estimator using the Matlab code on Wikipedia URL: N/A

Leyendas

Añade una explicación corta acerca de lo que representa este archivo

Elementos representados en este archivo

representa a

Historial del archivo

Haz clic sobre una fecha y hora para ver el archivo tal como apareció en ese momento.

Fecha y horaMiniaturaDimensionesUsuarioComentario
actual20:39 26 oct 2018Miniatura de la versión del 20:39 26 oct 2018733 × 561 (28 kB)JumpowTransferred from en.wikipedia: see original upload log above

La siguiente página usa este archivo:

Uso global del archivo

Las wikis siguientes utilizan este archivo:

Metadatos