Displaying 1 result from an estimated 1 matches for "thbad".
Did you mean:
thba
2009 Jun 22
1
How to make try to catch warnings in logistic glm
...s<-matrix(NA, nrow=nsim, ncol=6)
ow <- options("warn")
options(warn= 1)
outcome<-rep(c("G","B"), times=c(ngood, nbad))
outcome<-as.factor(outcome)
agood<-rep(c(0,1), times=c(27,1))
abad<-rep(c(0,1), times=c(14,7))
thgood<-rep(c(0,1), times=c(24,4))
thbad<-rep(c(0,1), times=c(12,9))
for (i in 1:nsim)
{
?tgood<-sample(agood, ngood)
?tbad<-sample(abad, nbad)
?hstop<-c(tgood, tbad)
?hstop<-as.factor(hstop)
?aldgood<-rnorm(ngood, mean=54, sd=8/0.675)
?aldbad<-rnorm(nbad, mean=64, sd=8/0.675)
?ald<-c(aldgood, aldbad)
?fatgood<...