[R] Proportion data GAM

[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 GAM with beta distribution from mgcv package.
> My dependent variable is a proportion continuously ranging from 0 to 1
> (whales density) and I have three co-variates Depth, Distance to Coast and
> Seabed Slope.
>
>>From what I read, beta distribution is the most appropriate for my response
> variable and not binomial.
> According to mgcv manual, it is possible to specify beta distribution on a
> GAM with the "betar" function, but I get the following error:
> could not find function "betar"
>
> My code is:
> library(mgcv)
> a2=gam(Density~s(DEPTH,k=4)+s(DISTCOAST_1,k=4)+s(SLOPE,k=4),
> family=betar(link="logit"),data=misti,gamma=1.4)
>
> The beta family is specified exactly as it is shown in the manual:
> bm <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),family=betar(link="logit"),data=dat)
>
> Does anyone know what it seems to be the problem?
> Thanks in advance,
> Rodrigo
>
>
> Rodrigo Tardin
>
> Shor Term Scholar - Duke Marine Lab. - Duke University
> Doutorando em Ecologia e Evolu��o - IBRAG - UERJ
> M.Sc em Biologia Animal - PPGBA - UFRRJ
> Especialista em Doc�ncia do Ensino Superior - IAVM
> Laborat�rio de Bioac�stica e Ecologia de Cet�ceos - UFRRJ/ IF/ DCA
>
> 	[[alternative HTML version deleted]]
>
>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>


-- 
Simon Wood, Mathematical Science, University of Bath BA2 7AY UK
+44 (0)1225 386603               http://people.bath.ac.uk/sw283