Output results to a single postscript document
Sorry, the message seems to have got botched. Here it is again:
Pele does open a PDF device (previously it was a postscript device). It
looks like what Pele is trying to do is "plot" the printed results of the
summary of a model + the AIC, together with the acf() and pacf() plots.
As Dieter indicated, Gavin Snow's textplot() function [in gplots package]
might help. This really is a job for Sweave, which Pele should have look at.
Pele: you could also look at Prof. Harrell's text.anova.Design function,
which you may be able to adapt to your needs. Even "simpler" would be to
write a function that uses text() [see: ?text] to plot the parts of your
model-summary and the AIC to an empty graph.
Regards, Mark.
Hello R users,
Below is the code and output of what I am trying to do. My goal is to
insert/print all items in the chart function into a pdf document. Only the
acf and pacf charts gets printed.
Again, thanks in advance for any help I can get!
####################################
options (scipen=999, digits=7)
library(lmtest)
library(dyn)
tdata <- ts(read.table("C:/test/data.csv" ,sep = ",",header=TRUE))
print(tdata)
res <- function(dep, indep, fn, env = parent.frame()) {
depnm <- deparse(substitute(dep))
indepnm <- deparse(substitute(indep))
fo <- sprintf("%s ~ lag(%s, -1) + %s", depnm, depnm, indepnm)
fo <- as.formula(fo, env = env)
mod <- do.call(dyn$lm, list(fo))
chart <- function(mod) {
pdf(file="results.pdf")
par(mfrow=c(2,1))
acf(residuals(mod))
pacf(residuals(mod))
summ<- summary(mod)
aic <- AIC(mod)
dev.off() }
output <- return(list(Charts=chart(mod))) }
results <- res(tdata[ ,"Y"], tdata[ ,"X_VARY"]) ; results
############output###############
options (scipen=999, digits=7)
library(lmtest)
library(dyn)
tdata <- ts(read.table("C:/test/data.csv" ,sep = ",",header=TRUE))
print(tdata)
Time Series: Start = 1 End = 20 Frequency = 1 Unit Y X_VARY K D_WAY 1 1 3 2 4 50 2 2 5 3 9 50 3 3 6 11 22 50 4 4 8 4 7 2 5 5 9 11 11 2 6 6 12 13 13 2 7 7 23 25 12 2 8 8 22 30 31 3 9 9 23 3 3 3 10 10 19 21 21 32 11 11 3 2 4 34 12 12 5 3 9 4 13 13 6 11 22 4 14 14 8 4 7 4 15 15 9 11 11 4 16 16 12 13 13 4 17 17 23 25 12 4 18 18 22 30 31 4 19 19 23 3 3 4 20 20 19 21 21 5
res <- function(dep, indep, fn, env = parent.frame()) {
+ depnm <- deparse(substitute(dep))
+ indepnm <- deparse(substitute(indep))
+ fo <- sprintf("%s ~ lag(%s, -1) + %s", depnm, depnm, indepnm)
+ fo <- as.formula(fo, env = env)
+ mod <- do.call(dyn$lm, list(fo))
+
+ chart <- function(mod) {
+ pdf(file="results.pdf")
+ par(mfrow=c(2,1))
+ acf(residuals(mod))
+ pacf(residuals(mod))
+ summ<- summary(mod)
+ aic<- AIC(mod)
+ dev.off() }
+
+ output <- return(list(Charts=chart(mod))) }
results <- res(tdata[ ,"Y"], tdata[ ,"X_VARY"]) ; results
$Charts null device http://www.nabble.com/file/p21888211/results.pdf results.pdf
Pele wrote:
Hello R users, I have been trying to output all my results (text, plots, etc) into the same postscript file as one document, but have been unable to...Can anyone help me improve my code below so that I can accomplish this? Currently I have to output them separately then piece them back together into one document..
[[elided Yahoo spam]]
View this message in context: http://www.nabble.com/Output-results-to-a-single-postscript-document-tp21884901p21889042.html Sent from the R help mailing list archive at Nabble.com.