When to use a GAM vs GLM
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty{ margin-bottom:0;
}
up vote
7
down vote
favorite
I realize this may be a potentially broad question, but I was wondering whether there are generalizable assumptions that indicate the use of a GAM (Generalized additive model) over a GLM (Generalized linear model)?
Someone recently told me that GAMs should only be used when I assume the data structure to be "additive", i.e. I expect additions of x to predict y.
Another person pointed out that a GAM does a different type of regression analysis than a GLM, and that a GLM is preferred when linearity can be assumed.
In the past I have been using a GAM for ecological data, e.g.:
- continuous timeseries
- when the data did not have a linear shape
- I had multiple x to predict my y that I thought to have some nonlinear interaction that I could visualize using "surface plots" together with a statistical test
I obviously don't have a great understanding of what a GAM does different than a GLM. I believe it's a valid statistical test, (and I see an increase in the use GAMs, at least in ecological journals), but I need to know better when it's use is indicated over other regression analyses.
regression generalized-linear-model gam
New contributor
add a comment |
up vote
7
down vote
favorite
I realize this may be a potentially broad question, but I was wondering whether there are generalizable assumptions that indicate the use of a GAM (Generalized additive model) over a GLM (Generalized linear model)?
Someone recently told me that GAMs should only be used when I assume the data structure to be "additive", i.e. I expect additions of x to predict y.
Another person pointed out that a GAM does a different type of regression analysis than a GLM, and that a GLM is preferred when linearity can be assumed.
In the past I have been using a GAM for ecological data, e.g.:
- continuous timeseries
- when the data did not have a linear shape
- I had multiple x to predict my y that I thought to have some nonlinear interaction that I could visualize using "surface plots" together with a statistical test
I obviously don't have a great understanding of what a GAM does different than a GLM. I believe it's a valid statistical test, (and I see an increase in the use GAMs, at least in ecological journals), but I need to know better when it's use is indicated over other regression analyses.
regression generalized-linear-model gam
New contributor
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
2
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday
add a comment |
up vote
7
down vote
favorite
up vote
7
down vote
favorite
I realize this may be a potentially broad question, but I was wondering whether there are generalizable assumptions that indicate the use of a GAM (Generalized additive model) over a GLM (Generalized linear model)?
Someone recently told me that GAMs should only be used when I assume the data structure to be "additive", i.e. I expect additions of x to predict y.
Another person pointed out that a GAM does a different type of regression analysis than a GLM, and that a GLM is preferred when linearity can be assumed.
In the past I have been using a GAM for ecological data, e.g.:
- continuous timeseries
- when the data did not have a linear shape
- I had multiple x to predict my y that I thought to have some nonlinear interaction that I could visualize using "surface plots" together with a statistical test
I obviously don't have a great understanding of what a GAM does different than a GLM. I believe it's a valid statistical test, (and I see an increase in the use GAMs, at least in ecological journals), but I need to know better when it's use is indicated over other regression analyses.
regression generalized-linear-model gam
New contributor
I realize this may be a potentially broad question, but I was wondering whether there are generalizable assumptions that indicate the use of a GAM (Generalized additive model) over a GLM (Generalized linear model)?
Someone recently told me that GAMs should only be used when I assume the data structure to be "additive", i.e. I expect additions of x to predict y.
Another person pointed out that a GAM does a different type of regression analysis than a GLM, and that a GLM is preferred when linearity can be assumed.
In the past I have been using a GAM for ecological data, e.g.:
- continuous timeseries
- when the data did not have a linear shape
- I had multiple x to predict my y that I thought to have some nonlinear interaction that I could visualize using "surface plots" together with a statistical test
I obviously don't have a great understanding of what a GAM does different than a GLM. I believe it's a valid statistical test, (and I see an increase in the use GAMs, at least in ecological journals), but I need to know better when it's use is indicated over other regression analyses.
regression generalized-linear-model gam
regression generalized-linear-model gam
New contributor
New contributor
edited yesterday
New contributor
asked yesterday
lueromat
386
386
New contributor
New contributor
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
2
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday
add a comment |
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
2
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
2
2
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday
add a comment |
2 Answers
2
active
oldest
votes
up vote
9
down vote
accepted
The main difference imho is that while "classical" forms of linear, or generalized linear, models assume a fixed linear or some other parametric form of the relationship between the dependent variable and the covariates, GAM do not assume a priori any specific form of this relationship, and can be used to reveal and estimate non-linear effects of the covariate on the dependent variable.
More in detail, while in (generalized) linear models the linear predictor is a weighted sum of the $n$ covariates, $sum_{i=1}^n beta_i x_i$, in GAMs this term is replaced by a sum of smooth function, e.g. $sum_{i=1}^n sum_{j=1}^q beta_i , s_j left( x_i right)$, where the $s_1(cdot),dots,s_q(cdot)$ are smooth basis functions (e.g. cubic splines) and $q$ is the basis dimension. By combining the basis functions GAMs can represent a large number of functional relationship (to do so they rely on the assumption that the true relationship is likely to be smooth, rather than wiggly). They are essentially an extension of GLMs, however they are designed in a way that makes them particularly useful for uncovering nonlinear effects of numerical covariates, and for doing so in an "automatic" fashion (from Hastie and Tibshirani original article, they have 'the advantage of being completely automatic, i.e. no "detective" work is needed on the part of the statistician').
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation inmgcv
does a lot of things you cannot do withglm
, but could have been done in that framework also ...
– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. withmgcv
)?
– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
|
show 1 more comment
up vote
7
down vote
I'd emphasize that GAMs are much more flexible than GLMs, and hence need more care in their use. With greater power comes greater responsibility.
You mention their use in ecology, which I have also noticed. I was in Costa Rica and saw some kind of study in a rainforest where some grad students had thrown some data into a GAM and accepted its crazy-complex smoothers because the software said so. It was pretty depressing, except for the humorous/admirable fact that they rigorously included a footnote that documented the fact that they'd used a GAM and the high-order smoothers that resulted.
You don't have to understand exactly how GAMs work to use them, but you really need to think about your data, the problem at hand, your software's automated selection of parameters like smoother orders, your choices (what smoothers you specify, interactions, if a smoother is justified, etc), and the plausibility of your results.
Do lots of plots and look at your smoothing curves. Do they go crazy in areas with little data? What happens when you specify a low-order smoother or remove smoothing entirely? Is a degree 7 smoother realistic for that variable, is it overfitting despite assurances that it's cross-validating its choices? Do you have enough data? Is it high-quality or noisy?
I like GAMS and think they're under-appreciated for data exploration. They're just super-flexible and if you allow yourself to science without rigor, they will take you farther into the statistical wilderness than simpler models like GLMs.
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how wellmgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!
– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
add a comment |
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
up vote
9
down vote
accepted
The main difference imho is that while "classical" forms of linear, or generalized linear, models assume a fixed linear or some other parametric form of the relationship between the dependent variable and the covariates, GAM do not assume a priori any specific form of this relationship, and can be used to reveal and estimate non-linear effects of the covariate on the dependent variable.
More in detail, while in (generalized) linear models the linear predictor is a weighted sum of the $n$ covariates, $sum_{i=1}^n beta_i x_i$, in GAMs this term is replaced by a sum of smooth function, e.g. $sum_{i=1}^n sum_{j=1}^q beta_i , s_j left( x_i right)$, where the $s_1(cdot),dots,s_q(cdot)$ are smooth basis functions (e.g. cubic splines) and $q$ is the basis dimension. By combining the basis functions GAMs can represent a large number of functional relationship (to do so they rely on the assumption that the true relationship is likely to be smooth, rather than wiggly). They are essentially an extension of GLMs, however they are designed in a way that makes them particularly useful for uncovering nonlinear effects of numerical covariates, and for doing so in an "automatic" fashion (from Hastie and Tibshirani original article, they have 'the advantage of being completely automatic, i.e. no "detective" work is needed on the part of the statistician').
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation inmgcv
does a lot of things you cannot do withglm
, but could have been done in that framework also ...
– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. withmgcv
)?
– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
|
show 1 more comment
up vote
9
down vote
accepted
The main difference imho is that while "classical" forms of linear, or generalized linear, models assume a fixed linear or some other parametric form of the relationship between the dependent variable and the covariates, GAM do not assume a priori any specific form of this relationship, and can be used to reveal and estimate non-linear effects of the covariate on the dependent variable.
More in detail, while in (generalized) linear models the linear predictor is a weighted sum of the $n$ covariates, $sum_{i=1}^n beta_i x_i$, in GAMs this term is replaced by a sum of smooth function, e.g. $sum_{i=1}^n sum_{j=1}^q beta_i , s_j left( x_i right)$, where the $s_1(cdot),dots,s_q(cdot)$ are smooth basis functions (e.g. cubic splines) and $q$ is the basis dimension. By combining the basis functions GAMs can represent a large number of functional relationship (to do so they rely on the assumption that the true relationship is likely to be smooth, rather than wiggly). They are essentially an extension of GLMs, however they are designed in a way that makes them particularly useful for uncovering nonlinear effects of numerical covariates, and for doing so in an "automatic" fashion (from Hastie and Tibshirani original article, they have 'the advantage of being completely automatic, i.e. no "detective" work is needed on the part of the statistician').
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation inmgcv
does a lot of things you cannot do withglm
, but could have been done in that framework also ...
– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. withmgcv
)?
– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
|
show 1 more comment
up vote
9
down vote
accepted
up vote
9
down vote
accepted
The main difference imho is that while "classical" forms of linear, or generalized linear, models assume a fixed linear or some other parametric form of the relationship between the dependent variable and the covariates, GAM do not assume a priori any specific form of this relationship, and can be used to reveal and estimate non-linear effects of the covariate on the dependent variable.
More in detail, while in (generalized) linear models the linear predictor is a weighted sum of the $n$ covariates, $sum_{i=1}^n beta_i x_i$, in GAMs this term is replaced by a sum of smooth function, e.g. $sum_{i=1}^n sum_{j=1}^q beta_i , s_j left( x_i right)$, where the $s_1(cdot),dots,s_q(cdot)$ are smooth basis functions (e.g. cubic splines) and $q$ is the basis dimension. By combining the basis functions GAMs can represent a large number of functional relationship (to do so they rely on the assumption that the true relationship is likely to be smooth, rather than wiggly). They are essentially an extension of GLMs, however they are designed in a way that makes them particularly useful for uncovering nonlinear effects of numerical covariates, and for doing so in an "automatic" fashion (from Hastie and Tibshirani original article, they have 'the advantage of being completely automatic, i.e. no "detective" work is needed on the part of the statistician').
The main difference imho is that while "classical" forms of linear, or generalized linear, models assume a fixed linear or some other parametric form of the relationship between the dependent variable and the covariates, GAM do not assume a priori any specific form of this relationship, and can be used to reveal and estimate non-linear effects of the covariate on the dependent variable.
More in detail, while in (generalized) linear models the linear predictor is a weighted sum of the $n$ covariates, $sum_{i=1}^n beta_i x_i$, in GAMs this term is replaced by a sum of smooth function, e.g. $sum_{i=1}^n sum_{j=1}^q beta_i , s_j left( x_i right)$, where the $s_1(cdot),dots,s_q(cdot)$ are smooth basis functions (e.g. cubic splines) and $q$ is the basis dimension. By combining the basis functions GAMs can represent a large number of functional relationship (to do so they rely on the assumption that the true relationship is likely to be smooth, rather than wiggly). They are essentially an extension of GLMs, however they are designed in a way that makes them particularly useful for uncovering nonlinear effects of numerical covariates, and for doing so in an "automatic" fashion (from Hastie and Tibshirani original article, they have 'the advantage of being completely automatic, i.e. no "detective" work is needed on the part of the statistician').
edited yesterday
answered yesterday
matteo
1,356513
1,356513
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation inmgcv
does a lot of things you cannot do withglm
, but could have been done in that framework also ...
– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. withmgcv
)?
– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
|
show 1 more comment
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation inmgcv
does a lot of things you cannot do withglm
, but could have been done in that framework also ...
– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. withmgcv
)?
– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
2
2
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation in
mgcv
does a lot of things you cannot do with glm
, but could have been done in that framework also ...– kjetil b halvorsen
yesterday
Well, but as said in comments, all of that can be done with glm's also ... I suspect the main difference is pragmatic. The R implementation in
mgcv
does a lot of things you cannot do with glm
, but could have been done in that framework also ...– kjetil b halvorsen
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
Yes, I agree with you, GAMs are an extension of GLMs. However the question was about when to use GAM and when to use GLM, and it seemed to me that the op meant "classical" forms of GLMs, which do not usually include a set of basis function as predictors and are not used to reveal/approximate unknown nonlinear relationship.
– matteo
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
thanks - this is helpful. and yes, I was talking about classic GLMs
– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. with
mgcv
)?– lueromat
yesterday
@ matteo just two more things: i) what exactly do you mean by "true relationship is likely to be smooth, rather than wiggly"? and ii) "particularly useful for uncovering nonlinear effects of numerical covariates" - how would one describe / quantify nonlinearity (e.g. with
mgcv
)?– lueromat
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
The true relationship might not actually be smooth, however GAMs typically control for model complexity by adding a "wiggliness" penalty during the process of likelihood maximization (usually implemented as a proportion of the integrated square of the second derivative of the estimated function). Nonlinear effects of numerical covariates means that the influence of a particular numerical variable on the dependent variable might, for example, not increases/decreases monotonically with the variable value, but have an unknown shape, e.g. with local maxima, minima, inflection points,...
– matteo
yesterday
|
show 1 more comment
up vote
7
down vote
I'd emphasize that GAMs are much more flexible than GLMs, and hence need more care in their use. With greater power comes greater responsibility.
You mention their use in ecology, which I have also noticed. I was in Costa Rica and saw some kind of study in a rainforest where some grad students had thrown some data into a GAM and accepted its crazy-complex smoothers because the software said so. It was pretty depressing, except for the humorous/admirable fact that they rigorously included a footnote that documented the fact that they'd used a GAM and the high-order smoothers that resulted.
You don't have to understand exactly how GAMs work to use them, but you really need to think about your data, the problem at hand, your software's automated selection of parameters like smoother orders, your choices (what smoothers you specify, interactions, if a smoother is justified, etc), and the plausibility of your results.
Do lots of plots and look at your smoothing curves. Do they go crazy in areas with little data? What happens when you specify a low-order smoother or remove smoothing entirely? Is a degree 7 smoother realistic for that variable, is it overfitting despite assurances that it's cross-validating its choices? Do you have enough data? Is it high-quality or noisy?
I like GAMS and think they're under-appreciated for data exploration. They're just super-flexible and if you allow yourself to science without rigor, they will take you farther into the statistical wilderness than simpler models like GLMs.
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how wellmgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!
– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
add a comment |
up vote
7
down vote
I'd emphasize that GAMs are much more flexible than GLMs, and hence need more care in their use. With greater power comes greater responsibility.
You mention their use in ecology, which I have also noticed. I was in Costa Rica and saw some kind of study in a rainforest where some grad students had thrown some data into a GAM and accepted its crazy-complex smoothers because the software said so. It was pretty depressing, except for the humorous/admirable fact that they rigorously included a footnote that documented the fact that they'd used a GAM and the high-order smoothers that resulted.
You don't have to understand exactly how GAMs work to use them, but you really need to think about your data, the problem at hand, your software's automated selection of parameters like smoother orders, your choices (what smoothers you specify, interactions, if a smoother is justified, etc), and the plausibility of your results.
Do lots of plots and look at your smoothing curves. Do they go crazy in areas with little data? What happens when you specify a low-order smoother or remove smoothing entirely? Is a degree 7 smoother realistic for that variable, is it overfitting despite assurances that it's cross-validating its choices? Do you have enough data? Is it high-quality or noisy?
I like GAMS and think they're under-appreciated for data exploration. They're just super-flexible and if you allow yourself to science without rigor, they will take you farther into the statistical wilderness than simpler models like GLMs.
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how wellmgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!
– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
add a comment |
up vote
7
down vote
up vote
7
down vote
I'd emphasize that GAMs are much more flexible than GLMs, and hence need more care in their use. With greater power comes greater responsibility.
You mention their use in ecology, which I have also noticed. I was in Costa Rica and saw some kind of study in a rainforest where some grad students had thrown some data into a GAM and accepted its crazy-complex smoothers because the software said so. It was pretty depressing, except for the humorous/admirable fact that they rigorously included a footnote that documented the fact that they'd used a GAM and the high-order smoothers that resulted.
You don't have to understand exactly how GAMs work to use them, but you really need to think about your data, the problem at hand, your software's automated selection of parameters like smoother orders, your choices (what smoothers you specify, interactions, if a smoother is justified, etc), and the plausibility of your results.
Do lots of plots and look at your smoothing curves. Do they go crazy in areas with little data? What happens when you specify a low-order smoother or remove smoothing entirely? Is a degree 7 smoother realistic for that variable, is it overfitting despite assurances that it's cross-validating its choices? Do you have enough data? Is it high-quality or noisy?
I like GAMS and think they're under-appreciated for data exploration. They're just super-flexible and if you allow yourself to science without rigor, they will take you farther into the statistical wilderness than simpler models like GLMs.
I'd emphasize that GAMs are much more flexible than GLMs, and hence need more care in their use. With greater power comes greater responsibility.
You mention their use in ecology, which I have also noticed. I was in Costa Rica and saw some kind of study in a rainforest where some grad students had thrown some data into a GAM and accepted its crazy-complex smoothers because the software said so. It was pretty depressing, except for the humorous/admirable fact that they rigorously included a footnote that documented the fact that they'd used a GAM and the high-order smoothers that resulted.
You don't have to understand exactly how GAMs work to use them, but you really need to think about your data, the problem at hand, your software's automated selection of parameters like smoother orders, your choices (what smoothers you specify, interactions, if a smoother is justified, etc), and the plausibility of your results.
Do lots of plots and look at your smoothing curves. Do they go crazy in areas with little data? What happens when you specify a low-order smoother or remove smoothing entirely? Is a degree 7 smoother realistic for that variable, is it overfitting despite assurances that it's cross-validating its choices? Do you have enough data? Is it high-quality or noisy?
I like GAMS and think they're under-appreciated for data exploration. They're just super-flexible and if you allow yourself to science without rigor, they will take you farther into the statistical wilderness than simpler models like GLMs.
edited yesterday
answered yesterday
Wayne
15.8k13873
15.8k13873
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how wellmgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!
– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
add a comment |
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how wellmgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!
– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how well
mgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!– lueromat
yesterday
I imagine that I most often do what those grad students did: throw my data in a gam and be dazzled by how well
mgcv
handles my data. I try to be parsimonious with my parameters, and I check how well the predicted values match my data. your comments are a good reminder to be a bit more rigorous - and maybe finally get simon woods book!– lueromat
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
Heck, I'll go so far as to use a smoother to explore a variable, and then either fix the degrees of freedom at a low value or eliminate the smooth and use, say, a squared term if the smoother was basically quadratic. A quadratic makes sense for an age effect, for example.
– Wayne
yesterday
add a comment |
lueromat is a new contributor. Be nice, and check out our Code of Conduct.
lueromat is a new contributor. Be nice, and check out our Code of Conduct.
lueromat is a new contributor. Be nice, and check out our Code of Conduct.
lueromat is a new contributor. Be nice, and check out our Code of Conduct.
Thanks for contributing an answer to Cross Validated!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstats.stackexchange.com%2fquestions%2f380426%2fwhen-to-use-a-gam-vs-glm%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
GAM's are used when the linear predictor depends linearly on unknown smooth functions of some predictor variables.
– user2974951
yesterday
The distinction is blurry as you can represent numeric covariables e.g. by a spline also in a GLM.
– Michael M
yesterday
2
While the distinction is blurry, gam's can represent interactions also the smae way as glm's so strict additivity is not needed, the big difference is in inference: gam's need special methods, since estimation is not done via projection, but via smoothing. What that does imply in practice, I don' t understand.
– kjetil b halvorsen
yesterday
GLM $subset$ GAM.
– usεr11852
yesterday