search

dlinzer / poLCA

Polytomous Variable Latent Class Analysis (R package)
https://dlinzer.github.io/poLCA/
48 stars 17 forks source link

add p-values for Chisq and Gsq #7

Open mikebarkmin opened 6 years ago

mikebarkmin commented 6 years ago

Calculate the p-values for Chisq and Gsq like:

C <- max(K.j) # number of categories
I <- J # number of items
df <- C^I - ret$npar - 1 # Degrees of freedom
Chisq.pvalue <- 1-pchisq(ret$Chisq,df)
Gsq.pvalue <- 1-pchisq(ret$Gsq,df)
anamgreco commented 2 years ago

Hello! Sorry but I am not able to compute my p-values from a poLCA object with this script... Is this an answer to a previous one? How can I adapt this to mine, please? Thanks!!

mikebarkmin commented 2 years ago

You can not use the lines after running poLCA. They are intended to run along. I have already modified poLCA to calculate the p-values (https://github.com/dlinzer/poLCA/pull/14), but the maintainers of this library seem to have abandoned it. So you probably must modify poLCA yourself.

anamgreco commented 2 years ago

Hello, again! Sorry for my delay, and thanks for your response and all your help. I managed to re-write the package by adding the three lines (183-185) showed here

However, now all p-values returned are the same for very different models (from 1 to 10 classes). Is that correct? I am sorry is the very first time I am using this indicator...

anamgreco commented 2 years ago

Hello, again! Sorry for my delay, and thanks for your response and all your help. I managed to re-write the package by adding the three lines (183-185) shown herehttps://github.com/dlinzer/poLCA/pull/14/commits/c0ed04aae02e81faa83e046034cb08a8b03092c7 However, now all p-values returned are the same for very different models (from 1 to 8 classes). Is that correct? I am sorry is the very first time I am using this indicator...

I can send database and code to check for reproducibility, but I do not want to ask you too much. Sorry again and thank you in advance!

Ana M. Greco- Dra. en Cerebro, cognición y conducta Grup de Recerca en Victimització Infantil i Adolescent (GReVIA) 722714704

De: Mike Barkmin @.> Enviado el: martes, 7 de junio de 2022 12:37 Para: dlinzer/poLCA @.> CC: Ana Martina Greco @.>; Comment @.> Asunto: Re: [dlinzer/poLCA] add p-values for Chisq and Gsq (#7)

You can not use the lines after running poLCA. They are intended to run along. I have already modified poLCA to calculate the p-values (#14https://github.com/dlinzer/poLCA/pull/14), but the maintainers of this library seem to have abandoned it. So you probably must modify poLCA yourself.

— Reply to this email directly, view it on GitHubhttps://github.com/dlinzer/poLCA/issues/7#issuecomment-1148493457, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AOAZ6CW3YVKFLEZQ5SOR2NLVN4Q4NANCNFSM4ESI6SFA. You are receiving this because you commented.Message ID: @.**@.>>

Aquest missatge, i els fitxers adjunts que hi pugui haver, pot contenir informació confidencial o protegida legalment i s’adreça exclusivament a la persona o entitat destinatària. Si no consteu com a destinatari final o no teniu l’encàrrec de rebre’l, no esteu autoritzat a llegir-lo, retenir-lo, modificar-lo, distribuir-lo, copiar-lo ni a revelar-ne el contingut. Si l’heu rebut per error, informeu-ne el remitent i elimineu del sistema tant el missatge com els fitxers adjunts que hi pugui haver.

Este mensaje, y los ficheros adjuntos que pueda incluir, puede contener información confidencial o legalmente protegida y está exclusivamente dirigido a la persona o entidad destinataria. Si usted no consta como destinatario final ni es la persona encargada de recibirlo, no está autorizado a leerlo, retenerlo, modificarlo, distribuirlo o copiarlo, ni a revelar su contenido. Si lo ha recibido por error, informe de ello al remitente y elimine del sistema tanto el mensaje como los ficheros adjuntos que pueda contener.

This email message and any attachments it carries may contain confidential or legally protected material and are intended solely for the individual or organization to whom they are addressed. If you are not the intended recipient of this message or the person responsible for processing it, then you are not authorized to read, save, modify, send, copy or disclose any part of it. If you have received the message by mistake, please inform the sender of this and eliminate the message and any attachments it carries from your account.

tbmpereira commented 1 year ago

@mikebarkmin

I wanted to thank you for your contribution to the poLCA package through the pull request you made to calculate Chisq and Gsq p-values. Your insights were very helpful!

I also wanted to mention that it would be great if the package maintainer could include p-values in the returned output. However, I wanted to point out that the most common method for estimating p-values in the case of LCA is through the parametric bootstrap method. While your suggested method is also valid, it may not be the most appropriate in this particular case.

Thank you again for your contribution to the package, and I look forward to seeing how it continues to evolve!