Vectorizing finding indexes

3 vues (au cours des 30 derniers jours)
Sargondjani
Sargondjani le 24 Juin 2012
hi, i have a vector X and a vector Y (both in ascending order). for every element in X i want the indexes of the values in Y between which this element lies. Actually, 1 index suffices (preferably the higher bound).
For example if have
Y=[0, 1, 2, 3, 4, 5]
X=[0.1, 2.5, 2.8, 4.1];
then I want to get as a result:
IND = [2 4 4 6]; %the higher bounds of the interval in Y in which the elements of X fall
I can do this with a for loop:
for ix=1:length(X);
IND(1,ix)=min(find(Y>X(ix)))
end
My question is whether it is possible to vectorize this, and how... Many thanks in advance!

Réponse acceptée

Walter Roberson
Walter Roberson le 24 Juin 2012
[counts, tIND] = histc( X, Y );
IND = tIND + 1; %to get the higher index
Note that this has questionable results in the case where an X is exactly equal to a Y: it will return the bin number of the next Y. You could, though, correct for this with
t = (X == Y(tIND));
IND(t) = IND(t) - 1;
It is also possible to correct for it at the time of histc(), but the code becomes notably more difficult to read.
  3 commentaires
Walter Roberson
Walter Roberson le 24 Juin 2012
While you are thinking about that, also think about what to do if an X is before or after all of the Y.
See the comments in http://www.mathworks.com/matlabcentral/answers/41814-interpolation-with-high-frequency-financial-data
Sargondjani
Sargondjani le 24 Juin 2012
that's not going to happen. i am doing value function iteration using Y as the grid (for the state variable), and i force my policy function, X, to be in/on the grid

Connectez-vous pour commenter.

Plus de réponses (3)

Jamie Rodgers
Jamie Rodgers le 24 Juin 2012
Try This: Your vector
Y=[0:1:1000];
X=[0.1, 2.5, 2.8, 4.1];
Vectorised code
Z1=arrayfun(@(a)intersect(Y,ceil(a)),X);
idx=arrayfun(@(x)find(Y==x),Z1);
  1 commentaire
Sargondjani
Sargondjani le 24 Juin 2012
this probably works best, because i would actually like to have index of lower bound and index of higher bound... this way i can get them both directly with an understandable method :-)

Connectez-vous pour commenter.


Andrei Bobrov
Andrei Bobrov le 24 Juin 2012
[idx,~] = find(bsxfun(@eq,ceil(X(:)'),Y(:)))
  1 commentaire
Sargondjani
Sargondjani le 24 Juin 2012
i really like this but if i change small thing i get results that i dont get!
if i remove the transpose (i would prefer the solution to have the same dimension as X and Y) then i get a weird result
and i get only ones if i remove the (:)
...so unfortunately this solution is probably out of my league

Connectez-vous pour commenter.


Sargondjani
Sargondjani le 24 Juin 2012
Hmmmmm, andrei and jamie: do your solutions only work when Y consists of integers??? In my example Y was just integers, but in fact Y is not integers...
when i tried your solutions with Y not being integers i got results that were not correct
  1 commentaire
Walter Roberson
Walter Roberson le 24 Juin 2012
You are correct: the expressions given by Jamie and Andrei as of the time of my writing this comment, only work for integer Y. My histc() based code does not depend on Y being integer.

Connectez-vous pour commenter.

Catégories

En savoir plus sur Logical dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by