在MATLAB中，给定一个关于值p的正态分布，如何才能看到p大于值p*的概率？

Question

我知道标题很模糊，老实说，我很难把我的问题用语言表达出来。我想说的是，给定一个初值，它是关于p的正态分布，我如何才能找到p>p*的概率，其中p*是一个临界值。

也就是说，我取p= 1g (我正在研究药物，所以它是1克)，并假设这是通过培养的细胞正常分布的。我知道细胞中含有超过一定数量的药物会杀死细胞，称之为p*。我怎么能说有多少细胞会因为药物的数量大于p*而死亡呢？

我希望这是合理的。

非常感谢。

Answer 1

你只是在寻找of或累积密度函数，它只是PDF或概率分布函数的积分。所以现在它取决于你是收集数据还是模拟数据。如果您正在收集数据，那么您将需要ecdf函数。如果您正在模拟数据，即您已经为您的分布选择了参数，那么您可以使用cdf函数。不管是哪种方式，这都是一个简单的问题，就是从这些函数的输出中读取您的观点。

编辑：

简单地更正一下上面的内容，CDF是P(X <= x)，您需要P(X >= x)，所以转换应该是很简单的。要么把你的数据翻过来，要么是1-CFD之类的。

Answer 2

您需要知道分布的平均值和标准差，然后使用normcdf (统计工具箱)：

mu = 1; %// your p value is the distribution mean, I gather
sigma = .2; %// you don't specify sigma in your question
x = 1.3; %// your p* value
result = 1-normcdf(x,mu,sigma);

如果您没有统计工具箱，可以使用erfc，这是一个标准的Matlab函数：

result = .5*erfc((x-mu)/sigma/sqrt(2));