Lucas,
Although I can't answer all your questions, I can give some suggestions.
I will assume you know how many years of data you
have. Suppose it is 2004 through 2009.
In each input file, the columns must have names,
so I will assume there is a column named "Month".
I am not sure if you want all of your output in
one file, or if you want one output file for each
year.
I am going to assume you want a separate output file for each year.
I am going to assume that when you do the
calculations for each month within each year, the
results
include two variables named "r1" and "pvalue" (as in your example).
Here is some example R code:
## make sure that R is working in the directory that contains the files
## this can be done from a menu in the R window,
if you are using R in Windows or Macintosh,
## or by using the setwd() command
setwd( "directory where my files are" )
years <- 2004:2009
for (yr in years) {
infile <- paste(yr,'.csv',sep='')
outfile <- paste(yr,'.txt',sep='')
clima <- read.csv2(infile,nrows=7)
## open the output file for the current year
sink(outfile)
for (mn in unique(clima$Month)) {
## do calculations for the current month
cat('Year: ', yr, ' Month:', mn, ' R1:', r1,' Pvalue:', pvalue,'\n')
## end of loop through months (in the current year)
}
## close the output file.
sink()
## end of loop through years
}
If you want all the output, including all the
years, in a single file, then put the sink()
commands outside the loop. Like this:
years <- 2004:2009
## open the output file for writing
sink("MyOutput.txt")
for (yr in years) {
infile <- paste(yr,'.csv',sep='')
clima <- read.csv(infile,nrows=7)
for (mn in unique(clima$Month)) {
## do calculations for the current month
## some results are in variables named r1 and pvalue
cat('Year: ', yr, ' Month:', mn, ' R1:', r1,' Pvalue:', pvalue,'\n')
## end of loop through months (in the current year)
}
## end of loop through years
}
## close the output file.
sink()
A few comments:
The sink() command must be given a file name, not a directory name. You had
sink("directore where my text file is")
You can try it first without the sink() commands,
and see the results on your computer screen. This
is a good way to make sure they are correct.
Since you already know your file names, 2004.csv,
2005.csv, and so on, it is not necessary to use
the dir() command to get them. And actually,
better to use list.files(), not dir().
If you do use the list.files() or dir() command,
try not to use it repeatedly inside a loop. It
should be sufficient to it once outside the loop,
for example:
myfiles <- list.files("directory where my files are")
nfiles <- length(myfiles)
for (year in 1:nfiles) {
filename <- myfiles[[year]] ## and I'm
pretty sure it should be myfiles[year], not
[[year]]
}
Also, the way you did it assumes that the dir()
command lists the files in the correct order, and
it also assumes that there are no other files in
the directory.
If you want to run R in a different directory
than the data files are in, you can change a few
things.
datadir <- "directory my files are in"
infiles <- list.files( datadir, full.names=TRUE)
nfiles <- length(infiles)
for ( i in 1:nfiles ) {
clima <- read.csv2(infiles[1], nrow=7)
sink( paste('outfile',i,'.txt',sep='')
## loop over months
## do calculations
## do the cat() commands
sink()
}
At 1:10 PM +0000 9/30/09, Lucas Sevilla García wrote:
Hi R community,
First of all, I want to thank everybody to share
their time solving R questions, You are great.
Ok, for my questions, I've been looking for a
solutions by myself, in forums but I'm just a
little bit desesperate so I hope somebody can
help me. I have built a code to read files from
a directory. These files are named by a year
(2004.csv, 2005.csv,...). When the code reads
first file (2004.csv), inside this file, there
is information about precipitation of every
months and I calculate different variables like
R square adjusted, p value or formula fit to the
data from linear regression. The code do more
things but for my question, to explain what I
need, that part of the code is enough. I want to
export to the text file, year and month apart of
some other variables, something like this:
Year: 2004 Month: January R1: 0.98 Pvalue: 0.03 ...
Year: 2004 Month: February R1:0.78 Pvalue:0.12 ...
I've seen that I can use order sink() and cat(),
so I would put those orders in my code, like
this:
nfiles<- length(dir("directory where my files are")) #Count file number
for(year in 1:nfiles) #Read first file
{
filename<-dir()[[year]]
#take first file and read filename, so if year
is 1, then filename will be 2004, is year is 2,
filename will be 2005,...
clima<-read.csv2(filename, nrows=7) #open 2004.csv
So, if want to export year to my text file I would do
for(year in 1:nfiles) #Read first file
{
sink("directore where my text file is")
filename<-dir()[[year]] #take first file and
read filename, so if year is 1, then filename will be 2004, is year is
2, filename will be 2005,...
cat(" Year: ",filename)
sink()
clima<-read.csv2(filename, nrows=7) #open 2004.csv
And in my text file would read
Year: 2004
Now, I want to the same to months. (I have built
a for loop to read months inside for loop to
read years). When I import a csv file I get
something like this
Januray February ....
1 3.0 4.1
2 1.4 3.7
3 0.2 1.5
4 6.7 4.1
.
.
.
I can use commands like clima$Januray or
clima[[1]] but I just get precipitation values.
However, I am not able to get the header of the
column. If I would able to do that I could do
the same as for years and export those headers
to my text file. Does anyone know how I could do
that? or does anyone know another way to do what
I need? Would anyone use sink() and cat()
commands to create a summary text like the one I
need to do?. Probably my for loop is not the
best, I am still a beginner with R, and probably
there are some better forms to express in R what
I need but I am working alone so there is nobody
in person to help me so I apologize for my
simple questions. Thanks in advance.
Lucas
_________________________________________________________________
[[alternative HTML version deleted]]
______________________________________________
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.
--
--------------------------------------
Don MacQueen
Environmental Protection Department
Lawrence Livermore National Laboratory
Livermore, CA, USA
925-423-1062
______________________________________________
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.
