Displaying 2 results from an estimated 2 matches for "xtild".
Did you mean:
xild
2013 Jan 18
0
problem that arises after using the new version of "BRugs"
...intKnots <- quantile(unique(x), seq(0, 1, length = (numIntKnots +
2))[-c(1, (numIntKnots + 2))])
}
numIntKnots <- length(intKnots)
allKnots <- c(rep(range.x[1], 4), intKnots, rep(range.x[2],
4))
K <- length(intKnots)
L <- 3 * (K + 8)
xtilde <- (rep(allKnots, each = 3)[-c(1, (L - 1), L)] + rep(allKnots,
each = 3)[-c(1, 2, L)])/2
wts <- rep(diff(allKnots), each = 3) * rep(c(1, 4, 1)/6,
K + 7)
Bdd <- spline.des(allKnots, xtilde, derivs = rep(2, length(xtilde)),
outer.ok = TRUE)$design
Omeg...
2012 Jul 09
1
R to winbugs interface
...mission denied
In addition in WinBUGS a window crops up which is named black box. and says
trap #060 and then gives some details. I am not being able to understand
what all hese error messages mean. I am also giving the R code and the code
for the model here :
R code :
regpractwo=function(Ytilde,D,Xtilde,a,b,n,FEV,Age,Smoke)
{
Xtildeinv=solve(Xtilde)
beta0=c(t(Xtildeinv%*%Ytilde))
C0=Xtildeinv%*%D%*%t(Xtildeinv)
C0inv=solve(C0)
r=dim(Xtildeinv)[1]
FEVdataBUGS=list("n","r","beta0","C0inv","a","b","FEV","Age","Smoke...