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@gmail.com Date: Tue, 23 Aug 2016 17:13:20 -0300 Subject: Re: [R-br] GOM.EM To: cruzkarol@hotmail.com; r-br@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@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@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.