Hi all! Maybe someone could help me with the following. I know this hasn't directly to do with ecology but I'm also using glm. I have a list of 16 genes and 10 samples. The samples are of two types, 4 Ctrl and 6 Diseased. If,
labelInd<-as.factor(c(rep("0",4),rep("1",6))) genes.glm<-glm(labelInd ~ ., family=binomial, data=mat) beeing "mat" the 10x16 matrix (without NAs), I got 17 values, first the intercept, 9 numerical values and "NA" for the last 7 genes. Does anybody you know why this is happening or how I can model using the 16 genes? I hope anyone could help me with this! Many thanks in advance, Paul [[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.