[R] source output differs from console output

[graham lawrence]

Dear R-help,
I would like to be able to run the following code sequence as a source 
routine.
If I paste it into R via the clipboard it works as expected, but if I source 
the code instead then the last 3 statements fail.
I've also tried writing to the file in place of the sink sequence, but that 
also hits a snag.  R 1.6, w98e2, dfr is a data frame containing the content 
of Julian Faraway's gala data frame.  I append the console log and 3 outputs 
to the file.
Thanks for your help.

qnum<-""
if(given[1]!="NA")qnum<-paste("\n",given[1],".",sep = "")
xfr<-dfr
if(given[2]!="NA")
{
x<-as.numeric(unlist(strsplit(given[2],",")))
xfr<-cbind(dfr[,x])
}
varlst<-names(xfr)
mr<-lm(eval(parse(text=paste(varlst[1],"~",paste(varlst[-1],collapse="+")))), 
data=xfr)
write(qnum,"math.txt",append = TRUE)
sink("math.txt",append = TRUE)
summary(mr)
sink()

The console log was as follows

>given
[1] "1"       "1,3,5,7"
>qnum<-""
>if(given[1]!="NA")qnum<-paste("\n",given[1],".",sep = "")
>xfr<-gala
>if(given[2]!="NA")
+  {
+  x<-as.numeric(unlist(strsplit(given[2],",")))
+  xfr<-cbind(gala[,x])
+  }
>varlst<-names(xfr)
>mr<-lm(eval(parse(text=paste(varlst[1],"~",paste(varlst[-1],collapse="+")))), 
>data=xfr)
>write(qnum,"math.txt",append = TRUE)
>sink("math.txt",append = TRUE)
>summary(mr)
>sink()
>write("now i run it as a source via mrfit.R","math.txt",append=TRUE)
>source("../sources/mrfit.R")
>write("now i run it with a write statement instead of the sink 
>sequence","math.txt",append=TRUE)
>qnum<-""
>if(given[1]!="NA")qnum<-paste("\n",given[1],".",sep = "")
>xfr<-dfr
>if(given[2]!="NA")
+  {
+  x<-as.numeric(unlist(strsplit(given[2],",")))
+  xfr<-cbind(dfr[,x])
+  }
>varlst<-names(xfr)
>mr<-lm(eval(parse(text=paste(varlst[1],"~",paste(varlst[-1],collapse="+")))), 
>data=xfr)
>write(qnum,"math.txt",append = TRUE)
>write(summary(mr),"math.txt",append = TRUE)
Error in cat(list(...), file, sep, fill, labels, append) :
        argument 1 not yet handled by cat
>l<-summary(mr)
>write(l,"math.txt",append = TRUE)
Error in cat(list(...), file, sep, fill, labels, append) :
        argument 1 not yet handled by cat
>l

Call:
lm(formula = eval(parse(text = paste(varlst[1], "~", paste(varlst[-1],
    collapse = "+")))), data = xfr)

Residuals:
    Min      1Q  Median      3Q     Max
-103.66  -53.31  -35.18   23.30  305.12

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept) 62.15560   22.98119   2.705 0.011901 *
Area         0.08465    0.02072   4.085 0.000375 ***
Nearest      0.37841    1.24235   0.305 0.763101
Adjacent    -0.01104    0.02072  -0.533 0.598893
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

Residual standard error: 94.45 on 26 degrees of freedom
Multiple R-Squared: 0.3914,     Adjusted R-squared: 0.3212
F-statistic: 5.573 on 3 and 26 DF,  p-value: 0.004334

The 3 outputs to the file math.txt follow

1.

Call:
lm(formula = eval(parse(text = paste(varlst[1], "~", paste(varlst[-1],
    collapse = "+")))), data = xfr)

Residuals:
    Min      1Q  Median      3Q     Max
-103.66  -53.31  -35.18   23.30  305.12

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept) 62.15560   22.98119   2.705 0.011901 *
Area         0.08465    0.02072   4.085 0.000375 ***
Nearest      0.37841    1.24235   0.305 0.763101
Adjacent    -0.01104    0.02072  -0.533 0.598893
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1

Residual standard error: 94.45 on 26 degrees of freedom
Multiple R-Squared: 0.3914,	Adjusted R-squared: 0.3212
F-statistic: 5.573 on 3 and 26 DF,  p-value: 0.004334

now i run it as a source via mrfit.R

1.
now i run it with a write statement instead of the sink sequence

1.






_________________________________________________________________
Internet access plans that fit your lifestyle -- join MSN. 


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._