威布爾檢查數據

Question

我目前正在使用與安裝日期和失效日期單位相對應的一些壽命數據。數據是現場數據，所以我確實有大量暫停（尚未出現故障的單元）。我想使用Scipy統計庫（將數據擬合到威布爾曲線並獲得分佈參數）對這些數據進行一些威布爾分析。我對Python和Scipy很陌生，所以我找不到一種方法將任何可用的威布爾分佈（dweibull，exponweibull，minweibull，maxweibull）中的暫停數據包含進來。有沒有簡單的方法來處理暫停？我不想重複這個輪子，但是我很難從我的數據中估計Weibull的參數。誰能幫我？

非常感謝！

Answer 1

如果我理解正確，那麼這需要使用刪失數據進行估計。

沒有一個scipy.stats.distribution會直接估計這種情況。你需要結合刪失觀測的非刪失和似然函數的似然函數。

您可以使用scipy.stats.distributions的pdf和cdf，或更好的sf，這兩個部分。然後，如果您也對參數估計的不確定性感興趣，則可以使用scipy優化來最小化負對數似然性，或者嘗試statsmodels中的GenericLikelihoodModel。