pdf_sigma_data_1d Subroutine

subroutine pdf_sigma_data_1d(self, data, h)

Arguments

Type IntentOptional Attributes Name
type(pdf_kernel) :: self
real(kind=8), dimension(:) :: data
real(kind=8) :: h

Calls

proc~~pdf_sigma_data_1d~2~~CallsGraph proc~pdf_sigma_data_1d~2 pdf_sigma_data_1d get_sd get_sd proc~pdf_sigma_data_1d~2->get_sd

Source Code

subroutine pdf_sigma_data_1d(self,data,h)
  type(pdf_kernel)     :: self
  real(8),dimension(:) :: data
  real(8)              :: h
  if(.not.self%status)stop "PDF_SET_SIGMA: PDF not allocated"
  h = get_sd(data)
  h = (4d0/3d0/size(data))**(1d0/5d0)*h !Silverman's rule of thumb.
end subroutine pdf_sigma_data_1d