Domanda

I would like to generate 6 plots in one layout, but my code generates only the last plot. This is my code:

i<-c(500,1500,2500)
for(n in i)
{
  par(mfrow = c( 6, 1 ))
  window<-333
  Filename2<-paste("/home/blabla/Bilder/6Fenster",n,".png",sep="")

  png(filename=Filename2,width = 800,height = 600)  

  plot(datalist[n:(n+window-1),1],type="l",col="blue",xlab=n+window,ylab="values")

  x.specR<-spectrum(datalist[n:(n+window-1),1],plot=FALSE)
  plot((1:180),x.specR$spec,type="l",xlab="Frequenz",ylab="Frequenz")  

  hfft<-hanning.window(333)*datalist[n:(n+window-1)]
  x.specHF<-spectrum(hfft,plot=FALSE)
  plot((1:180),x.specHF$spec,type="l",xlab="Frequenz",ylab="Frequenz")

  m<-acf(datalist[n:(n+window-1),1],lag.max=333,plot=FALSE)$acf[-1]
  plot(m,type="l",ylim=c(-1,1),axes=FALSE, frame.plot=TRUE)
  axis(1,at=333/(1:333),labels=rev(333/(1:333)),axTicks(4))
  axis(2, at=c(-1,1))

  x.specA<-spectrum(m,plot=FALSE)
  plot((1:180),x.specA$spec,type="l",xlab="Frequenz",ylab="Frequenz")

  ha<-hanning.window(333)*m
  x.specHA<-spectrum(ha,plot=FALSE)
  plot((1:180),x.specHA$spec,type="l",xlab="Frequenz",ylab="Frequenz")

  dev.off()
}
È stato utile?

Soluzione

The par function should be called after you open the graphics device, which in this case is done with png. So you simply need to move that line.

Autorizzato sotto: CC-BY-SA insieme a attribuzione
Non affiliato a StackOverflow
scroll top