[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