Hi,


I am running a negative binomial model using Gamlss and when I try to include 
random effect, I get the following message:



Warning messages:
1: In vcov.gamlss(object, "all") :
  addive terms exists in the mu formula standard errors for the linear terms 
maybe are not appropriate
2: In vcov.gamlss(object, "all") :
  addive terms exists in the  sigma  formula standard errors for the linear 
terms maybe are not appropriate
3: In summary.gamlss(model1) :
  summary: vcov has failed, option qr is used instead

I still get the output for the fixed effects but not for the random effects. 
Can you recommend on a solution for this problem?



Below is the code that I wrote:



install.packages("gamlss")
library(utils)
library(MASS)
library (gamlss)

Data_Reduced = read.table(file.choose(), header=TRUE, sep="\t")
model1<-gamlss(Events~ factor(Solo_Time) + factor(Fleet)+ Father_SS + Mother_SS 
+ random(factor(Driver_ID)) + offset(logDuration), 
sigma.fo=~random(factor(Driver_ID)), data = na.omit(Data_Reduced), family=NBI, 
na.action = "na.omit")
summary(model1)


Thanks in advance.

Best regards,
Haneen


--------------------------------------------------------------------------------------
Haneen Farah, Ph.D.
Department of Transport Sciences
School of Architecture & Build Environment
KTH- Royal Institute of Technology
Teknikringen 72, SE-100 44 Stockholm
Tel: +46-8-7907977
fha...@kth.se<mailto:fha...@kth.se>
http://home.abe.kth.se/~fhanin/


        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to