Skip to content
Back to formatted view

Raw Message

Message-ID: <0A66BC2F-9829-4AE0-8673-F67F0DC2F7A8@comcast.net>
Date: 2011-10-11T17:13:06Z
From: David Winsemius
Subject: Help to write to a file
In-Reply-To: <CAHKvgGspnr7hjvW9_g7vY53dHPu_RXTVFcvux0OHWU+yd8bDUA@mail.gmail.com>

On Oct 11, 2011, at 1:03 PM, Sergio Ren? Araujo Enciso wrote:

> Dear all:
>
> I am having some problems to use the function "sink()". Basically I  
> am doing
> a loop over two files which contain unit-root variables. Then on a  
> loop, I
> extract every i element of both files to create an object called z.  
> If z
> meets some requirements, then I perform a unit root test (ADF test),
> otherwise not. As this process is repeated several times, for each i  
> I want
> to get the summary of the ADF test on a common file. For that I use  
> the
> function "sink()". My code runs fine, but I do not get anything  
> written on
> the text file where my results are supposed to be saved. The code is  
> below
>
> setwd("C:\\Users\\Sergio Ren?\\Dropbox\\R")
>
> library("urca")
>
> P1<-read.csv("2R_EQ_P_R1_500.csv")
> P2<-read.csv("2R_EQ_P_R2_500.csv")
>
> d<-(1:1000)
> sink ("ADF_results_b_1.txt")
>
> for (i in seq(d))
> {
> z.1<-P1[i]*-1-P2[i]*-1
> if (all(z.1<=0)) {r=1} else {if (all(z.1>=0)) {r=1} else {r=2}}
> if (r==1) {ADF<-ur.df(ts(z.1), lags=1, type='drift')}
> if (r==1) {summary(ADF)}

You may need to print() that summary-object inside the for-function.

 > summary(a)
    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
    1.00    1.75    2.50    2.50    3.25    4.00
 > sink("test.txt")
 > for(i in 1) summary(a)
 > sink()    # No test.txt file created
 > sink("test2.txt")
 > for(i in 1) print( summary(a) )
 > sink() # The expected file created

This relates to the FAQ about similar puzzling behavior with plotting  
lattice , grid or ggplot objects.

> }
> sink()
>
> Any suggestion of what I might be doing wroong?
>
> best regards,
>
> Sergio Ren?

David Winsemius, MD
West Hartford, CT