2013-04-17 3 views
1

Я использую следующий код в OpenBUGS выполнить анализ:OpenBUGS - переменный не определен

model 
{ 
for(i in 1:467) 
{ 
probit(p[i])<-gamma0+gamma1*drug[i]+gamma2*CD41[i] 
R[i]~dbern(p[i]) 
junk[i]<-ID[i] 
} 
gamma0~dnorm(0,.0001) 
gamma1~dnorm(0,.0001) 
gamma2~dnorm(0,.0001) 
} 
ID[] drug[] CD41[] R[] 
1 0 114 NA 
2 1 40 NA 
3 1 12 0 
4 0 15 0 
.... 
END 

И я получаю следующее сообщение об ошибке: Переменная CD41 [] не определена. Я не уверен, как решить эту проблему, поэтому любая помощь будет принята с благодарностью.

ответ

2

Вы уверены, что правильно передали векторную переменную CD41 в свой список данных?

Если вы используете R, попробуйте позвонить CD41<-as.vector(CD41), прежде чем передавать его в OpenBUGS. Помните, что 1 матрица столбца и вектор - разные вещи как для BUG, ​​так и для R.

+0

Я не уверен, что вы имеете в виду. Формат, который я использовал, обычно заключается в том, как я передавал большие переменные векторы в OpenBUGS в прошлом. Я использую программное обеспечение OpenBUGS, а не пакет R. Изменить: если я не могу понять это, я попытаюсь реализовать его в R-пакете позже в качестве альтернативы и посмотреть, исправляет ли он это. У меня сейчас нет времени на это. –

+0

Я вижу, но сначала тщательно проверьте, не пропустили ли вы что-либо набрали, передав данные в BUGS или что-то в этом роде, может быть ошибка. –

+0

Таким образом, я нашел ошибку, у меня было NA в моих данных для зависимой переменной. Спасибо за вашу помощь. –

Смежные вопросы