R help - Oct 2007 - Batch-Plot Histograms

Hi everyone,

 

anybody’s got an idea why the following script doesn’t produce batch-histograms?

I’m using Windows XP and R 2.5.1.

 

Here’s the script:

 

matrix<-read.csv("C:\\Stadtwerke_mit_Umlage.csv",
header=TRUE,sep=";",dec=".")

Stadtwerke<-colnames(matrix)

Bereich_blau<-66.67

Schleifen<-dim(matrix)[2]

Plan<-c(35.46,1.75,63.49,13.52,39.38,21.83)

Median<-apply(matrix,2,median)

 

breaks<-round((0.01)*dim(matrix[1]))

 

for(i in 1:Schleifen){

  hist(matrix[,i],main=paste("Risiko-Chancenprofil der Stadtwerke ",

  Stadtwerke[i]),col="blue",xlab="EBIT in Mio. €",
ylab="Häufigkeit",breaks[1])

  abline(v=Median[i],lwd=2,col="grey",lty=3)

  axis(side=1,at=Median[i],labels="Median",col="grey")

  hist_data<-hist(matrix[,i],plot=F,breaks[1])

  Quantil_unten<-quantile(matrix[,i],((100-Bereich_blau)/2)/100)

  Quantil_oben<-quantile(matrix[,i], (100-(100-Bereich_blau)/2)/100)

 

  #Plotten der roten unterer Teil

    for(j in 1:length(hist_data$breaks)){

      if (hist_data$breaks[j]<Quantil_unten)

      hv<-j

    }

  for(k in 1:hv){

         
rect(hist_data$breaks[k],0,hist_data$breaks[k+1],hist_data$counts[k],col="red")

        }

 

    for(l in 1:length(hist_data$breaks)){

      if (hist_data$breaks[l]<Quantil_oben)

      hv2<-l+1

    }

    for(m in hv2:length(hist_data$breaks)){

         
rect(hist_data$breaks[m],0,hist_data$breaks[m+1],hist_data$counts[m],col="red")

        }

  dev.copy(device =
postscript,file=paste("C:\\",Stadtwerke[i],".ps"))

  dev.off()

}

 

It seems, that the error is in line „dev.copy“… but I don’t understand why. It
doesn’t matter if I’m using device = jpeg or postscript. Maybe I must add a
waiting time of 2 seconds or something like that?

 

Thanks a lot,

Thomas


	[[alternative HTML version deleted]]

why don't you just output to postscript instead of using dev.copy:

postscript("yourfile.ps")
.......your script......
dev.off()

On 10/12/07, Thomas Schwander <thomas_schwander at web.de>
wrote:
> Hi everyone,
>
>
>
> anybody's got an idea why the following script doesn't produce
batch-histograms?
>
> I'm using Windows XP and R 2.5.1.
>
>
>
> Here's the script:
>
>
>
> matrix<-read.csv("C:\\Stadtwerke_mit_Umlage.csv",
header=TRUE,sep=";",dec=".")
>
> Stadtwerke<-colnames(matrix)
>
> Bereich_blau<-66.67
>
> Schleifen<-dim(matrix)[2]
>
> Plan<-c(35.46,1.75,63.49,13.52,39.38,21.83)
>
> Median<-apply(matrix,2,median)
>
>
>
> breaks<-round((0.01)*dim(matrix[1]))
>
>
>
> for(i in 1:Schleifen){
>
>  hist(matrix[,i],main=paste("Risiko-Chancenprofil der Stadtwerke
",
>
>  Stadtwerke[i]),col="blue",xlab="EBIT in Mio. ?",
ylab="H?ufigkeit",breaks[1])
>
>  abline(v=Median[i],lwd=2,col="grey",lty=3)
>
>  axis(side=1,at=Median[i],labels="Median",col="grey")
>
>  hist_data<-hist(matrix[,i],plot=F,breaks[1])
>
>  Quantil_unten<-quantile(matrix[,i],((100-Bereich_blau)/2)/100)
>
>  Quantil_oben<-quantile(matrix[,i], (100-(100-Bereich_blau)/2)/100)
>
>
>
>  #Plotten der roten unterer Teil
>
>    for(j in 1:length(hist_data$breaks)){
>
>      if (hist_data$breaks[j]<Quantil_unten)
>
>      hv<-j
>
>    }
>
>  for(k in 1:hv){
>
>         
rect(hist_data$breaks[k],0,hist_data$breaks[k+1],hist_data$counts[k],col="red")
>
>        }
>
>
>
>    for(l in 1:length(hist_data$breaks)){
>
>      if (hist_data$breaks[l]<Quantil_oben)
>
>      hv2<-l+1
>
>    }
>
>    for(m in hv2:length(hist_data$breaks)){
>
>         
rect(hist_data$breaks[m],0,hist_data$breaks[m+1],hist_data$counts[m],col="red")
>
>        }
>
>  dev.copy(device =
postscript,file=paste("C:\\",Stadtwerke[i],".ps"))
>
>  dev.off()
>
> }
>
>
>
> It seems, that the error is in line ?dev.copy"? but I don't
understand why. It doesn't matter if I'm using device = jpeg or
postscript. Maybe I must add a waiting time of 2 seconds or something like that?
>
>
>
> Thanks a lot,
>
> Thomas
>
>
>        [[alternative HTML version deleted]]
>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>

-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem you are trying to solve?