## ----------------------------------------------------------------------------
##   R scripts for Haerdle/Mueller/Sperlich/Werwatz: "Nonparametric and
##   Semiparametric Modelling", Springer Series in Statistics, 2004
## ----------------------------------------------------------------------------
##   Script      SPMhisdiffbin
## ----------------------------------------------------------------------------
##   Description  computes histograms with different binwidths
## ----------------------------------------------------------------------------
##   Author       Marlene Mueller, 2008/11/21
## ----------------------------------------------------------------------------

x <- read.csv("D:\\stockres.txt")
x <- unlist(x)

breaks <- function(x,x0,h){
  b <- floor((min(x)-x0)/h) : ceiling((max(x)-x0)/h)
  b <- b*h+x0
  return(b)
}

par(mfrow=c(2,2))
hist(x,freq=FALSE,breaks=breaks(x,x0=0,h=0.007),
     main="x0=0, h=0.007",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0,h=0.02 ),
     main="x0=0, h=0.02 ",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0,h=0.05 ),
     main="x0=0, h=0.05 ",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0,h=0.1  ),
     main="x0=0, h=0.1  ",xlab="x",ylab="fh", ylim=c(0,17))
par(mfrow=c(1,1))