mcmc

1열

1답변

Metropolis-Hastings 방법을 사용하여 Matlab의 적분 값을 계산하십시오.

Metropolis-Hasting의 방법을 사용하여 Matlab의 적분을 평가하는 데 약간의 문제가 있습니다. 적분은 0에서 무한대까지의 e^(x^-2)입니다. 내가 작성한 코드는 오류가 발생하지 않지만, 1) 내가하고 싶은 일을하고 있는지 잘 모르겠다. 2) 내가 원하는대로해도 2) '추출' '코드가 나는이 문제 같은 느낌 clc clear all

0열

1답변

pymc3의 다 변수 선형 회귀

나는 최근에 emcee 만 독점적으로 사용한 후 pymc3을 학습하기 시작했으며 일부 개념적 문제가 발생합니다. 제 7 장 Hogg's Fitting a model to data으로 연습하고 있습니다. 이것은 임의의 2 차원 불확실성을 갖는 직선에 대한 mcmc 피트를 포함한다. 나는 이것을 emcee에서 아주 쉽게 수행했지만, pymc이 나에게 몇 가지

0열

1답변

MCMCglmm 이항 모델은

이전에 나는 R 패키지 MCMCglmm와 이항 모델을 추정 할 수 있습니다. 모델은 절편과 경사면을 고정 된 부분과 무작위 부분으로 통합해야한다. 허용되는 prior을 어떻게 지정해야합니까? (하지만, 훨씬 더 복잡한 환경에서, here is a similar question 참고.) 데이터는 다음과 같은 형태가 있다고 가정 : y x cluster

2열

1답변

MCMCglmm의 "형질"의 의미

this post과 마찬가지로 MCMCglmm이라는 표기법, 특히 trait의 의미로 고민하고 있습니다. 나에게 오류 메시지 For error structures involving catgeorical data with more than 2 categories pleasue use trait:units or variance.function(trait):un

1열

1답변

MCMCglmm 베이 즈 팩터

MCMCglmm 패키지를 사용하여 데이터를 피팅하는 베이지안 모델을 사용하려고합니다. 그러나 두 모델을 비교하기 위해 Bayes Factor를 계산하는 방법에 대한 참조는 찾을 수 없습니다. 의견이나 제안이 있으십니까?

1열

1답변

MCMCglmm에서 중도 절단 데이터 모델링을위한 적절한 데이터 구조는 무엇입니까

내 행동 (Cort 및 Pred)이 개별 동작에 미치는 영향을 추정하기 위해 혼합 모델을 만들고 싶습니다. 내 데이터가 우연히 검열됨에 따라 (주어진 시간 내에 피난처에서 출현하지 않는 개인은 최대 점수 인 1200 점을 얻고 있습니다.), 나는 'cengaussian'가족과 함께 MCMCglmm 패키지를 사용하려고합니다. 내 원래의 데이터 구조는 다음과

2열

1답변

내가 여기 R에서 작업하고

메트로 폴리스 - 헤이스팅스 알고리즘 (MCMC)를 사용하여 프로 빗 모델을 시뮬레이션하는 것은 내 코드입니다 : 나는 문제가 추측 n <- 1000 trueB0 <- 0 trueB1 <- 0.1 prior.mean0 <- 0 prior.sd0 <- 10 prior.mean1 <- 0 prior.sd1 <- 10 x <- rnorm(n, 0,

1열

1답변

병렬로 여러 베이지안 체인을 실행 중 R : 3 노드에서 오류가 발생했습니다. 첫 번째 오류 : 잘못된 치수 수

여러 코어에서 베이지안 통계 함수의 복제 체인을 실행하려고합니다. 코어 당 1 개의 체인. 아래의 mcmc 함수는 단일 mcmc 체인을 실행하는 스크립트입니다. 내 생각으로는 mcmc 함수를 세 번 실행할 수 있습니다. 각 인스턴스는 별도의 코어에 있습니다. 수정하려고 시도한 몇 가지 예를 발견했지만 적절하게 실행하지 못했습니다. 다음과 같은 오류가 발생합

0열

1답변

JAGS에서 사후 분포 해석

이산 분포의 사후 확률을 요약하면 둘 이상의 값이 제공됩니다. 내가 어디로 잘못 가고 있니? 이가 들쭉날쭉 의해 발생되는 후방이다

0열

1답변

베타 회귀 분석을위한 'zoib'R 패키지를 사용할 때 'ts'함수의 오류

R의 베타 회귀 분석을 수행하기 위해 'zoib'패키지를 사용하고 있습니다. 41 페이지의 예제를 복제하려고합니다. 는 R 저널에 발표 된 패키지 저자 : 내가하려고하면 오류가 점점 오전, library(zoib) data("GasolineYield", package="zoib") GasolineYield$batch <- as.factor(Gasolin