[英]Error when generating pdf using script in R
我正在使用R循環遍歷數據框的列,並繪制結果分析圖。 腳本運行時我沒有收到任何錯誤,但它生成了一個無法打開的pdf。
如果我運行腳本的內容,它可以正常工作。 我想知道循環的速度是否有問題,所以我試圖強迫它暫停。 這似乎沒有什么區別。 我對人們有任何建議感興趣,而且我對R也很陌生,所以我也歡迎如何改進方法的建議。 謝謝。
for (i in 2:22) {
# Organise data
pop_den_z = subset(pop_den, pop_den[i] != "0") # Remove zeros
y = pop_den_z[,i] # Get y col
x = pop_den_z[,1] # get x col
y = log(y) # Log transform
# Regression
lm.0 = lm(formula = y ~ x) # make linear model
inter = summary(lm.0)$coefficients[1,1] # Get intercept
slop = summary(lm.0)$coefficients[2,1] # Get slope
# Write to File
a = c(i, inter, slop)
write(a, file = "C:/pop_den_coef.txt", ncolumns = 3, append = TRUE, sep = ",")
## Setup pdf
string = paste("C:/LEED/results/Images/R_graphs/Pop_den", paste(i-2), "City.pdf")
pdf(string, height = 6, width = 9)
p <- qplot(
x, y,
xlab = "Radius [km]",
ylab = "Population Density [log(people/km)]",
xlim = x_range,
main = "Analysis of Cities"
)
# geom_abline(intercept,slope)
p + geom_abline(intercept = inter, slope = slop, colour = "red", size = 1)
Sys.sleep(5)
### close the PDF file
dev.off()
}
這條線應該是
print(p + geom_abline(intercept = inter, slope = slop, colour = "red", size = 1))
在pdf設備中,ggplot(和lattice)僅在顯式打印時寫入文件。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.