Re: [R-es] Interpretar intercepto de modelo ZINB

[Freddy Omar López Quintero]

Hola.

Esta es una pregunta interesante porque, en mi caso, me hace preguntarme
hace cuánto tiempo no interpreto realmente un intercepto, no solo este tipo
de modelos, sino en cualquier modelo regresivo. ¡Hace mucho tiempo le
presto solo atención atención los coeficientes de las variables!

Ahora, sin tener la respuesta cabal a tu pregunta, yo creo que debes apelar
al génesis del modelo. El modelo 'inflado' es una mezcla de distribuciones
donde la parte que es 'inflada' puede/debe tener una estructura distinta a
los otros componentes de la mezcla. Así, el 0 es explicado tanto por una
parte binaria y otra parte binomial negativa y los interceptos, desde mi
punto de vista, solo contribuyen un poco a la mezcla total dentro de su
propio componente. Dicho eso, yo apostaría por interpretar cada intercepto
como usualmente lo haces, sin olvidar la función de enlace y el componente
al que pertenece y comentando que son parte de un modelo de mezclas (ZINB,
en este caso). Es mi apuesta, no sé si haya una respuesta oficial.

He husmeado asimismo un poco los libros de J. Hilbe (QEPD) y no hace mucho
énfasis en los interceptos sino más bien en las otras variables
explicativas una vez el modelo ha sido estimado advirtiendo que la
interpretación es más compleja que con otro tipo de modelos (Hilbe, 2012,
pág. 376).

Ojalá algo sirva.

¡Salud!

2018-04-18 16:21 GMT-03:00 :

> Buenas tardes,
>
> ¿Cómo interpretarías el intercepto que da R en un modelo de ceros
> inflados? Por un lado en la parte de conteo tengo un intercepto de -4.2 y
> en la parte de ceros de 102, ambos salen significativos (***). ¿Qué me
> dirían?
>
> Gracias
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] r python interfaz

[Freddy Omar López Quintero]

La herramienta seguramente es buena; pero a mí particularmente no deja de
sorprenderme el cisma que desde RStudio se está provocando en el sentido
que ya casi ni usamos R sino R + X, donde X es cualquiera cosa que produzca
RStudio, aún cuando esto duplique/ignore/arrebate/pase por alto/opaque el
trabajo de mucha gente antes (incluyendo, sospecho yo, el mismo R base).
Así es la evolución de las cosas, me contesto.

2018-04-04 19:40 GMT-03:00 Javier Marcuzzi 
:

> Estimados
>
> Recuerdo que hace unos años en esta lista con uno de los Carlos se había
> creado un hilo sobre R y python, hoy leo una noticia, por lo menos para mí,
> que podría entrar en esa situación entre lo mejor de dos alternativas.
>
> Lógicamente, es nuevo para mí, no puedo opinar nada más que compartir la
> información, probar, usar, y cada uno puede utilizar o rechazar la
> herramienta.
>
> Miren http://blog.rstudio.com/2018/03/26/reticulate-r-interface-to-python/
>
> Javier Rubén Marcuzzi
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Experiencias de uso de algoritmo SMOTE y relacionados.-

[Freddy Omar López Quintero]

Hola gente,

¡Espero que estén muy bien!

Consultas: ¿Podrían compartir su experiencia (verbalmente, no hablo del
código) utilizando el algoritmo SMOTE o relacionados para atacar el poblema
de desbalance?¿Tuvieron éxito?¿Qué magnitud tenía el desbalance?

Esto es para corroborar/desmentir mi propia experiencia en este inciso de
la vida...

¡Gracias!

-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Goal seek

[Freddy Omar López Quintero]

Hola.

De acuerdo a lo que comentas, lo que dijo don Javier y una búsqueda de
internet, parece que sí, con un loop o un optimize se puede alcanzar lo que
describes. Aquí un enlace que parece atacar tu inquietud:

https://stackoverflow.com/questions/28575395/back-solving-a-function-or-goal-seek-in-r

No obstante, siempre un ejemplo específico sería mejor para tratar.

Salud.

2018-03-27 17:18 GMT-03:00 Juan Carlos Cardenas Oliveros <
juan.carde...@pucp.pe>:

>
> Hola Freddy, respecto a tu consulta.
>
> Por ejemplo, supongamos que necesita pedir prestado dinero. Usted sabe
> cuánto dinero quiere, cuánto tiempo quiere que pague el préstamo y cuánto
> puede pagar cada mes. Puede usar Goal Seek para determinar qué tasa de
> interés necesitará asegurar para cumplir con su objetivo de préstamo.
> Goal Seek solo funciona con un valor de entrada variable. Si quieres
> aceptar más de un valor de entrada; por ejemplo, tanto el monto del
> préstamo como el monto del pago mensual de un préstamo, utiliza el
> complemento Solver.
>
> Saludos
>
>
> El 27 de marzo de 2018, 14:42, Freddy Omar López Quintero <
> freddy.lopez.quint...@gmail.com> escribió:
>
>> Hola.
>>
>> ¿Y qué es *el* goal seek?
>>
>> ...
>>
>> 2018-03-27 16:39 GMT-03:00 Juan Carlos Cardenas Oliveros <
>> juan.carde...@pucp.pe>:
>>
>>> Hola
>>>
>>> es la primera vez que escribo por acá.
>>>
>>> Alguno de ustedes tiene el goal seek  hecho en R?
>>>
>>> Gracias
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ___
>>> R-help-es mailing list
>>> R-help-es@r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>>
>>
>>
>>
>> --
>> «...homines autem hominum causa esse generatos...»
>>
>> Cicero
>>
>>
>


-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Goal seek

[Freddy Omar López Quintero]

Hola.

¿Y qué es *el* goal seek?

...

2018-03-27 16:39 GMT-03:00 Juan Carlos Cardenas Oliveros <
juan.carde...@pucp.pe>:

> Hola
>
> es la primera vez que escribo por acá.
>
> Alguno de ustedes tiene el goal seek  hecho en R?
>
> Gracias
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] hacer un paquete en R

[Freddy Omar López Quintero]

Hola.

Algunas fuentes son:

   1. https://cran.r-project.org/doc/contrib/Leisch-CreatingPackages.pdf
   2.
   
https://support.rstudio.com/hc/en-us/articles/200486488-Developing-Packages-with-RStudio
   3. https://www.youtube.com/watch?v=9PyQlbAEujY
   4. https://hilaryparker.com/2014/04/29/writing-an-r-package-from-scratch/

ordenadas según el top de búsqueda para mi google.

Saludos.


2018-03-15 13:51 GMT-03:00 Javier Valdes Cantallopts (DGA) <
javier.val...@mop.gov.cl>:

> Hola a todos:
>
>
>
> Tengo un modelo que logré construir en R. Sin embargo, me gustaría que
> dicho modelo pusiese funcionar como un paquete de R, en el sentido que
> todas los pasos y fórmulas estén contenidas en un solo “paquete” y no como
> está desplegada ahora, en donde se observan casi sesenta líneas de código.
> Con el paquete sería más limpio el despliegue del modelo.
>
> Alguien ha construido un paquete en R, es muy complejo? Es lo que necesito?
>
> Saludos a todos.
>
>
>
> [image: Descripción: FIRMA3]
>
>
>
> --
>
> CONFIDENCIALIDAD: La información contenida en este mensaje y/o en los
> archivos adjuntos es de carácter confidencial o privilegiada y está
> destinada al uso exclusivo del emisor y/o de la persona o entidad a quien
> va dirigida. Si usted no es el destinatario, cualquier almacenamiento,
> divulgación, distribución o copia de esta información está estrictamente
> prohibido y sancionado por la ley. Si recibió este mensaje por error, por
> favor infórmenos inmediatamente respondiendo este mismo mensaje y borre
> todos los archivos adjuntos. Gracias.
>
> CONFIDENTIAL NOTE: The information transmitted in this message and/or
> attachments is confidential and/or privileged and is intented only for use
> of the person or entity to whom it is addressed. If you are not the
> intended recipient, any retention, dissemination, distribution or copy of
> this information is strictly prohibited and sanctioned by law. If you
> received this message in error, please reply us this same message and
> delete this message and all attachments. Thank you.
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
«...homines autem hominum causa esse generatos...»

Cicero
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Alternativa a rrefine

[Freddy Omar López Quintero]

¡Hola!

Espero que estén muy bien.

Consulta: ¿alguno tiene experiencia (y que la pueda compartir) utilizando
alguna alternativa a rrefine (
https://cran.r-project.org/web/packages/rrefine/vignettes/rrefine-vignette.html
)?

He utilizado rrefine y open refine (http://openrefine.org/) en general para
corregir conjuntos de datos chicos, pero ya cuando los datos son demasiados
para ocupar un csv, no me parece la mejor opción y por eso es mi inquietud.

¡Muchas gracias!

-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Filtrado de variables

[Freddy Omar López Quintero]

Es cierto, una vez más entendí todo mal. Mis disculpas.

La solución con SQL no va nada corta, me parece:

# filtrar el mes en el que todos los campos de la columna RESULTADO son NA

sqldf("
select *
from DF
where mes in
(
select a.mes from
(
select mes, count(*) as total_NA
from
(select *
from DF
where resultado=' NA ')
group by mes
) a,
(
select mes, count(*) as total_MES
from DF
group by mes
) b

where total_na=total_mes and a.mes=b.mes

)
")

¡Salud!

2018-02-22 9:01 GMT-03:00 Carlos Ortega :

> Nope..,
>
> Sí, puede ser "sqldf" pero de esta forma no te muestra los meses que
> tienen todos sus meses con "NA"..
> Tienes que hacer primero un conteo de meses y luego un "left join" con el
> conteo de "NAs" por mes y quedarte con el mes que tiene los dos valores
> iguales (número de meses y número de NAs)...
>
> Estaba probando la opción, pero no he terminado de completarla...
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El 22 de febrero de 2018, 12:53, Freddy Omar López Quintero <
> freddy.vat...@gmail.com> escribió:
>
>> El jue, 22-02-2018 a las 09:52 +0100, Carlos Ortega escribió:
>>
>> Aquí tienes un par de formas...
>>
>>
>> ... y aún una más, con el siempre fiel lenguaje SQL:
>>
>>
>> library(sqldf)
>>
>> # corriendo antes lo que hizo don Carlos pero dejándolo como un data.frame 
>> tradicional:
>>
>> DF<-as.data.frame(DF)
>>
>> sqldf("select * from DF where resultado=' NA '")
>>
>>
>> ¡Saludos!
>>
>>
>> --
>>
>> «...homines autem hominum causa esse generatos...»
>>
>> Cicero
>>
>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>



-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Filtrado de variables

[Freddy Omar López Quintero]

El jue, 22-02-2018 a las 09:52 +0100, Carlos Ortega escribió:
> Aquí tienes un par de formas...

... y aún una más, con el siempre fiel lenguaje SQL:

library(sqldf)# corriendo antes lo que hizo don Carlos pero dejándolo como un 
data.frame tradicional:DF<-as.data.frame(DF)sqldf("select * from DF where 
resultado=' NA '")
¡Saludos!

-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] SOBRESCRIBIR EN LA PRIMERA LÍNEA DE TEXTOS

[Freddy Omar López Quintero]

Hola.

El lun, 19-02-2018 a las 16:37 +0100, Carlos Córcoles escribió:
> No se si me podrías dar  desde esta lista alguna sugerencia...

Tomando un subconjunto del clásico conjunto de datos mtcars y
convirtiéndolo a json para ejemplificar, puedes hacer lo siguiente:

library(jsonlite)
line<-toJSON(head(mtcars), pretty=T)write(line,file="map_wgs91.js")
line<-toJSON(head(mtcars), pretty=T)write(line,
file="map_wgs91.js")fConn <- file('map_wgs91.js', 'r+')Lines <-
readLines(fConn)writeLines(c("var statesData =\n", Lines), con =
fConn)close(fConn) 

que lo tomé, básicamente, de aquí:

http://r.789695.n4.nabble.com/Appending-strings-at-the-beginning-of-a-t
ext-file-td901370.html

Ojalá sea de utilidad.


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] help (Crear bucle para generar un dataframe a partir de datos .mdb)

[Freddy Omar López Quintero]

El jue, 15-02-2018 a las 16:18 -0300, Matias Parra escribió:
> el problema que veo me aparece, es que previo al sqlQuery, debo
> llamar la base de datos con  odbcConnectAccess ("nombre archivo"). si
> lo hago base por base me demoraría bastante.

Si RODBC no te funciona como esperas y tienes los .mdb, podrías
utilizar también la librería Hmisc, con algo como:

library(Hmisc)mdb.get(paste(año, ".mbd", sep=''), tables=paste("tabla",
s, "_", año, sep=''))
que quizá es un poco más limpio **pero** realmente no sé si vaya a
hacer el trabajo más rápido. Supongo que esto depende más de las bases
de datos de access y su tamaño que de la importación propiamente (¿?).

Ojalá sea de utilidad.

¡Salud!


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] help (Crear bucle para generar un dataframe a partir de datos .mdb)

[Freddy Omar López Quintero]

¡Hola!

El jue, 15-02-2018 a las 12:43 -0300, Matias Parra escribió:
> library(RODBC)
> setwd("C:/Users/M/Documents/R data/")
> base<-odbcConnectAccess("2014")
> datos1<-sqlQuery(base ,"SELECT * FROM TABLA1_2014")
> datos2<-sqlQuery(base ,"SELECT * FROM TABLA2_2014")
> datos3<-sqlQuery(base , "SELECT * FROM TABLA3_2014")
> base1<-Reduce(function(...) merge (...,by="ID", all=T),
> list(datos1,datos2,datos3))

Utilizando la función paste() y algún bucle por aquí y por allá podrías
reutilizar tu mismo código. 

Por ejemplo:

sqlQuery(base, paste(paste("select * from tabla", s, sep=''), "_", año,
sep=''))
variando s de 1 a 3 y año de 1995 a 2014, haría el trabajo de
seleccionar lo que necesitas para esas tablas y períodos. Luego
restaría que la información la acumules o guardes en una lista o donde 
sea más conveniente.

Ojalá sea de ayuda.

¡Salud!


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Ayuda

[Freddy Omar López Quintero]

Holap.

El mié, 14-02-2018 a las 13:41 +0100, Francisco Javier Ibáñez López
escribió:
> Tengo una base de datos con 28 sujetos, a los cuales se les mide un 
> montón de elementos (tengo 27 variables dependientes) en tres
> tiempos 
> distintos (un total de 84 mediciones por variable). Busco comparar
> unas 
> variables con otras.
> 
> ¿Alguien me ayuda con el modelo estadístico a seguir para poder
> analizarlos?

Además de los modelos individuales longitudinales que puedes ajustar
para cada variable dependiente, puedes explorar:

 * Modelos de medidas repetidas para datos multivariantes. Aquí se
evalúa cómo evolucionan p variables respuesta sobre n individuos en t
momentos diferentes. Son modelos intermedios entre MANOVAs y modelos de
curvas de crecimiento. La §6.9 del clásico Johnson & Wichern (2007)
está dedicado a esto y R cuenta con la función manova() en la base.
 * Parecido a lo anterior, pero cuando la respuesta no es normal
(continua) hay modelos generalizados multivariantes para alcanzar lo
mismo, aunque no he visto un modelo de estos manejando tantas variables
respuestas (27) como en tu caso. Hubo un paquete de nombre sabreR (http
://www.sabre.lancs.ac.uk/sabreRuse_intro.html) que servía para estos
menesteres. Como estoy desactualizado y estoy viendo que ya no existe
en CRAN, alguien más seguramente te dará el nombre de un paquete vivo
para hacer esto.

Ahora, un posible problema que puedas tener, es que tienes más o menos
la misma cantidad de variables que de individuos, y cualquier matriz de
covarianza que calcules, si la puedes calcular, será inestable y poco
convincente. Cuenta la leyenda que regularizando (modelos lasso y
relacionados) puedes sacarle provecho. Si lo haces, podrías
compartirnos luego cómo jeje.

Ojalá sirva de algo.

¡Salud!


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] ANOVA no da el mismo resultado

[Freddy Omar López Quintero]

Hola.
El mar, 13-02-2018 a las 12:32 +0100, Yesica Pallavicini Fernandez
escribió:
> 1) en anova2 quiero testar el modelo
> 
> Y~T,
> Y~A
> Y~B
> Y~TxA
> Y~AxB
> Y~AxBxT
> ¿Está bien planteada la función? anova2<-lm(d$Y~d$T*d$A*d$B)
> 
> 2) En la salida de ambas anovas, el resultado es diferente para la
> misma
> variable (T). Es mas, en la anova 1 el resultado es no significativo
> y en
> anova2 si lo es ¿A que se debe esto?
> 
> 3) Estos análisis los he repetido en SAS ¿ Por que puede dar
> diferente
> resultado en un programa y en otro?

Si tu objetivo es el anova, yo creo que es mejor utilizar la función
aov().

A continuación un ejemplo utilizando los datos npk:

> summary(aov(yield ~ N*P*K, npk))
Df Sum Sq Mean Sq F value Pr(>F)  
N1  189.3  189.28   6.161 0.0245 *
P18.48.40   0.273 0.6082  
K1   95.2   95.20   3.099 0.0975 .
N:P  1   21.3   21.28   0.693 0.4175  
N:K  1   33.1   33.14   1.078 0.3145  
P:K  10.50.48   0.016 0.9019  
N:P:K1   37.0   37.00   1.204 0.2887  
Residuals   16  491.6   30.72 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

que produce lo mismo que:

proc anova data=npk;
class n p k;
model yield=n|p|k ;
run;

SourceDFAnova SSMean SquareF ValuePr > 
FN1189.2816667189.28166676.160.0245P18.40166678.40166670.270.6082N*P121.281666721.28166670.690.4175K195.201666795.20166673.100.0975N*K133.13533.1351.080.3145P*K10.48166670.48166670.020.9019N*P*K137.001666737.00166671.200.2887

Ojalá sea de utilidad.

¡Salud!

-- 
«...homines autem hominum causa esse generatos...»
 
Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] RStudio, Spss, SAS, Stata, txt

[Freddy Omar López Quintero]

¡La magia del tidyverse!

2018-02-12 12:08 GMT-03:00 Javier Marcuzzi 
:

> Estimados
>
> Hoy instale RStudio, y por sorpresa o desconocimiento de mi parte, veo que
> hay la posibilidad de importar datos excel, spss, sas, stata, texto,
> recuerdo consultas en la lista al respecto, posiblemente mucho se
> soluciones con las herramientas propuestas por RStudio.
>
> Javier Rubén Marcuzzi
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Avro.-

[Freddy Omar López Quintero]

Hola.

Espero que estén bien.

Les cuento. Tengo acceso a unos datos .avro en un s3 de amazon pero estoy
desconcertado con la forma de leerlos desde sparkR (*no sparklyr*) sin
tener que descargar el archivo.

Como antecedente, descargando el archivo, he utilizado felizmente las
instrucciones:

Sys.setenv(SPARK_HOME="/opt/apache-spark/")

library(SparkR, lib.loc = c(file.path(Sys.getenv("SPARK_HOME"), "R",
"lib")))

sparkR.session(master = "local[*]", sparkConfig = list(spark.driver.memory
= "2g"), sparkPackages = c("com.databricks:spark-avro_2.11:4.0.0"

dd <- read.df(path="/home/vate01/Descargas/TR_TRN_1.avro", source =
'com.databricks.spark.avro')

Pero quisiera evitar el paso de la descarga, teniendo las credenciales de
la cuenta de amazon.

¡Gracias!

-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] error en función ggadjustedcurves, paquete survminer

[Freddy Omar López Quintero]

Hola.

Sin tener una muestra de los datos es un poco difícil saber bien la causa,
pero parece que la forma en la que tienes los datos no es coherete para
ggadjustedcurves. Veo en la documentación de esta función:

*data* a dataset for predictions. If not supplied then data will be
> extracted from the fit object.
>

lo que me hace sospechar que deberías probar sin pasarle este argumento a
la función y/o ver si siendo tipo data.frame funciona. Talvez tienes los
datos de otra forma (¿array, matrix, ..?).

Saludos.

2018-01-27 11:50 GMT-03:00 Patricio Suárez Gil :

> Tengo un modelo de regresión de Cox y quiero obtener el plot ajustado por
> una covariable (sexo) con la función ‘ggadjustedcurves’, pero me da el
> siguiente error:
>
> > cox2 <- coxph(os ~ imc_25 + sexo.1, data = datos)
> > cox2
> Call:
> coxph(formula = os ~ imc_25 + sexo.1, data = datos)
>
>   coef exp(coef) se(coef) z p
> imc_25  -0.621 0.5370.299 -2.08 0.038
> sexo.1M  0.714 2.0420.387  1.84 0.065
>
> Likelihood ratio test=6.23  on 2 df, p=0.0444
> n= 76, number of events= 54
> > ggadjustedcurves(cox2, variable = datos$sexo.1, data = datos)
> Error in `[.data.frame`(data, , variable) : undefined columns selected
>
> Agradezco cualquier feedback.
>
> Saludos,
>
>
> Patricio
>
>
> Patricio Suárez Gil
> Unidad de Investigación Área V-Gijón
> Planta 5ª Impar
> Hospital Universitario de Cabueñes
> C/Prado, 395
> 33394 Gijón (Asturias)
> Tfno: 985 185 000 (Ext. 85715)
> @uinvest_psg
> unidadinvestigacion.ar...@sespa.es
> ESPAÑA
>
>
>
>
>
>
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es




-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] ejecución de script si un archivo es colocado en un directorio

[Freddy Omar López Quintero]

Hola.
El mar, 23-01-2018 a las 19:44 -0500, patricio fuenmayor escribió:
> He pensado ejecutar el llamado del script por intervalos de tiempo y
> validar si el archivo se encuentra en el directorio y ejecutarlo

Aunque yo no lo he usado, creo que esta utilidad:

https://github.com/ar-/incron

describe tu situación. Aquí un tutorial:

https://www.howtoforge.com/tutorial/trigger-commands-on-file-or-directo
ry-changes-with-incron/

Saludos.

-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Sobre paquete slam.-

[Freddy Omar López Quintero]

¡A diferencia de la mayoría de las veces, encontré la respuesta pronto!

A los que les suceda lo mismo (y posiblemente para mí mismo en el futuro),
mirar las versiones de gcc y de gfortran:

[ec2-user@ip ~]$ gcc --version
gcc (GCC) 7.2.1 20170915 (Red Hat 7.2.1-2)
Copyright (C) 2017 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

[ec2-user@ip ~]$ gfortran --version
GNU Fortran (GCC) 6.4.1 20170727 (Red Hat 6.4.1-1)
Copyright (C) 2017 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

Al cambiarlos a la misma (6.4.1), todo funcionó.

¡Gracias!

2018-01-23 14:55 GMT-03:00 Freddy Omar López Quintero <
freddy.lopez.quint...@gmail.com>:

> ¡Hola!
>
> Espero que estén bien.
>
> Estoy intentando instalar el paquete slam en la variante de Centos de
> amazon y obtengo el siguiente problema:
>
> gcc -m64 -std=gnu99 -shared -L/usr/lib64/R/lib -L/usr/local/lib64 -o
> slam.so apply.o dll.o grouped.o sparse.o util.o -L/usr/lib64/R/lib -lRblas
> -lgfortran -lm -lquadmath -L/usr/lib64/R/lib -lR
> /usr/bin/ld: cannot find -lgfortran
> /usr/bin/ld: cannot find -lquadmath
>
> ¿Alguien sabría por dónde estaría la solución a este incidente? Hace unos
> meses no tenía ningún problema (gcc y gcc-fortran están instalados).
>
> ¡Gracias de antemano!
>
> --
> «...homines autem hominum causa esse generatos...»
>
> Cicero
>
>


-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Sobre paquete slam.-

[Freddy Omar López Quintero]

¡Hola!

Espero que estén bien.

Estoy intentando instalar el paquete slam en la variante de Centos de
amazon y obtengo el siguiente problema:

gcc -m64 -std=gnu99 -shared -L/usr/lib64/R/lib -L/usr/local/lib64 -o
slam.so apply.o dll.o grouped.o sparse.o util.o -L/usr/lib64/R/lib -lRblas
-lgfortran -lm -lquadmath -L/usr/lib64/R/lib -lR
/usr/bin/ld: cannot find -lgfortran
/usr/bin/ld: cannot find -lquadmath

¿Alguien sabría por dónde estaría la solución a este incidente? Hace unos
meses no tenía ningún problema (gcc y gcc-fortran están instalados).

¡Gracias de antemano!

-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] STATA base de datos

[Freddy Omar López Quintero]

El mié, 06-12-2017 a las 12:55 +0100, Antonio Rodriguez Andres
escribió:
> me sale problema
> de memoria.

Pregunta posiblemente tonta: ¿tienes suficiente memoria para procesar
un archivo de tales dimensiones? Puede que ni aún cambiando la manera
de leer el archivo realmente lo puedas procesar.
> Entonces, una solucion es intentar pasar de STATA a CSV y luego usar
> el
> comando fread, y finalmente cargar los datos como RData

¿Has probado el paquete haven y específicamente su función read_dta?
Parece que es bastante más eficiente que la de foreign y preserva más
características que trae el formato de stata.

¡Ojalá algo sirva!

¡Salud!


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Búsqueda de palabras en una variable de R

[Freddy Omar López Quintero]

El mar, 28-11-2017 a las 03:42 +0100, miriam.alz...@unavarra.es
escribió:
> Tengo un vector de 40 palabras (marca) y necesito saber si en una de
> las
> variables del data.frame (datos) se incluye alguna de esas 40
> palabras. Si
> se incluye alguna de ellas, me gustaría crear una variable dummy
> siendo 1
> que incluye alguna palabra y 0 que no incluye.
> 
> ¿Qué paquete me recomendáis? ¿Cuál sería el comando a ejecutar?

Lo que describes luce como minería de texto y lo que parece que quieres
es una porción de la matriz que llaman Term-Document Matrix. El paquete
por excelencia para estos menesteres es tm:

https://cran.r-project.org/web/packages/tm/

que tiene su buena viñeta

https://cran.r-project.org/web/packages/tm/vignettes/tm.pdf

Ojalá sirva.

Saludos.


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Tibble o data.table?

[Freddy Omar López Quintero]

El mié, 22-11-2017 a las 20:51 +, Jesús Para Fernández escribió:
> Me parece que queda mas limpio, pero no se porque no le pillo la
> gracia

No estás solo en el mundo: si bien muchas cosas facilitan la vida
utilizando elementos de tidyverse, en su conjunto, yo no lo compro,
tampoco. Son como un dialecto en el idioma R.
¡Salud!


-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Es todo python

[Freddy Omar López Quintero]

> El 11-11-2017, a las 08:51, Jesús Para Fernández 
>  escribió:
> 
> Acabo de entrar en kaggle, hacia tiempo que no lo hacia, y veo con 'estupor' 
> que practicamente la gente trabaja con python.
> 
> Es cierto que la mayor parte de competiciones son usadas con redes 
> neuronales, donde python si tiene ventaja frente a R, pero me da miedo ver 
> como python para temas de machine learning parece estar ganando la partida...

¿Qué es «machine learning», qué es la estadística?¿Qué es python, qué es R?

Yo creo que R fue concebido como un software para cómputo estadístico, con un 
poderoso motor gráfico. Sería raro que en un concurso estadístico, R no 
resultara victorioso. Los estadísticos somos su audiencia por excelencia, pero 
ahora hablan de ciencia de datos, que me parece aún un término un poquín difuso.

python es un software de propósito general, popular entre computistas, 
informáticos, ingenieros y, como no, algunos estadísticos. No me alarma que, en 
grupos multidisciplinarios, un programa de propósito general resulte más útil 
que uno cuyo objetivo es particular.

En mi percepción, la mayoría de las personas que participan en kaggle y 
similares, no son estadísticos de formación (cuestión que no invalida nada, por 
si acaso) y por eso no han estado tan expuestos a R y utilizan otras 
herramientas de mayor visibilidad, como python. Quién sabe cuánta razón tendré, 
pero no quiero ni pensar en el porcentaje que utilizan otros excelentes 
programas estadísticos, como stata o el mismísimo SAS, en este tipo de 
escenarios.

En el mundo real, alejados de competiciones y ejercicios académicos, python 
funciona mejor para aglutinar los resultados que se van obteniendo en las 
distintas áreas. Esa ha sido mi experiencia.

¡Salud!
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Filtrar datos con una excepción

[Freddy Omar López Quintero]

El lun, 23-10-2017 a las 03:54 +0200, miriam.alz...@unavarra.es
escribió:
> Gracias...sí pero creo que no me has entendido...El código tampoco es
> el
> que necesito. Necesito todas las observaciones de Datos, excepto
> aquellas
> para las que se cumpla conjuntamente que evollucionsi=0 y
> evolucionno>0.

Ok, ha quedado claro que yo no he entendido, pero esa misma lógica que
mencionas (que debe ser muy fácil) es la que debes incluir en tu
filtro, sea cual fuere ¿podrías intentarlo, tal vez?. Reemplaza los
iguales, desiguales, mayores o menores en él y se acabó.



-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Filtrar datos con una excepción

[Freddy Omar López Quintero]

> El 22-10-2017, a las 22:41, miriam.alz...@unavarra.es escribió:
> 
> Gracias por la respuesta. Quizá me he explicado mal.

Oh no, he leído mal. 

Solo cambia la condición de == a !=:

filter(Datos, evolucionsi!=0, evolucionno != 1)

No olvides leer las ayudas que ofrecen las mismas funciones. Por ejemplo acá: 
http://dplyr.tidyverse.org/reference/filter.html 
.

¡Salud!


«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Filtrar datos con una excepción

[Freddy Omar López Quintero]

El lun, 23-10-2017 a las 02:00 +0200, miriam.alz...@unavarra.es
escribió:
> Necesitaría tener todos los datos excepto cuando se de
> simultáneamente
> que: Datos$evolucionsi=0 y Datos$evolucionno=0.

Podrías usar, entre muchas opciones, la función filter de la librería
dplyr:

filter(Datos, evolucionsi==0, evolucionno == 1)

Saludos.

-- 
«...homines autem hominum causa esse generatos...»

Cicero
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Estudios de no inferioridad en R

[Freddy Omar López Quintero]

> El 18-10-2017, a las 11:15, José Martín Arévalo  escribió:
> 
> ver si tres tratamientos quirúrgicos son equivalentes en el tratamiento de 
> una patología

¿Te refieres a un análisis de varianza? 

En este caso, la referencia inmediata es la función aov. Debes definir tu 
variable respuesta (numérica) y tus tratamientos en un factor (vector con 
categorías) en un data.frame y usar la sintaxis aproximada:

aov(respuesta~tratamientos, data=datos)

La ayuda de la función (?aov) te dará algunos ejemplos.

Pero si no te estás refiriendo a un análisis de varianza, por favor háznoslo 
saber para ser más específicos. Una minúscula muestra de tus datos también 
sería prudente.

¡Salud!
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] ggplot ayuda con argumento

[Freddy Omar López Quintero]

Hola.

Plagiando parte de lo que está en:

>
>
> https://stackoverflow.com/questions/12018499/how-to-put-labels-over-geom-bar-for-each-bar-in-r-with-ggplot2
>

se puede usar algo como:

dat <- read.table(text = "sample Types Number
> sample1 A   3641
> sample2 A   3119
> sample1 B   15815
> sample2 B   12334
> sample1 C   2706
> sample2 C   3147", header=TRUE)
>
> ggplot(data=dat, aes(x=Types, y=Number, fill=sample)) +  geom_bar(position
> = "dodge", stat="identity") + geom_text(aes(label=Number),
> position=position_dodge(width=0.9), vjust=-0.25)
>

Esto fue hecho para una tabla que ya tenía agrupadas las cantidades. En
cualquier caso, las modificaciones que debes hacer parece que son remover
los elementos del dodge y acomodar el vjust.

Otro ejemplo hay en

https://stackoverflow.com/questions/32986651/display-frequency-and-count-on-a-bar-chart
>
​
Ojalá este plagio casi absoluto sirva de algo.

¡Salud!

2017-10-11 22:19 GMT-03:00 Kelly Johanna Monroy Malagon <
kjmonr...@unal.edu.co>:

> Buenas Noches
>
> Estoy trabajando con R, con ggplot
> y deseo insertar el valor en la parte superior de cada uno de los
> rectángulos del gráfico, pero no se que argumento debo implementar,,
> cualquier ayuda o concejo sera bien recibido
>
> Muchas Gracias
>
>
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...homines autem hominum causa esse generatos...»

Cicero

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Subtotales de matriz

[Freddy Omar López Quintero]

Hola.

Calcular los subtotales cada 50 filas se puede hacer de varias formas. Una
es esta:

ko<-read.csv('Tdist.csv',h=F)
> ko<-data.frame(aux=rep(1:ceiling(nrow(ko)/50),each=50,length=nrow(ko)),ko)
> aggregate(.~aux,data=ko,sum)
>

donde aux es una variable auxiliar que da cuenta de la cantidad de grupos
de 50 filas que hay. El último grupo es menor y tiene solo 9 filas.

Este código es un poco rudimentario. Casi seguro con data.table o dplyr (de
los que soy ajeno) se puede hacer algo más grato.

Ojalá sirva, sin embargo.

¡Salud!

2017-09-22 1:37 GMT-03:00 Manuel Máquez :

> Estimado Freddy:
> Te agradezco la rapidez de tu respuesta.
> Trato de obtener los subtotales de cada 50 filas de Tdist.csv, es por ello
> que adjunte dicha tabla; y el objeto de haber enviado el otro archivo, era
> poner un ejemplo de los resultados que se obtendrían, al efectuar los
> subtotales de la misma Tdist.csv, pues es parte de ésta.
> Estoy a tus órdenes para cualquier otra aclaración.
> Atentamente;
>
> *M​**ANOLO MÁRQUEZ P.*
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es




-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Subtotales de matriz

[Freddy Omar López Quintero]

Hola.

Yo creo que no se entiende bien lo que quieres hacer. ¿Los subtotales que
mencionas son cada 50 filas o cada cuánto exactamente? O, específicamente
¿de acuerdo a cuál variable se crean o se deben crear estos subtotales?

Salud.

2017-09-21 21:51 GMT-03:00 Manuel Máquez :

> Buenas noches Colegas:
> He estando buscando la manera de encontrar los subtotales, de una tabla
> que tiene algo más de 750 registros, para ello uso un bucle for, pero
> siempre me da el total, a pesar de que dentro del bucle está el obtener la
> suma parcial en un registro de otra matriz.
> He tratado con colSum, applay(ko, 2, sum) y da sólo el gran total; también
> con  mr <- rbind(ko, 2, sum) pero en este caso el sistema me indica "cannot
> coerce type 'builtin' to vector of type 'list'"
> Adjunto un archivo llamado, Ejem.csv en donde vacié el contenido de los 50
> primeros registros, el registro 52 tiene los totales que se deberían
> asentar en el primero de la matriz resultante; y en los registros 55 a 104
> están los 50 siguientes, para terminar con el 106 que contiene el total de
> esos 50 anteriores.
> Es de hacer notar que no me importaría si en lugar de salir lo que indica
> el registro 106, se tuviera la suma de 52 y 106.
> Si hubiera alguna duda, adjunto también Tdist.csv que es la totalidad de
> los datos.
> Ojalá que alguno de Uds. me pudiera orientar, ya, anticipo las gracias más
> cumplidas.
>
> *MANOLO MÁRQUEZ P.*
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Pseudo R2 en modelos ZINB

[Freddy Omar López Quintero]

Hola.

El lun, 18-09-2017 a las 23:53 +0200, miriam.alz...@unavarra.es
escribió:
> Sin embargo, me gustaria saber si existen otras medidas que digan por si
> solas si el modelo es bueno o no, estilo un pseudo-R2.

El reciente trabajo

Jacob Martin & Daniel B. Hall (2016): R² measures for zero-
inflatedregression models for count data with excess zeros, Journal
of Statistical Computation and Simulation, DOI:
10.1080/00949655.2016.1186166

ofrece una variante. Puedes mirar su ecuación (5) para la idea inicial
y la ecuación (7) para el caso de la binomial negativa. Solo necesitas
tener algunos valores de la función de verosimilitud y calcular su
cociente.

Ojalá sea de utilidad.
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Evaluar el valor p de cada factor en un modelo ZINP

[Freddy Omar López Quintero]

​Hola Daniel.

2017-09-03 18:18 GMT-03:00 Daniel Chan Espinoza 

> Me gustaría saber como evaluar la significancia de cada factor, utilizando
> la prueba likelihood ratio (LRT), basada en la distribución de Chi2 en un
> modelo con inflación de ceros .
>

​Mucha gente te ayudará casi con seguridad en este foro, pero es importante
que nos muestres qué has hecho para saber por dónde comenzar o saber dónde
estás estancado.

Por ejemplo ¿estás usando una librería? En caso afirmativo, seguramente te
arroje las verosimilitudes y esto te pueda servir más fácilmente para ir
probando qué predictoras te son útiles y cuáles no. Si NO estás usando una
librería, la lógica prevalece pero debes calcular las verosimilitudes de
forma artesanal y proceder como se sugiere.

Lo otro, es que, como siempre, es posible que cada uno de nosotros
interprete diferente, y nadie quiere eso.

¡Salud!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Chi cuadrado ajustado

[Freddy Omar López Quintero]

Hola.

El mié, 30-08-2017 a las 16:51 -0300, Sebastian Gadea escribió:
> A la que me refiero es a la que se utiliza para probar la independencia de
> dos variables entre sí, mediante la presentación de los datos en tablas de
> contingencia. Pero no se si se puede realiza esa prueba Chi-cuadrado
> ajustando por otra variable.

Yo creo que te refieres a una prueba de independencia condicional (es
decir, ver si un par -al menos- de variables son independientes
condicionadas por una tercera). Una manera de hacerlo es con la prueba
de Cochran–Mantel–Haenszel, en en principio sigue una hipergeométrica
pero gracias a la magia que nos regalan las aproximaciones asintóticas,
podemos decir que se distribuye chi-cuadrado (que es lo que estás
preguntando exactamente).

Si he interpretado correctamente, entonces puedes echarle un ojo a la
función mantelhaen.test, que está en el paquete base de R.

Ojalá sirva de algo.

¡Salud!
[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Calcular probabilidad binomial negativa

[Freddy Omar López Quintero]

Hola.

No he revisado las expresiones de binomial negativa (lo mejor que
puedes/debes hacer, es escribir con claridad el código donde obtienes el
error para que no estemos adivinando) pero entiendo que te refieres al
combinatorio de 756 en 117 que sería con choose(755,117).

combn se usa para listar las posibles combinaciones de los elementos (que
con semejante cantidad no los podrás ver en ninguna parte).

¡Salud!

2017-08-29 3:02 GMT-03:00 Manuel Máquez :

> Estoy tratando de calcular la probabilidad de referencia con los siguientes
> datos:
> Eventos: 756
> Número de éxitos: 118
> Probabilidad de éxito: 118/756
> Probabilidad de fracaso: 319/378
> Hasta aquí llevo:  6.240113e -143
> El problema aparece cuando pongo el número de combinaciones es decir
> 755C117 porque RStudio me dice:
> Error in matrix(r, nrow = len.r, ncol = count) :
>   valor de 'ncol' no válido (demasiado grande o NA)
> Además: Warning message:
> In combn(756, 117) : NAs introduced by coercion to integer range
> La respuesta que obtuve por Libre Office es: 6.23540E -03
> Me podrían hacer favor de ayudarme indicando el procedimiento correcto.
> Anticipo las gracias más cumplidas.
> *MANOLO MÁRQUEZ P.*
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema con shiny en Fedora

[Freddy Omar López Quintero]

​Hola Guillermo.
​

"Reinstalling packages with the new version of R" dice "The following code
> should be run as root after upgrading R on Linux". Esto podría ser una
> solución pero ciertamente no sé cómo hacerlo porque en la terminal:
> $ sudo update.packages(ask = FALSE, checkBuilt = TRUE)
> no se puede ejecutar.
>

​¿Servirá escribir desde la terminal

sudo R​
>

​ingresar la contraseña y luego actualizar con

update.packages(ask = FALSE, checkBuilt = TRUE)​


?

​Haznos saber.

¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] impute TS y resultados en consola.

[Freddy Omar López Quintero]

Hola.

Yo no he usado el paquete antes y es altamente probable que esté
entendiendo mal la pregunta pero intentaré responder.

Me parece que el argumento principal de las funciones del paquete son un
vector o una serie de tiempo. Dependiendo de cuál función te estés
refiriendo pudiera haber un comportamiento distinto (quizá habrías de
aclararlo), pero con la mayoría que probé creo que solo reemplazando el
nombre con el resultado de la función misma es suficiente.

Usando el ejemplo de na.mean():

x <- ts(c(2,3,4,5,6,NA,7,8))
>

no sé si baste con:

x <- na.mean(x)
>

y así para las demás funciones. Si lo que quieres es ver específicamente
cuáles fueron los valores imputados, con

na.mean(x)[is.na(x)]
>

los tendrías (antes de hacer la imputación, claro está).

Ojalá sirva de algo.

¡Salud!

2017-07-21 23:07 GMT+02:00 Javier Valdes Cantallopts (DGA) <
javier.val...@mop.gov.cl>:

> Hola estimados.
>
>
>
> Alguien ha utilizado alguna vez el paquete IMPUTE_TS?
>
>
>
> La verdad es que es bastante recomendable, sin embargo hay un tema que no
> he podido resolver.
>
>
>
> El programa te entrega los resultados *en la consola*…sin embargo, la
> idea es rellenar los (*N.A) del dataset ORIGINAL*.
>
>
>
> La pregunta es: Cómo rellenar los (*N.A)* *del dataset* *original* si
> tengo los resultados en la consola? O como llevar los *resultados de la
> consola a un formato que me permita hacerlo*?
>
> Me explico bien?
>
>
>
> Saludos a todos.
>
>
>
> [image: Descripción: FIRMA3]
>
>
>
> --
>
> CONFIDENCIALIDAD: La información contenida en este mensaje y/o en los
> archivos adjuntos es de carácter confidencial o privilegiada y está
> destinada al uso exclusivo del emisor y/o de la persona o entidad a quien
> va dirigida. Si usted no es el destinatario, cualquier almacenamiento,
> divulgación, distribución o copia de esta información está estrictamente
> prohibido y sancionado por la ley. Si recibió este mensaje por error, por
> favor infórmenos inmediatamente respondiendo este mismo mensaje y borre
> todos los archivos adjuntos. Gracias.
>
> CONFIDENTIAL NOTE: The information transmitted in this message and/or
> attachments is confidential and/or privileged and is intented only for use
> of the person or entity to whom it is addressed. If you are not the
> intended recipient, any retention, dissemination, distribution or copy of
> this information is strictly prohibited and sanctioned by law. If you
> received this message in error, please reply us this same message and
> delete this message and all attachments. Thank you.
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Leer archivos con read.csv

[Freddy Omar López Quintero]

> La idea de Javier Nieto me parece buena, try cath es muy útil, a esto le
> sumaría print, como para ver en pantalla el nombre del archivo y un Si o
> No, mientras se cargan los datos, aunque esto es solo gusto personal.
>

​De eso se trata ¿no?: que la persona que consulta obtenga una o varias
respuestas con posibles soluciones para que pueda experimentar con ellas y
adecuarlas a su situación.​


Saludos.

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Leer archivos con read.csv

[Freddy Omar López Quintero]

Hola.

Una manera puede ser tomar una lista de todos los archivos de un directorio
y trabajar con ellos, puesto que ciertamente existen. ¿Cómo? Supongamos que
tienes una carpeta llena de archivos csv. Entonces, puedes guardar el
nombre de tus archivos en un objeto con

lista_archivos<-dir()
>

y esto te dará todos los csv de esa carpeta (si tienes otros archivos
también estarán allí, pero insistamos por ahora en suponer que son solo
csvs...).

Luego puedes hacer algo como

read.csv(lista_archivos[s])
>

dentro del bucle que ya estabas utilizando, para s entre 1 y el tamaño de
lista_archivos.

Finalmente, si en tu carpeta no hay solo csvs sino otros archivos, puedes
usar alguna variante de

grepl('.csv', lista_archivos[s])
>

para quedarte solo con aquellos archivos que cumplan con ser csv, por
ejemplo. Esto ya depende mucho de la entropía que haya en tus carpetas jeje.

Ojala sirva.

¡Salud!

2017-07-13 13:18 GMT-04:00 Javier Nieto :

> Hola, no he visto tu código pero puedes utilizar la función try.
>
>
> Saludos
>
> 
> De: R-help-es  en nombre de WILMER
> CONTRERAS SEPULVEDA 
> Enviado: jueves, 13 de julio de 2017 12:00:59 p. m.
> Para: r-help-es@r-project.org
> Asunto: [R-es] Leer archivos con read.csv
>
> Buen día para todos.
>
> Tengo un pequeño problema en el quisiera que me colaboraran.
>
> Estos días estoy trabajando en la lectura de una gran cantidad de archivos
> de extención .csv, se trata de una lista de documentos de 11 estaciones
> meteorologicas. Cada estación esta clasificada por meses y cada mes por
> dias, asi que tengo un archivo cada día. Por ejemplo tengo un archivo .csv
> para el día uno del mes de enero, otro archivo para el día dos del mes de
> enero, otro para el día tres del mes de enero y así sucesivamente hasta
> completar todos los meses.
>
> Mi problema es, que hay algunos archivos faltantes. Me refiero a que por
> ejemplo el día 5 del mes de febrero hace falta.
>
> Al tratar de leer los archivos FALTANTES con read.csv() me arroja el
> siguiente error:
>
> Error in file(file, "rt") : cannot open the connection
> In addition: Warning message:
> In file(file, "rt") :
>   cannot open file 'C:/R/Zonas rurales/Corponor/2017/7. PUNTO DE
> CONTROL/Febrero/27-02-2017.csv': No such file or directory
>
> Lo que es logico, debido a que no estan.
>
> Sin embargo al arrojarme este error se rompe el ciclo for y no me permite
> continuar trabajando.
>
> Mi pregunta es:
>
> *¿Existe alguna forma de indicarle al comando read.csv() que evite los
> archivos faltantes? *
>
> Muchas Gracias de antemano.
>
> *Wilmer Contreras Sepulveda.*
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> [[alternative HTML version deleted]]
>
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Favor ayuda

[Freddy Omar López Quintero]

2017-07-02 11:23 GMT-04:00 Gustavo Adolfo Redondo <
gustavo.redo...@ustadistancia.edu.co>:

> Por favor con que instrucción pudo reconfigurar esto, para que siempre
> aparezca en la parte inferior izquierda y no en una alterna ?
>

​Parece un problema de vieja data:

https://support.rstudio.com/hc/en-us/community/posts/207640217-Graphic-devices-are-opened-outside-RStudio
>

de esos que suelen resolverse reinstalando/actualizando Rstudio...

¿Sobre qué sistema estás y qué versión de R y RStudio manejas?

¡Salud!

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Ordenar vector

[Freddy Omar López Quintero]

2017-06-30 5:14 GMT-04:00 Isidro Hidalgo Arellano :

> x[order(x,  decreasing = TRUE)]


​También con la función sort:

sort(x, decreasing=T)
>

¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema en lectura de datos. Memoria insuficiente

[Freddy Omar López Quintero]

2017-06-27 13:12 GMT-04:00 Antonio Rodriguez Andres <
antoniorodriguezandre...@gmail.com>:

> Que archivo has intentado leer el del VWS o el de la link que te envié?


​Ambos. El del enlace que me enviaste es bastante más pequeño (me tendrías
que decir qué variables seleccionar para que sean exactamente iguales -o
enviárnoslo en un enlace drive o dropbox o algo así-). ​

¡
​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema en lectura de datos. Memoria insuficiente

[Freddy Omar López Quintero]

2017-06-27 11:42 GMT-04:00 Antonio Rodriguez Andres <
antoniorodriguezandre...@gmail.com>:

> Estaba usando el ESS pero me he dado cuenta que si tiene csv, yo lo bajé
> de aquí en formato SPSS pero me da el error de falta de memoria
>
>
> http://www.europeansocialsurvey.org/downloadwizard/
>

​De acuerdo. He usado la fuente que envías y otra que hallé en

http://www.worldvaluessurvey.org/WVSDocumentationWVL.jsp
>

​para descargar algunos datos y tuve inconvenientes para importarlos con la
función de foreign en su versión stata:

Error in `levels<-`(`*tmp*`, value = if (nl == nL) as.character(labels)
> else paste0(labels,  :
>   factor level [31] is duplicated
>

y con el formato de SPSS:

Error in read.spss() : error reading system-file header
>

pero no al usar haven, que los leyó aparentemente bien. Sí noté que con el
archivo .sav tarda bastante más en lanzar el error.

Naturalmente el archivo que es bastante grande (de casi 570MB) requerirá
mucho más tiempo para la importación (yo reinicié el ordenador para dejarle
el trabajo único de la importación).

Aventuro finalmente a decir que el problema es únicamente por la capacidad
del ordenador que utilizas y el esfuerzo que debe hacer la función para
convertir los datos. ¿Pudiste leer los datos en csv? (¡sería una 'solución'
feliz!)

Por otro lado, he leído que en Windows se puede usar la función
memory.limit() para modificar la cantidad de memoria para utilizar (
https://stackoverflow.com/questions/1395229/increasing-or-
decreasing-the-memory-available-to-r-processes).

¿Sirven de algo estas divagaciones? jeje.

¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema en lectura de datos. Memoria insuficiente

[Freddy Omar López Quintero]

2017-06-27 11:34 GMT-04:00 Antonio Rodriguez Andres <
antoniorodriguezandre...@gmail.com>:

> Si de hecho he utilizado el read_spss por el tema de las labels. Pero
> obtengo el mismo mensaje de problema de memoria.


​Vale. He buscado "WVS_Longitudinal_1981_2014_spss_v2015_04_18.sav"​ en
google pero mi búsqueda ha sido infructuosa. ¿Es posible que envíes un
enlace con los datos para probar? A ver si descubrimos alguna fisura por
dónde atacar.


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema en lectura de datos. Memoria insuficiente

[Freddy Omar López Quintero]

2017-06-27 11:10 GMT-04:00 Antonio Rodriguez Andres <
antoniorodriguezandre...@gmail.com>:

> El
> survey solo tiene los datos en formato SPSS o en formato de STATA (.dta).
>

​En el pasado yo he utilizado read_spss y read_dta ​del paquete haven y me
han funcionado mejor que las funciones de foreign, sobre todo por las
etiquetas que tanto SPSS y stata ocupan, que al leerse con las clásicas
funciones desaparecen. ¿Has probado estas funciones?

¡​Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problema con Histograma con porcentajes usando ggplot

[Freddy Omar López Quintero]

2017-06-18 23:28 GMT-04:00 Antonio Rodriguez Andres <
antoniorodriguezandre...@gmail.com>:

> Me puede recomendar algún libro donde poder empezar.


​De ggplot2, sin duda el libro de su (¿cismático?) creador es una
referencia obligada:

Wickham: ggplot2 Elegant Graphics for Data Analysis, 2016​
>

y uno muy bueno de

Chang: R Graphics Cookbook, 2012
>

¡
​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Incluir versión

[Freddy Omar López Quintero]

¡​Hola!​

2017-06-18 1:08 GMT-04:00 Manuel Máquez :

> Tengo la inquietud de poner el número de versión al usar write.csv de una
> tabla y no sé donde puedo encontrar la manera de hacerlo.
>

​Yo no entiendo bien. ¿El número de versión de R, de write.csv, de la
tabla..?​¿De qué exactamente?

¡​Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Regresión ponderada

[Freddy Omar López Quintero]

2017-06-14 15:35 GMT-04:00 Gerardo Gold :

> Para el caso de datos agrupados, donde ya se tiene la desviación estándar o
> varianza calculada ajustar un modelo ponderado es fácil, pero no se como
> hacerlo con datos sin agrupar.
>

​Si no entiendo mal, el núcleo del problema son las varianzas de la
relación. Una opción sería incluir entonces una estructura para la
varianza. Tomando arbitrariamente el conjunto de datos swiss como ejemplo,
podría ser algo como:

library(nlme)
> gls(Fertility~Education, data=swiss, weights=varIdent(form=~1|Education))​


​y/o con el argumento correlation de la misma función.​

​Quizás sirva.

¡Salud!​



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] errror al determinar puntos óptimos de corte (librería: OptimalCutpoints)

[Freddy Omar López Quintero]

2017-06-10 11:51 GMT-04:00 Fernando Sanchez via R-help-es <
r-help-es@r-project.org>:

> library(OptimalCutpoints)prediccion<-c(0.49165923,0.52759793,0.30213400,0.
> 33468349,0.14979703,0.47401846,0.52216404,0.42018794,0.92168073,0.
> 76893929,0.83362668,0.38251162,0.70803701,0.49165923,0.94462558)
> real<-c(0,1,0,0,0,0,1,1,1,1,1,0,1,0,1)datos_OPTIMO<-cbind(prediccion,real)
> cutpoint1 <- optimal.cutpoints(X = "prediccion", status =
> "real",tag.healthy = 1, methods = "Youden", data =
> datos_OPTIMO,categorical.cov =NULL, pop.prev = NULL,control =
> control.cutpoints(), ci.fit = TRUE)
>

​Creo que el detalle está en que tus datos son una matriz y no un
data.frame.​ Añadiendo:

datos_OPTIMO<-data.frame(datos_OPTIMO)
>

a tu código, obtengo:

​> cutpoint1
>
Call:
> optimal.cutpoints.default(X = "prediccion", status = "real",
> tag.healthy = 1, methods = "Youden", data = datos_OPTIMO,
> categorical.cov = NULL, pop.prev = NULL, control =
> control.cutpoints(),
> ci.fit = TRUE)
>
> Optimal cutoffs:
>   Youden
> 1 0.1498
>

Ignoro si es lo que se espera de respuesta.

¡Salud!

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Color en líneas (ggplot2)

[Freddy Omar López Quintero]

2017-06-08 12:54 GMT-04:00 Javier Marcuzzi 
:

> ¿me hice entender?


​No.

Para salir del escollo lo convertiré a gráficos​ base y continuaré con mi
vida 😂...

​Au revoir.​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Series de tiempo

[Freddy Omar López Quintero]

2017-06-09 16:06 GMT-04:00 WILMER CONTRERAS SEPULVEDA 
:

> Si, tengo un analisis en mente, se trata de la proyección a 15 años de la
> demanda de enrgia electrica para mi universidad.


​Es que allí está el asunto don Wilmer: hay muchas maneras de hacer una
proyección. La más sencilla sería usar alguna forma del tiempo como
predictora(s) y la función lm para regresar a su respuesta. Otra manera
sería hacer un suavizado exponencial y allí el paquete forecast es uno de
los que lleva la batuta. Si quisiera hacer un clásico arima, podría
utilizar la función arima de R base o, para un arma, la función​ arma de la
librería tseries. También podría estar interesado, dependiendo de sus
datos, en incorporar la variabilidad de la varianza en un modelo (g)arch,
que también está en tseries...

R cuenta con un view de series de tiempo aquí:

https://cran.r-project.org/web/views/TimeSeries.html
>

Allí encontrará más ideas de las que yo pueda listar.

¡Salud!

PS: he puesto en copia a la lista, por si acaso alguien más tiene esta
misma inquietud.


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Series de tiempo

[Freddy Omar López Quintero]

​Hola Wilmer,​

2017-06-09 9:34 GMT-04:00 WILMER CONTRERAS SEPULVEDA :

> quisiera saber que paquete me recomienda


​Antes, ¿tienes algún tipo de análisis en mente? R trae de por sí ya varias
utilidades (lm, ts, ar, arima, etc.) pero sin un objetivo claro es un poco
difícil sugerir​.

¡
​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Color en líneas (ggplot2)

[Freddy Omar López Quintero]

2017-06-08 9:41 GMT-04:00 Javier Marcuzzi :

> R Graphics Cookbook, y hay unos capítulos que entiendo que tratarían su
> problema


¡​Gracias don Javier!​

Efectivamente, he revisado el libro, especialmente la sección de nombre *Adding
Fitted Lines from an Existing Model*, pero no se muestra cómo añadir
leyendas a las líneas añadidas (¡sospechoso!).

ggplot espera que la mayoría de lo que se vaya a dibujar esté dentro del
mismo conjunto de datos para poder añadirle una columna identificadora al
mismo que sirva para generar la leyenda, pero acá se van añadiendo las
líneas una a una, de forma que no hay una variable por la cual agrupar... Y
si añado el color dentro del aes, entonces se generan dos leyendas
independientes, una con el tipo de línea, y otra con el color... :-(

¡Mil gracias!
​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Color en líneas (ggplot2)

[Freddy Omar López Quintero]

¡Gracias!

Pero al hacerlo R se cae y aparece un "R session aborted"... :-\

¿Habrá alguna manera manual?

2017-06-07 17:53 GMT-04:00 Carlos Ortega :

> Hola,
>
> Para este tipo de ajustes, vaya u otros sobre "ggplot2", este asistente
> ayuda mucho a modificar y adaptar la mayor parte de los elementos que
> forman un gráfico "ggplot"...
>
> https://github.com/calligross/ggthemeassist
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
>
> El 7 de junio de 2017, 23:17, Freddy Omar López Quintero <
> freddy.lopez.quint...@gmail.com> escribió:
>
>> ¡Hola!
>>
>> Espero que estén muy bien.
>>
>> Le he dado muchas vueltas a este asunto pero ggplot lleva ganada la
>> batalla.
>>
>> Les consulto. Tengo un conjunto de datos longitudinales y se hicieron
>> varios ajustes sobre ellos.
>>
>> Ahora, intento colocar las predicciones sobre el mismo gráfico con el
>> código adjunto. Mi problema exacto es que quisiera que la leyenda, además
>> del tipo de línea, pudiera mostrar el color de la línea. En este momento,
>> solo los individuos tienen líneas grises continuas, sin embargo la leyenda
>> solo muestra la línea continua, y negra.
>>
>> ¿Alguien podría ayudarme?
>>
>> ¡Gracias!
>>
>> --
>> «Pídeles sus títulos a los que te persiguen, pregúntales
>> cuándo nacieron, diles que te demuestren su existencia.»
>>
>> Rafael Cadenas
>>
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] OPINIONES comparacion personal de R en linux y windows

[Freddy Omar López Quintero]

2017-06-07 10:33 GMT-04:00 Javier Valdes Cantallopts (DGA) <
javier.val...@mop.gov.cl>:

> No sé qué opinan?


​Con R propiamente, yo creo que no encuentro mayor diferencia puesto que
uso RStudio y en general es la misma experiencia de usuario (con conjuntos
de datos bastante modestos), pero las diferencias me vienen cuando hay
herramientas o paquetes que hay que compilar o quiero trabajar más rápido
con la consola (tools>shell). Desde windows ignoro muchas cosas,
simplemente, y no les puedo sacar provecho. ¿Pero es esto un problema de R
o del usuario (yo)? :D

¡Saludos!​



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Color en líneas (ggplot2)

[Freddy Omar López Quintero]

¡Hola!

Espero que estén muy bien.

Le he dado muchas vueltas a este asunto pero ggplot lleva ganada la batalla.

Les consulto. Tengo un conjunto de datos longitudinales y se hicieron
varios ajustes sobre ellos.

Ahora, intento colocar las predicciones sobre el mismo gráfico con el
código adjunto. Mi problema exacto es que quisiera que la leyenda, además
del tipo de línea, pudiera mostrar el color de la línea. En este momento,
solo los individuos tienen líneas grises continuas, sin embargo la leyenda
solo muestra la línea continua, y negra.

¿Alguien podría ayudarme?

¡Gracias!

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas
parametros<-structure(list(X = 1:51, X.1 = c(93.6443045794884, 0.274514923404001, 
-0.618328826511992, 6.35004789645277, 80.4486680704678, 0.398671210849383, 
-0.36311742015548, 1.03735833192979, 123.433079354653, 0.0102818971461322, 
0.0737926758635394, 0.667956051190706, 0.166023046221436, 0.169835680441718, 
92.8727273392054, 0.280548440675895, -0.591788496784862, 6.34069803543858, 
0.00396918686734812, 0.003873777667181, 0.00397871442819185, 
8.62002235763456e-05, 0.000177982071249404, -7.20746737162845e-05, 
83.6179497099678, 0.37456999758485, -0.395741042334919, 1.35275234352414, 
9.19666589730974, 0.049465436765756, 0.262280809304406, 90.4949487031123, 
0.309798683863988, -0.455118957698306, 4.59410797144556, 6.70214183242127, 
2.55241318343321, 0.0040969271690877, 0.00406057805801237, 0.00408729052989319, 
-0.0001337437288817, -9.02974739023956e-05, -4.79605394640869e-05, 
84.1324991529296, 0.371924990227874, -0.394364833252805, 12.2942829860544, 
15.1762245140945, 9.39329446699581, 15.2777604868252, 1.12559789523233
), X2.5. = c(89.8075680676307, 0.24062495144807, -0.77721816380053, 
6.10328321127024, 82.0970864037386, 0.41628213216484, -0.316459851584057, 
1.10925276972414, 151.676569972767, 0.0136445644723838, 0.0964128077066584, 
0.566110617222703, 0.0719498743654596, 0.288082253061202, 88.7480021052549, 
0.248299724724991, -0.738083589646706, 6.13608196789711, 0.0033131765251182, 
0.00331992605912848, 0.00330510660938651, -0.000385787681769791, 
-0.000279920354928013, -0.000648710689208116, 82.1887422364474, 
0.361429662565472, -0.444308809159175, 1.26215061383911, 8.34267614185542, 
0.0373843251863912, 0.231315025808771, 88.4278512805936, 0.28043065708031, 
-0.576980216992324, 3.72759905121779, 2.51202639984708, 2.50514689364638, 
0.00352687109447984, 0.00343525744154969, 0.00338207649947237, 
-0.000660750333682392, -0.000468516920705751, -0.000454316005673436, 
82.6462425607854, 0.359068069451452, -0.441167465031998, 8.62796498220036, 
10.1227010680911, 5.12805200272888, 9.26381073259622, 1.01759087592246
), X97.5. = c(98.3877041960296, 0.309426923468285, -0.477640818649688, 
6.61776093280683, 83.7455047370094, 0.433893053480296, -0.269802283012633, 
1.18612987361049, 186.382629348511, 0.0181069832731281, 0.12596683046256, 
0.447134504827384, 0.0136480105954702, 0.411703462075963, 97.0714933511228, 
0.3231018418536, -0.440032765543431, 6.5910793206374, 0.00458902518252994, 
0.00459529042640684, 0.00459033307804952, 0.000433596593495416, 
0.000469433582168282, 0.000280674306287311, 84.9503686450073, 
0.388274472592534, -0.348087133287464, 1.45030802655994, 10.2002714211677, 
0.0609950356940555, 0.299171819309348, 93.548966629373, 0.333696770103566, 
-0.363375587398017, 5.0879499410823, 9.76329673573814, 2.7431077428306, 
0.00474466476593635, 0.0046929951241291, 0.00478980004233669, 
0.000449783679712311, 0.000413896440866447, 0.000443348858707374, 
85.4751963528418, 0.384064677811932, -0.348277751176412, 17.4449169639356, 
22.4575183191046, 16.3379448906165, 24.2774760534726, 1.23218199233448
)), .Names = c("X", "X.1", "X2.5.", "X97.5."), row.names = c(NA, 
-51L), class = "data.frame")

dd<-structure(list(x = c(1L, 2L, 3L, 4L, 5L, 3L, 4L, 5L, 6L, 1L, 
2L, 3L, 4L, 5L, 6L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 
5L, 6L, 1L, 4L, 5L, 1L, 2L, 3L, 5L, 3L, 4L, 5L, 4L, 1L, 2L, 3L, 
5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 6L, 1L, 2L, 3L, 4L, 
5L, 6L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 6L, 1L, 2L, 2L, 3L, 
4L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 
4L, 5L, 6L, 7L, 1L, 2L, 4L, 1L, 3L, 4L), y = c(28.032936, 47.78408, 
59.070448, 66.124428, 71.767612, 59.6954225, 66.54971, 73.4039975, 
78.8874275, 30.3693767567568, 47.0061689189189, 61.0834545945946, 
66.2024675675676, 71.3214805405405, 76.4404935135135, 31.292405333, 
44.719304, 55.460822933, 66.202341867, 74.258481067, 
79.629240533, 82.314620267, 32.530764556962, 44.8980734177215, 
57.265382278481, 65.9224984810127, 70.8694220253165, 75.8163455696203, 
29.5625964864865, 72.2403994594595, 77.4134664864865, 44.697832238806, 
60.9025955223881, 72.6878779104478, 80.0536794029851, 60.0298032876712, 
67.7313161643836, 74.1492435616438, 63.8927616901408, 41.41994, 
53.44556, 66.9743

Re: [R-es] Ingresar datos dentro de R sin Excel

[Freddy Omar López Quintero]

2017-04-28 10:35 GMT-03:00 WILMER CONTRERAS SEPULVEDA 
:

> Quisiera saber si existe alguna función o libreria en R que simule la
> ventana de Excel para ingresar datos. No estoy hablando de importar datos
> de excel a R, No.
> Quisiera que dentro de R existiera el entorno para ingresar los datos
> rapidamente y no pasar por la molestia de ingresar los datos a Excel y
> luego tener que importarlos a R.
>

​Hola Wilmer.​

​Creo que en general no se opera de la manera que describes. No debemos
pasar por excel para trabajar con los datos. ​Javier te sugiere utilizar
archivos csv que creo es bastante expedito.

Sin embargo, una manera nativa, aunque rudimentaria, de emular una especie
de hoja de cálculo es la función data.entry. Algo como:

> a_rellenar<-data.frame(matrix(ncol=2,nrow=10))
> > data.entry(a_rellenar)
>

Ojalá te sirva.


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Intervalos de confianza de la varianza de los residuos enunmodelo no lineal.-

[Freddy Omar López Quintero]

2017-04-01 16:06 GMT-03:00 :

> busque un documento donde hay algo de variancias y una librería de R,


¡​Gracias Javier!​¡Muy buen documento! Lo tendré bajo la manga para cuando
deba enfrentarme a un modelo jerárquico con varianza heterogénea.

También agradezco a Rubén por el tip sobre el uso de bootstrap para
recuperar el parámetro.

Por ahora me contentaré entonces con usar teoría asintótica clásica
(Galland, 1987; Seber & Wild, 2003) y aproximaré el intervalo con:

resumen<-summary(fm1DNase1)
> c(resumen$df[2]*resumen$sigma^2/qchisq(.975,resumen$df[2],lower.tail=T),
> resumen$df[2]*resumen$sigma^2/qchisq(.975,resumen$df[2],lower.tail=F))
>

puesto que sigma2 sigue una chi cuadrado con n-p grados de libertad.

¡​Gracias!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Intervalos de confianza de la varianza de los residuos en unmodelo no lineal.-

[Freddy Omar López Quintero]

​​Gracias Javier,​

2017-04-01 12:07 GMT-03:00 :

> Mi duda es la siguiente, la varianza residual en su modelo, es homogénea o
> heterogénea,


​La varianza es homogénea, común a todas las observaciones.

Digamos que el modelo es el siguiente,

y=f(x, betas)+e

con f alguna función no lineal cuyos parámetros son betas y con e~N(0,
sigma2). Uno suele estar centrado en obtener información sobre los betas y
marginar un poco la inferencia sobre el parámetro sigma2 (incluso en la
regresión lineal). Es sobre este que me pregunto si existirá alguna función
que devuelva su intervalo de confianza.

Ojalá haya respondido la inquietud.

¡Gracias!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Intervalos de confianza de la varianza de los residuos en un modelo no lineal.-

[Freddy Omar López Quintero]

Hola Rubén y Carlos,

¡Muchas gracias por sus respuestas!

Me parece sin embargo que ellas van encaminadas a obtener intervalos de
confianza para la predicción.

En la salida del ejemplo de la ayuda, se obtiene:

Residual standard error: 0.01919 on 13 degrees of freedom
>

Ese sería (si no he entendido mal 😌) el estimado de la desviación de los
errores del modelo y es para este (o su cuadrado) que quisiera un intervalo
de confianza.

¡Muchas gracias!


2017-04-01 8:48 GMT-03:00 Carlos Ortega :

> Hola,
>
> Esto te puede ayudar:
>
> http://stackoverflow.com/questions/32459480/r-confidence-bands-for-
> exponential-model-nls-in-basic-graphics
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El 1 de abril de 2017, 4:21, Freddy Omar López Quintero <
> freddy.lopez.quint...@gmail.com> escribió:
>
>> Hola amigos,
>>
>> Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en
>> el ejemplo de la ayuda:
>>
>> DNase1 <- subset(DNase, Run == 1)
>> fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)
>> summary(fm1DNase1)
>>
>> Aquí se está modelando la densidad óptica de un ensayo relacionada de
>> forma
>> no lineal (logística) con (el logaritmo) de la concentración de una
>> proteína.
>>
>> ¿Sabe alguien si existe una función que sirva para calcular un intervalo
>> de
>> confianza para la varianza de los errores?
>>
>> ¡Muchas gracias!
>>
>>
>> --
>> «Pídeles sus títulos a los que te persiguen, pregúntales
>> cuándo nacieron, diles que te demuestren su existencia.»
>>
>> Rafael Cadenas
>>
>> [[alternative HTML version deleted]]
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Intervalos de confianza de la varianza de los residuos en un modelo no lineal.-

[Freddy Omar López Quintero]

Hola amigos,

Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en
el ejemplo de la ayuda:

DNase1 <- subset(DNase, Run == 1)
fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)
summary(fm1DNase1)

Aquí se está modelando la densidad óptica de un ensayo relacionada de forma
no lineal (logística) con (el logaritmo) de la concentración de una
proteína.

¿Sabe alguien si existe una función que sirva para calcular un intervalo de
confianza para la varianza de los errores?

¡Muchas gracias!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] names

[Freddy Omar López Quintero]

​Hola.​

2017-03-29 7:37 GMT-03:00 Mauricio Monsalvo :

> Tengo un set de datos importados desde SPSS (que usa etiquetas).


¿
​Has probado con read_sav() del paquete haven? Suele importar mejor que las
funciones de foreign.

Cuéntanos.​



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Alternativa a RStudio

[Freddy Omar López Quintero]

2017-03-26 9:50 GMT-03:00 Manuel Spínola :

> Y qué tal Jupyter Notebooks para R


Particularmente yo lo he intentado algunas veces, pero no me ha parecido
muy... fluido. Para preparar código que deba ser distribuido de forma
pedagógica, me parece excelente, pero para desarrollo y trabajo 'real' (por
decirlo de alguna forma), no me ha parecido cómodo.

​Mención aparte merece la experiencia específica con python. Muchos
usuarios python minimizan el hecho que a veces no es transparente mantener
python en el sistema, dejando todo a terceros (como anaconda)​ sin
mencionar esa incómoda convivencia entre python 2 y python 3. No han sido
buenas mis experiencias compartiendo notebooks cuando quienes lo quieren
ejecutar no tienen mucha experiencia previa y falta una librería aquí y
otra allá...

Bueno, esa ha sido mi experiencia. Otros usuarios te comentarán seguramente
las suyas :)

¡Salud!



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] GLM con clusters

[Freddy Omar López Quintero]

2017-03-22 13:48 GMT-03:00 Sebastian Gadea :

> O sea, quitan el supuesto de independencia entre observaciones de un mismo
> grupo, en este lo indica la variable red. Pero ésto no afecta la estimación
> de los coeficientes, si los intervalos de confianza.
>

Podrías utilizar el paquete pglm (glm para datos de panel).

Para ello debes preparar previamente tus datos utilizando pdata.frame con
el id red. Luego en model del pglm indicas 'pooling' (además de la familia
binomial). Algo como:

pglm(y ~ x1+x2, data=tus_datos, family = binomial, model = "pooling")

​Ojalá sea de utilidad.

¡
​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] GLM con clusters

[Freddy Omar López Quintero]

​Hola.​

2017-03-22 12:33 GMT-03:00 Sebastian Gadea :

> quiero realizar una regresión logística,con la función GLM,
> pero ademas quiero designarle a las observaciones clusters.
>

​Si no entiendo mal, quieres estimar un modelo marginal; entonces la opción
inmediata (para mí) está en el paquete gee. Su función del mismo nombre
tiene un ejemplo para datos binarios.

Si el modelo que quieres no es marginal sino un modelo mixto más general,
podrías utilizar (entre muchísimas opciones) el paquete lme4,
específicamente la función glmer.

Cabe sin embargo la posibilidad que haya entendido todo mal y debas darnos
un poco más de información para intentar ayudarte mejor.

¡Salud!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Alternativa a RStudio

[Freddy Omar López Quintero]

​Hola Javier,​

Alguno utiliza una alternativa a RStudio


​En particular, yo me paseo entre RStudio (sobre todo por demostraciones
rápidas que debo hacer) y la consola junto al siempre fiel gedit, de
Gnome.​

Este último par es suficiente para mí. En Windows talvez yo lo intentaría
con NotePad++. Husmeando en la web encuentro el proyecto NppToR:
https://sourceforge.net/projects/npptor/

¡Si lo pruebas y no es mucho abuso, no dejes de darnos una
retroalimentación de él!

¡​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] R - Markdown

[Freddy Omar López Quintero]

2017-03-21 21:50 GMT-03:00 Genaro Contreras :

> Si no la tienes puedes instalar miktek que es una distribución
> de látex para Windows
>

​Y si tienes instalado miktex, debes especificar el lugar donde se
encuentra pdflatex en tu ordenador, que RStudio no lo encuentra. Aunque
puede variar un poco según la versión del Windows que tengas, este enlace
(hecho para Windows 7) puede serte de utilidad:

http://stackoverflow.com/questions/27004849/object-
pdflatex-not-found-on-windows-7

¡Salud!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Recta de regresión en un Scatterplot

[Freddy Omar López Quintero]

2017-03-13 17:00 GMT-03:00 MARIA RUBIO APARICIO :

> Creo que abline solo vale para regresion simple. Alguna otra sugerencia?


Ah, había entendido mal. Quizás esto sea más apropiado:

> plot(x=d_asims, y=d_curts, xlab="Skewness", ylab="Kurtosis")
> > lines(sort(d_asims), predict(model)[order(d_asims)], col=2)
>

¡Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Recta de regresión en un Scatterplot

[Freddy Omar López Quintero]

>
>
> 2017-03-13 13:39 GMT-03:00 Maria Rubio :

> > model=lm(d_curts ~ d_asims + I(d_asims^2))
>  > plot(d_asims,d_curts,xlab="Skewness", ylab="Kurtosis")
>

​Podrías añadir:

>abline(model)​
>

​Seguramente te arroje un warning porque tienes una variable más
(I(d_asims^2)) y el dibujo se hace sobre dos ejes.​

¡
​Salud!​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Función Hoy()

[Freddy Omar López Quintero]

Hola.

Quizás te funcione alguna variación de date().

¡Salud!

2017-03-13 15:08 GMT-03:00 Andres Hirigoyen :

> ¿Existe en R una función Hoy() como en excel? De forma que cada vez que
> corra un srcipt las funciones que usen como dato la fecha del día se
> actualicen automáticamente.
>
> Gracias
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] conversion semanas epidemiologicas a calendario

[Freddy Omar López Quintero]

​Hola.​

2017-02-23 12:49 GMT-03:00 Luis Alfonso LOPEZ ALVAREZ :

> tengo unas series de tiempo
> por semanas epidemiologicas. Alguien tiene un script que permita
> convertirlas/ajustarlas a dias calendario
>

​Una _función_ que parece hacer esto es ​epiweekToDate() de la
librería EpiWeek.

​¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Rutina de cálculo

[Freddy Omar López Quintero]

​Hola.​

2017-02-21 15:08 GMT-03:00 Andres Hirigoyen :

> ¿Cómo puedo automatizar esto para no tener que hacer de forma manual  un
> subset por Región?
>

​Así es un poco difícil ayudar más porque no sabemos cómo es realmente la
estructura de datos o con qué peculiaridades puedes toparte; pero en
general librerías como dplyr (y otras del Wickham Universe :-D) ​son muy
útiles para este tipo de menesteres. Échale un ojo a las 'grouped
operations':

https://cran.rstudio.com/web/packages/dplyr/vignettes/introduction.html

¡Salud!



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Sobre gráficos

[Freddy Omar López Quintero]

2017-02-20 11:10 GMT-03:00 hibiki :

> desgraciadamente no me abre, da un error, es que aki hay muchas
> restricciones absurdas que impiden el acceso incluso a sitios inofensivos
> pero bueno eso no viene al caso, ya tengo el link y en cuanto encuentre la
> via lo reviso


​

Los temas que no vienen al caso también son importantes. Échale un vistazo
al navegador tor

https://www.torproject.org/projects/torbrowser.html.en


​¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Geany y RKWard con R

[Freddy Omar López Quintero]

2017-02-17 15:41 GMT-03:00 Santiago Barbini :

> ¿alguno sabe como se puede configurar correctamente Geany o RKWard para
> enviar el script/s a la consola de R?
>

​Yo no, aunque no creo que sea culpa de las actualizaciones de ubuntu que
haya dejado de funcionar sino del plugin que sospecho que utilizabas al
cual no le han hecho mantenimiento hace varios años... En fin ¿es
imprescindible Geany o RKWard o puedes probar con otros editores de texto?
En este último caso, podrías echarle un ojo a sublime. Aquí un micro
tutorial de cómo echarlo a andar:

http://hagutierrezro.blogspot.cl/2016/10/sublime-text-3-
alternative-to-rstudio.html

¡Salud!​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Ayuda con Date

[Freddy Omar López Quintero]

​Hola.​

2017-02-16 9:34 GMT-03:00 Mauricio Monsalvo :

> ¿Podrían por favor ayudarme a correr la sintaxis correcta?


Esta es otra opción:

​datos<-read.csv('datos.csv', sep=';', row.names=1)
> ​​datos$pprfecbaja<-​ replace(as.Date(datos$pprfecbaja),
> datos$pprfecbaja=="0001-01-01", NA)


​Saludos.​


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Eliminación de filas en data frame según versión del fichero de origen

[Freddy Omar López Quintero]

​Hola.​

2017-01-23 14:27 GMT-03:00 Rubén Coca :

> El caso es que para un mismo id y date debo quedarme con la observación que
> tenga la versión más alta (descartando el resto).
>

Si es válido apoyarse en SQL, yo usaría algo como:

> library(sqldf)

> sqldf("select rowid, id, max(value) from df group by id")

  rowid   id max(value)

1 4 0001   19.57054

2 5 0002   14.70713

3 7 0003   19.34788


donde se reporta el número de la fila por si se quiere hacer otro join.​

​Saludos.​

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] ¿Está R perdiendo la batalla?

[Freddy Omar López Quintero]

¡​Hola!​

2017-01-15 10:54 GMT-03:00 Jesús Para Fernández <
j.para.fernan...@hotmail.com>:

> Hace poco me puse al tema del big data y la verdad es que een este campo
> la decisión parece clara. Python ahora mismo está un paso por delante de R,
> aunque sparklyR puede igualar la contienda.
>
> Pero lo que me empieza a preocupar es que parece que si nos alejamos del
> Big Data, Python tb está ganando adeptos a pasos agigantados. ¿Está
> perdiendo R la batalla?
>

¡​Este es un tema fascinante!​

Yo me pregunto: de forma corporativa ¿ha tenido R un verdadero espacio? En
mi modesta opinión, creo que no ha asustado realmente a una gran
corporación, donde SAS, por ejemplo, ya tiene el mercado hecho. Me aventuro
a pensar que R tiene espacio en compañías más jóvenes pero ¿para grandes
volúmenes de datos ha realmente competido?¿ha sido su target?

​Python es un programa general que en sí mismo puede incluir de todo, como
estadística, optimización, manejo de todo tipo de procesos, etc. etc. y por
eso es muy simpático a mucha gente.​

Mi muy personal opinión es que me asusta más no saber SAS que no saber
Python 😂.

​En el apartado académico, por otro lado, me parece, R es hegemónico para
hacer estadística (sobre todo por estadísticos).​..

¡Salud!

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Big data con R

[Freddy Omar López Quintero]

2017-01-04 13:09 GMT-03:00 Francisco Rodríguez :

> casi todo se reduce a aplicar un algoritmo sobre una tabla de 70 GB
> directamente


En muchos entornos, tristemente, esto es así. Un ejemplo muy difundido
conque yo me he topado es que hay personas que al aplicar estos algoritmos
directamente es que no reparan en la naturaleza de las variables
(categóricas, fechas, etc.)​ quedando cuestiones como que la probabilidad
de éxito de la compra de un producto, depende de la marca (brand) 5.5
(¡cuando las marcas eran categóricas, pero nadie lo notó!). Sin embargo, la
crossvalidación, los nodos, las cajas, todo, arrojó ese como mejor
resultado.

Esto ya sucedía muchos años antes de que el big data, machine learning o el
data scientism estuvieran de moda, y se veía mucho cuando las personas
utiliza(ba)n softwares potentísimos como SPSS, en los que se corre el
riesgo de hacer clic en los botones equivocados y errar todo el trabajo.

No sé a quién se pueda responsabilizar de este tipo de situaciones. ¿La
ansiedad de obtener resultados rápidamente, tal vez? En clases formales en
general se hace énfasis en diferenciar para qué sirven los procedimientos,
cuándo y cómo aplicarlos.

​En fin, respondiendo a

Parece que la tendencia es hacia el uso de entornos hadoop y spark. ¿qué
> opináis al respecto? ¿Es una tendencia pasajera?


no, no creo que sea algo pasajero, al menos pensando como un corporativo
(que no lo soy 😂). Es algo, en mi opinión, que ha venido para quedarse.
Puede ser que hadoop y spark propiamente pasen de moda (no tan) pronto,
pero, computacionalmente, el presente y futuro ya es el cálculo distribuido
y los esfuerzos se irán poniendo en hallar formas de implementar lo
existente y lo por descubrir en nodos, hilos, clusters, nubes...

¡En fin, esto da para mucho debate!

¡Salud!


-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Suma condicionada.-

[Freddy Omar López Quintero]

Hola comunidad,

Les consulto mi duda. Quisiera generar N números aleatorios tal que su suma
esté predeterminada a un valor k. Sé que con la distribución Dirichlet, por
ejemplo, se podrían obtener números aleatorios tal que su suma es la
unidad. Por ejemplo

MCMCpack::rdirichlet(1,rep(10,7))

genera 7 números tal que su suma es 1 (el parámetro rep(10,7) es cualquier
cosa).

Pero me pregunto si se podrá realizar algo parecido variando la suma (no
necesariamente 1) y que los valores no solo sean positivos (como en la
Dirichlet).

¡Gracias de antemano!

-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Identificar porción de una imagen.-

[Freddy Omar López Quintero]

Gracias Carlos y Javier.

No. No es reconocimiento de imágenes.

Por ahora me he contentado con la (muy rudimentaria) solución:

image(array(findInterval(as.numeric(lennon), c(0, 150,  200), left.open =
F, all.inside = T), dim = c(256, 256)), col=grey(seq(0,1,,256)) )

donde el c(0, 150, 200) lo encontré al tanteo. Nada elegante, pero me dice
que los pixeles de la frente (y nariz) tienen valores sobre 150.

¡Mil gracias!

On Thu, Sep 1, 2016 at 5:31 PM,  wrote:

> Estimado Freddy
>
>
>
> Un artículo reciente que se me paso comentar en el correo anterior
>
>
>
> https://www.r-bloggers.com/image-recognition-tutorial-in-
> r-using-deep-convolutional-neural-networks-mxnet-package/
>
>
>
>
>
> Javier Rubén Marcuzzi
>
>
>
> *De: *Freddy Omar López Quintero 
> *Enviado: *jueves, 1 de septiembre de 2016 16:10
> *Para: *r-help-es@r-project.org
> *Asunto: *[R-es] Identificar porción de una imagen.-
>
>
>
> ¡Hola muchachos!
>
>
>
> Tengo un problema que no he podido resolver.
>
>
>
> Les cuento: tengo una imagen y quisiera conocer cuáles son los valores de
>
> los pixeles en un área más o menos delimitada. Yo creo que mi problema
>
> principal es que no sé cómo podría hacer la selección de esta área más o
>
> menos delimitada para luego conocer el valor en cada pixel :-\
>
>
>
> Un ejemplo:
>
>
>
> library(fields)
>
> data(lennon)
>
> image(lennon, col=grey(seq(0,1,,256)) )
>
>
>
> Para esa imagen, digamos, estoy interesado en conocer los (valores de los)
>
> pixeles de la frente del señor Lennon.
>
>
>
> ¿Alguien tendrá una idea de cómo podría proceder?
>
>
>
> ¡Gracias!
>
>
>
> --
>
> «...my role is to be on the bottom of things.»
>
>
>
> Donald Knuth
>
>
>
> [[alternative HTML version deleted]]
>
>
>
> ___
>
> R-help-es mailing list
>
> R-help-es@r-project.org
>
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>



-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Identificar porción de una imagen.-

[Freddy Omar López Quintero]

¡Hola muchachos!

Tengo un problema que no he podido resolver.

Les cuento: tengo una imagen y quisiera conocer cuáles son los valores de
los pixeles en un área más o menos delimitada. Yo creo que mi problema
principal es que no sé cómo podría hacer la selección de esta área más o
menos delimitada para luego conocer el valor en cada pixel :-\

Un ejemplo:

library(fields)
data(lennon)
image(lennon, col=grey(seq(0,1,,256)) )

Para esa imagen, digamos, estoy interesado en conocer los (valores de los)
pixeles de la frente del señor Lennon.

¿Alguien tendrá una idea de cómo podría proceder?

¡Gracias!

-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Muestrear de una normal multivariante.-

[Freddy Omar López Quintero]

2016-07-17 17:49 GMT-04:00 Carlos J. Gil Bellosta :

> Si son valores arbitrarios los que quieres colocar en Sigma_UZ, hazlos
> más pequeños (p.e., dividiéndolos por 10).
>

​¡Gracias!¡Funciona perfectamente!

Si quisiera establecer escenarios como: covarianzas bajas, covarianzas
moderadas y altas ¿tiene algún consejo que pueda seguir? (imagino que en
este punto solo cuento con covarianzas bajas).

¡Gracias y disculpen el abuso!


-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Caracteres especiales en gráficos ggplot

[Freddy Omar López Quintero]

Hola.

¿Podrías darnos un ejemplo de lo que te sucede para intentar reproducirlo?
Así es más fácil que te ofrezcan algunas pistas de la solución.

¡Hasta luego!


2016-07-17 2:38 GMT-04:00 Leslie Vargas Vásquez :

> Hola, no logro que los gráficos con ggplot reconozcan los caracteres
> especiales como tildes y ñ conoces si se requiere de un paquete especial
> para  trabajar estos gráficos con MacBook.
>
> Gracias!
>
> Leslie Vargas
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Muestrear de una normal multivariante.-

[Freddy Omar López Quintero]

¡Hola a todos!

Estoy intentando muestrear de una normal multivariante donde hay dos grupos
de variables que deben tener una relación "manipulable" entre sí pero
ignoro cómo hacerlo.

Les cuento, he intentado lo siguiente:

# covarianzas del primer grupo de variables:
Sigma_U <- matrix(c(.25, .2, .2, .25), ncol=2)

# covarianzas del segundo grupo de variables:
Sigma_W <- diag(2)

# covarianzas _arbitrarias_ entre los dos grupos de variables
Sigma_UZ <- matrix(rnorm(4), nrow=2)

# consolidación de las covarianzas anteriores:
Sigma<-rbind(
cbind(Sigma_U, Sigma_UZ),
cbind(t(Sigma_UZ), Sigma_W)
)

# muestreo:
MASS::mvrnorm(1, mu=rep(0, 4), Sigma=Sigma)

​De donde recibo:

Error in mvrnorm(1, mu = rep(0, 4), Sigma = Sigma) :
  'Sigma' is not positive definite


El error (creo yo) está la generación de esas covarianzas arbitrarias para
que esta matriz consolidada, Sigma, sea definida positiva. Mi necesidad es
que las matrices ​Sigma_U y Sigma_W sean las que he definido pero saber o
poder ubicar las covarianzas de Sigma_UZ para que no arroje error.

¿Alguien sabe cómo podría hacer?¿Qué pasos debería seguir?

¡Gracias!

-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Error en "optim" modelo APARCH sstd

[Freddy Omar López Quintero]

​Hola.​

2016-05-25 13:06 GMT-04:00 Mª Ángeles Navarro :

> Mi problema es que he ejecutado el modelo,con el mismo código, en
> diferentes series (de mismo tamaño muestral) obteniendo resultados
> adecuados,  pero en concreto con una de las series me da el siguiente error:
> Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = excess) :
>  non-finite value supplied by optimIn addition: Warning messages:1: In
> sqrt(diag(fit$cvar)) : NaNs produced2: In sqrt(diag(fit$cvar)) : NaNs
> produced3: In sqrt(diag(fit$cvar)) : NaNs produced4: In
> sqrt(diag(fit$cvar)) : NaNs produced5: In sqrt(diag(fit$cvar)) : NaNs
> produced6: In sqrt(diag(fit$cvar)) : NaNs produced
>


Entiendo que has aplicado el mismo modelo a cada una de las series, pero
esto no tiene por qué ser así. El modelo que intentas estimar no está
funcionando numéricamente para este conjunto de datos. ¿Has probado otro
orden para esa serie en especial?

​Salud.​


-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Latex y gráfico

[Freddy Omar López Quintero]

​Hola.​

2016-05-13 0:22 GMT-03:00 Javier Marcuzzi :

> Probé de su sugerencia el ejemplo 3.4, pero me falla, copio y pego todo el
> código junto al error
>
>
>
> > library(tikzDevice) tikz(
> ​​
> ’figs/simpleEx.tex’,width=3.5,height=3.5) plot(1,main=’Hello World!’) 
> dev.off()
>
> Error: unexpected symbol in "library(tikzDevice) tikz"
>
>
​Si realmente copiaste y pegaste desde el pdf, sospecho que el error tuvo
que ver con el símbolo del apóstrofo que está antes de figs/sigmplEx.tex (o
con todos ellos). Prueba por favor reemplazando todos los apóstrofos con
unos tipeados por ti mismo. Ojalá funcione.

¡Salud!​



-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Latex y gráfico

[Freddy Omar López Quintero]

Hola Javier,

Yo creo que TikZ es el mejor camino. Échale un ojo a al paquete tikzDevice:

https://cran.r-project.org/web/packages/tikzDevice/vignettes/tikzDevice.pdf

y alucinarás con lo que te ofrece :D...

Au revoir!

2016-05-09 19:13 GMT-03:00 Javier Marcuzzi 
:

> Estimados
>
> Hay varias alternativas para utilizar latex dentro de R, concretamente por
> ejemplo una ecuación matemática dentro de un gráfico.
>
> Pero para un caso se me ocurre utilizar líneas que conecten partes de
> gráficos, por ejemplo, supongamos que hay dos barras con tres o cuatro
> divisiones, digamos, 25% cada una, de la primera desde una de estas me
> gustaría realizar una flecha, hacia la segunda barra.
>
> En latex hay algunas alternativas como xy-pic, podría utilizar esto dentro
> de un archivo Rnw, colocando el gráfico en una matriz de latex, pero
> posiblemente alguna persona escribió al respecto, y me gustaría poder leer
> la experiencia de otros para no cometer los mismos errores, u obtener
> soluciones alternativas a las que estoy pensando.
>
> Gracias
>
> Javier Rubén Marcuzzi
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«...my role is to be on the bottom of things.»

Donald Knuth

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Curso R por Microsoft

[Freddy Omar López Quintero]

Tantas manos ya pondrán este caldo morado... Me espanta tanto microsoft
tras R.

Saludos a todos.

2015-08-13 21:04 GMT-03:00 Javier Rubén Marcuzzi <
javier.ruben.marcu...@gmail.com>:

> Estimados
>
> Vi publicado un curso sobre R donde Microsoft está interesado, desconozco
> si nos es de utilidad, pero creo que es bueno que una empresa con mucho
> capital económico como tecnológico invierta en R. Habrá que esperar a
> futuro para ver los resultados y si son de nuestro agrado, pero por ahora
> podemos compartir el curso.
>
> http://blog.datacamp.com/microsoft-launches-first-free-online-r-course-edx/
>
> Javier Rubén Marcuzzi
> Técnico en Industrias Lácteas
> Veterinario
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Duda interpolación (package ' gstat ')

[Freddy Omar López Quintero]

Hola Marcos,

¿El problema persiste si pruebas con un subconjunto de los datos?

Saludos.

2015-08-06 12:34 GMT-03:00 Marcos Bermejo :

> Sale plano sí.
>
> Ya se que sin tener los datos y el código es un poco difícil, pero es que
> mis datos ocupan mucho, es imposible.
> Seguiré mirando por internet.
>
> Muchas gracias Rubén.
>
> Un saludo,
>
>
> > To: r-help-es@r-project.org
> > From: rubenfca...@gmail.com
> > Date: Thu, 6 Aug 2015 14:21:47 +0200
> > Subject: Re: [R-es] Duda interpolación (package ' gstat ')
> >
> > Hola Marcos,
> >
> >  Parece que el problema es con el ajuste del variograma (sale
> > plano?), sin más información no se exactamente que puede estar pasando...
> >
> >  Si me envías el código completo y los datos lo miro con más detalle
> > (e incluso te doy una alternativa no paramétrica con el paquete npsp).
> >
> >  Un saludo, Rubén.
> >
> >
> >
> > El 04/08/2015 a las 11:24, Marcos Bermejo escribió:
> > > Hola,
> > >
> > ># Hacemos el KED. Ver funci�n "krige()":
> > >  KED.rad <- krige(
> > >formula=pluvPcp~layer,  # covariable ->
> radar
> > >locations=lluvia.rad.pluv.spdf,
> > >newdata=radarGrid,  # podr�a ser
> cualquier objeto Spatial
> > >model=v.fit,# modelo de
> semivariograma.
> > >maxdist=Inf
> > >  )
> > >
> > > Esta es la funci�n que me interpola los datos de lluvia. El error que
> me da es:
> > >
> > > "solve.c", line 88: singular matrix in function Usolve()
> > >
> > > "lufactor.c", line 208: singular matrix in function m_inverse()
> > > Error in predict.gstat(g, newdata = newdata, block = block, nsim =
> nsim,  :
> > >m_inverse
> > > In addition: Warning message:
> > > In fit.variogram(vg.aux, model = vgm(psill = 0.1, model = "Gau",  :
> > >Warning: singular model in variogram fit
> > >
> > >
> > > Mi funci�n del variograma es :
> > > v.fit <- fit.variogram(vg.aux, model=vgm(psill=0.15, model='Gau',
> range=5000,
> > > nugget=0.05))
> > >
> > > �Alguien me podr�a ayudar?
> > >
> > > Gracias de antemano.
> > >
> > > Un saludo,
> > >
> > > [[alternative HTML version deleted]]
> > >
> > >
> > >
> > > ___
> > > R-help-es mailing list
> > > R-help-es@r-project.org
> > > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
> >
> >   [[alternative HTML version deleted]]
> >
> > ___
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] duda

[Freddy Omar López Quintero]

2015-06-22 12:38 GMT-03:00 jbetancourt :

> Gracias funcionó
>
> si , efectivamente , así se eliminan, pero es que son casos diarios de una
> enfermedad y no qusiera que el espacio se eliminara ?como tratarlo de
> manera que se evalue el día sin tener en cuenta el valor, como un valor
> perdido?
>
>

​Puedes sustituir los valores con:

x[x==0]<-NA​

​Pero no creo que sea lo más conveniente porque si no hubo casos esos esos
días, el 0 tiene sentido. Diferente es un valor faltante donde no tienes
información.

Opinando un poco más de la cuenta, pienso que quizás un método para
variables de conteo (¿Poisson, binomial negativa?) sea más conveniente para
lo que has descrito.​

Saludos.

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] serie y ceros

[Freddy Omar López Quintero]

​Holap.

es decir que no se tengan en
>  cuenta para el análisis de la serie.
>

​Reemplaza x con:​


x<-​x[x!=0]​
​
¡​Salud!​

2015-06-22 11:04 GMT-03:00 jbetancourt :

>  Estimados
>
>
> En el presente script quisiera que me mostraran la manera de que los valores 
> cero que aparecen en la serie de casos diarios no se lean, es decir que no se 
> tengan en
>  cuenta para el análisis de la serie. Adjunto la serie
>
>
>
> rm(list = ls())
>  x<-scan('D:/Public/Documents/R/EPICALC/e.csv')
> x
>  library(forecast)
> fit <- auto.arima(x)  #que no tenga en cuenta los valores cero
> plot(forecast(fit,h=15))
>
> Saludos
> José
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] funciones en estructura data.table

[Freddy Omar López Quintero]

​​No limpió, no fijó y no dio esplendor nunca...


2015-06-17 21:26 GMT-03:00 Patricio Fuenmayor Viteri <
patricio.fuenma...@outlook.com>:

>
>
> -Mensaje original-
>
> De: "Patricio Fuenmayor Viteri" 
> Enviado: 17 de junio de 2015 6:37 PM
> Para: r-help...@r-ptoject.org, patricio.fuenma...@outlook.com
> Asunto: funciones en estructura data.table
>
> hola...
> luego de haber migrado la mayoria de mis procesos a data.table. aun me
> kedan ciertos detalles por afinar.. eate es un ejemplo de eso..
> tengo una tabla de emails y necesito hacer un analisis de dominios... para
> separar el nombre de usuario el dominio uso la funcion str_split del
> paquete stringr... pero el resultado me queda en una sola columna y
> separada por una coma... mi pregunta es ... como puedo hacer ke los
> multiples resultados de una funcion ... se guarden en diferentes columnas
> de mi data.table
> email01[,.(str_split(email))]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Problemas al cargar Rcomander en consola de Rstudio

[Freddy Omar López Quintero]

Holap.

Loading required package: splinesLoading required package:
> RcmdrMiscLoading required package: carLoading required package:
> sandwichError : .onLoad failed in loadNamespace() for 'Rcmdr',
> details:
>   call: structure(.External(.C_dotTclObjv, objv), class = "tclObj")
>   error: [tcl] invalid command name "tk_messageBox".
> In addition: Warning message:In fun(libname, pkgname) : couldn't
> connect to display ":0"Error: package or namespace load failed for
> ‘Rcmdr’


Quizás peco de ingenuo pero: ¿has verificado que car y sandwich estén
instalados?

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Regresión logística

[Freddy Omar López Quintero]

​Holap.​

ran out of iterations and failed to converge


​Prueba aumentando el número de iteraciones, con el argumento maxit:

​GLM <- bigglm(In.hospital_death ~ GCS + BUN, data = DatosGLM, family =
> binomial(logit), maxit=1000)​


​Salud.​

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Serie temporal interrumpida del tipo AirPassengers

[Freddy Omar López Quintero]

2015-06-12 18:42 GMT-03:00 Dalios Castellano Marrero <
dcastell...@aqualogy.net>:

> Necesito detectar si existe o no un cambio de tendencia y si dicho cambio
> es significativo, para una serie temporal del tipo AirPassengers, en la que
> a partir de un determinado momento se ha hecho una campaña (supongamos que
> una promoción de vuelos).
>

​Esta descripción encaja con lo que le llaman análisis de intervención. El
procedimiento no es inmediato y se debe hacer una búsqueda cuidadosa del
modelo adecuado. La librería TSA (del libro de ​Jonathan Cryer and Kung-Sik
Chan) tiene la función arimax() que aborda este problema. Hay una guía aquí:

http://econometricsense.blogspot.com/2012/01/time-series-intervention-analysis-wih-r.html

También se podrían analizar cambios estructurales. Te recomiendo que
revises algunos de los trabajos de Achim Zeileis al respecto, por ejemplo:

http://cran.r-project.org/web/packages/strucchange/vignettes/strucchange-intro.pdf


> Para ello he pensado varios métodos:
>
> Usar la descomposición espectral de la muestra [decompose(AirPassengers)]
> y luego una Regresión Joinpoint para la serie de tendencia,
> pero desconozco como hacer eso con R.
>
> Y por otro lado, he leído algo sobre la librería segmented para analizar
> cambios de tendencias a lo largo de series temporales interrumpidas,
> pero también desconozco como realizar dicho procedimiento en R.
>

De estos puntos estoy seguro que habrá alguien en el grupo que conozca más
:)

¡​Salud!​

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] pregunta

[Freddy Omar López Quintero]

​Holap.​

x <- c 
(1,3,1,1,1,1,2,1,1,2,1,1,2,1,1,1,2,1,2,2,1,1,1,3,1,1,2,2,2,2,1,1,1,2,1,1,1,1,1,2,2,5,1,1,1,1,2,1,1,1,2,4,1,2,1,3,1,1,1,1,3,1,2,1,1,3,1,3,3,1,2,1,2,2,2,3,1,2,1,1,1,1,2,1,1,1,1,1,1)
>
>
> names(x) <- format(index(zooreg(1:89, start =as.Date("2015-01-01"))), 
> '%Y-%m-%d')
>

​Puedes usar:

x<-scan('ruta/a/tu/archivo/a.txt')​

y luego el código que ya utilizas.

​¡Salud!​

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Sobre optimización con Rcpp.-

[Freddy Omar López Quintero]

2015-06-09 14:36 GMT-03:00 :

>  En lo personal nunca compile algo escrito por mi utilizando Rcpp, pero si
> algo de C++, hay compiladores (creo que ahora deben ser casi todos), que
> tienen una opción de autovectorización, posiblemente una alternativa es
> buscar una optimización utilizando vectores. Otra puede ser la programación
> en paralelo. Posiblemente hoy en día hay muchas alternativas pero requieren
> tener conocimientos informáticos elevados (por lo menos para mí).


​Gracias amigo Rubén. Estas soluciones por ahora van más allá de lo que
ando buscando, creo. Mi problema es la maximización de una función, como lo
hace la función optim​ de R, pero me intriga saber si hay una manera sabia
o establecida para usar Rcpp en estos menesteres.

¡Salud!

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] librerias para programar una interfaz grafica en R ?

[Freddy Omar López Quintero]

Hola Eric,

Michael Lawrence y John Verzani tienen un libro relativamente reciente al
respecto:

Lawrence & Verzani: Programming Graphical User Interfaces in R, (2012).

Ojalá te sirva.

Saludos.

2015-06-09 17:07 GMT-03:00 eric :

> Estimados, he estado mirando en la internet si es posible programar con R
> una interfaz grafica para pedir el ingreso de datos por ejemplo, pero
> invariablemente obtengo paginas que hablan de interfaces graficas (rkward,
> rstudio, r-commander, etc) para R o de las librerias para hacer graficos
> (ggpolt, lattice, etc) ... alguien me puede orientar donde buscar o si
> estas existen ?
>
> Saludos y gracias,
>
> Eric.
>
>
>
>
> --
> Forest Engineer
> Master in Environmental and Natural Resource Economics
> Ph.D. student in Sciences of Natural Resources at La Frontera University
> Member in AguaDeTemu2030, citizen movement for Temuco with green city
> standards for living
>
> Nota: Las tildes se han omitido para asegurar compatibilidad con algunos
> lectores de correo.
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Sobre optimización con Rcpp.-

[Freddy Omar López Quintero]

Hola compañeros de R,

Antes he utilizado Rcpp (y armadillo) para lograr velocidad en procesos
llenos de bucles y de álgebra lineal; pero no me había visto en la
necesidad de maximizar (optimizar) ninguna función (una verosimilitud, en
este caso).

¿Alguien tiene conocimiento de cuál es la forma más conveniente para
optimizar funciones dentro de un programa en Rcpp?

¡Muchas gracias!

-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] comportamiento de data.table al hacer calculos por grupos

[Freddy Omar López Quintero]

Holap.

Reproduje tu ejemplo y ciertamente las filas se duplican innecesariamente.
No supe cómo corregir este comportamiento. Si no se pudiera arreglar esto,
por favor considera la librería sqldf y obtendrás el resultado correcto:

>sqldf("select sol, dia, con, avg(media) as ave from dd group by sol, dia,
con")
   sol dia con   ave
1  con   1   0 -29.37000
2  con   2   0 -31.65000
3  con   3   0 -28.25000
...

donde dd es tu conjunto de datos.

Salud.

2015-05-14 20:00 GMT-03:00 eric :

> Estimada comunidad tengo un problema del que no encuentro datos que me
> ayuden mucho en la web.
>
> Estoy haciendo calculos por grupos con data,table. Tengo un archivo
> (zp.res) con tres columnas que clasifican los datos (sol, con, dia)  y una
> columna de datos numericos (media), de la siguiente forma:
>
>  sol con dia media
>   1: con   0   1 -22.6
>   2: con   0   1 -36.6
>   3: con   0   1 -35.6
>
> y quiero calcular el promedio de "media" (la col 4) agrupando con las
> variables sol,con,dia. Lo hago de la siguiente forma:
>
> med <- zp.res[, mean(media), by="sol,dia,con"]
>
> cuando reviso "med" esta todo bien, se han calculado las medias y el
> objeto tiene solo las filas que corresponden a los promedios con sus
> respectivas columnas sol,con,dia que los identifican. Pero como notaran por
> el codigo, la columna con el resultado no tiene un nombre asi es que R la
> bautiza como V1. Quise ponele un nombre y agregue este al codigo de la
> siguiente forma:
>
> med <- zp.res[, "ave":=mean(media), by="sol,dia,con"]
>
> el problema es que ahora el objeto med tiene el mismo numero de filas que
> zp.res y repite el promedio para cada dato dentro del grupo obtenido con
> by="sol,dia,con". Esto no me sirve pues tengo que graficar los promedios ...
>
> m puede que mi explicacion sea algo confusa, espero que me entiendan.
>
> Encontre que luego puedo renombrar la columna, pero no lo quiero hacer,
> pues pienso que deberia ser estandar poder ponerle el nombre y que se
> construya bien el objeto con los promedios inmediatamente. Ademas el
> promedio es solo uno de los calculos que debo hacer y los otros tambien
> quedan con el nombre V1 en la mismo data.table.
>
> Alguna idea de como hacerlo ?
>
> Adjunto archivo con datos.
>
> Saludos y gracias, eric.
>
>
>
>
>
> --
> Forest Engineer
> Master in Environmental and Natural Resource Economics
> Ph.D. student in Sciences of Natural Resources at La Frontera University
> Member in AguaDeTemu2030, citizen movement for Temuco with green city
> standards for living
>
> Nota: Las tildes se han omitido para asegurar compatibilidad con algunos
> lectores de correo.
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es