How to calculate residence time?

Question

Akash Pandya 2019-3-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451041-how-to-calculate-residence-time

回答： Image Analyst 2019-3-19

Hi,

I have a time column in picoseconds and 442 distance columns measured in nanometres. I know how to extract the time points that meet the distance threshold I have set.

x = xlsread('10_1.xlsx')
x = x(x(:,2)<=0.5,:);

I want to calculate the total residence time for each of the columns. How can I do this?

Akash

2 个评论
显示无隐藏无

dpb 2019-3-19

Dunno...what's the working definition of "residence time" here?

Akash Pandya 2019-3-19

So the residence time refers to the time a small molecule is <=0.5 nanometers away from a protein. I want to find the total time where the distance (nm) meets this criteria.

Akash

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2019-3-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451041-how-to-calculate-residence-time#answer_366184

在 MATLAB Online 中打开

Use regionprops() and ask for PixelIdxList. Then use PixelIdxList(end) - PixelIdxList(1) for each region found. Something like

props = regionprops(vec, 'PixelIdxList');  % Vec is a 1-d list of binary (boolean/logical) data of whether x is above 0.5 or not.
% Now find the delta time for all regions found in this vector, vec.
for k = 1 : length(props)
    startingIndex = props(k).PixelIdxList(1);
    endingIndex   = props(k).PixelIdxList(end);
    deltaTime(k) = x(endingIndex) - x(startingIndex); % x is the time vector
end