weibull

    2熱度

    2回答

    我想在給定間隔內生成一些威布爾隨機數。例如,來自威布爾分佈的具有形狀2和在區間(0,10)中的比例30的20個隨機數。 rweibull R中的函數產生具有給定形狀和比例值的威布爾分佈的隨機數。有人可以提出一種方法嗎?先謝謝你。

    1熱度

    1回答

    我有一些簡單的事件數據時間,沒有協變量。我試圖將威布爾分佈擬合到它。所以我有以下代碼。一切都很好,直到我加載我的首字母縮寫。它說「這個鏈包含未初始化的變量」。但我不明白。我認爲Weibull dist只有2個參數,我已經指定了它們。您能否提一些建議?謝謝! model { for(i in 1 : N) { t[i] ~ dweib(r, mu)I(t.cen[

    3熱度

    2回答

    我想估計3p威布爾分佈的尺度,形狀和閾值參數。 我到目前爲止已經做的是以下幾點: 指的這個帖子,Fitting a 3 parameter Weibull distribution in R 我用過的功能 EPS = sqrt(.Machine$double.eps) # "epsilon" for very small numbers llik.weibull <- function(sha

    2熱度

    1回答

    使用SciPy,我試圖從this question重現威布爾擬合。 import numpy as np from scipy.stats import genextreme import matplotlib.pyplot as plt data=np.array([37.50,46.79,48.30,46.04,43.40,39.25,38.49,49.51,40.38,36.98,4

    3熱度

    1回答

    我有一個數據集,我想適合已知的概率分佈。目的是在數據生成器中使用擬合的PDF - 這樣我就可以從已知的(合適的)PDF中採樣數據。數據將用於模擬目的。目前我只是從正態分佈抽樣,與實際數據不一致,因此模擬結果不準確。 我第一次想用下面的方法: Fitting empirical distribution to theoretical ones with Scipy (Python)? 我首先想到的是

    0熱度

    1回答

    我嘗試用optim()函數估計三個參數a,b0和b1。但我總是得到錯誤: 優化錯誤(par = c(1,1,1),fn = logweibull,method =「L-BFGS-B」,: L-BFGS-B需要有限值'fn 「 t<-c(6,6,6,6,7,9,10,10,11,13,16,17,19,20,22,23,25,32,32,34,35,1,1,2,2,3,4,4,5,5,8,8,8,8

    2熱度

    2回答

    我想創建一個QQ圖來測試,如果我的數據可以由Weibull分佈使用命令 建模使用數據 x =c(3.367, 0.769,0.8,1,1.2) qqplot(x,'weibull') 我不斷收到帶有與錯誤 "In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion" ,並找不出原因。這是否意味着我無法將We

    0熱度

    1回答

    我想通過在R 2與一個給定的審查矢量施加到下面的數據來估計威布爾分佈的極大似然參數: 數據= 9 2 11 49 7 5 3 36 30 6 62 5 3 29 29 1 13 1 24 11 9 4 7 15 11 15 1 1 1 1 1 2 6 12 12 28 14 14 57 17 4 2 3 6 21 6 16 19 28 18 19 9 59 12 3 27 8 26 19 47 6

    2熱度

    1回答

    我試圖找出一個隨機變量事件超過特定值的概率,即pr(x> a),其中a是某個常數,通常遠高於x的平均值,並且x不是任何標準的高斯分佈。所以我想要擬合一些其他的概率密度函數,並把x的PDF從a到inf的積分。由於這是尖峯建模的問題,我認爲這是一個極值分析問題,並且發現威布爾分佈可能是合適的。 關於極值分佈,威布爾分佈有一個非常「不易實現」的積分,因此我想我可以從Scipy得到pdf,並做一個黎曼和。

    3熱度

    1回答

    我正在研究一個涉及組合兩個威布爾分佈並因此創建雙峯值曲線的項目。然後,我的目標是使用這個預測。我在網上搜索,我似乎無法找到任何東西,或者如果R有一個函數,允許我結合兩個Weibulls。 下面顯示了我用來創建兩個Weibull分佈的代碼,我希望這兩個分佈組合起來構成一個概率密度函數。 curve(dweibull(x, scale=30.59898985, shape=2.27136646),fr