ОписаниеKernel density estimation, comparison between rule of thumb and solve-the-equation bandwidth.png
У этого изображения отсутствует описание и может отсутствовать другая информация.
Медиафайл обязан иметь описание (сводку) для информирования других о содержании, авторе, источнике и дате, если это возможно. Если вы знаете такую информацию или имеете к ней доступ, пожалуйста, добавьте её на страницу описания файла.
Этот график желательно воссоздать или аккуратно преобразовать в векторный формат SVG. Это даёт несколько преимуществ, прочитать о которых подробнее вы можете на странице Commons:Media for cleanup. Если вам уже сейчас доступна векторная версия данного изображения, загрузите её, пожалуйста, а затем замените этот шаблон на следующий: {{Vector version available|Имя загруженного файла.svg}}.
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
Лицензирование
User:Kernel estimator, владелец авторских прав на это произведение, добровольно публикует его на условиях следующей лицензии:
Лица, связанные с работой над этим произведением, решили передать данное произведение в общественное достояние, отказавшись от всех прав на произведение по всему миру в рамках закона об авторских правах (а также связанных и смежных прав), в той степени, которую допускает закон. Вы можете копировать, изменять, распространять, исполнять данное произведение в любых целях, в том числе в коммерческих, без получения на это разрешения автора.
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse
Файл содержит дополнительные данные, обычно добавляемые цифровыми камерами или сканерами. Если файл после создания редактировался, то некоторые параметры могут не соответствовать текущему изображению.