[R-sig-ME] Random Effect Tobit

Dimitris Rizopoulos d@r|zopou|o@ @end|ng |rom er@@mu@mc@n|
Wed Sep 15 03:36:04 CEST 2021


Yes, it seems to me that the non-adaptive GH is used according to procedure described in Section 3.4 here: https://cran.r-project.org/web/packages/censReg/vignettes/censReg.pdf

The adaptive GH updates the position of the abscissas and the weights during the optimization to achieve a better approximation of the log-likelihood.



��
Dimitris Rizopoulos
Professor of Biostatistics
Erasmus University Medical Center
The Netherlands
________________________________
From: Ben Bolker <bbolker using gmail.com>
Sent: Wednesday, September 15, 2021 3:28:43 AM
To: Dimitris Rizopoulos <d.rizopoulos using erasmusmc.nl>; r-sig-mixed-models using r-project.org <r-sig-mixed-models using r-project.org>
Subject: Re: [R-sig-ME] Random Effect Tobit

    Can you explain a bit more? Is this the distinction that is
sometimes made between Gauss-Hermite quadrature and *adaptive* GHQ, i.e.
whether the expansion is centered at the population-level value or at
zero?  (If so, I would *definitely* recommend GLMMadaptive, non-adaptive
GHQ can be terrible in some cases ...)

   cheers
    Ben Bolker

On 9/14/21 9:26 PM, Dimitris Rizopoulos wrote:
> Indeed, but the censReg does not seem to update the abscissas and
> weights during the optimization procedure.
>
> Best,
> Dimitris
>
> ��
> Dimitris Rizopoulos
> Professor of Biostatistics
> Erasmus University Medical Center
> The Netherlands
> ------------------------------------------------------------------------
> *From:* R-sig-mixed-models <r-sig-mixed-models-bounces using r-project.org> on
> behalf of Ben Bolker <bbolker using gmail.com>
> *Sent:* Wednesday, September 15, 2021 2:08:13 AM
> *To:* r-sig-mixed-models using r-project.org <r-sig-mixed-models using r-project.org>
> *Subject:* Re: [R-sig-ME] Random Effect Tobit
> This CrossValidated answer:
> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstats.stackexchange.com%2Fquestions%2F544511%2Fzero-inflated-gaussian-for-weights-below-zero-recorded-as-0%3Fnoredirect%3D1%23comment999636_544511&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=c4FixGlxg5IY%2BOsgfi2CfQsBANPImB5n6k8eT7Ordag%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstats.stackexchange.com%2Fquestions%2F544511%2Fzero-inflated-gaussian-for-weights-below-zero-recorded-as-0%3Fnoredirect%3D1%23comment999636_544511&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=c4FixGlxg5IY%2BOsgfi2CfQsBANPImB5n6k8eT7Ordag%3D&reserved=0>
>
>    says that the censReg package can handle random effects Tobit models
> via Gauss-Hermite Quadrature (GLMMadaptive is great, but it's nice to
> have options)
>
> On 9/14/21 3:32 AM, Dimitris Rizopoulos wrote:
>> You can use the GLMMadaptive package (https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdrizopoulos.github.io%2FGLMMadaptive%2F&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=vPMiaxS3qhMVxCTCCk%2F25Rofp80fuyVxvc4BQTuPaQ0%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdrizopoulos.github.io%2FGLMMadaptive%2F&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=vPMiaxS3qhMVxCTCCk%2F25Rofp80fuyVxvc4BQTuPaQ0%3D&reserved=0>)
> and the censored.normal() family object. An example is given here:
> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdrizopoulos.github.io%2FJMbayes2%2Farticles%2FNon_Gaussian_Mixed_Models.html%23censored-linear-mixed-models-1&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=EHNedoTKhPZCOb20NCgnUMwovT%2FNQlqJsfWL3nK%2BzV0%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdrizopoulos.github.io%2FJMbayes2%2Farticles%2FNon_Gaussian_Mixed_Models.html%23censored-linear-mixed-models-1&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=EHNedoTKhPZCOb20NCgnUMwovT%2FNQlqJsfWL3nK%2BzV0%3D&reserved=0>
>>
>>
>> Best,
>> Dimitris
>>
>> -----Original Message-----
>> From: R-sig-mixed-models <r-sig-mixed-models-bounces using r-project.org> On Behalf Of Hedyeh Ahmadi
>> Sent: Tuesday, September 14, 2021 12:25 AM
>> To: r-sig-mixed-models using r-project.org
>> Subject: [R-sig-ME] Random Effect Tobit
>>
>> Hello all,
>> I was wonder if anyone is aware of an R package with random effect Tobit. Any help would be greatly appreciated.
>>
>> Best,
>>
>> Hedyeh Ahmadi, Ph.D.
>> Statistician
>> Keck School of Medicine
>> Department of Preventive Medicine
>> University of Southern California
>>
>> Postdoctoral Scholar
>> Institute for Interdisciplinary Salivary Bioscience Research (IISBR) University of California, Irvine
>>
>> LinkedIn
>> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.linkedin.com%2Fin%2Fhedyeh-ahmadi&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=sgV%2BBBAHfwTRqn9fTiOYsGsIgkvsOnNLaAJ9hAziMv8%3D&reserved=0<https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.linkedin.com%2Fin%2Fhedyeh-ahmadi&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=sgV%2BBBAHfwTRqn9fTiOYsGsIgkvsOnNLaAJ9hAziMv8%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.linkedin.com%2Fin%2Fhedyeh-ahmadi&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=sgV%2BBBAHfwTRqn9fTiOYsGsIgkvsOnNLaAJ9hAziMv8%3D&reserved=0<https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.linkedin.com%2Fin%2Fhedyeh-ahmadi&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=sgV%2BBBAHfwTRqn9fTiOYsGsIgkvsOnNLaAJ9hAziMv8%3D&reserved=0>>
>>
>>
>>        [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-sig-mixed-models using r-project.org mailing list
>> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=JeXCvNShZ9xfCFlX7qBV%2BlR7RguHfsgEAoSPHBJzI%2FY%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=JeXCvNShZ9xfCFlX7qBV%2BlR7RguHfsgEAoSPHBJzI%2FY%3D&reserved=0>
>>
>> _______________________________________________
>> R-sig-mixed-models using r-project.org mailing list
>> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273318983%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=JeXCvNShZ9xfCFlX7qBV%2BlR7RguHfsgEAoSPHBJzI%2FY%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273328938%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Bla%2FRWD6uYWMRaG6dLGKyoOVzMuWrXH%2B4HzftJmbLk4%3D&reserved=0>
>>
>
> --
> Dr. Benjamin Bolker
> Professor, Mathematics & Statistics and Biology, McMaster University
> Director, School of Computational Science and Engineering
> Graduate chair, Mathematics & Statistics
>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273328938%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Bla%2FRWD6uYWMRaG6dLGKyoOVzMuWrXH%2B4HzftJmbLk4%3D&reserved=0
> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-mixed-models&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C596d927815de4130a05708d977e8293c%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C1%7C637672661273328938%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Bla%2FRWD6uYWMRaG6dLGKyoOVzMuWrXH%2B4HzftJmbLk4%3D&reserved=0>

--
Dr. Benjamin Bolker
Professor, Mathematics & Statistics and Biology, McMaster University
Director, School of Computational Science and Engineering
Graduate chair, Mathematics & Statistics

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list