Hi Everyone, I am having trouble running glmmLasso.
An abbreviated version of my dataset is here: https://drive.google.com/open?id=0B_LliPDGUoZbVVFQS2VOV3hGN3c Activity is a measure of brain activity, Novelty and Valence are categorical variables coding the type of stimulus used to elicit the response, ROI is a categorical variable coding three regions of the brain that we have sampled this activity from, and STAIt is a continuous measure representing degree of a specific personality trait of the subjects. Subject is an ID number for the individuals the data was sampled from. Before glmmLasso I am running: KNov$Subject <- factor(KNov$Subject) to ensure the subject ID is not treated as a continuous variable. If I run: glm1 <- glmmLasso(Activity~as.factor(Novelty) + as.factor(Valence) + STAIt + as.factor(ROI) + as.factor(Valence):as.factor(ROI), list(Subject=~1), data = KNov, lambda=10) summary(glm1) I don't get any warning messages, but the output contains b estimates only, no SE or p-values. If I try to include a 3-way interaction, such as: glm2 <- glmmLasso(Activity~as.factor(Novelty) + as.factor(Valence) + STAIt + as.factor(ROI) + as.factor(Novelty):as.factor(Valence):as.factor(ROI), list(Subject=~1), data = Nov7T, lambda=10) summary(glm2) I get the warnings: Warning messages: 1: In split.default((1:ncol(X))[-inotpen.which], ipen) : data length is not a multiple of split variable 2: In lambda_vec * sqrt(block2) : longer object length is not a multiple of shorter object length And again, I do get parameter estimates, and no SE or p-values. If I include my continuous variable in any interaction, such as: glm3 <- glmmLasso(Activity~as.factor(Novelty) + as.factor(Valence) + STAIt + as.factor(ROI) + as.factor(Valence):as.factor(ROI) + as.factor(Novelty):STAIt, list(Subject=~1), data = Nov7T, lambda=10) summary(glm3) I get the error message: Error in rep(control$index[i], length.fac) : invalid 'times' argument and no output. If anyone has an input as to (1) why I am not getting SE or p-values in my outputs (2) the meaning of there warnings I get when I include a 3-way variable, and if they are something to worry about, how to fix them and (3) how to fix the error message I get when I include my continuous factor in an interatction, I would be very appreciative. Thanks! Walker ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.