[R-br] GOM.EM
Anna Karoline R. da Cruz
cruzkarol em hotmail.com
Quarta Agosto 24 21:54:57 BRT 2016
Olá Julimar!
Obrigada pelo retorno.
Então, seria interessante conhecer essa versão. Mas pelo que entendi, ainda não tem um pacote no R...
Como faço para usá-lo?
Abs.
Anna.
From: julimarsp.r em gmail.com
Date: Tue, 23 Aug 2016 17:13:20 -0300
Subject: Re: [R-br] GOM.EM
To: cruzkarol em hotmail.com; r-br em listas.c3sl.ufpr.br
Olá Anna!
Eu até tenho conhecimento da existência dessa função gom.em() contida na package "sirt", que foi construída pelo alemão Alexander Robitzsch.
Contudo, eu utilizo uma modelagem um pouquinho diferente. A função gom.em() utiliza-se de modelos mistura e tem sua forma baseada no modelo preconizado pela Erosheva [1].
O modelo que utilizo baseia-se no modelo preconizado por WOODBURY & CLIVE [2] [3] [4]. É uma versão que eu mesmo escrevi do software GOM3.EXE [5]. O software GOM3.EXE foi desenvolvido em ambiente MS-DOS para estimar os modelos Grade of Membership (GoM). Este software, que possui código fonte aberto e livre, foi escrito originalmente por Woodbury e Clive em 1975, posteriormente modificado por Peter Charpentier e Burt Singer no Departamento de Epidemiologia e Saúde Pública da Escola de Medicina da Universidade de Yale – EUA, (1996).
Neste último caso, a interpretação dos lambdas em GOM3.EXE e gom3b [6] são idênticas, pois pertecem ao mesmo tipo de modelo. Nesse contexto, a solução desenvolvida por André Junqueira Caetano e eu na PUC-Minas possui a mesma interpretação de GOM3.EXE. Eu diria mesmo que "GoMRcpp.R" é um GOM3.EXE "bombado".
A Daniela Petruzalek e eu estamos conversando para transformar o R-Script "GoMRcpp.R" em uma package do R. Melhorando a documentação e o código também.
Se essa forma de interpretação, provida pelo "GoMRcpp.R", lhe for útil terei prazer em lhe ajudar a compreender seus dados.
[1] Erosheva, E. A., Fienberg, S. E., & Joutard, C. (2007). Describing disability through individuallevel mixture models for multivariate binary data. Annals of Applied Statistics, 1, 502-537.
[2] WOODBURY, M. A.; CLIVE, J. Clinical pure types as a fuzzy partition. Journal of Cybernetics, 4, p. 111-121, 1974.
[3] WOODBURY, M. A.; CLIVE, J.; GARSON JR., A. Mathematical typology: a grade of membership technique for obtaining disease definition. Computers and Biomedical Research, 11:277-298, 1978.
[4] MANTON, Kenneth G.; WOOBURY, Max A.; TOLLEY, H. Dennis. Statistical applications using fuzzy sets. New York: John Wiley & Sons, 1994.
[5] CHARPENTIER, Peter; SINGER, Burt. GOM3.EXE. Yale, Departamento de Epidemiologia e Saúde Pública da Escola de Medicina da Universidade de Yale – EUA, 1996. Disponível para download em: <http://lib.stat.cmu.edu/DOS/general/.index.html>. Acesso em: 17 fev. 2011.
[6] LAGE, Rafael Kelles Vieira. gom3b. Belo Horizonte: s/n, 2005. 1 cd-rom.
[7] PINTO, Julimar Santos; CAETANO, André Junqueira. A Heterogeneidade da Vulnerabilidade Social das Juventudes: Uma Perspectiva Empírica Através do Método Grade of Membership. Revista Mediações, Londrina, v. 18, n. 1, 2013. ISSN: 1414-0543. Disponível para download em: <https://github.com/julimarsp-into-r>. Acesso em: 23 de ago. 2016.
Warm Regards,
Julimar Pinto
Em 22/08/2016 14:45, "Anna Karoline R. da Cruz via R-br" <r-br em listas.c3sl.ufpr.br> escreveu:
Boa tarde, pessoal!
Estou usando o "GoM.em" para um conjunto de dados.
Minha dúvida é: como descrevo os perfis a partir do lambda que ele fornece?
OBS: tenho algumas as variáveis com mais de duas categorias.
> mod1$lambda
[,1] [,2] [,3]
[1,] 0.125000000 0.375000000 0.62500000
[2,] 0.999000000 0.758576148 0.94791392
[3,] 0.999000000 0.999000000 0.99900000
[4,] 0.999000000 0.192048304 0.99900000
[5,] 0.958999554 0.239871163 0.99900000
[6,] 0.999000000 0.582693490 0.60513053
[7,] 0.760374047 0.999000000 0.99900000
[8,] 0.001000000 0.046672897 0.90641452
[9,] 0.475789685 0.883268404 0.10884461
[10,] 0.355654569 0.232261520 0.05983177
[11,] 0.236404268 0.001038477 0.12928587
[12,] 0.009139838 0.364151602 0.28363626
[13,] 0.001000000 0.604041999 0.99900000
[14,] 0.268329215 0.339902463 0.99900000
[15,] 0.001000000 0.001000000 0.78822352
[16,] 0.031297098 0.338268504 0.99900000
[17,] 0.001000000 0.001000000 0.26163060
[18,] 0.999000000 0.999000000 0.99900000
[19,] 0.001000000 0.001000000 0.13074344
[20,] 0.999000000 0.999000000 0.99900000
[21,] 0.999000000 0.999000000 0.99900000
[22,] 0.999000000 0.999000000 0.99900000
[23,] 0.001000000 0.001000000 0.00100000
[24,] 0.768230734 0.999000000 0.47083018
[25,] 0.001000000 0.001000000 0.39344518
[26,] 0.532563903 0.999000000 0.00100000
[27,] 0.001000000 0.001000000 0.00100000
[28,] 0.289341873 0.999000000 0.21041057
[29,] 0.001000000 0.001000000 0.13074344
[30,] 0.001000000 0.790611117 0.00100000
[31,] 0.001000000 0.001000000 0.00100000
[32,] 0.531285647 0.999000000 0.00100000
[33,] 0.001000000 0.001000000 0.26163060
[34,] 0.283878725 0.999000000 0.00100000
Desde já agradeço a ajuda!
Abs.
Anna.
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160825/b8909ee4/attachment.html>
Mais detalhes sobre a lista de discussão R-br