Matlab中的样本t检验

Question

在 Matlab 中，我必须生成两个具有分布的高斯随机样本：

1. N(100,5)
2. N(115,15)

然后运行模拟 1000 次，并计算从 2 到 50 的每个样本量在 1000 次中被拒绝的次数。

最后，我需要根据样本量绘制H0 的拒绝结果。

h = zeros(1000,1);
k = 0;

for i = 1:1000
    r1 = ();
    r2 = ();
    for j= 2:50
        r1(j-1)=normrnd(100,5,[1,j]);
        r2(j-1)=normrnd(110,15,[1,j]);
    end
    h(i)=ttest2(r1,r2)
    if h(i)==1
        k=k+1;
    end
    percent_of_rejection_Null=(k/1000*100);
end
%plot(h vs sample_size)

有人可以帮我更正上面的代码吗？

Answer 1

您的代码存在一些问题：

1. ri 是一个可变大小的向量，因此您需要将其附加到一个可以存储不同大小元素的结构 -> 元胞数组。
2. percent_of_rejection_Null 需要是一个向量来存储所有的值。
3. 您没有具体说明您要查找的图类型，我展示了条形图，但很容易调整为您需要的图。
4. h 的计算有问题，它的向量每次迭代都是 50，因此需要计算 1 的总和。

这段代码是：

h = zeros(1000,1);
percent_of_rejection_Null = zeros(1000,1);
k = 0;

for i = 1:1000
    r1 = {};
    r2 = {};
    for j= 2:50
        r1{j-1}=normrnd(100,5,[1,j]);
        r2{j-1}=normrnd(110,15,[1,j]);
    end
    h = cellfun(@ttest2, r1, r2);
    k = sum(h==1);
    percent_of_rejection_Null(i) = (k/1000*100);
end
figure;
plot(percent_of_rejection_Null)

结果图是：

如果输出不是您想要的，请编辑您的问题以便我进行调整。