search for: binomialff

...replace=TRUE) table(x,y) ## y ## 0 1 ##x 0 10 17 ## 1 11 12 Now, it seems to me that if I use vglm to fit cbind(x,y) it simply produces separate fits for x and y, without taking account of their association. Thus, first (some lines of the output removed): summary(vglm(cbind(x,y)~1,family=binomialff(mv=TRUE))) ##Coefficients: ## Value Std. Error t value ##(Intercept):1 -0.16034 0.28375 -0.56508 ##(Intercept):2 0.32277 0.28652 1.12651 ## ##Number of linear predictors: 2 ##Names of linear predictors: logit(E[x]), logit(E[y]) Now if I shuffle say y: y1<-sample(y,50...

...els. 1. When fitting smoothed parameters (e.g. s(X, df=2)) changing the degrees-of-freedom has no effect on the level of smoothing (e.g. number of knots for the spline). This appears to be the case even for the example on the vgam help page. data(hunua) fit2 = vgam(agaaus ~ s(altitude, df=2), binomialff, hunua) > coef(fit2) (Intercept) s(altitude, df = 2) -1.1661259280 0.0003932463 fit10 = vgam(agaaus ~ s(altitude, df=10), binomialff, hunua) > coef(fit10) (Intercept) s(altitude, df = 10) -1.1661259280 0.0003932463 2. There may be a bug w...

...trying to figure out make this seemingly simple modification to the code. If anyone has the patience to step me through the process, or refer me to the relevant information, I would be very grateful. 2) Using the loglog function in VGAM: I tried fitting the model: m<- vglm(Y~X, family=binomialff(link="loglog"), data = d) but I get the following error? Error in lm.fit(X_vlm, z_vlm, ...) : NA/NaN/Inf in foreign function call (arg 4) In addition: Warning message:In log(log(theta)) : NaNs produced?. If I run the same model using cloglog as the link, I get a result that looks...

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil --

Help needed with lm function: Dear R's, Could anyone tell me how to replace the link function (probit logit, loglog etc.) in lm with an abitrary user-defined function? The task is to perform ML Estimation of betas for a dichotome target variable. Maybe there is already a package for this (I did not find one). Any hints or a code excerpt would be welcome! Thank you -Jeff jeff.pr2 (at)

...nction betageometric Beta-geometric Distribution Family Function betaprime The Beta-Prime Distribution bilogistic4 Bivariate Logistic Distribution Family Function binom2.or Bivariate Logistic Regression binom2.rho Bivariate Probit Model binomialff Binomial Family Function bisa Birnbaum-Saunders Distribution Family Function brat Bradley Terry Model bratt Bradley Terry Model With Ties calibrate.qrrvglm Calibration for CQO, UQO and CAO models cao Fitt...