от
Я бы хотел модель соотношения риска для данных в опросе, но у меня возникли трудности разработки проблемы pymc3. Коэффициент риска-отношение вероятности наступления определенного исхода дается экспозиция, чтобы вероятность без воздействия. Так как я понимаю, соотношение риск: (заболеваемости в облученной группе) / (заболеваемость в группе notExposed) Я пробовал постановка задачи как показано в следующем примере кода:
import pymc3

exposed_num_observations = 25
exposed_observed = 11
notExposed_num_observations = 853
notExposed_observed = 202

with pm.Model() as rr_mdoel:
    # Prior
    p_exposed = pm.Uniform("p_exposed", 0 ,1)
    p_notExposed = pm.Uniform("p_notExposed", 0 ,1)
    sd = pm.Uniform("sd", 0,5)

    # Likelihood
    exposed = pm.Binomial("exposed", 
                      p=p_exposed, 
                      n=exposed_num_observations,  
                      observed=exposed_observed)
    notExposed = pm.Binomial("notExposed", 
                             p=p_notExposed, 
                             n=notExposed_num_observations,  
                             observed=notExposed_observed)

    risk_ratio = pm.Normal("risk_ratio", mu=(p_exposed/p_notExposed), sd=sd)

    # Inference Run/ Markov chain Monte Carlo 
    trace = pm.sample(10000, chains=2, tune=1000)
Я ожидал, что
risk_ratio
всегда быть выше 0, и я ожидал, что
sd
будет соответствовать некоторое значение. Результаты (как показано на
az.plot_posterior(trace, round_to=3)
например) показывает негативные значения для
risk_ratio
и
sd
по-прежнему однородным.

Ваш ответ

Отображаемое имя (по желанию):
Конфиденциальность: Ваш электронный адрес будет использоваться только для отправки уведомлений.
Анти-спам проверка:
Чтобы избежать проверки в будущем, пожалуйста войдите или зарегистрируйтесь.
Добро пожаловать на сайт ByNets, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...