16

在 Stata 中编程时,我经常发现自己在编程中使用循环索引。例如,我将遍历变量nominalprice 和realprice 的列表:

local list = "nominalprice realprice"
foreach i of local list {
  summarize `i'
  twoway (scatter `i' time)
  graph export "C:\TimePlot-`i'.png"
}

这将绘制名义价格和实际价格的时间序列,并导出一个名为 TimePlot-nominalprice.png 的图表和另一个名为 TimePlot-realprice.png 的图表。

在 R 中,我想出的做同样事情的方法是:

clist <- c("nominalprice", "realprice")
for (i in clist) {
  e <- paste("png(\"c:/TimePlot-",i,".png\")", sep="")
  eval(parse(text=e))
  plot(time, eval(parse(text=i)))
  dev.off() 
}

这个 R 代码对我来说看起来不直观和混乱,我还没有找到在 R 中做这种事情的好方法。也许我只是没有以正确的方式思考这个问题?你能建议一种更好的使用字符串循环的方法吗?

4

4 回答 4

20

正如其他人所暗示的那样,如果您有一个包含名为nominalpriceand的列的数据框,这会更容易realprice。如果你不这样做,你总是可以使用get. 你根本不需要parse这里。

clist <- c("nominalprice", "realprice")
for (i in clist) {
   png(paste("c:/TimePlot-",i,".png"), sep="")
   plot(time, get(i))
   dev.off() 
}
于 2009-11-02T20:48:40.887 回答
2

如果您的主要问题是需要键入 eval(parse(text=i)) 而不是“i”,则可以创建一个更易于使用的函数来评估字符串中的表达式:

e = function(expr) eval(parse(text=expr))

然后R示例可以简化为:

clist <- c("nominalprice", "realprice")
for (i in clist) {
  png(paste("c:/TimePlot-", i, ".png", sep=""))
  plot(time, e(i))
  dev.off() 
}
于 2009-11-02T20:37:48.287 回答
1

我看不出您的原始解决方案有什么特别错误,只是我不知道您为什么使用 eval() 函数。这对我来说似乎没有必要。

您还可以使用应用功能,例如 lapply。这是一个工作示例。我将虚拟数据创建为zoo()时间序列(这不是必需的,但是因为您无论如何都在使用时间序列数据):

# x <- some time series data
time <- as.Date("2003-02-01") + c(1, 3, 7, 9, 14) - 1
x <- zoo(data.frame(nominalprice=rnorm(5),realprice=rnorm(5)), time)
lapply(c("nominalprice", "realprice"), function(c.name, x) { 
  png(paste("c:/TimePlot-", c.name, ".png", sep=""))
  plot(x[,c.name], main=c.name)
  dev.off()
}, x=x)
于 2009-11-02T20:31:15.563 回答
1

使用 ggplot2 并重塑:

library(ggplot2)
library(reshape)
df <- data.frame(nominalprice=rexp(10), time=1:10)
df <- transform(df, realprice=nominalprice*runif(10,.9,1.1))
dfm <- melt(df, id.var=c("time"))
qplot(time, value, facets=~variable, data=dfm)
于 2009-11-02T20:48:10.793 回答