https://nd.psychstat.org/ 2026-04-08T03:14:41+00:00 https://nd.psychstat.org/ https://nd.psychstat.org/_media/wiki/dokuwiki.svg text/html 2016-01-24T14:48:26+00:00 Anonymous (anonymous@undisclosed.example.com) https://nd.psychstat.org/lab/zzhang/consultant/winter test<-read.table('/Users/zzhang4/Desktop/test.txt') test<-t(test) i1<-seq(1,34,2) i2<-seq(2,34,2) dset<-cbind(test[i1,1], test[i2,1],test[i1,2], test[i2,2],test[i1,3], test[i2,3]) dset1<-cbind(rbind(dset[,1:2], dset[,3:4],dset[,5:6]), rep(1:3, each=17), rep(c(rep(0,9),rep(1,8)),3)) dset1<-cbind(dset1[,2], dset1[,1]-dset1[,2], dset1[,3:4]) m1<-glm(dset1[, 1:2]~factor(dset1[,3])+dset1[,4], family=binomial("logit"))