[R] Trying to "pretty up" output from R job

Paul E Johnson pauljohn at ukans.edu
Sun Aug 6 22:47:19 CEST 2000


Running R 1.1 on RedHat Linux 6.2.

I need to write a shell script that goes through a bunch of directories
of simulation output, creating summary files that have the mean and
standard deviation of the variables found in the data files in each
directory.  I've got the R code doing almost the right thing. It reads
in data, then gets the mean and standard deviation for the numeric
variables, puts those results together, and then grabs the directory
name from the system and jams that in to the results as well. Here is
what I have:

data<-read.table("lastline.txt",header=T,as.is = TRUE)
indices<-1:dim(data)[2]
indices<-na.omit(ifelse(indices*sapply(data,is.numeric),indices,NA))
mean<-sapply(data[,indices],mean)
sd<-sapply(data[,indices],sd)
dir<-system("pwd",TRUE)
newOutput<-rbind(mean,sd)
transNewOutput<-t(newOutput)
transNewOutput<-rbind(dir,transNewOutput)
write.table(transNewOutput,file="Atestoutput",quote=FALSE, sep="\t")

The output file looks like this:

mean    sd
dir     /home/pauljohn/swarm/src        /home/pauljohn/swarm/src
T       732.1   148.618520536522
Run     49.5    29.011491975882
Seed    499366074.13    296739559.579976
IntAct  43.85   5.592879672432
Change  0       0

Now the questions:
1. I want 3 columns of output, but here the first line only has two
items.  I can't see any way to make it put out a lable for the first
column, like "Variable" or some such.

2. Can you tell me how to truncate the numerical printout so it only
shows two numbers after the decimal point?

3. Is my way of keeping track of the directory completely ridiculous? 
How else would you do it?

4. My next chore is to find how to make this thing run through a series
of directories, summarizing each.  Would you recommend I do this with a
shell script or something within R itself?
-- 
Paul E. Johnson                       email: pauljohn at ukans.edu
Dept. of Political Science            http://lark.cc.ukans.edu/~pauljohn
University of Kansas                  Office: (785) 864-9086
Lawrence, Kansas 66045                FAX: (785) 864-5700
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list