search for: y_u

Displaying 1 result from an estimated 1 matches for "y_u".

Did you mean: __u
2010 Dec 30
1
Different results in glm() probit model using vector vs. two-column matrix response
...<- c(8,8,8,8,8,8,8,8,8,8) r <- c(0,1,3,8,8,0,0,0,4,5) p <- r/n d <- data.frame(group, dose, ldose, n, r, p) SF <- cbind(success=d$r, failure=d$n - d$r) #80 record set-up dose2<-c(7,8,9,10,11) doserep<-sort(rep(dose2,8)) x<-c(doserep,doserep) log10x<-log10(x) y_U<-c(rep(0,8), 1, rep(0, 7), 1, 1, 1, rep(0,5), rep(1, 16)) y_M<-c(rep(0,24), rep(1,4), rep(0,4), rep(1,5), rep(0,3)) y<-c(y_U, y_M) trt<-c(rep(1, 40), rep(0, 40)) # print x & y's for both SF y ldose log10x # analysis with 10 records and 80 records f1 <- glm(SF ~ group + ldos...