Opened 14 years ago
Last modified 14 years ago
#1050 accepted task
Funcion GLM ( General Linear Modelo) de R con el API de TOL
Reported by: | jlaybar | Owned by: | Víctor de Buen Remiro |
---|---|---|---|
Priority: | lowest | Milestone: | TOL Packages |
Component: | Math | Version: | |
Severity: | minor | Keywords: | |
Cc: |
Description
Hola
Os envio la función GLM ( General Linear Models) de R escrita con el API de tol
Puede realizar modelos Logit y Probit con la familia Binomial
" family = binomial(link=probit)\n"
" family = binomial(link=logit)\n"
Para el resto de familias hay que modificar el archivo Rglm.tol y cambiar el case donde pone:
Real mod = Case ( model == "Logit" , 1,
model == "Probit", 2,
1 , 0
);
" family = gaussian(link=identity)\n"
" family = Gamma(link =inverse)\n"
" family = inverse.gaussian(link = 1/mu2)\n"
" family = poisson(link = log)\n"
Un saludo
Jaime
Attachments (1)
Change History (4)
Changed 14 years ago by
comment:1 Changed 14 years ago by
Priority: | highest → lowest |
---|---|
Severity: | blocker → minor |
comment:2 Changed 14 years ago by
Status: | new → accepted |
---|
comment:3 Changed 14 years ago by
Component: | ASCII data files → Math |
---|---|
Milestone: | → TOL Packages |
Type: | defect → task |
Buen trabajo Jaime.
Justo acabamos de crear el paquete QltvRespModel que va en esa dirección pero siempre viene bien contrastar diferentes soluciones y además el GLM nos da más generalidad por el momento.
En favor de QltvRespModel hay que decir que permite introducir también restricciones de desigualdad lineal e información a priori, y también que proporciona estimación máximo verosímil y simulación bayesiana, pero de momento sólo admite Logit y Probit, aunque está preparado para generalizarse.
Si te parece podemos crear un paquete con lo que has hecho para que todo el mundo pueda utilizarlo de una forma más cómoda. Yo te puedo guiar y ayudar a diseñarlo.