pc = read.table('../../Data/nodal.dat', header=FALSE) colnames(pc) = c('xray', 'stage', 'grade', 'age', 'acid', 'node') ##print(pc) resp = as.vector(pc$node) index = seq(1, dim(pc)[1]) lr = glm(resp ~ index, family=binomial, data=pc) print(summary(lr))