由于我是winbugs的新手,我需要在R中进行贝叶斯分析.为此,我需要将下面的winbugs代码转换为R代码:
model{
# model’s likelihood
for (i in 1:n){
time[i] ~ dnorm( mu[i], tau ) # stochastic componenent
# link and linear predictor
mu[i] <- beta0 + beta1 * cases[i] + beta2 * distance[i]
}
# prior distributions
tau ~ dgamma( 0.01, 0.01 )
beta0 ~ dnorm( 0.0, 1.0E-4)
beta1 ~ dnorm( 0.0, 1.0E-4)
beta2 ~ dnorm( 0.0, 1.0E-4)
# definition of sigma
s2<-1/tau
s <-sqrt(s2)
# calculation of the sample variance
for (i in 1:n){ c.time[i]<-time[i]-mean(time[]) }
sy2 <- …比方说我有一个矩阵,X我需要在R中使用它.我使用colnames命令命名每一列.但是,当我键入列的名称时,没有任何内容出现.为了说明,我使用了如下代码:
colnames(X) <- c("column1","column2")
当我键入X,column1并column2出现在列的顶部.但是,当我输入column1或者column2,它们无法找到.
有谁知道需要做什么？