Define new variables to each instance of a string of repetitive values

1 次查看(过去 30 天)
I would like to create a variable for each instance where I have a series of uninterrupted, repeated values. For example, one of my data columns is a behavioral code timeseries where 0 is when an animal is on the surface and 1 is when the animal is underwater (see below). I'd like to obtain the indices of the first dive (in this case (5:8)), the second dive (in this case (16:19)), the third dive, etc. so that I can plot each dive v. time elapsed, make these separate variables, and run analyses on each dive.
dive = [0
0
0
0
1
1
1
1
0
0
0
0
0
0
0
1
1
1
1
0
0
0]
Thanks in advance for your help.
  3 个评论
Jessica Kendall-Bar
I just found Jan's "RunLength" function, which seems like it would be a good solution, but I am a little confused on how to implement it.

请先登录,再进行评论。

采纳的回答

Matt Macaulay
Matt Macaulay 2018-6-28
编辑:Matt Macaulay 2018-6-28
The RunLength function on file exchange would be great for this. You could use it like this on the dive variable in the question:
[B, N, IB] = RunLength(dive);
% IB is an array of the dive start times (indexes)
IB = IB(B>0);
% N is an array of the dive lengths
N = N(B>0);
% Plot the dive start times against the dive lengths
plot(IB, N, 'kx')
xlabel('Start Time')
ylabel('Dive length')
  3 个评论
Stephen23
Stephen23 2018-6-28
"I've never downloaded and used a function before"
  • Click the big blue Download button.
  • Unzip it onto the MATLAB Search Path (e.g. the current directory).
If the submission contains MATLAB code then it is ready to use. If it contains some mex code (as this submssion does) then this must be compiled for your specific PC: Jan Simon normally writes some self-compiling code, that automagically compiles the first time it is called.
Jessica Kendall-Bar
Thanks, I had downloaded it, but hadn't included it in the current directory. I just tried it and it worked!

请先登录,再进行评论。

更多回答(1 个)

Walter Roberson
Walter Roberson 2018-6-28
mask = [false, dive.', false];
starts = strfind(mask, [0 1]);
stops = strfind(mask, [1 0]) - 1;
divetimes = arrayfun(@(startidx, stopidx) times(startidx:stopidx), 'uniform', 0);
Now starts is a vector of the indices at which each dive started, and stops is a vector of corresponding indices at which each dive stopped. It is assumed that before and after the recorded information that they are above water.
divetimes() will be a cell array of times extracted from your times vector -- for example you might use this if your recordings are not at uniform time intervals.

类别

Help CenterFile Exchange 中查找有关 Just for fun 的更多信息

产品


版本

R2017a

Community Treasure Hunt

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

Start Hunting!

Translated by