Arnaud -
A reproducible example would be nice, but something
like this might work, assuming your data is in a file
called 'traj.txt':
dat = readLines('traj.txt')
triallines = grep('trial',dat)
starts = triallines + 2
ends = c(triallines[-1] - 1,length(dat))
thetrials = sub('# *trial *(\\d+)','\\1',dat[triallines],perl=TRUE)
pieces = mapply(function(start,end,trial)
data.frame(trial,read.table(textConnection(dat[start:end]))),
starts,ends,thetrials,SIMPLIFY=FALSE)
answer = do.call(rbind,pieces)
names(answer) = c('trial','t','X','Y')
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
spec...@stat.berkeley.edu
On Mon, 15 Mar 2010, arnaud chozo wrote:
Hi,
I'm trying to read some trajectory files (text files) which have the form:
# trial n
# t X Y
0 1 2
0.2 1 3
0.4 1.2 4
...
# trial n+1
# t X Y
0 1 2
0.2 1.3 3.3
0.4 1.5 5
...
...
where the symbol # means that the line is not a numeric value, but I still
need to read the number of the trial (n=1, 2, etc.).
I wonder if there is an easy way to read this file in order to plot the
trajectory for each trial.
Thanks,
Arnaud
[[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.
______________________________________________
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.
