[R-br] GOM.EM
Julimar Pinto
julimarsp.r em gmail.com
Terça Agosto 23 17:13:20 BRT 2016
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/20160823/98fd2375/attachment.html>
Mais detalhes sobre a lista de discussão R-br