Foro de discusión y soporte para usuarios de R

Un foro abierto a todos los temas relacionados con R que quiera.
Nota: Una url estable a este foro es https://knuth.uca.es/R-foro.

Análisis de datos extremos con R

de Francisco Serrano -

Buenas! Me presento, me llamo Francisco Serrano y estoy haciendo mi trabajo de fin de grado con análisis de datos extremos. Con ayuda de R y el paquete in2extRemes (interfaz del paquete extRemes) estoy tratando de comparar diferentes distribuciones de extremos, pero estoy muy verde en el software y no se por donde tirar en estos puntos:

- Estoy tratando de representar las distribuciones de Gumbel, Frechet y Weibull. Varias representaciones para cada distribución con diferentes variables en un mismo grafico. El problema está en que no se por donde tirar con Frechet y Weibull. Tengo el script con el que he logrado el de Gumbel por si es útil.

- El análisis de extremos con el paquete in2extRemes es muy cómodo para analizar Pareto, Gumbel o la GEV. Sin embargo, no puedo analizar Frechet o Weibull (tengo atravesadas estas dos distribuciones ) 

¿Alguien puede echarme una mano? Muchísimas gracias de antemano


PD: Soy nuevo en el foro y creo que no ha salido este tema. Disculpad si he metido la pata.

Modelo de Simulación Complejo en R

de Javier de la Cruz -

Hola buenas, lo primero de todo, soy nuevo en el foro así que me presento:

Me llamo Javier de la Cruz y ahora mismo me encuentro realizando mi Proyecto de Fin de Grado en Ingeniería Mecánica sobre la idoneidad de la ejecución de cálculos computacionales intensos en entornos HPC (High Peformance Computing).

Se trata de realizar simulaciones con software que requiera elevada potencia de cáculo (como software FEM o CFD) para comparar el rendimiento de diferentes máquinas.

Ahora mismo quiero hacer una simulación pesada utilizando el software estadístico R, y necesito que por favor alguien me envíe un archivo que pueda cargar en R y que tarde unas 4 horas aproximadamente en simularse, es decir, pesado. 

Si alguien me lo envía o me facilita un lugar donde descargarme modelos de este tipo para simular en R, se lo agradeceré enormemente.

Un saludo y muchas gracias por adelantado.


Paquete "Rfacebook"

de David Elías -

Buen día a todos.

Quisiera saber si hay alguien que haya usado efectivamente el paquete Rfacebook.

Me anoté como desarrollador, pero no logro realizar estadísticas.
Se conecta efectivamente R al sitio de facebook, pero no me otorga información para realizar estadísticas.

Un saludo.
David Elías

Error usando el paquete dsm

de Milagros Antun -

Hola, intento utilizar el paquete dsm,(library(Distance); library(dsm)) , siguiendo el Anexo de Miller (http://onlinelibrary.wiley.com/store/10.1111/2041-210X.12105/asset/supinfo/mee312105-sup-0001-AppendixS1.pdf?v=1&s=ced953b57365e5eb5753f0ad76dcc02c26918736 ). 

 
Trabajo con tres dataframes, cuyas str se detallan:

1) segdata:
 data.frame':193 obs. of  17 variables:
 $ Sample.Lab    : int  1 2 3 4 5 6 7 8 9 10 ...
 $ Transect.Label: Factor w/ 56 levels "1","100","101",..: 36 36 36 36 36 20 56 52 52 52 ...
 $ Effort        : int  1800 1800 1800 1800 1800 1800 1800 1800 1800 1800 ...
 $ x             : num  4443636 4437817 4442085 4440564 4439117 ...
 $ y             : num  5267395 5271579 5268309 5269266 5270337 ...
 $ ID_ESTRATO    : int  3 2 3 2 2 2 2 2 4 2 ...
 $ NDVI2010      : num  1813 1816 1804 1807 1816 ...
 $ NDVI2011      : num  2007 1943 1935 1894 1893 ...
 $ NDVI2012      : num  1705 1736 1686 1691 1729 ...
 $ NDVI2013      : num  2206 2305 2145 2211 2279 ...
 $ PROM_NDVI     : num  2218 2313 2148 2206 2275 ...
 $ DIST_PUEST    : num  959 455 2652 3194 1394 ...
 $ DIST_CUADR    : num  1482.1 137.5 549.9 62.9 514.8 ...
 $ DIST_MOLIN    : num  794 5022 2519 4156 5715 ...
 $ X_4326        : num  -63.7 -63.8 -63.7 -63.7 -63.7 ...
 $ Y_4326        : num  -42.7 -42.7 -42.7 -42.7 -42.7 ...
 $ O.KM2_2015    : num  64.1 34.6 43.4 44.4 46.6 ...
 
2) obsdata:
'data.frame':399 obs. of  6 variables:
 $ Especie.    : Factor w/ 1 level "Oveja": 1 1 1 1 1 1 1 1 1 1 ...
 $ size        : int  3 1 5 18 6 2 6 3 5 2 ...
 $ distance    : int  210 178 65 210 250 37 72 350 380 320 ...
 $ object      : int  1 2 5 7 8 13 14 20 30 31 ...
 $ Sample.Label: int  26 26 30 30 30 29 28 27 31 31 ...
 $ Effort      : num  1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 ... 
 
3)disdata
'data.frame':399 obs. of  7 variables:
 $ x       : num  4418278 4418667 4421229 4421308 4421308 ...
 $ y       : num  5299140 5298846 5295963 5295805 5295805 ...
 $ Especie.: Factor w/ 1 level "Oveja": 1 1 1 1 1 1 1 1 1 1 ...
 $ size    : int  3 1 5 18 6 2 6 3 5 2 ...
 $ distance: int  210 178 65 210 250 37 72 350 380 320 ...
 $ object  : int  1 2 5 7 8 13 14 20 30 31 ...
 $ Effort  : num  1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 ...
 
 
Luego de ajustar una función de detección con el script: 
 
hr.model <-ds(distdata,truncation="10%",transect ="line",dht.group= FALSE,key ="hr", convert.units = 1,adjustment =NULL) 

Intento ajustar mis datos a un modelo aplicando gams por medio del paquete dsm, corriendo el siguiente script:

mod1<-dsm(count ~ s(x,y, k=6), ddf.obj=hr.model, segdata, obsdata, engine = "gam",convert.units = 1, family = quasipoisson(link = "log"),group = FALSE,gamma = 1.4,control = list(keepData = TRUE),availability = 1,segment.area = NULL, weights = NULL)

Aquí es cuando tengo inconvenientes, ya que me sale el error:

Error in fix.by(by.x, x) : 'by' must specify a uniquely valid column

Alguien podría ayudarme a resolverlo?? Desde ya muchas gracias!!!

Problema con el paquete dsm en R

de Milagros Antun -

Hola, he intentado usar el paquete dsm, siguiendo los ejemplos de http://cran.r-project.org/web/packages/dsm/dsm.pdf y el tutorial de Miller http://onlinelibrary.wiley.com/store/10.1111/2041-210X.12105/asset/supinfo/mee312105-sup-0001-AppendixS1.pdf?v=1&s=ced953b57365e5eb5753f0ad76dcc02c26918736.

Mi problema surge cuando intetno correr el básico script :

mod1<-dsm(count ~ s(x,y, k=6), ddf.obj=hr.model, segdata, obsdata, engine = "gam",convert.units = 1, family = quasipoisson(link = "log"),)

Me aparece el siguiente warning:

Error in fix.by(by.x, x) : 'by' must specify a uniquely valid column

Alguien sabe que error puedo estar cometiendo?

Muchas Gracias

Liberada la versión 3.2.0 del paquete R-UCA

de Manuel Muñoz Márquez -

Buenas:

Quiero informarles de que el pasado día 10 de junio se liberó la versión 3.2.0 del paquete R-UCA basada en la versión 3.2.0 de R que incluye importantes novedades. Además incluye la versión 2.0-4 de RcmdrPlugin.UCA con opciones para realizar test de aleatoriedad para variables numéricas y factores de dos niveles.

Puede descargar la nueva versión del paquete desde https://knuth.uca.es/R/R-UCA

Si ha instalado el paquete le agradeceríamos contestase a la encuesta sobre el paquete R-UCA

Saludos

AYUDA CON PROGRAMA

de Hugo Alberto Ruiz Cancino -

Hola buenas tardes, mi problema es el siguiente, tengo este codigo

##Media de alfa
> limit=10000000
> for(i in 1:n)
+ {
+ if(alfa.m[i]>limit)alfa.m=alfa.m[-i]
+ }

Error en if (alfa.m[i] > limit) alfa.m = alfa.m[-i] :
valor ausente donde TRUE/FALSE es necesario

Y estoy usando ese codigo para eliminar unos valores de infinito que me salen, pero no se porque me sale el error en if. Alguien podria ayudarme. Gracias

tipos de apuntamiento en Rcmdr

de Bárbara Lopez -

Buenas tardes, me piden que calcule la curtosis de una variable en funcion de otra, aquí van las variables:

mujer    20
mujer    3
hombre 19
hombre 10
mujer    14
hombre 1
mujer   11
hombre  2
mujer    13
mujer    9
hombre 1
mujer    16
mujer    10
hombre  6
hombre  1

Rcmdr te permite calcular el apuntamiento, pero te pregunta de qué tipo es, si tipo 1, 2, y 3. ¿Cómo puedp saber qué tipo es el que tengo q usar?

Gracias!

Trabajo con R

de alberto del rio -

Hola buenas soy estudiante de economia en la complutense y tengo que hacer un trabajo con R, la verdad es que no tengo mucha idea y me gustaria saber, que no se si es aqui, si alguien me podria ayudar con algun comando, interpretacion etc..muchas gracias

Error en R

de Marta García Flores -

Al intentar usar la función "attach" con un archivo csv (tabla.csv) me da el siguiente error:

número mágico de archivo de restauración inválido (el archivo puede estar dañado)--ningún dato recargado

Además: Warning message:

file 'tabla.csv' has magic number '                               '

Use of save versions prior to 2 is deprecated.

 

¿A qué puede deberse el error y cómo lo soluciono? Gracias, es urgente