Questions tagged [quasi-likelihood]
In GLMs, quasi-likelihood estimation is a way to allow over- or under-dispersion by choosing an appropriate variance function.
151
questions
0
votes
0
answers
43
views
Approximating profile likelihood confidence intervals when I only have a score function and not a likelihood
I'm working on a modeling problem where I can define a score function that looks a lot like a binomial likelihood, but the model isn't really binomial. I'd like to use profile likelihood to estimate ...
- 332
2
votes
1
answer
33
views
Likelihood ratio test for model specification with boundary Null
I am interested in understanding the asymptotic distribution of Likelihood ratio (LR) test statistic for model specification. I am focusing on the case in which the null hypothesis is of the form (i.e....
- 21
2
votes
2
answers
216
views
Poisson model for non-integer
I have a GLM with (quasi)poisson family.
My dataset has 3 variables:
rate_data
rate_benchmark
X
So fitting the model:
...
- 23
0
votes
1
answer
33
views
Why does using proportion DV (vs. count) change the results?
I'm running a regression model using a count variable as the DV and I'm trying to understand what changes when using it as a count outcome vs. a proportion.
REF: https://ademos.people.uic.edu/...
- 331
0
votes
0
answers
33
views
How can I estimate the Relative Risk using quasi-Poisson regression?
In studies that relate pollution to hospitalizations for certain diseases, the authors usually present the results in the form of relative risk for a 95% confidence interval.
If I have the regression ...
- 21
2
votes
1
answer
27
views
Can I apply quasi-Poisson regression on discontinuous data?
I'm trying to relate the concentrations of certain air pollutants to the incidence of a disease, but I don't have continuous data on that pollution. I have daily data on cases and deaths from the ...
- 21
1
vote
0
answers
54
views
Interpretation of p-values in glm() and sjPlot::tab_model() for quasi-poisson regression
I am running a quasi-poisson regression predicting a count outcome from 4 variables of interest.
I understand that these estimates need to be exponentiated to correctly interpret the results. I used ...
- 21
0
votes
0
answers
20
views
Why not characterise quasi-Poissons by an actual density?
In the discussions I have seen about quasi-Poisson regressin methods, I glean that the method is to merely assume (in the standard notation) $E(Y_i) = \phi \text{Var}(Y_i)$ given the explanatory ...
- 101
1
vote
1
answer
42
views
Comparing two quasipoisson glm models
Have these results but not sure how to interpret them to pick the best fitting model, I think the high p value suggests I should go with model 1? Also in terms of the order that I put the models in, ...
user358155
1
vote
1
answer
65
views
Bnomial (Logit) regression for proportion/percentage data
I have run a binomial (logit) regression on some proportion data as the dependent variable in an Interrupted Time Seies Analysis [see below]:
...
- 41
2
votes
0
answers
76
views
What are the assumptions of fractional response regression model?
I employed fractional response regression proposed by Wooldridge and Papke for my study as the dependent variable is a proportion (remaining/whole) and is between 0 to 1 (including 0 and 1). I want to ...
- 143
2
votes
0
answers
39
views
Is Quasi-Likelihood As "Strong" As The Standard Likelihood?
I have been trying to learn more about the differences between Quasi Likelihood compared to the standard Likelihood. When learning about this, the following points come to mind:
The Quasi-Likelihood ...
- 7,048
1
vote
1
answer
59
views
WHy is the over dispersion in this poisson and quasi-poisson the same?
I have a zero inflated count data, on which I have run a poisson and quasi poisson reg using glm().
The output from a poisson model is as follows:
...
- 15
1
vote
1
answer
41
views
Terminology question regarding a certain "partial maximum likelihood" which approximates the marginal likelihood
Suppose that we have a model with many parameters, which we'll partition into two subvectors called $\theta$ and $\lambda$. In this situation, $\lambda$ corresponds to those parameters that are really ...
- 207
1
vote
0
answers
201
views
Fitting a regression model on bounded data
Typically, in a regression equation, the dependent variable is unbounded i.e. follows normal distribution. But sometimes it may happen that dependent variable is bounded i.e. dependent variable is ...
- 841
1
vote
0
answers
26
views
How to interprete the outcome for quasi-binomial models
I am currently trying to understand the quasi-binomial model and have a rather basic question, which I seem not to be able to see an answer. The outcomes in binomial models are binary. So, 0 vs 1, car ...
- 73
1
vote
0
answers
217
views
Serious Coding Error in QIC function in geepack? [closed]
I believe the QIC function in geepack has a significant error. The function appears to incorrectly specify the independence model, which is needed to calculate QIC. The function will therefore often ...
- 51
0
votes
0
answers
25
views
Can the offset variable be smaller than the outcome variable in a Poisson model?
I am modeling a number of bags of local crop harvested from a village in a given harvest over time.
The village is divided into 48 communities and I have data on the total number of bags of crops from ...
- 13
0
votes
1
answer
63
views
quasi-likelihood estimates of beta
working through Peter McCullagh's glm book and having a hard time with understanding quasi-likelihood. I'm working on this question below and I think I need to find the quasi-likelihood estimates and ...
0
votes
1
answer
139
views
quasi likelihood for ungrouped binary data
I read in one of the textbooks that for ungrouped binary data the dispersion parameter should always be $\phi = 1$.
Do you know why it is the case?
- 641
0
votes
0
answers
48
views
Find the optimal linear combination of the following estimating function / quasi-score
Consider a random variable $Y_1 \sim Bin(n,\theta)$ and $Y_2 \sim Bin(n,\theta^2)$ where $Y_1$ is independent from $Y_2$.
Consider the residuals $R_1 = Y_1 - n\theta$ and $R_2 = Y_2 - n\theta^2$. Find ...
0
votes
2
answers
555
views
Appropriate Regression Model for Proportions and Rate Data
I have a problem where my dependent variable is given as a click-through rate and thus bounded [0,1]. While I have the traffic for each sample (a combination of design factors) and could reconstruct a ...
- 1
2
votes
0
answers
307
views
Comparing performance of Quasi-binomial model and Beta-binomial model
I read some books in biostatistics about fitting binary date with Beta-Binomial regression model and Quasi-Binomial regression model. It proposes a setting:
Setting:
Assuming we have a sequence of ...
- 31
2
votes
1
answer
97
views
Correct GLM or NLS to model exponential model with response variable with positive and negative values
I have been struggling to find the right way to model this dataset, this is a Data Frame with the dataset:
...
- 343
1
vote
1
answer
85
views
How to calculate model weight for GLM-quasi poisson model
I am running several models with the quasi-Poisson family, I have calculated QAIC for each model but I wanted to know the weight of each individual model. I tried AICcmodavg but it did not work.
Is ...
- 21
0
votes
1
answer
80
views
Quasi-likelihood can't be generated by any valid probability distribution
I am learning about quasi-Poisson and i'm stuck at the concept of quasi-likelihood function.
In wikipedia, it is said that:
The term quasi-likelihood function was introduced by Robert Wedderburn in ...
2
votes
1
answer
85
views
Mean-variance relationship in the quasi-likelihood
I have some questions regarding the quasi-likelihood model of GLM:
I understand that one reason to use quasi-likelihood in GLM is over-dispersion. This seems to justify using the quasi-Poisson, or ...
- 2,870
2
votes
1
answer
147
views
How to compute the gradient for a GARCH with the package rugarch in R
I am estimating a GARCH(1,1) with external regressors and the package rugarch allows me to do it easily.
However, to compute QMLE robust standard errors, I need the ...
2
votes
1
answer
457
views
How to compute the sandwich variance ML estimator in R
I'm currently estimating a DCC-type model by maximum likelihood. Im using the command solnp and it return an object where I can compute the Hessian H evaluated at ...
0
votes
1
answer
133
views
Test for significant differences for data between 0 and 1
I have to test for significant differences between scenarios. Data consist of the length of a segment divided by the total length of the network. They are distributed between 0 (never equal to 0) and ...
0
votes
1
answer
189
views
Output in two-part Fractional Response Model (FRM) package incorrect?
I was experimenting with the Fractional Response Model (FRM) package, and decided to replicate the results using the base GLM package to better understand the theory. I am able to replicate the ...
- 113
3
votes
1
answer
265
views
Quasi-likelihood function
I got stuck in the derivation of the quasi-likelihood function. Namely, given an i.i.d sample $\{Y_i,X_i \}_{i=1}^n$ with $n$ the sample size, let the conditional mean and variance functions be ...
- 132
0
votes
0
answers
55
views
Comparing the marginal effect of a GLM to the OLS estimates
My question is, whether there is any way to (somewhat) compare the marginal effect of a GLM estimate to an OLS estimate. As in, "since the OLS and GLM results are very similar, I will favour OLS ...
- 359
1
vote
1
answer
705
views
Comparing the marginal effects of glm output to polr output
I have a dependent variable that is technically ordinal, so I ran a ordered probit model (polr).
However, an ordered probit model does not produce any residuals ...
- 359
1
vote
0
answers
135
views
How flexible is Stata's ivpois? Could I use it for a (quasi) binomial distribution?
According to this post on statalist, Stata's ivpois (an instrumental variable approach) is pretty flexible, with very little assumptions.
The problem mentioned in ...
- 359
1
vote
0
answers
451
views
How to do a Control Function (CF) / Two Stage Residual Inclusion (2SRI) with an ordinal dependent variable in the first stage and a glm in the second
I am trying to use a Control Function (CF) / Two Stage Residual Inclusion (2SRI) approach, because the modeled relationship that I am trying to estimate is non-linear (my dependent variable has a ...
- 359
1
vote
2
answers
382
views
Getting the (Stata) margins from fractional regression (=glm with family quasibinomial) for an ordinal variable in R
I first found this really nice Stata video on fractional regression (the dependent variable is a proportion including 0 and 1). I am especially interested in how he applies the margin approach to ...
- 359
3
votes
2
answers
639
views
What is the interpretation of a glm coefficient on a dependent variable that has a % interpretation
I have a dependent variable that takes on values between 0 and 1, including 0 and 1.
The variable signifies a proportion (0 = nothing, 1 = all). I am running a model of the type:
...
- 359
0
votes
1
answer
101
views
Changes in significance between Poisson & quasipoisson glm
I am fairly new to GLMs, and am currently practising and testing with an insurance dataset, after many tries, I am modeling the frequency (counting model of the number of claims) and I have several ...
0
votes
0
answers
25
views
A problem of likelihood function in a dynamic setting?
I'm having a problem regarding perhaps conditional maximum likelihood problem, but I'm not sure.
Suppose time horizon we consider is $T=4$, our goal is to minimize the loss function
$$
\sum_{t=2}^T L(...
- 101
1
vote
0
answers
75
views
Relationship between quasi-likelihood and "GLM-conjugate" models?
Suppose we have a response variable that represents proportions with a poorly defined denominator. Two ways to handle this (1) a quasibinomial model, which assumes only that the variance is ...
- 38k
0
votes
0
answers
32
views
Is there any sense of applying cross validation to quasibinomial variable?
We can apply without any doubts cross validation for linear models as well as for binary models.
For linear models for example we can output RMSE and MAE and Accuracy for binary models.
But I have ...
- 442
1
vote
1
answer
130
views
Comparing model efficiency
I hope you all don't mind me asking this question.
I have two models :
general linear mixed effects model
...
- 477
4
votes
0
answers
197
views
Why is my quasibinomial GLM estimator biased - Monte Carlo simulation
I'm playing with some Monte Carlo simulations to get an idea of the properties of some linear and non-linear models.
The linear OLS model in my case is specified as:
$Y_t = \beta_0 + \beta_1x+ \...
- 93
1
vote
0
answers
44
views
If my data is overdispersed when binomial but not when Poisson does that alone mean I should use Poisson GLM?
Top is Poisson, bottom is Binomial
Source of variation is different people observed over a week long period. 118 people to be exact.
want to make predictions about the proportion of days/week one ...
- 31
1
vote
0
answers
81
views
Avoiding data dredging with quasibinomial GLMs
Because my data is a bit overdispersed I am using quasibinomial GLM to analyse them however this means I cannot use AICs to compare my models.
I am therefore using drop1 and update functions to do ...
- 31
1
vote
0
answers
69
views
Equivalence between prediction interval and hypothesis test on Poisson predicted values with 2/3 power transformation
I fitted a quasi-Poisson regression model to predict the value of a certain counting variable.
Once the predicted value $\widehat{\mu_{0}}$ was obtained, I calculated the upper limit of its prediction ...
- 93
1
vote
0
answers
68
views
Overdispersion in logistic model
I'm relatively a newbie in R, and I've been trying to make a silly example of logistic regression to predict, according to Age and Sex whether someone dies of corona or not. I'm from Colombia, so my ...
0
votes
1
answer
888
views
Interpreting coefficents on quasibinomial model
This might seem like a pretty basic question but I've scoured seemingly everywhere and can't get a definitive answer.
I have a response parameter "rr" bound by 0-1 which is essentially the ...
1
vote
1
answer
536
views
