search for: data_reduc

...mmary: 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.om...

I have data that looks like this: > df1 group id 1 red A 2 red B 3 red C 4 blue D 5 blue E 6 blue F I want a list of the groups containing vectors with the ids. I am avoiding subset(), as it is only recommended for interactive use. Here's what I have so far: df1 <- data.frame(group=c("red", "red", "red", "blue",