[R] R에서 WINBUGS를 불러서 짜는데 프로그램이 ERROR가 뜹니다. 고수님들 답변좀 부탁드립니다.

Bytarra1

먼저 데이터는 아래의 내용인데, 이파일을 “d:/work/count.txt”로 저장했습니다.

===========================
id at.bat hit salary
A 6 5 8000
B 20 12 7500
C 5 3 7300
D 28 14 7000
E 8 4 6500
F 45 21 6000
G 39 16 5500
H 5 2 5500
I 81 26 4500
J 83 21 4000
K 8 2 3200
L 24 6 3000
M 5 1 3000
N 25 5 2500
O 41 8 2300
P 13 0 2000
Q 5 0 2000

===========================

그리고 아래는 winbugs에서 사용했던 model 입니다.

이부분의 속성은 “D:/work/baseball.txt” 파일로 저장했습니다.

====================================================

model
{
for(i in 1:K){
  y[i]~dbin(p[i],n[i])
  u[i]~dnorm(0, tau)
    logit(p[i])<- beta0 + x[i]*beta1 + u[i]
}
beta0~dnorm(0.0, 1.0E-6 )
beta1~dnorm(0.0, 1.0E-6 )
tau~dgamma(0.001, 0.001)
sigma <- 1/sqrt(tau)
}

이제 부터 R 파일 입니다.

==========================================

require(R2WinBUGS)   # R2WinBUGS 패키지 인스톨 한후
baseball <- read.table("D:/work/count.txt",header=T)  #데이터를 불러들였습니다.
baseball
J <- nrow(baseball)
at.bat <- baseball$at.bat
hit <- baseball$hit
salary <- baseball$salary
data <- list("J", "at.bat", "hit", "salary")  
# 데이터를 다음과 같은 목록으로 재배치 하였습니다

inits =function() {list(beta0=0,beta1=0,tau=1)}
#초기값을 주고~~~이부분이 정확히 맞는지 의심스러움??

parameters <- c("beta0","beta1","tau")      # parameters 지정해주고
baseball.sim <- bugs(data,inits,model.file="D:/work/baseball.txt",
                parameters=c(“beta0″,”beta1″,”tau”),n.chains=3,n.iter=1000 )  
# bugs를 이용해 model을  불러들었는데.error가 뜹니다.

print(baseball.sim)    # bugs가 못아가지 않아서 이부분도 실행안됨!!

==============================================================================

R 고수님들 무한 리플과 답변 부탁드립니다..

Legacy document_srl: 81074 / Legacy URL: http://www.statedu.com/QnA/81074

Similar Posts