## ----------------------------------------------------------------------------
##   R scripts for Haerdle/Mueller/Sperlich/Werwatz: "Nonparametric and
##   Semiparametric Modelling", Springer Series in Statistics, 2004
## ----------------------------------------------------------------------------
##   Script       SPMhisdiffori
## ----------------------------------------------------------------------------
##   Description  computes histograms with different origins
## ----------------------------------------------------------------------------
##   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.04),
     main="x0=0, h=0.04",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0.01,h=0.04),
     main="x0=0.01, h=0.04",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0.02,h=0.04),
     main="x0=0.02, h=0.04",xlab="x",ylab="fh", ylim=c(0,17))
hist(x,freq=FALSE,breaks=breaks(x,x0=0.03,h=0.04),
     main="x0=0.03, h=0.04",xlab="x",ylab="fh", ylim=c(0,17))
par(mfrow=c(1,1))