R 中具有 Beta 分布的 EM 算法
[英]EM algorithm with Beta Distribution in R
我想編寫 R 代碼來構建 dirichlet 混合 model。我用於 beta 分布的對數似然如下： (,)=(−1)ln¯+(−1)ln(1−)¯+lnΓ(+)−lnΓ()−lnΓ() 我需要幫助來初始化參數（alpha、beta 和為代碼中使用的 gamma 分布設置額外參數）。 這是我想要 ...
如何遍歷多列
[英]How to loop through many columns
我在 pandas dataframe 中有大約 88 列。 我正在嘗試應用一個為每列計算單個值的公式。 如何切換出每列的名稱，然后從等式中構建一個新的單行 dataframe？ 下面是方程（線性混合模型），它為每列產生一個值。 下面是列名稱的示例 ['OBJECTID'，'Benthic_Mo'， ...
點符號不適用於 MixtureExpectation？
[英]Dot-notation not working with MixtureExpectation?
這是我的第一篇文章，所以對於任何格式問題，我深表歉意。 我正在嘗試計算 Julia 中一組數字的期望值，假設概率分布是兩個 Beta 分布的混合。 使用下面的代碼給出了下面看到的錯誤 MethodError: 沒有方法匹配 *(::MixtureExpectation{Vector{Iterab ...
R 中的混合 copula
[英]mixture copula in R
我想用混合copula進行可靠性分析，現在，在朋友的幫助下，我已經完成了'RVMs_fitted'。現在我想進行概率積分變換（PIT），但是RVINEPIT的function不能使用，因為 RVINEPIT(data，RVM)，這個 RVM 不是 RVINEMATRIX 這是我的代碼： ...
設置 mixfitEM 圖的軸限制
[英]Set Axis Limits of mixfitEM plot
我想為mixR包的mixfit函數的輸出圖的 x 軸設置限制。 此輸出屬於mixfitEM類。 可重現的示例：首先我模擬了對數法線的混合。 然后我嘗試使用 mixR 中的mixR mixfit()將兩個對數法線擬合到該數據中。 在第一次嘗試中，我嘗試使用base包中的plot()函數； 但pl ...
在 R 中運行 JAGS model 時出現“無法為偏差設置跟蹤監視器沒有觀察到的隨機節點”錯誤
[英]"Failed to set trace monitor for deviance There are no observed stochastic nodes" error when running JAGS model in R
我正在嘗試讓 N-mixture model 在 JAGS 中運行，以獲得我當前所在的統計信息 class。但是，我不斷收到錯誤消息“無法為偏差設置跟蹤監視器，沒有觀察到的隨機節點”嘗試運行 model。目前我只是試圖獲得一個非常基本的 model，沒有運行協變量，以確保我的所有格式都正確，但仍然無 ...
如何在 R 中使用 qgcomp() 創建帶有 for 循環的矩陣
[英]How to create a matrix with for loop using qgcomp() in R
我總是使用 for 循環來創建線性和邏輯回歸 output 的矩陣，但是使用 qgcomp() 這樣做有困難。 如果有人有經驗或建議，我將不勝感激。 我更喜歡這種方法而不是函數，因為只需單擊一下。 輸出的矩陣僅包含“NA”，即使我刪除了 try catch 錯誤，它也會指出錯誤：$ 運算符對原子向量 ...
高斯混合：最小值必須小於或等於最大值
[英]Gaussian Mixture: minvalue must be less than or equal to maxvalue
在我的探索性分析中，我目前正在使用高斯混合 Model 來排除異常值，方法是為每部電話 model（總共 4 個獨特的模型）繪制輪廓 plot。 我正在使用這 2 個變量（ r_max ， b_max ）來檢測一個點是否是異常值。 這是我的數據集（180 行，4 列）： 這是我創建輪廓 plot 的 ...
高斯混合模型 - 獲得給定概率值的輪廓：Matlab
[英]Gaussian mixture model - get the contour of given probability value : Matlab
我需要確定適合數據的 GMM 的 99% 概率輪廓。 按照這個示例，我希望能夠指定要繪制的輪廓以及它們的 x、y。 因此，在下圖中，我希望能夠用黑色虛線“k”繪制 99%。 知道如何做到這一點嗎？ ...
Tensorflow 值錯誤：尺寸必須相等：LSTM+MDN
[英]Tensorflow ValueError: Dimensions must be equal: LSTM+MDN
我正在嘗試基於此實現（ https://www.katnoria.com/mdn/ ）使用 LSTM + 混合密度網絡進行下一個單詞預測 model。 輸入：300維詞向量*窗口大小（5）和21維數組（c）表示文檔的主題分布，用於訓練隱藏的初始狀態。 Output：混合系數*num_gaussi ...
pymc3：截斷正常混合
[英]pymc3: Truncated Normal mixture
我想在 pymc3 中混合兩個 TruncatedNormal 分布。 我正在嘗試修改這段文檔。 參見示例 #2 的泊松。 我得到的錯誤是： 將混合模型應用於截斷法線的正確方法是什么？ ...
基於新均值和方差的正確樣本數
[英]Correct number of samples, based on new mean and variance
我有一組樣本，從特定的高斯混合物中采樣。 然后我更新混合參數（均值、方差、權重）。 現在，我想通過新參數更正（移動）這些樣本。 我不想從新混合物中隨機采樣，只需移動舊樣本以匹配新參數。 我怎樣才能做到這一點 ...
通過在 R 中的 mclust package 分析 BIC 的可重復性導致潛在的 class 分析
[英]Reproducibility of BIC results in latent class analysis via mclust package in R
我正在使用 mclust package 在 R 中進行潛在的 class 分析。 我向 select 提供了 BIC 的代碼，最適合 model 或最佳類數，但注意到該解決方案不可重現，因為它給出了不一致的 Z78E6221F6393D1356681IDB398F14CE6DZ 任何時候重新運行相 ...
來自 TensorFlow 概率混合密度網絡的混合參數
[英]Mixture parameters from a TensorFlow Probability mixture density network
如何從使用 TensorFlow Probability 創建的混合密度網絡中獲取混合參數？ 我正在嘗試學習一些關於混合密度網絡的知識，並在此處的 TensorFlow Probability 文檔中遇到了一個示例。 順便說一句，我是這個東西的初學者。 請參閱下面的完整代碼，以上面的示例為起點。 我 ...
用 Julia 繪制分布的混合
[英]Plotting mixtures of distributions with Julia
我想用 Julia 繪制兩個 1d 高斯分布的混合。 我不確定最好的方法是什么。 我正在嘗試使用Distributions.jl並特別針對兩者使用d1 = Normal(0.0, 1.0)和d2 = Normal(1.0, 1.8)定義兩個高斯使用MixtureModel(Normal[ Norma ...
繞過 'BLAS/LAPACK 例程'DGEBAL' 在 CNmixt fn (ContaminatedMixt pkg) 中給出錯誤代碼 -3'
[英]Bypassing 'BLAS/LAPACK routine 'DGEBAL' gave error code -3' in CNmixt fn (ContaminatedMixt pkg)
我正在使用ContaminatedMixt package 、 CNmixt function對多維數據進行聚類。 我的數據維度范圍從 3 到 9。大多數數據都可以毫無問題地進行聚類，但我一直在努力對 7D 數據進行聚類。 抱歉無法包含獨立的代表，但我知道錯誤可靠地發生在我自己的數據中。 您可以通 ...
傳遞給 fitgmdist 的有效起始條件是什么？
[英]What would be valid starting conditions to pass to fitgmdist?
我正在嘗試使用八度音階統計 package 中的 fitgmdist function。 當我讓它使用默認的 k++ 方法來查找初始參數時，它就起作用了。 然而，結果並不一致，有時只是平面錯誤。 這就是為什么我希望能夠通過查看直方圖來傳遞均值的初始值。 但由於某種原因，function 不會接受它 ...
正態混合分布
[英]Normal Mixture Distribution
我正在嘗試創建一個 qqplot 並針對 25% N(μ=0,σ=4) 和 75% N(μ=4,σ=2) 的正態混合分布運行 KS 測試。 我怎樣才能使我的 qqplot 和 KS 測試適應這個分布？ 我不認為我的 abline 是正確的，我的 KS 測試並沒有真正正確地反映分布。 任何幫助，將不勝 ...
如何從石榴高斯混合物中獲取均值和協方差值 model
[英]How to get Mean and Covariance value from pomegranate Gaussian Mixture model
在 scikit 中學習高斯混合 model 我們可以通過以下方式獲得均值和協方差 但在石榴高斯混合的情況下 model 表示沒有名為“covariances_”和“means_”的屬性非常感謝您的寶貴時間。 ...