Displaying 1 result from an estimated 1 matches for "doserep".
2010 Dec 30
1
Different results in glm() probit model using vector vs. two-column matrix response
...)
dose <- rep(c(7, 8, 9, 10, 11), 2)
ldose <- log10(dose)
n <- 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...