On Mar 28, 2010, at 12:51 PM, Muhammad Rahiz wrote:
Dear Jannis,
Thanks for the tip. It works but I'd like to improve on the way I
did it.
x <- array(1:50,dim=c(10,10)) # data
h1 <- c("ncols 180") # header, line 1
h2 <- c("nrows 290") # header, line 2
h3 <- c("xllcorner -200000") # header, line 3
h4 <- c("yllcorner -200000") # header, line 4
h5 <- c("cellsize 5000") # header, line 5
h6 <- c("NODATA_value NA") # header, line 6
h0 <- rbind(h1,h2,h3,h4,h5,h6)
write(ho,"file.txt")
write.table(x,"file.txt",row.names=FALSE,col.names=FALSE,append=TRUE)
The header should display as;
ncols 180
nrows 290
xllcorner -200000
yllcorner -200000
cellsize 5000
NODATA_value NA
Executing that code appears to deliver what you requested. So rather
than asking us to read you mind, can you explain what is missing?
ncols 180
nrows 290
xllcorner -200000
yllcorner -200000
cellsize 5000
NODATA_value NA
1 11 21 31 41 1 11 21 31 41
2 12 22 32 42 2 12 22 32 42
3 13 23 33 43 3 13 23 33 43
4 14 24 34 44 4 14 24 34 44
5 15 25 35 45 5 15 25 35 45
6 16 26 36 46 6 16 26 36 46
7 17 27 37 47 7 17 27 37 47
8 18 28 38 48 8 18 28 38 48
9 19 29 39 49 9 19 29 39 49
10 20 30 40 50 10 20 30 40 50
--
David.
Jannis wrote:
Hallo,
have a look at the documentation of the different output functions
(write(),cat(),write.table()). Most of them have an argument like
"append". If you set this to TRUE, your command adds the lines you
are
writing into the file to the lines already existent in that file.
So if
you run this write() command first with your additional data (header
data) and then with the "real" data and append=TRUE, your information
ends up in the same file.
HTH
Jannis
Muhammad Rahiz schrieb:
Hello all,
I'd like some advise on this. When I read my files, I pass the
argument, skip=6, to skip 6 lines of header information. After
performing the desired calculations, I have the output. Now I want
to copy the 6 lines of skipped information to the output.
What I've been doing so far is to use a shell script (i.e. the cat
function) to copy the header information to the output file. This
means I need 2 additional files, namely a file containing the
header information only, and the shell script to run the the
function.
Is there a direct way to perform this in R?
The following forms the main component of my R script;
x <- as.matrix(read.table("file.txt",skip=6))
y <- x*2
write.table(y,"output.txt",row.names=FALSE,col.names=FALSE)
Thanks
Muhammad
______________________________________________
R-help@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@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@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.
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@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.
