Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Is there a may to run in Lmer a log-gaussian link function? #102

Open
PuddleJumper2018 opened this issue Mar 15, 2022 · 3 comments
Open

Comments

@PuddleJumper2018
Copy link

I'm currently using pymer4 for quite a big dataset with over 4000 samples and (obviously) have distribution issues. Is there a way to specify model link function as log-gaussian? I believe in R glmer() takes it as gaussian(link = "log")

@ejolly
Copy link
Owner

ejolly commented Mar 15, 2022

Unfortunately currently not. We only support the following families with their default link functions:

implemented_fams = [
            "gaussian",
            "binomial",
            "gamma",
            "inverse_gaussian",
            "poisson",
        ]

I can look into adding an update that supports using a non-default link function.

@DaviYokogawa
Copy link

@ejolly where can I change the code to implement this?
My idea is use the families like the control, it's possible?

@ejolly
Copy link
Owner

ejolly commented Mar 29, 2023

@DaviYokogawa Sorry for the delay. If you're interested in those changes I would check out lines 424 and 440 of Lmer.py. Those are the lines that run the model similar to calling lmer or glmer in R and should be able to handle the same arguments. Though you'll want to check out the rpy2 docs for how to structure those arguments.

@ejolly ejolly added this to Roadmap Sep 8, 2023
@ejolly ejolly moved this to Long Term maybe in Roadmap Sep 8, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
Status: Long Term maybe
Development

No branches or pull requests

3 participants