search for: betabin

Hi, can anyone help me fit betabinomial model to the following dataset where each iD is a cluster in itself , if i use package aod 's betabinom model it gives an estimate of zero to phi(the correlation coeficient ) and if i fix it to the anova type estimate obtained from icc( in package aod) then it says system is ex...

...R 2.5.0alpha and noticed new warnings related to the extractor "$" when using contributed packages. For instance: > library(RODBC) Warning message: $ operator is not valid for atomic vectors, returning NULL > library(aod) Package aod, version 1.1-18 > data(orob2) > m1 <- betabin(cbind(y, n-y) ~ 1, random = ~ 1, data = orob2) > head(residuals(m1)) 1 2 3 4 5 6 -1.4629774 -0.7897877 -1.3754283 0.1029910 -0.3556985 1.4030026 Warning message: $ operator not defined for this S4 class, returning NULL In these two cases,...

...cattered in different packages. The proposed functions should be considered as complements to more sophisticated methods such as generalized estimating equations (GEE) or generalized linear mixed effect models (GLMM). Index (edited): betabin Beta-Binomial Model for Proportions donner Test of Proportion Homogeneity using Donner's Adjustment icc Intra-Cluster Correlation negbin Negative-Binomial Model for Counts quasibin Quas...

...cattered in different packages. The proposed functions should be considered as complements to more sophisticated methods such as generalized estimating equations (GEE) or generalized linear mixed effect models (GLMM). Index (edited): betabin Beta-Binomial Model for Proportions donner Test of Proportion Homogeneity using Donner's Adjustment icc Intra-Cluster Correlation negbin Negative-Binomial Model for Counts quasibin Quas...

I have a problem, would really appreciate if anyone could help me. I am trying to use bootstrap to resample a dataframe n then use that sample of dataframe to calcuate parameter estimates using betabinomial model. I am kind of lost coz the bootstrap function requires us to give statistic (A function which when applied to data returns a vector containing the statistic(s) of interest) . my dataframe looks like : y<-c(1,0,2,2,1,0,1) n<-rep(2,length(y)) z1<-rep(1,...

hi, a) i have something like: ecdfgrp1<-ecdf(subset(mydata,TMT_GRP==1)$Y); ecdfgrp2<-ecdf(subset(mydata,TMT_GRP==2)$Y); how can i plot the difference between these 2 step functions? i could begin with ecdfrefl<-function(x){ecdfgrp2(x)-ecdfgrp1(x);} ... what next? b) if i have a vector with repeated numeric values how can i get the subset without repeated values .e.g (0,4,0,2,2)

I am getting "Chi-squared approximation may be incorrect in: chisq.test(x)" with the data bleow. Frequency distribution of number of male offspring in families of size 5. Number of Male Offspring N 0 518 1 2245 2 4621 3 4753 4 2476 5

...b2 Blood Group System ABO The ABO Blood Group System acat Ordinal Regression with Adjacent Categories Probabilities benini Benini Distribution Family Function betaII Beta Distribution of the Second Kind betabin.ab Beta-binomial Distribution Family Function betabinomial Beta-binomial Distribution Family Function betaff The 2-parameter Beta Distribution Family Function betageometric Beta-geometric Distribution Family Function betaprime The Bet...

Hi all- I fit a zero-inflated Poisson model to model bycatch rates using an offset term for effort. I need to apply the fitted model to a datasets of varying levels of effort to predict the associated levels of bycatch. I am seeking assistance as to the correct way to code this. Thanks in advance! Laura [[alternative HTML version deleted]]

I am trying to model data in which subjects are followed through time to determine if they fall, or do not fall. Some of the subjects fall once, some fall several times. Follow-up time varies from subject to subject. I know how to model time to the first fall (e.g. Cox Proportional Hazards, Kaplan-Meir analyses, etc.) but I am not sure how I can model the data if I include the data for those