绘制置信区间信号的平均值答案

【问题标题】：plot mean value of a signal for confidence interval绘制置信区间信号的平均值
【发布时间】：2019-01-18 04:48:46
【问题描述】：

我在 y 轴（偏航率）上有一个信号，在 x 轴（时间）上有一个时间信号。我试图找到这个信号的置信区间。但我无法绘制信号的平均值。我认为因为平均值有 2 个向量，但信号大小为 5001。我该如何解决？我感谢您的回答！

N = size(yawrate,1);  % Number of ‘Experiments’ In Data Set

yMean = mean(yawrate);  % Mean Of All Experiments At Each Value Of ‘time’`

ySEM = std(yawrate)/sqrt(N);  % Compute ‘Standard Error Of The Mean’ Of All Experiments At Each Value Of ‘time’

CI95 = tinv([0.025  0.975], N-1);  % t-score Calculate 95% Probability Intervals Of t-Distribution

yCI95 = bsxfun(@times, ySEM, CI95(:));  % Calculate 95% Confidence Intervals Of All Experiments At Each Value Of ‘x’         

figure

plot(time, yMean)   % Plot Mean Of All Experiments

hold on

plot(time, yCI95+yMean)  % Plot 95% Confidence Intervals Of All Experiments

hold off

grid

【问题讨论】：

yawrate 和 time 的大小是多少？ “平均值有 2 个向量”是什么意思？
非常感谢您的回答。偏航率和时间的大小是 5001。我的意思是平均值的大小显示 1 1。
代码是为yawrate 编写的，它是一个NxM 矩阵，具有M 时间步和N 实验。在这种情况下，mean(yawrate) 是一个具有M 值的向量。 time 也需要有 M 值。如果你的N 为 1，则这段代码没有用。
不，我的 N 是 5001。
不，你的 M 是 5001，你的 N 是 1。

标签： matlab plot mean confidence-interval

【解决方案1】：

我会发表评论，但我没有足够的声誉。由于我在回答，我还将使用 cmets 中已有的内容。关于 cmets 中已经指出的事情，这个问题有点不清楚：对于某些 M，size(yawrate) 是 [1,N] 或 [N,M]？同样的问题代表time。如果我从您的回答和您的代码中正确理解，您假设 yawrate 是 NxM 和 yMean 是大小 1xM。在这种情况下，您也需要时间成为1xM。这是 Cris Luengo 写的。

另外，yCI95+yMean 意味着将包含置信区间的2x1 向量添加到1xM 值或yMean。这会给你一个2xM 数组。这在每个 Matlab 版本中都不起作用，我认为它需要第一个是 2x1 而不是 1x2，所以你可能需要这样做

plot(time, yCI95(1)+yMean)
plot(time, yCI95(2)+yMean)

再次，假设我理解您的代码。检查您所做的尺寸假设。

【讨论】：

谢谢马科斯。其实我在问，你问我什么:) 我不明白置信区间。我有信号和信号图，就像这张照片一样。大小（偏航率）和大小（时间）是 [5001 1]，但平均值是数字和大小 [1 1]。如果我绘制（时间，平均值），它肯定是一条线。在这种情况下，我如何创建置信区间的上下边界？这就是重点，我不明白。

【解决方案2】：

yMean 似乎是一个单一的值，所以当你绘制它时，你会为每个时间值得到一个点。如果您的采样时间过长，您可以通过在绘图函数中使用 t=0:0.01:time(end) 而不是时间来增加线条的可见性。

【讨论】：