我可以在命令行上尝试使用4x4绘图,但在下面的代码中不能使用for -loop和png元素。我很沮丧问题是什么,我现在已经试着理解这个时间而且无法前进。
将下面的代码复制到paper.Rnw -file后,只需$ R CMD Sweave paper.Rnw
,$ pdflatex paper.tex
和$ evince paper.pdf
,您将获得1x1的绘图而不是4x4的绘图。为什么呢?
$ mkdir Pictures
$ mkdir Data
$ cd Data
$ wget https://noppa.aalto.fi/noppa/kurssi/mat-2.3128/harjoitustyot/Mat-2_3128_data_2.xls
$ cd ..
$ cat paper.Rnw
\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage{amsfonts}
\usepackage{amsmath}
\usepackage{graphicx}
\begin{document}
<<echo=FALSE>>=
library(gdata)
f<-read.xls('./Data/Mat-2_3128_data_2.xls')
for(i in 1:ncol(f))
{
png(paste('./Pictures/CO',i,'.png',sep =''))
par(mfrow=c(2,2))
plot(pacf(na.omit(f[,i])), main=paste('Pacf', i))
plot(acf(na.omit(f[,i])), main=paste('Acf', i))
plot(na.omit(f[,i]), main=paste('Data', i), type='l')
}
@
\subsection{Time-serie C1}
\includegraphics{./Pictures/CO1.png}
\subsection{Time-serie C2}
\includegraphics{./Pictures/CO2.png}
\subsection{Time-serie C3}
\includegraphics{./Pictures/CO3.png}
\subsection{Time-serie C4}
\includegraphics{./Pictures/CO4.png}
\subsection{Time-serie C5}
\includegraphics{./Pictures/CO5.png}
\end{document}
答案 0 :(得分:5)
函数acf
和pacf
的参数plot
默认为TRUE
。因此,不需要绘图语句(否则您将获得acf
的两个图和pacf
的两个图。请注意,我还添加了dev.off()
声明。
这应该有效:
for(i in 1:ncol(f))
{
png(paste('./Pictures/CO',i,'.png',sep =''))
par(mfrow=c(2,2))
pacf(na.omit(f[,i]), main=paste('Pacf', i))
acf(na.omit(f[,i]), main=paste('Acf', i))
plot(na.omit(f[,i]), main=paste('Data', i), type='l')
dev.off()
}
@