1,764 questions
1
vote
1
answer
33
views
Bayesian ordinal logistic model using rstan
Is my syntax below for an ordinal logistic model correct? I get an error message which I don't understand.I got the code from a published paper illustrating a graded response model (ordinal logistic, ...
- 247
0
votes
0
answers
38
views
Can't install PyCuba
I am unable to use PyCuba, a Python interface for the multidimensional-integration library Cuba, on macOS 15.5. PyMultiNest, which requires PyCuba, is common for Bayesian analysis in the scientific ...
- 21
1
vote
0
answers
64
views
How can I fix model parametrization for truncated normal distribution in pymc?
I am trying to setup a Bayesian metanalysis of report rates of events. My target estimate is the probability that participants will report X in each study. Every participant provides ~20-40 response, ...
- 11
1
vote
0
answers
34
views
JAGS model insensitive to data
I am new to JAGS and working on an assignment for a class. I have been trying to get my model to work for the past few hours and something just isn't clicking, and I hope someone here may be able to ...
- 11
2
votes
1
answer
114
views
Integrating a density to find the upper limit to give a certain probability
I have a posterior density function and I am trying to integrate its posterior cumulative distribution function (CDF) such that I want to find the upper bound of integration so that the probability is ...
- 345
0
votes
0
answers
16
views
Scaling Data and Prior in Bayesian while using ADVI
Assume I have data X with normal likelihood and Prior Mean C as Laplace prior as inputs to the model.
For easier convergence, I multiply each observation of Data by C and then use Prior Mean 1. And I ...
- 1
1
vote
0
answers
99
views
How to estimate a reliable average using bootstrapping with weighted data and measurement uncertainties?
I have a small dataset where each data point (y_i) comes with its own measurement uncertainty (e_i). I want to estimate a reliable average value for (y) while properly considering the uncertainties (...
- 233
0
votes
1
answer
127
views
Prior distribution for zero-inflated Poisson MCMCglmm?
I'm quite new to Bayesian statistics and not sure how to determine the dimensions of my priors. I am using code from Ben Bolker to run quite a big MCMCglmm with 52 response variables, 7 fixed effects, ...
0
votes
0
answers
70
views
How to resolved the import error with scipy when using keras_tuner?
I have download the tensorflow nightly version 2.19.0.dev20241219 (cpu version) and download the version keras_tuner 1.4.7 and keras 3.7.0. I also check my system has the scipy version 1.15.0. ...
0
votes
0
answers
29
views
In R2JAGS, how can I code an interaction between two variables within a survival model?
My model has three independent variables: snow (a categorical variable of whether or not snow was present, which is arranged as a matrix by time and individual), sex (also categorical), and HFL (a ...
0
votes
0
answers
29
views
Alternative to forecast sequential regressions
Seeking advice on how I could forecast the below better. Ultimately I am looking to build a 5-year rent growth forecast, but as part of the forecasting process, I will need other variables.
Currently, ...
- 313
1
vote
0
answers
12
views
Pogit Model: Correlated or Orthogonal Covariates?
I’m working with the following Pogit model:
Y ~ Pois(lambdap)
log(lambda) = alpa_0 + alpha_1x_1 + alpha_2x_2
logit(p) = gamma_0+gamma_1t_1 + gamma_2*t_2
I’m conducting a simulation study with two ...
- 11
1
vote
1
answer
144
views
Use Bayesian PyMC linear model on out-of-sample data
I am trying to fit a linear model to data using Bayesian inference technique. For this, I thought of using PyMC. Naturally, after training a model, I want to test its performance on new data and that'...
- 1,267
0
votes
1
answer
59
views
Predictions of coupled PyMC linear bayesian models for new inputs
I would like to model two observables, O_2 and O_3 using bayesian linear regression: O_2 ~ m_2 * O_1 + q_2 and O_3 ~ m_3 * O_2 + q_3. I would like to couple them via O_2. I am using the PyMC framework ...
- 113
0
votes
1
answer
207
views
Bayesian optimization in OpenMDAO?
I have a model built in OpenMDAO that I am interested in trying with Bayesian optimization. I have found some work in this area for older versions of OpenMDAO, but nothing for the current OpenMDAO ...
- 2,605
0
votes
1
answer
91
views
emcee refuses to properly explore function
I have a 1D function that looks like this:
i.e. a pronounced drop towards a stable value around 0. The function is written as:
((1. / np.sqrt(1. + x ** 2)) - (1. / np.sqrt(1. + C ** 2))) ** 2
where ...
- 43k
0
votes
1
answer
78
views
BRMS: on what scale are the estimated coefficients of a model with a log link function?
I am using the brms package to fit a model (brms). However, I am unsure on how to interpret the resulting coefficients/effects. Specifically, on what scale are the coefficients if the link function is ...
- 43
3
votes
0
answers
30
views
SMAC with MFFacade returns None
I'm using SMAC (Sequential Model-based Algorithm Configuration) to optimize hyperparameters for training a model (specifically, a PPO agent) using Stable Baselines3. During the training phase, the ...
- 31
0
votes
1
answer
64
views
GPBoost and Bayesian optimization
I am trying to create a parameter grid for a gpboosted model. Currently, I am using the grid tune search method, but I have been reading on Bayesian Optimization in the package rBayesianOptimization. ...
0
votes
1
answer
389
views
PyMC use JAX to compute sample_posterior_predictive with GPU
I'm able to use my GPU to speed up my PyMC model sampling with m.sampling_jax.sample_numpyro_nuts, and it works well in my Linux environment. However, when sampling the posterior_predictive to extend ...
- 127
0
votes
1
answer
30
views
How do I fit a distribution to more than one site at a time in rstan?
There is probably a really easy solution to this, but I’m coming up short. I have code in rstan where I fit a Pareto distribution to data, and I’m looking to find out what the alpha parameter is for ...
- 1
1
vote
0
answers
104
views
How to present findings from the mcmc_areas() and conditional_effects using r brms and bayesplot packages?
I am fairly new to bayesian statistics but have really fallen in love with what it can do and how it presents findings compared to frequentist statistics especially for categorical or ordinal datasets....
1
vote
1
answer
29
views
Simulation: Use of Integrate() in conjunction with custom functions
In attempting to create a simple simulation for psycholinguistics, I have encountered the following error:
> expected_success_rate(170)
> Error in if (height \>= threshold) { : the condition ...
- 11
1
vote
1
answer
713
views
YOLOv8: Optimising for map with confidence and iou in prediction
I'm trying to figure out what the best conf and iou is for the model.pred.
from ultralytics import YOLO
import pandas as pd
import numpy as np
df= pd.DataFrame()
# Load a model
for i in range(1,105):
...
- 868
0
votes
0
answers
114
views
How to do a Bayesian Repeated Measures ANOVA analysis with missing cells in R?
Suppose I am trying to measure people's Happiness depending on Time of Day (Morning/Afternoon/Evening) and Ice Cream Flavor (Chocolate/Strawberry/Vanilla) through a 1-7 scale. Participants are divided ...
