[R] Output results to a single postscript document

[Mark Difford]

Hi David, Pele:


David Winsemius wrote:
> 
>>> I don't see anywhere that you opened a pdf device. When I try :
> 
>>> pdf("test.pdf") and then run your code I get what looks like the  
>>> desired output sitting in my working directory:
> 
> 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.
> 
> 
> On Feb 7, 2009, at 8:16 AM, Pele wrote:
> 
>>
>> 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
>>
>>
>> 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]]
>>
>>
>>
>> http://www.nabble.com/file/p21888211/data.csv data.csv
>> -- 
>> View this message in context:
>> http://www.nabble.com/Output-results-to-a-single-postscript-document-tp21884901p21888211.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
> 
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
> 
> 

-- 
View this message in context: http://www.nabble.com/Output-results-to-a-single-postscript-document-tp21884901p21889020.html
Sent from the R help mailing list archive at Nabble.com.