Displaying 2 results from an estimated 2 matches for "model_poisson".
2012 Oct 03
1
Errors when saving output from WinBUGS to R
...1, x2=x2, x3=x3, N=N)
inits<- function(){
list(beta0=0, beta1=0, beta2=0, beta3=0)}
parameters <- c("beta0", "beta1", "beta2", "beta3")
#call winbugs using bugs() function
model <- bugs(data, inits, parameters,
model.file="C:/Jinmodel/test_codes/model_poisson.txt", n.chains=1,
n.iter=30, n.burnin=0,n.thin=1, codaPkg=TRUE, DIC=TRUE, debug=TRUE,
bugs.directory="C:/Users/Yiyi/Desktop/winbugs14/WinBUGS14")
###########################################
I can get pretty good traceplots & parameter estimate outputs in WinBUGS.
However, when...
2010 Jun 02
1
Problems using gamlss to model zero-inflated and overdispersed count data: "the global deviance is increasing"
...al type I and type II NBI NBII and zero-inflated negative binomial ZINBI) while using cubic splines with a df=3. The best fitting model will then be choses on the basis of its AIC.
However, these models frequently fail to converge, and I am not sure why, and what to do about it. For example:
> model_Poisson <- gamlss(duck ~ cs(HHCDI200,df=3) + cs(HHCDI1000,df=3) + cs(HHHDI200,df=3) + cs(HHHDI1000,df=3) + cs(LFAP200,df=3),data=data,family= PO)
GAMLSS-RS iteration 1: Global Deviance = 1350.623
GAMLSS-RS iteration 2: Global Deviance = 1350.623
> model_ZIPoisson <- gamlss(duck ~ cs(HHCDI200,df...