[R] R vglm new family writing: mix Poisson/multinomial

From: Benoit J <benoit_jottreau_at_yahoo.fr>
Date: Mon, 7 Jan 2008 01:28:08 -0800 (PST)

Hi dear R users,

1)
I would like to know if there is a simple way to define a vglm family which would be a mix of poisson variables and bernoulli variables (0/1 response) for idea this would be invoked like this:

 vglm(...,family=mixpoissonmultinom(npoisson,n01response))

where the n's give the number of each type of response.

2)

and a simpler question : How to use constraints in rrvglm?

for example, i want to regress X1 X2 X3 X4 (poisson variables with log-link) on two factor-type variables N1, N2 but i would like to use X1 and X2 [or something like 1(X1>X2)] in the estimation of X3 and X4 (but not in X1 X2 estimation of course).

i can do it with simple vglm (with the constraints matrices) but with rr-vglm i can't.

Excuse me if my questions are silly and for my bad english.

Benoit J.
R-fan

PS: the final purpose would be to regress something like X1,X2,X3,X4 1(X1>X2) 1(X3>X4) on N1, N2 (for the part depending on X1,X2) and N1,N2,X1,X2 [or N1,N2,1(X1>X2) ] for the responses depending on X3,X4.

X3,X4 could be interpreted like second period variables and X1,X2 like first period variables.

-- 
View this message in context: http://www.nabble.com/R-vglm-new-family-writing%3A-mix-Poisson-multinomial-tp14660710p14660710.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Mon 07 Jan 2008 - 09:31:28 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 07 Jan 2008 - 10:30:05 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive