[R-es] Paquete dismo, cálculo coeficiente de variación
Marcelino de la Cruz Rot
m@rcelino@del@cruz @ending from urjc@e@
Mie Jun 20 14:21:09 CEST 2018
Hola Jaume:
Lo más rápido para ver el código, al ser un método S4 es escribir:
findMethods(biovars)
La descripción de la función (con las definiciones y comentarios si los
tuviese) los tendrías que ver rebuscando en el código fuente del
paquete. Generalmente en la carpeta "R" del paquete con el código
fuente comprimido (dismo_1.1-4.tar.gz) que puedes encontrar en CRAN
(https://cran.r-project.org/src/contrib/dismo_1.1-4.tar.gz) o en este
caso, directamente desde su página de Github
(https://github.com/cran/dismo/blob/master/R/biovars.R).
Un saludo,
Marcelino
El 19/06/2018 a las 23:22, Jaume Tormo escribió:
> Hola Jorge y Marcelino,
>
> Muchas gracias a los dos.
> Para futuras dudas ¿Cómo puedo mirar el código y la definición de una
> función?
>
> Gracias de nuevo.
>
> Jaume.
>
> El 19 de junio de 2018, 13:16, Jorge Virto <jorge.virto using ehu.es
> <mailto:jorge.virto using ehu.es>> escribió:
>
> Hola,
>
> en la misma definición de la función:
>
> # P15. Precipitation Seasonality(Coefficient of Variation)
>
> # the "1 +" is to avoid strange CVs for areas where mean rainfaill
> is < 1)
>
> p[,15] <- apply(prec+1, 1, cv)
>
> Un saludo,
>
> Jorge
>
> On Martes, 19 de Junio de 2018 13:07:27 Marcelino de la Cruz Rot
> escribió:
>
> > Hola Jaume:
>
> >
>
> > Si miras el código de biovars() verás que la variable bio15 (el
>
> > coeficiente de variación de la precipitación) la obtiene sumando
>
> > previamente 1 a "prec":
>
> >
>
> > p[, 15] <- apply(prec + 1, 1, cv)
>
> >
>
> >
>
> > Es decir, en tu caso,
>
> >
>
> > > cv(prec+1)
>
> >
>
> > [1] 109.9637
>
> >
>
> > Lo que coincide con lo proporcionado por biovars, mientras que
>
> >
>
> > > cv(prec)
>
> >
>
> > [1] 112.5923
>
> >
>
> > coincide con tu cálculo manual.
>
> >
>
> > En la página web de WorldClim, donde definieron originalmente esta
>
> > variable no he encontrado una explicación de este sumatorio, pero
>
> > tampoco he buscado mucho...Tal vez en alguno de los artículos en
> los que
>
> > se describe esta base de datos den razón de este sumatorio.
>
> >
>
> > Saludos,
>
> >
>
> > Marcelino
>
> >
>
> > El 19/06/2018 a las 11:43, Jaume Tormo escribió:
>
> > > Estimados erreros,
>
> > >
>
> > > Estoy intentando entender como calcula el paquete dismo (
>
> > > https://cran.r-project.org/web/packages/dismo/index.html
> <https://cran.r-project.org/web/packages/dismo/index.html>) un
> coeficiente
>
> > > de
>
> > > variación. Os pongo un ejemplo:
>
> > > tmin <- c(10,12,14,16,18,20,22,21,19,17,15,12) # temperatura
> mínima media
>
> > > mensual de un año
>
> > > tmax <- tmin + 5 # temperatura máxima media mensual de un año
>
> > > prec <- c(0,2,10,30,80,160,80,20,40,60,20,0) #precipitación
> media mensual
>
> > > de un año
>
> > > biovars(prec, tmin, tmax) #este comando calcula una serie de
> variables
>
> > > relevantes para la distribución de especies.
>
> > >
>
> > > #El resultado es:
>
> > > bio1 bio2 bio3 bio4 bio5 bio6 bio7 bio8 bio9 bio10
>
> > >
>
> > > bio11
>
> > > [1,] 18.83333 5 29.41176 384.5501 27 10 17 22.5 13.83333 23.5
>
> > > 13.83333
>
> > >
>
> > > bio12 bio13 bio14 *bio15* bio16 bio17 bio18 bio19
>
> > >
>
> > > [1,] 502 160 0 *109.9637* 320 2 260 2
>
> > >
>
> > > La que a mi me interesa es bio15 el coeficiente de variación de la
>
> > >
>
> > > precipitación, pero no me da lo que yo esperaba que me diera:
>
> > >> (sd(prec)/mean(prec))*100
>
> > >
>
> > > [1] *112.5923*
>
> > >
>
> > > ¿Hay otra forma de calcular el Coeficiente de variación? Me estoy
>
> > > perdiendo
>
> > > algo. En el manual del paquete (
>
> > > https://cran.r-project.org/web/packages/dismo/dismo.pdf
> <https://cran.r-project.org/web/packages/dismo/dismo.pdf>), no
> explica como
>
> > > calcula el indice de variación ¿Hay alguna forma de rebuscar
> en dentro de
>
> > > biovars para saber que está haciendo?
>
> > >
>
> > > Muchas gracias.
>
> --
>
> --
>
> Jorge Virto Moreno
>
> Dpt. Economia Aplicada III
>
> Facultad de Economía y Empresa UPV-EHU
>
> Avda. Lehendakari Agirre, 83
>
> 48015 Bilbao
>
> Spain
>
> Tel: 94 601 3851
>
> Fax: 94 601 3754
>
> eman ta zabal zazu
>
> _ _
>
> | \___ |___\-\_
>
> | ___] __ | |
>
> | [_ __ [_ |_| |
>
> |__ _] [_ |___] /
>
> | [_ | ______/
>
> |___ \__| |
>
> | |
>
> |______|
>
> Universidad del País Vasco
>
> Euskal Herriko Unibertsitatea
>
>
>
>
> --
> Jaume Tormo.
> https://es.linkedin.com/in/jaumetormo
> https://acercad.wordpress.com/
--
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España
Más información sobre la lista de distribución R-help-es