2014-02-28 49 views
1

我知道標題是模糊的,誠實地說,我很難將我的問題放在單詞上。我想說,給定一個初始值,這是一個關於p的正態分佈,我怎麼能找到p> p *的概率,其中p *是一個臨界值。在MATLAB中,給定關於值p的正態分佈,我怎麼能看到p大於p *的概率?

I.e.我拿p = 1克(我正在用藥物,所以它的1克),並假設這通常通過細胞在文化中分配。我知道在細胞中有一定量的藥物會殺死細胞,稱之爲p *。我怎麼能說,有多少細胞會因爲其中的藥物量大於p *而死亡?

我希望這是有道理的。

非常感謝。

回答

0

您只是在尋找CDF或累積密度函數,它只是PDF或概率分佈函數的積分。那麼現在取決於您是否收集數據或模擬數據。如果你正在收集數據,那麼你會想要ecdf函數。如果您正在模擬數據,即您已爲您的分配選擇了參數,那麼您可以使用cdf函數。無論哪種方式,這將是一個簡單的事情,從這些函數的輸出讀取你的觀點。

編輯:

只是一個快速修正上述情況,CDF是P(X < = x)和你想P(X> = X),所以它應該是微不足道的轉換。要麼翻轉你的數據,要麼可能是1-CFD或類似的東西。

+0

完美!自從我做任何數據已經太久了,不敢相信我忘記了CDF的。 – user3365435

0

你需要知道的均值和您的分佈的標準偏差,然後用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)); 
相關問題