поиск источников в R-сценарии для MCM C (Байесовская статистика) Анализ в R - PullRequest
Новая
       95

поиск источников в R-сценарии для MCM C (Байесовская статистика) Анализ в R

0 голосов
/ 06 апреля 2020

Я пытаюсь запустить Байесовский анализ линейной регрессии. Я пытаюсь найти файл Jags.R, но продолжаю сталкиваться с ошибками, которые не могу устранить.

Вот основной код 

myData = read.csv(file="HeightWeightData.csv" )
xName = "height" ; yName = "weight"
source("Jags.R")

# Run the MCMC process:
mcmcCoda = genMCMC( data=myData , xName=xName , yName=yName )

# Display posterior distribution:
plotMCMC( mcmcCoda , data=myData , xName=xName , yName=yName ,
compValBeta1=0.0 , ropeBeta1=c(-0.5,0.5) ,
pairsPlot=TRUE )

Исходный файл Jags.R выглядит следующим образом: 

 model{
 for (i in 1:Ntotal) 
  {
  zy[i] ~ dt( zbeta0 + zbeta1 * zx[i] , 1/zsigma^2 , nu )
  } 
  zbeta0 ~ dnorm( 0 , 1/(20)^2 )
  zbeta1 ~ dnorm( 0 , 1/(20)^2 )
  zsigma ~ dunif( 1.0E-4 , 1.0E+4 )
  nu ~ dexp(1/40)
  }

Я продолжаю сталкиваться с тремя ошибками в исходных зазубринах. Файл R. 

Error: unexpected '{' in "model{"
Error in Ntotal : object 'Ntotal' not found
Error: unexpected '}' in "}"
...