similar to: problem in reading trajectory file

Hi all, I need a little help for construct an state's map on R. The first problem is to get the data. I have a datafile of longitude and latitude in the follow format: trajectory latitude longtude T -22.045618 -51.287056 T -22.067078 -51.265888 T -22.067039 -51.207249 T -22.059690 -48.089695 T -22.075529 -48.074608 T -22.072460 -48.044472 T -22.062767 -48.298473 T -22.077349

Hello. I have 5 measurement points, my dependent variable is ordinal (0 - 3), and I want to visualize my data. I'm pretty new to R. What I want is to find out whether people with different baseline covariates have different trajectories, so I want a plot with the means trajectory of my dependent variable (the individual points do not make a lot of sense in ordinal data) on each measurement

dear helpers, I wonder if there is a way to read a molecular dynamic trajectory file ( binary file) produced by CHARMM into R. Something like that in matlab. Actually this will save tremendous effort in post processing. best regards karim

Hello, I'm attempting to create a smooth, 3D plot of a trajectory (rather than the cloud or wireframe functions). I would rather the individual data points not be visible. I've had no luck finding this on the graphics or help pages. Thank you in advance. Chris Some example data, just in case: 25 32 40 12 25 32 2 12 25 2 2 12 20 2 2 1 20 2 6 1 20 5 6 1 5 5 6

Hi there, I want to make trajectory plots for data as follows: ID time y 1 1 1.4 1 2 2.0 1 3 2.5 2 1.5 2.3 2 4 4.5 2 5.5 1.6 2 6 2.0 ... That is, I will plot a growth curve for each subject ID, with y in the y axis, and time in the x axis. I would like to have all growth curves in the same plot. Is there

Hi all, I am trying to plot the fitted trajectories for each individual from an individual growth model (fit with a linear mixed effects model in lme). How can I plot each person's trajectory in the *same* panel, along with the mean-level trajectory? Below is an image of a plot similar to what I'm trying to create (from: http://jpepsy.oxfordjournals.org/content/31/10/1002/F6.large.jpg):

>If you are looking at radioactive decay maybe but how often do >you actually see exponential KM curves in real life? Exponential curves are rare. But proportional hazards does not imply exponential. > A trial design could in fact try to get all the control sample to "event" at the same time if enough was known about prognostic factors and natural trajectory You are a

I'm trying to understand how to plot individual growth curve trajectories, with the overall mean trajectory superimposed (preferably in a slightly thicker line, maybe in black) over the individual trajectories. Using the sleepstudy data in lme4, here is the code I have so far: library(lme4) library(lattice) xyplot(Reaction ~ Days, data = sleepstudy, group = Subject, type = 'l')

Hi all, I have a string vector like that: x=c("1\t\t", "2", "3\t\t\t") I need to remove all the occurrences of "\t", in order to obtain: x = "1" "2" "3" I'm trying to use the function substring2, and it works for each component, for example: substring2(x[1], "\t") ="" gives x = "1"

Hi all, I have a main file main.R in which I include some other R files. For example, in main.R I have: source("functions.R"). When I modify the file "functions.R", I'd like R to take into account the changes and to reload the file functions.R when I run main.R Is it possible? Thanks, Arnaud [[alternative HTML version deleted]]

Hi, I need to select a substring from the filename of a file in a list. I can find all the filenames in a directory by list.files(fn_dir) Each filename has the format (the following is an example because field lengths are variable) XXX184_YYY_ZZZ.dat XXX: 2 or 3 alphabetic characters 184: 3 numeric characters ... I'd like to select the substring XXX which has variable length (from each

Hi all, I'd want to plot a segment from a line specified by slope and intercept. I want to plot this line between two limits, x1 and x2, without imposing these limits to the hole plot as it is the case with scale_x_continuous(limits=c(x1,x2)) or with xlim and ylim. Any idea? Arnaud Chozo [[alternative HTML version deleted]]

Dear all, Fernando Colchero, Owen Jones and Maren Rebke, Max Plank Institute for Demographic Research, present BaSTA -Bayesian Survival Trajectory Analysis. Fernando, Owen and Maren have put together this beautiful video exploring research on ageing and and how to deal with incomplete data. To install BaSTA http://basta.r-forge.r-project.org BaSTA users mailing list,

Hi all, I'd like to use the Hmisc::summarize function, but it uses a function (FUN) of a single vector argument to create the statistical summaries. Consider an easy case: I'd like to compute the correlation between two variables in my dataframe, grouped according to other variables in the same dataframe. For exemple, consider the following dataframe D: V1 V2 V3 A 1 -1 A 1

I am trying to simulate the trajectory of the pension assets of one person. In C-like syntax, it looks like this: daily.wage.growth = 1.001 # deterministic contribution.rate = 0.08 # deterministic 8% Wage = 10 # initial Asset = 0 # initial for (10,000 days) { Asset += contribution.rate * Wage

Hi all, I'd want to plot a table of numbers such that the values are represented by gray level. Is there an easy way to do that using ggplot2? Thanks, Arnaud [[alternative HTML version deleted]]

Hello, Does anyone know how to make "movies" in R by making a sequence of plots? I'd like to animate a long trajectory for exploratory purposes only, without creating a bunch of image files and then using another program to string them together. In Splus I would do this using double.buffer() to eliminate the flickering caused by replotting. For instance, with a 2-D trajectory

Hello, I'm using R code that includes a residual permutation that was written as a supplement to the paper: Turner et al. 2010. A general hypothesis-testing framework for stable isotopes ratios in ecological studies. Ecology 91:2227-2233. The supplemental code is available at: http://www.esapubs.org/archive/ecol/E091/157/suppl-1.htm When I execute the function, no warnings are given

Hi the list, I am defining S4 objet. Is it possbile to define a method that change the slot of an object without using <- ? My object contain a numeric and a matrix. At some point, I would like to impute the missing value in the matrix. So I would like to use something like : ----------------- setClass("MyObj",representation(time="numeric",traj="matrix")) a

Hello, I'm currently trying to find a method to interpolate or smooth data that represent a trajectory in space. For example, I have an ordered (=time) set of (x,y) tuples which constitute a path in a 2D space. Is there a way using R to interpolate between these points in a way similar to spline interpolation so that I get a smooth path in space? Greetings, Dieter -- Dieter Vanderelst