Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.
up vote 4 down vote favorite

I have a plot of a GAM which models deforestation severity.

Component smooth functions of a GAM

The explanatory variables are transformed, with log10, inverse, and square root. Is it meaningful to interpret a gradient like this: 'Deforestation severity increases as log10(Distance to forest edge) increases' or 'Deforestation severity increases as distance to forest edge increases' for plot e (note, the lower the value of the y axis the greater the deforestation severity). If not, how could I plot the component smooth functions untransformed?

share|improve this question
1  
Do you mean plot e, not plot c, when you refer to "distance to forest edge increases"? – Peter Ellis Sep 2 '12 at 20:13
1  
Your question about "how could I plot the component smooth functions untransformed" is a question for your particular software. It looks a bit like R but could you confirm what you are using. There's no conceptual difficulty - you just reverse the transformation for the purpose of the plot. Arguably this is not a statistical question, but because it is related to interpretation however I think it is still on-topic for CV. – Peter Ellis Sep 2 '12 at 20:18
Yes, I did mean e, thanks. Yes, this modelling is done using R. The difficulty of reversing the transformations is that these plots are produced using the plot.gam function in the mgcv package. It uses the variables from the model to plot the smooth components. If I were to model the untransformed data and then plot that, would the smoothes be comparable to the transformed data? – gisol Sep 2 '12 at 21:25
At the moment, I have described things like this: "This trend is reflected in Figure 4.5e, where the log10 of distance to forest edge shows an initially flat relationship with deforestation severity, until it drops sharply, suggesting core forest away from the fringes of the forest patch is susceptible to severe deforestation." So I acknowledge the transformation, and mentally untransform the data to interpret the trends. – gisol Sep 2 '12 at 21:26
2  
Hi, sorry I don't have time to do a proper answer but no, if you modelled the untransformed data and then plotted that it would not give you the same results. You need to somehow extract the numeric results from plot.gam and then reverse transform them. – Peter Ellis Sep 2 '12 at 22:12
show 1 more comment

1 Answer

up vote 2 down vote

You can transform the graph by changing trans in plot.gam(..., trans = I, ...) to something else instead of I, for example the inverses of the functions you mentioned.

HTH.

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.