[R] Output results to a single postscript document

[David Winsemius]

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:

test.pdf is attached below but my Mac email client sometimes does  
unexpected (to me anyway) things to graphics attachments.
-- 
David Winsemius
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.pdf
Type: application/pdf
Size: 163 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090207/533935ff/attachment-0002.pdf>
-------------- next part --------------




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.