On Aug 19, 2012, at 4:34 PM, White, William Patrick wrote:
Also it occurred to me that my initial explanation was not
explicitly clear as to what the desired output is. What I am trying
to get is a moving absolute deviation between the two sets of numbers.
The phrase "a moving absolute deviation" admits of several
interpretations. I suggest you post the correct answer for some simple
cases or that you be more mathematical in your description (as is
suggested in the Posting Guide.)
set.seed(123)
X <- sample(-5:5, 10)
Y <- sample(-5:5, 10);
> X
[1] -2 2 5 4 1 -5 -3 3 -4 0
> Y
[1] 5 -1 1 4 -5 0 -4 3 -2 2
> abs( tail(X,9) - head(Y,9) )
[1] 3 6 3 3 0 3 7 7 2
So this is c( abs(X[2] -Y[1]), abs( X[3]-Y[2], .....)
This is not to be confused with the mean absolute deviation, or the
median absolute deviation which are both something different and not
what i am after.
Again. Not a clear description (of what you do do not want), given
that the problem involves two vectors.
________________________________________
From: David Winsemius [dwinsem...@comcast.net]
Sent: Sunday, August 19, 2012 4:03 PM
To: David Winsemius
Cc: White, William Patrick; r-help@r-project.org
Subject: Re: [R] moving distance between two sets of data
On Aug 19, 2012, at 12:58 PM, David Winsemius wrote:
On Aug 19, 2012, at 12:04 PM, White, William Patrick wrote:
On the surface this seems pretty simple, but I flummoxed. I have
two sets of numbers they bounce around zero, positive one and
negative one. They have a relationship between them, where one
diverges away from the other. I want create a second set of numbers
that tracks that divergence.
#Lets make some data like mine, kinda
Firstset <- runif(100, min = -1 , max =1)
Secondset <- runif(100, min = -1 , max =1)
#So something like:
Divergence <- abs (Firstset - Secondset)
#but this doesn't work because when Firstset is at .5 and Secondset
is at -.25 it returns .25 instead of .75
abs( .5 - (-.25) ) should NOT return .25 so you need to produce a
better example or point to specifics in the example you offered. If
what you wanting what you are getting, then use set.seed(123) and
refer to specific values.
I meant to write: "If you are not getting what you are wanting .... "
abs( .5 - (-.25) )
[1] 0.75
--
David.
#another possibility is:
Divergence <- abs (Firstset) - abs (Secondset)
#but when Firstset is at .5 and Secondset is at -.5 it returns 0
instead of 1
#It seems like there is a better way to do this. Any ideas?
______________________________________________
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.
David Winsemius, MD
Alameda, CA, USA
______________________________________________
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.
David Winsemius, MD
Alameda, CA, USA
David Winsemius, MD
Alameda, CA, USA
______________________________________________
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.
