add p-values for Chisq and Gsq

[mikebarkmin]

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]

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]

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 (#14), but the maintainers of this library seem to have abandoned it. So you probably must modify poLCA yourself.

[anamgreco]

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]

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 here<c0ed04a> 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 (#14<#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 GitHub<#7 (comment)>, or unsubscribe<https://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.

[marcelo-pereira-pcst]

@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!