[R] Get the difference of values to their own median value

Felix Wave felix-wave at vr-web.de
Wed May 2 09:44:16 CEST 2007


Hello,
I've got a matrix (mail end) with the colnames x, y, z. In this matrix
are different measurements. x and y are risign coordinates.

With the following line I got the median value of z for all "x" AND "y" witch
are the same (not every measurment in my list hast the same number of 
"x" and "y" values. Sometimes lines are missing.
>MEDIAN  <- na.omit( aggregate(INPUT[,3], by=list(INPUT[,2],INPUT[,1] ), FUN=median ) )


To see the failin of my measurments I want to get the difference between
every measurment and the detected median values.
The aim would be to have a list with different colums (for every measurment)
with the difference between measurment 1,2,3, ..., x and the median.

Has anybody an idea?

I though I could split the measurment list. The beginning of every measurment 
can be found with
>START     <- grep(" 0 0.0 0.00000000", INPUT) 
But not every measurment has the same length and I don't have always the
same number of measurments.



I hope you can help me.

Thank's a lot.
Felix




###############
## My R Code ##
###############
INPUT           <- readLines(dat.dat)
INPUT           <- gsub("^ ", "", INPUT)
INPUT           <- t( sapply( strsplit(INPUT, split=" "), as.numeric ) )
colnames(INPUT) <- c("x", "y", "z" )






#############
## dat.dat ##
#############
29 4.5 1.505713
29 4.6 1.580402
29 4.7 1.656875
29 4.8 1.735054
30 0 0
30 0.1 0.00096108
30 0.2 0.00323831
29 4.5 1.495148
29 4.6 1.568961
29 4.7 1.644467
30 0 0
30 0.1 0.00093699
30 0.2 0.00319411
30 0.3 0.00676619"



More information about the R-help mailing list