pymc3

3熱度

1回答

我想用我的PyMC3 LR模型來獲得預測變量y的值的80％HPD範圍，因爲新數據可用。因此，外推y的值的可信分佈值爲x的新值不在我的原始數據集中。型號： with pm.Model() as model_tlr: alpha = pm.Normal('alpha', mu=0, sd=10) beta = pm.Normal('beta', mu=0, sd=10)

0熱度

1回答

在Windows 10上使用PYMC3 - theano無法導入名稱'floatX'

我努力讓PYMC3在Windows上正確安裝。我已經嘗試過使用Anaconda軟件包，通過conda install -c conda-forge pymc3以及virtualenv只使用pip。它似乎安裝正常，但運行import pymc3時出現以下錯誤，但失敗。研究表明，可能存在一些錯過的依賴關係。 >>> import pymc3 WARNING (theano.sandbox.cuda

0熱度

1回答

Pymc3多維變量方程

我想在pymc3中使用多維變量來描述共享一個公因子的4個變量的系統。也許我正在考慮這個錯誤，但是當我設置模型時，我不斷從pymc3中得到有關輸入尺寸不匹配的錯誤。我嘗試了幾個變化無濟於事。我如何使用pymc3的形狀參數來描述一個具有共享觀察變量的結構的4個觀察變量的系統？ import pymc3 import numpy # Simulate some data shared_obs =

1熱度

1回答

PyMC3模型不產生接近

我已經使用PyMC3以下在下面的示例中所提供的步驟中創建的MCMC模型所觀察到的數據的示例數據： https://pymc-devs.github.io/pymc3/notebooks/LKJ.html 結果，我得到類似於以下圖片：模擬樣本跟隨觀察數據具有相同的模式。但是，X和Y的範圍要小得多。有誰知道我可以如何獲得涵蓋與觀察數據類似的數字範圍的樣本？

3熱度

1回答

隨機效應的多級迴歸中的變量解釋

我有一個看起來像下面的數據集（顯示的前5行）。 CPA是來自不同廣告航班實驗（治療）的觀察結果。航班分層次地分組在活動中。 campaign_uid flight_uid treatment CPA 0 0C2o4hHDSN 0FBU5oULvg control -50.757370 1 0C2o4hHDSN 0FhOqhtsl9 control 10.963426 2 0C2o4hHD

1熱度

1回答

在gpu上的pymc3上的ADVI

以下示例： http://twiecki.github.io/blog/2016/07/05/bayesian-deep-learning/，可以使用每個變量的float32值在gpu上訓練網絡。但是，奇怪的是，如果有更多的圖層被傳遞到神經網絡（例如每層有800個神經元的5層），或者當使用較少的圖層將神經元的數量設置爲較高的值時，我收到如下錯誤消息： 0%| | 1/50000 [00:00<6:

2熱度

1回答

PYMC3：NUTS難以從分層零充氣伽馬模型採樣

我試圖從a paper from Richard McElreath複製數據分析，其中他將數據與分層零膨脹Gamma模型擬合。這些數據是關於狩獵返回約15000狩獵之旅從約150二十多年的獵人。由於許多狩獵旅行的回報率爲零，因此該模型假定每次旅行具有零回報率的概率爲pi，而1 - pi爲正數回報率的概率爲alpha和beta之間的Gamma分佈。預測變量是年齡，該模型使用年齡多項式（最高爲3）來

2熱度

1回答

簡單貝葉斯網絡通過蒙特卡羅馬爾可夫鏈移植到PyMC3

我正在將Simple Bayesian Network via Monte Carlo Markov Chain的示例從PyMC2移植到PyMC3，它的工作原理。結果可以在以下的gist文件中找到關於GitHub的pymc3_rain_sprinkler_grass_simple_bayesian_network.py. 我想通過提供證據，例如延長原始示例草地是溼的，然後讓PyMC3給我答案，如

1熱度

2回答

pymc3分層模型與多個觀察，不計算MCMC期間的可能性？

我試圖用pymc3和運行陷入困境的地方看起來像我的先驗正在取樣，但我可能是不適合的分層模型。我曾嘗試複製描述爲here的方法，該方法適用於多重觀察。我正在使用未標記的包中的R中找到的數據集鳧，加載爲熊貓數據幀。您可以從我的GitHub頁面here獲取我正在使用的CSV。簡而言之，每行表示最多三個沿着單個樣線收集的觀測值（y1，y2和y3列）。我嘗試刪除任何y列或協變量（海拔，森林覆蓋率，調查長度

2熱度

2回答

定製Theano運算數值積分

我試圖寫一個自定義Theano歐普其數字集成了兩個值之間的函數。 Op是一種PyMC3的定製可能性，涉及一些積分的數值計算。我不能簡單地使用@as_op裝飾器，因爲我需要使用HMC來執行MCMC步驟。任何幫助將不勝感激，因爲這個問題似乎已經走到了好幾次，但一直沒有得到解決（例如https://stackoverflow.com/questions/36853015/using-theano-wit