Displaying 1 result from an estimated 1 matches for "statusagga".
2006 Feb 01
1
glm-logistic on discrete-time methods with individual and aggregated data
...4), rep(1:2,2), rep(1:3 , 2), rep(1:4,
2))
statusINDa <- c(rep(0 ,16), rep(0:1,2), rep(c(0,0,1), 2),
rep(c(0,0,0,1),2))
datiINDa <- as.data.frame(cbind(timesINDa, statusINDa))
datiINDa$timesINDa <- as.factor(datiINDa$timesINDa)
timesAGGa <- c( 1:4, 1:2, 1:3, 1:4)
statusAGGa <- c(rep(0,4), 0:1, c(0,0,1), c(0,0,0,1))
weightAGGa <- c(rep(4,4), rep(2,2), rep(2,3), rep(2,4))
datiAGGa <- as.data.frame(cbind(timesAGGa, statusAGGa, weightAGGa))
datiAGGa$timesAGGa <- as.factor(datiAGGa$timesAGGa)
coef(glm(statusINDa ~ timesINDa, family=binomial, data=dati...