Re: [R] VGAM package : Frechet distribution - 2 parameter estimation

Hello, a day or two ago I submitted VGAM 0.9-5 to CRAN, which has myriads of changes to family functions---their names, their arguments, and their order thereof. Especially regarding family functions for discrete and continuous distributions. In a nutshell, I found lots of inconsistencies while

[R] Testing proportional odds assumption in R

Hi, M1 and M2 are extreme in that all or none of the variables have parallel lines on the logit scale. One can try fitting a partial POM, which remains fraught (but not as much as M2) because if the lines intersect for a particular variable where the data lie then there will be numerical problem

Re: [R] VGAM constraints-related puzzle

Edward Wallace gmail.com> writes: > > Hello R users, > I have a puzzle with the VGAM package, on my first excursion into > generalized additive models, in that this very nice package seems to > want to do either more or less than what I want. > > Precisely, I have a 4-component outcome, y, and

Re: [R] constraint matrices in vglm (VGAM package)

Hi Lin, try multi2 = vglm(case123con ~ SNP_A1+SNP_A2+age, multinomial(parallel = TRUE ~ SNP_A1+SNP_A2 - 1), work.analy) or multi3 = vglm(case123con ~ SNP_A1+SNP_A2+age, multinomial(parallel = FALSE ~ 1 + age), work.analy) After your fit,

Re: [R] VGAM and vglm

Hi Ted, yes, vglm(..., family=binomialff(mv=TRUE)) does treat each response as independent. You can see that with coef(fit, matrix=TRUE) where fit is the object. If you want to model two dependent binary responses then you can try the VGAM family functions binom2.or(), loglinb2(), binom2.rho().