Outpt with zeros when should be values

Question

Krispy Scripts 2016-12-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/315002-outpt-with-zeros-when-should-be-values

评论： Krispy Scripts 2016-12-6

intervals.mat

I have a matrix of two columns. The first column is cue times, as in times an animal was cued to do something. The second column is times that the animal responded. If the animal correctly responded then the animal will respond within a range greater than .15 seconds of cue and less then 5 seconds from cue. If so this is a correct behavior and needs to go into one matrix(correcttrials), if not then this is an error behavior and needs to go in a second matrix(errortrials).

The original Matrix (I have also attached it):

intervals = 
4358          735.1798
0970  1442.0839
2985  1682.2252
41885  1961.58675
95985  2507.14965
20085  2681.83055
4032  2682.0705
0439  2777.6710
62435  3498.9947
4852  0
38675  0
0074  0
6089  0
9099  0
73055  0
63135  0
3324  0
97285  0
3738  0
87465    0

This is the code I have written, but I am not sure why the correct output has zeros in it and the error output gets all zeros.

Here is my code:

data=intervals;
cue=data(:,1);
response=data(:,2);
 for k=1:length(cue)
    for i=1:length(response)
        if reward(i)>= (cue(k)+.15) && reward(i)<= (cue(k)+5)
            correcttrials(i)=response(i); 
        else
            errortrials(i)=response(i);
        end
    end
end

Here are the Matrices I am getting out:

Correct Trials:

correcttrials=
735.1798  1442.0839  1682.2252  1961.58675  735.1798  0  2682.0705  1442.0839  0  1682.2252  1961.58675  0  0  0  2682.0705

Error Trials:

errortrials = 
0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0

I am intending to get:

Correct Trials:

correcttrials = 
735.1798  1442.0839  1682.2252  1961.58675  2682.0705

Error Trials:

errortrials = 
2507.14965  2681.83055  2777.6710  3498.9947

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Image Analyst 2016-12-1

We can't run this code because you have not defined the reward array.

请先登录，再进行评论。

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

Answer 1

Andrei Bobrov 2016-12-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/315002-outpt-with-zeros-when-should-be-values#answer_245601

在 MATLAB Online 中打开

intervals = [...
4358          735.1798
0970  1442.0839
2985  1682.2252
41885  1961.58675
95985  2507.14965
20085  2681.83055
4032  2682.0705
0439  2777.6710
62435  3498.9947
4852  0
38675  0
0074  0
6089  0
9099  0
73055  0
63135  0
3324  0
97285  0
3738  0
87465    0];
d = intervals;
dn0 = d(d(:,2)~=0,2)';
da = abs(d(:,1) - dn0);
t = any(da > .15 & da < 15);
correcttrials = dn0(t)
errortrials = dn0(~t)

2 个评论
显示无隐藏无

Krispy Scripts 2016-12-2

Thank you!

Krispy Scripts 2016-12-6

If I wanted to instead get the values in column one that meet the criteria instead of column 2 could that be easily done?

请先登录，再进行评论。

Answer 2

Image Analyst 2016-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/315002-outpt-with-zeros-when-should-be-values#answer_245590

在 MATLAB Online 中打开

Try this vectorized way:

load('intervals.mat');
cueTimes = intervals(:,1); % Extract the first column.
responseTimes = intervals(:,2); % Extract the second column.
% Get the time differences.
timeDifferences = responseTimes - cueTimes
% Convert from milliseconds to seconds. ??????????
timeDifferences = timeDifferences / 1000;
correctTrialsRows = (timeDifferences > 0.15) & (timeDifferences < 5)
correctTrials = responseTimes(correctTrialsRows)
errorTrials = responseTimes(~correctTrialsRows)

None of the things was correct so I figured the times might be in milliseconds whereas your criteria was in seconds, so I converted. I'm still not sure if this gives you what you want though.

2 个评论
显示无隐藏无

Krispy Scripts 2016-12-2

It does not need to be converted into milliseconds. It also cannot be assumed that the value in one column corresponds to another. So thats why I did not look to get a time difference.

My method was to take each value in the first column, for your code cueTimes, and read through the entire second column. If any value is within the cutTime window, greater than .15 and less than 5 (seconds), than put that second column responeTime in correctTrials. If it does, not then it should go in the errorTrials.

Does that help at all? I am not sure about how to write for this though..

Image Analyst 2016-12-2

编辑：Image Analyst 2016-12-2

Then none of your time differences are in the range 0.15 to 5 seconds. Let's take a typical start/stop pair - the pair in row #5. Start/cue time of 734.95985 seconds, and response (stop time) at 2507.14965 seconds. That's 1771 seconds. And all the rest are like that. Not one of them is in the range 0.15 to 5 - they're all hundreds.

请先登录，再进行评论。

Outpt with zeros when should be values

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

2 个评论
显示无隐藏无

更多回答（1 个）

2 个评论
显示无隐藏无

另请参阅

类别

标签

Community Treasure Hunt

Outpt with zeros when should be values

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

2 个评论 显示 无隐藏 无

更多回答（1 个）

2 个评论 显示 无隐藏 无

另请参阅

类别

标签

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

2 个评论
显示无隐藏无

2 个评论
显示无隐藏无