Displaying 1 result from an estimated 1 matches for "log_v".
Did you mean:
log_p
2007 Nov 12
4
a repetition of simulation
...1, .6)
x[i] <- rbinom(1, 1, .95)
y[i] <- z[i]*x[i]
if (y[i]==1) aps[i] <- rnorm(1,mean=13.4, sd=7.09) else aps[i] <-
rnorm(1,mean=12.67, sd=6.82)
if (y[i]==1) tiss[i] <- rnorm(1,mean=20.731,sd=9.751) else tiss[i] <-
rnorm(1,mean=18.531,sd=9.499)
}
v <- data.frame(y, aps, tiss)
log_v <- glm(y~., family=binomial, data=v)
summary(log_v)
I want to do a repetition of this 600 times (I want to have 600 logistic
models), and see all the coefficients of the covariates aps & tiss.
Thanks in advance,
Sigalit.
[[alternative HTML version deleted]]