[R] Proportion data GAM

2014-10-15 Thread Rodrigo Tardin
Hi all, I am not sure if this is the right place for this question or if there is one more specific. Anyway, I hope somebody can help me. I am trying to run a GAM with beta distribution from mgcv package. My dependent variable is a proportion continuously ranging from 0 to 1 (whales density) and

Re: [R] Proportion data GAM

2014-10-15 Thread Simon Wood
Can you give the result of typing sessionInfo() in the session where this happens, please? On 15/10/14 16:48, Rodrigo Tardin wrote: Hi all, I am not sure if this is the right place for this question or if there is one more specific. Anyway, I hope somebody can help me. I am trying to run a

Re: [R] Proportion data GAM

2014-10-15 Thread Simon Wood
Rodrigo, OK, it looks as if your mgcv help files/manual are somehow out of sync with the package version you have loaded. 'betar' is only available from mgcv 1.8. If you update to the current mgcv from CRAN then this problem should be solved. best, Simon ps. beta regression is only

Re: [R] Proportion data GAM

2014-10-15 Thread Rodrigo Tardin
Hi Simon, The result of sessionInfo() is: R version 3.1.0 (2014-04-10) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base

Re: [R] Proportion data GAM

2014-10-15 Thread Rodrigo Tardin
Thanks Simon! That worked! I did not constrain my k as you suggested, but when I saw my results, my degrees of freedom are not larger than 1, the REML is negative and all covariates are not significant (what it does not make sense). Is there something wrong? Here's the results of summary (a2)