You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Wrote the maintainer of VGAM who has fixed the problem where the CI for N is not printed.
The weight problem is deeper in several packages and likely not fixable.
Try the following code fragment
data_rodli <- data.frame(cap_hist=c("11""10","01"),
t1 = c(1,1,0),
t2 = c(1,0,1),
freq =c(57, 52, 120))
need to expand capture history by frequency
rodli.vgam.expand <- plyr::adply(rodli.vgam, 1, function(x){
x[ rep(1, x$freq),]
})
head(rodli.vgam.expand)
rodli.vgam.m.t <- vglm(cbind(t1,t2) ~ 1,
posbernoulli.t, data = rodli.vgam.expand)
summary(rodli.vgam.m.t)
try with weights
rodli.vgam.m.t2 <- vglm(cbind(t1,t2) ~ 1,
posbernoulli.t, data = rodli.vgam, weights=freq)
summary(rodli.vgam.m.t)
The text was updated successfully, but these errors were encountered: