Re: [R-es] realizar ANOVAs en Loop

[Jorge I Velez]

Hola Yésica,

A lo mejor hay otras formas, pero esta funciona bien:

R> u <- unique(datos$iteraccion)
R> resultado <- lapply(u, function(i){
  aov(valor ~ Grupo, data = subset(datos, iteraccion == i))
})
R> names(resultado) <- u

Para acceder a la tabla ANOVA para iteraccion T simplemente haces

> resultado[['T']]
Call:
   aov(formula = valor ~ Grupo, data = subset(datos, iteraccion ==
i))

Terms:
   Grupo Residuals
Sum of Squares   0.70326  62.83939
Deg. of Freedom213

Residual standard error: 2.19859
Estimated effects may be unbalanced

Si quieres obtener los estadísticos correspondientes basta con hacer

R> summary(resultado[['T']])
Df Sum Sq Mean Sq F value Pr(>F)
Grupo2   0.70   0.352   0.073   0.93
Residuals   13  62.84   4.834

Espero sea de utilidad.

Saludos,
Jorge.-


On Tue, Jul 18, 2023 at 1:25 AM Yesica Pallavicini Fernandez <
yesipa...@gmail.com> wrote:

> Buenos días y gracias de antemano por vuestra ayuda.
>
> Necesito realizar una serie de ANOVAS en loop.
> Os adjunto unos datos ficticios en este email.
> Dichos datos tienen 3 variables:
> 1)Valor: corresponde a la variable dependiente y es numérica
> 2) Grupo: Corresponde a la variable independiente y es u factor
> 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los
> ANOVAs con las variables anteriores y es un factor.
>
> Abajo os pego el código con el que he estado trabajando pero que no lo
> tengo bien, porque no puedo acceder a los resultados.
> Os agradeceria mucho si:
> -Me podéis ayudar a mejorar este código para que funcione
> -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop
> porque no he dado con los blogs adecuados.
> -Que me digáis cuál es vuestro libro/web de estadística favorito para
> profundizar en la matemática detrás de los análisis.
>
> Mil gracias compañeros
> Yésica
>
> library(agricolae)
> library(readxl)
> datosa_fict <- read_excel("datosa-fict.xlsx")
> #Cambiar nombre a la base de datos
> datos=datosa_fict
> #Copiar la variable sobre la que hacer el loop
>
> iter=datos$iteraccion
> #Crear el loop
> for(i in iter) {
>   res=aov(valor~Grupo,data=datos)
>
> }
> #Salvar los resultados
> ANOVA(res)
> ___
> 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

Re: [R-es] Importar gráficos Word

[Jorge I Velez]

Hola Javier,
Usa Quarto (https://quarto.org/). En 5 mins tienes un documento funcional
de Word con todo lo que necesitas y sin tantos pasos.
Saludos,
Jorge.-

El El mar, 27 de jun. de 2023 a la(s) 10:17 p. m., Javier Gómez Gonzalez <
zaraga...@gmail.com> escribió:

> Hola a todos.
>
> Quisiera saber cual es la mejor forma de insertar gráficos hechos con
> ggplot2 en Word. Ya sé que lo mejor es usar R-Markdown, pero todavía no he
> empezado con él. Lo que hago es guardar los gráficos con ggsave en formato
> png con una resolución de 600 y luego los inserto en un documento de Word.
> Suelo utilizar una relación ancho alto 16:9 ¿Alguna sugerencia respecto a
> la resolución y a la relación ancho alto a usar en los gráficos?
>
>  Javier Gómez González
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Analisis multivariante con Primer

[Jorge I Velez]

Javier,
Por qué no ir más allá, usando R?
Mira
https://juba.github.io/explor/
Saludos,
Jorge.-

El El mié, 26 de abr. de 2023 a la(s) 9:55 a. m., Javier Gómez Gonzalez <
zaraga...@gmail.com> escribió:

> Hola amigos:
>
> Quisiera saber si alguna opción en R del software primer para análisis de
> datos multivariantes en ecología
>
> Gracias
>
>
> Javier
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] (sin asunto)

[Jorge I Velez]

Manuel,

Si los archivos están en workig directory, basta con que hagas

files <-  list.files(pattern = ".xlsx")
list_df <- lapply(files, function(f){
  data <- read_excel(f)
  data
}

Para acceder al primer archivo basta con hacer

list_df[[1]]

Si los archivos tiene el _mismo_ nombre en las columnas, puedes unirlos
haciendo

datos_completos <- do.call(rbind, list_df)

Espero sea de utilidad.

Saludos,
Jorge.-



On Tue, Feb 14, 2023 at 12:08 AM Manuel Mendoza 
wrote:

> Buenos días, tengo un conjunto de 10 bases de datos en excel cuyos nombres
> van de BIO1 a BIO10. Quiero hacer un loop en el que me coja una cada vez.
> He probado lo que os pongo abajo (y otras cosas parecidas) y obviamente no
> funciona. ¿Conocéis alguna forma?
> Gracias como siempre,
> Manuel
>
> BIOs<-
> ("BIO1”,”BIO2”,”BIO3","BIO4”,”BIO5”,”BIO6”,”BIO7”,”BIO8","BIO9”,”BIO10”)
> for (j in 1:length(BIOs)) {
>   data <- read_excel("BIOs[j].xlsx")
> ...
>
> [[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

Re: [R-es] minera de texto 2

[Jorge I Velez]

José,

Puedes usar

discurso <- scan('ruta/al/archivo/de/texto.txt', what = 'character')

Luego puedes tomar el objeto "discurso" y procesarlo adecuadamente.

Saludos,
Jorge.-


On Wed, Oct 26, 2022 at 12:26 PM Jose Betancourt Bethencourt <
betans...@gmail.com> wrote:

> -Estimados
>
> Como lograr en mineria de texto en R que se lea un archivo de txt que
> tiene un discurso escrito, no en columnas, sino texto integro
>
> saludos
>
> José
>
> ___
> 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

Re: [R-es] error en matriz de correlación

[Jorge I Velez]

Ricardo,

No sé cuál es el objetivo del análisis de correlación y entiendo que la
respuesta de Marcelino responde a lo que preguntas.  Sin embargo, yendo un
poco más allá, me pregunto la validez de lo que quieres hacer, sobretodo
con la estructura de datos que tienes a la mano.

Al parecer tienes medidas repetidas --- tiene sentido *esa *matriz de
correlación que quieres calcular?   Quizás esa sea la primera pregunta que
debes resolver.

Saludos,
Jorge.-




On Mon, Sep 26, 2022 at 7:12 AM Ricardo Alva  wrote:

> Muchas gracias por la aclaración, ya me estaba preocupando jajajaja.
>
>
>
> Enviado desde mi Galaxy
>
>
>
>  Mensaje original 
> De: Marcelino de la Cruz Rot 
> Fecha: 26/9/2022 2:19 a. m. (GMT-05:00)
> A: r-help-es@r-project.org
> Asunto: Re: [R-es] error en matriz de correlación
>
> Hola Ricardo:
>
> Es correcto lo que dices de declarar las variables según su tipo. El
> "problema" con cor() es que explícitamente requiere que las variables
> sean numéricas. Sin embargo, cuando calcula el coeficiente de
> correlación de Kendall (o el de Spearman) lo hace en función del *orden*
> de los valores de esa variable numérica.
> Es decir, no es R el que sabe que Kendall es la prueba que corresponde a
> tus datos, sino tú como usuario, y por eso te pide que incluyas el
> argumento method = "kendall" cuando los datos son ordinales (o cuando
> son numéricos pero por su distribución no puedes confiar en la
> aproximación de Pearson).
>
> Un saludo,
> Marcelino
>
>
> El 26/09/2022 a las 9:00, Ricardo Alva escribió:
> > Amigos gracias por sus respuestas. Disculpen mi ignorancia, pero no se
> > supone que antes de que empezar a trabajar con el data frame, no se
> > debe de declarar el tipo de variable que le corresponde a cada una y
> > así el sistema reconozca si la variable es numérica, ordinal,
> > dicotómoca o categorica? Por ejemplo en este caso, si le pongo
> > numérica a todas las variables, el sistema como va saber que kendall
> > es la prueba que le corresponde a mis datos dado que son ordinales, o
> > este tema de numéricas es sólo una excepción para que corra el algoritmo?
> >
> > Siempre trabajo declarando previamente las variables y no he tenido
> > inconvenientes para ejecutar por ejemplo kmedias, kmodas, kprototipe,
> > PCA, FAMD, o lo estoy haciendo mal?
> >
> > Agradesco sus comentarios.
> >
> >
> >
> >
> >
> > Enviado desde mi Galaxy
> >
> >
> >
> >  Mensaje original 
> > De: Marcelino de la Cruz Rot 
> > Fecha: 26/9/2022 1:06 a. m. (GMT-05:00)
> > A: r-help-es@r-project.org
> > Asunto: Re: [R-es] error en matriz de correlación
> >
> > Hola Ricardo:
> >
> > Por precisar lo que dice Emilio, la conversión debe hacerse por
> > columnas, por ejemplo:
> >
> > a <- apply(a, 2, as.numeric)
> >
> > Un saludo,
> > Marcelino
> >
> >
> > El 26/09/2022 a las 7:54, Emilio L. Cano escribió:
> > > Hola Ricardo,
> > >
> > > Tienes que convertir las variables a numéricas, como indica el error
> > y la ayuda de la función. Como tus valores están ordenados,
> > simplemente convierte con as.numeric() y ya te funcionará.
> > >
> > > Un saludo,
> > >
> > > Emilio L. Cano
> > > http://emilio.lcano.com
> > >
> > >
> > >
> > >
> > >> El 26 sept 2022, a las 3:48, Ricardo Alva 
> > escribió:
> > >>
> > >> Amigos buen día.
> > >> Estoy intentando generar la matriz de correlaciones por cualquiera
> > de los 3 métodos de la función cor, pero siempre me arroja el mismo
> error.
> > >>
> > >> cor(a,method = 'kendall')
> > >> Error in cor(a, method = "kendall") : 'x' must be numeric
> > >>
> > >> a es un data frame con 61 filas, compuesto de 5 variables, todas
> > ordinales y cargadas desde un archivo de excel que mantiene la
> > siguiente forma:
> > >> agencia|personeria|constitucion|producto|actividad
> > >> 1 1 2 1 1
> > >> 2 1 5 1 1
> > >> 2 2 41 5
> > >> 2 2 24 1
> > >>
> > >> He tratado consignado todas las variables como ordinales y luego
> > también como factor, pero me sigue saliendo el mismo error, por mas
> > que elijo cualquier método de correlación.
> > >>
> > >> #Convirtiendo todas las variables a factor:
> > >> str(a)
> > >> 'data.frame': 61 obs. of  5 variables:
> > >> $ agencia : Factor w/ 3 levels "1","2","3": 1 2 2 2 2 2 1 2 2 2
> ...
> > >> $ personeria  : Factor w/ 2 levels "1","2": 1 1 1 1 1 2 1 2 2 2 ...
> > >> $ constitucion: Factor w/ 5 levels "1","2","3","4",..: 2 2 2 2 5 2
> > 1 3 4 5 ...
> > >> $ producto: Factor w/ 2 levels "1","4": 1 1 1 1 1 1 1 1 1 1 ...
> > >> $ actividad   : Factor w/ 5 levels "1","2","3","4",..: 1 1 1 1 1 1
> > 5 1 3 2 ...
> > >>
> > >> #Convirtiendo todas las variables a ordinales.
> > >> str(a)
> > >> 'data.frame': 61 obs. of  5 variables:
> > >> $ agencia : Ord.factor w/ 3 levels "1"<"2"<"3": 1 2 2 2 2 2 1 2
> > 2 2 ...
> > >> $ personeria  : Ord.factor w/ 2 levels "1"<"2": 1 1 1 1 1 2 1 2 2 2
> ...
> > >> $ constitucion: Ord.factor w/ 5 levels "1"<"2"<"3"<"4"<..: 2 2 2 2
> > 5 2 1 3 4 5 ...
> > >> $ producto

Re: [R-es] error en matriz de correlación

[Jorge I Velez]

Hola Carlos,
Qué quieres hacer exactamente?
Lo que intentas tiene poco sentido.
A lo mejor si nos indicas qué estás buscando sea más fácil ayudarte.
Saludos,
Jorge.-

El El dom, 25 de sep. de 2022 a la(s) 8:48 p. m., Ricardo Alva <
kalo_a...@hotmail.com> escribió:

> Amigos buen día.
> Estoy intentando generar la matriz de correlaciones por cualquiera de los
> 3 métodos de la función cor, pero siempre me arroja el mismo error.
>
> cor(a,method = 'kendall')
> Error in cor(a, method = "kendall") : 'x' must be numeric
>
> a es un data frame con 61 filas, compuesto de 5 variables, todas ordinales
> y cargadas desde un archivo de excel que mantiene la siguiente forma:
> agencia|personeria|constitucion|producto|actividad
> 1 1 2 1 1
> 2 1 5 1 1
> 2 2 41 5
> 2 2 24 1
>
> He tratado consignado todas las variables como ordinales y luego también
> como factor, pero me sigue saliendo el mismo error, por mas que elijo
> cualquier método de correlación.
>
> #Convirtiendo todas las variables a factor:
>  str(a)
> 'data.frame': 61 obs. of  5 variables:
>  $ agencia : Factor w/ 3 levels "1","2","3": 1 2 2 2 2 2 1 2 2 2 ...
>  $ personeria  : Factor w/ 2 levels "1","2": 1 1 1 1 1 2 1 2 2 2 ...
>  $ constitucion: Factor w/ 5 levels "1","2","3","4",..: 2 2 2 2 5 2 1 3 4
> 5 ...
>  $ producto: Factor w/ 2 levels "1","4": 1 1 1 1 1 1 1 1 1 1 ...
>  $ actividad   : Factor w/ 5 levels "1","2","3","4",..: 1 1 1 1 1 1 5 1 3
> 2 ...
>
> #Convirtiendo todas las variables a ordinales.
> str(a)
> 'data.frame': 61 obs. of  5 variables:
>  $ agencia : Ord.factor w/ 3 levels "1"<"2"<"3": 1 2 2 2 2 2 1 2 2 2
> ...
>  $ personeria  : Ord.factor w/ 2 levels "1"<"2": 1 1 1 1 1 2 1 2 2 2 ...
>  $ constitucion: Ord.factor w/ 5 levels "1"<"2"<"3"<"4"<..: 2 2 2 2 5 2 1
> 3 4 5 ...
>  $ producto: Ord.factor w/ 2 levels "1"<"4": 1 1 1 1 1 1 1 1 1 1 ...
>  $ actividad   : Ord.factor  w/ 5 levels "1"<"2"<"3"<"4",..: 1 1 1 1 1 1 5
> 1 3 2 ...
>
> Por favor alguien que me pueda decir que estoy haciendo mal, o es que las
> variables en el archivo de excel deben de estar en texto?, lo cual me
> parecería algo raro.
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] RMD

[Jorge I Velez]

Hola José,

Has intentado agregando

lang: "es-ES"

en el YAML?

Espero sea de utilidad.

Saludos,
Jorge.-


On Thu, Sep 22, 2022 at 5:53 AM Jose Betancourt Bethencourt <
betans...@gmail.com> wrote:

> Estimados cuando hago un Rmd , al salir en html los acentos no salen
> sino otros signos no deseados. ?Como puedo evitar eso?
>
> saludos
> --
> Dr. Jose A. Betancourt Bethencourt
> Universidad de Ciencias Medicas Carlos j. Finlay
>
> ___
> 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

Re: [R-es] (sin asunto)

[Jorge I Velez]

Manuel,

Si entiendo bien, algo como esto te servirá:

# máximo por fila
apply(probs, 1, max, na.rm = TRUE)

# probabilidades promedio por columna
colMeans(probs, na.rm = TRUE)

# probabilidad promedio para todos
mean(c(probs), na.rm = TRUE)

Saludos,
Jorge.-

El El vie, 26 de mar. de 2021 a la(s) 1:55 p. m., Manuel Mendoza <
mmend...@fulbrightmail.org> escribió:

> Muy buenas, tengo una matriz (probs) de 100 x 2, que son probabilidades.
> > head(probs, 3)
> 1   2
> [1,] 0.8282016   0.1717984
> [2,] 0.1288460   0.8711540
> [3,] 0.8830735   0.1169265
>
>   A partir de ella obtengo un vector de 100 elementos que incluye el valor
> máximo de los dos. Utilizo
>
>   Maxprob<-c()
> for(j in 1:100){
>   Maxprob[[j]] <- c(max(probs[j,]))
> }
>
> > head(Maxprob,3)
> [[1]]
> [1] 0.8282016
> [[2]]
> [1] 0.871154
> [[3]]
> [1] 0.8830735
>
> Esto lo hice hace años y calculaba la media de todos los valores de
> Maxprob con mean(Maxprob), y funcionaba, puesto que lo utilicé a menudo
> durante años, con distintas dfs. Pero ahora me dice: Warning message: In
> mean.default(Maxprob) : argument is not numeric or logical: returning NA
>
> He probado a calcular los valores máximos con:
>  Maxprob<-pmax(c(Prob[,1:2])) y hace bien la media con mean(Maxprob), pero
> el resto del código sale distinto que antes. De hecho, ya no sale bien.
> Llevo un buen rato dándole vueltas pero no doy con la clave.
> A ver si alguien pudiera darme aunque sea una pista.
>
> Gracias, como siempre,
> Manuel
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] eliminar filas que cumplen cierta condición.

[Jorge I Velez]

Manuel,

Mira

?subset
?grep
?gsub

En casos como este, no está de más hacer

names(with(data, table(Clase)))

y verificar si existen espacios o no.

Felicidades,
Jorge.-




El El jue, 25 de mar. de 2021 a la(s) 12:12 p. m., Manuel Mendoza <
mmend...@fulbrightmail.org> escribió:

> Buenas tardes, hay algo que me tiene desconcertado (prueba de mi
> desconocimiento):
> ¿cómo puede obtener estas dos cosas de una misma df con 1484 filas?
>
> > table(data$Clase)
>  CYT  ERL  EXC  ME1  ME2  ME3  MIT  NUC  POX  VAC
>  4635   35  44  51163   244  429 20 30
>
> > table(data$Clase == "ERL")
> FALSE
>  1484
>
> Llevo un rato tratando eliminar esas 5 muestras con Clase=ERL usando:
>
> data <- data[data$Clase != "ERL",] o también,   data <-  data$Clase !=
> "ERL"
> y no hacen nada. hice entonces la segunda tabla y vi que ¿no están?
>
> Esto otro:
> > table(data$Clase != "ERL")
> me da:
> TRUE
> 1484
> Claro,
>
> Y probé con filter(), pero tampoco las identificaba. Imagino que estoy
> asumiendo algo que no es.
> 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
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Cambiar datos de un DF

[Jorge I Velez]

Mira la ayuda de %in%  Quizás es lo que buscas.   —JIV

El El dom, 13 de sep. de 2020 a la(s) 7:38 a. m., Samura . <
tontit...@hotmail.com> escribió:

> Hola a tod@s
>
>
>
> ¿Cómo puedo cambiar varios datos por otros en un data frame?
>
>
>
> Por ejemplo
>
>
>
> col1<-c("uno","dos",3,4,"cinco",6,"siete",8,9,"diez")
>
> col2<-c(1,2,3,4,5,6,7,8,9,10)
>
>
>
> df<-data.frame(col1,col2)
>
> df$col1<-as.character(df$col1)
>
> df
>
>
>
> col1 col2
>
> uno1
>
> dos2
>
>   33
>
>   44
>
>   cinco5
>
>   66
>
>   siete7
>
>   88
>
>   99
>
> diez   10
>
>
>
> Ahora quiero cambiar por ejemplo los números por letras (o cualquier otro
> dato en cualquier otro ejemplo)
>
>
>
> df$col1[df$col1==3]<-"tres"
>
> df$col1[df$col1==4]<-"cuatro"
>
> df$col1[df$col1==6]<-"seis"
>
> df$col1[df$col1==8]<-"ocho"
>
> df$col1[df$col1==9]<-"nueve"
>
>
>
> ¿No hay otra forma de ponerlo para no ir uno a uno?
>
>
>
> he probado con
>
>
>
> df$col1[df$col1==c(3,4,6)]<-c("tres","cuatro","seis")
>
>
>
> pero nada, pq creo que tendría q ponerlos todos, solo quiero poner los que
> quiero cambiar.
>
>
>
> 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
>
> --
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Cálculo - intervalo de confianza - modelo nls - predict

[Jorge I Velez]

Marcelino dio en el blanco!
@María:  Para más detalles sobre cómo realizar predicción en modelos nls
desde base, dale una mirada a la ayuda de ?predict.nls
Felicidades,
Jorge.-



On Tue, Sep 1, 2020 at 7:43 AM Marcelino de la Cruz Rot <
marcelino.delac...@urjc.es> wrote:

> Hola:
>
> Igual es porque has escrito:
>
> predict(model,data.frame(t = new_juliano),interval =
> "conficende",level = 0.95)
>
>
> en vez de:
>
> predict(model,data.frame(t = new_juliano),interval =
> "confidence",level = 0.95)
>
>
> Un saludo,
>
> Marcelino
>
> El 01/09/2020 a las 12:57, María Ángeles Onieva escribió:
> > Buenas tardes,
> >
> > Quisiera obtener el intervalo de confianza (y también intervalos de
> > predicción) para los valores predichos en un modelo nls.
> > ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el
> valor
> > listado -además de en el gráfico-) o por bootstrap?
> >
> > Os copio el código del ajuste del modelo y predicción para los 3 días
> > siguientes:
> >
> > *#Ajuste del modelo*
> >
> > model = nls(formula = N~K*exp(-log(K/N0)*exp(-a*(t-t0))),
> >  data = datos,
> >  start = list(K=300, a = 0.25))
> >
> > *#Predicción para 3 días*
> >
> >   new_juliano =
> >
> c(juliano,juliano[(length(juliano))]+1,juliano[(length(juliano))]+2,juliano[(length(juliano))]+3)
> >casos_predichos = predict(model,data.frame(t = new_juliano),interval =
> > "conficende",level = 0.95)
> >
> > Teóricamente debería devolver los intervalos con esto último, sin
> embargo,
> > no los obtengo.
> >
> > Muchas gracias de antemano.
> > Un cordial saludo,
> >
>
> --
> Marcelino de la Cruz Rot
> Depto. de Biología y Geología
> Física y Química Inorgánica
> Universidad Rey Juan Carlos
> Móstoles España
>
> ___
> 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

Re: [R-es] State-of-the-art NLP models from R

[Jorge I Velez]

Diego,
Existe una razón (de peso) para no hacer este análisis directamente con
Python?
Saludos,
Jorge.-

El El lun, 3 de ago. de 2020 a la(s) 11:58 a. m., Diego Martín <
ako.siste...@gmail.com> escribió:

> Estimado Carlos:
>
>  Hago lo que me dices y creo que poniéndote el
> resultado seré más expresivo y claro, que tratando de explicarlo.
>
> (r-reticulate) MacBook-Pro-de-Diego:~ Diego$ python
> Could not find platform independent libraries 
> Could not find platform dependent libraries 
> Consider setting $PYTHONHOME to [:]
> Python path configuration:
>   PYTHONHOME = (not set)
>   PYTHONPATH = (not set)
>   program name = '/anaconda3/bin/python'
>   isolated = 0
>   environment = 1
>   user site = 1
>   import site = 1
>   sys._base_executable = '/anaconda3/bin/python'
>   sys.base_prefix =
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol'
>   sys.base_exec_prefix =
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol'
>   sys.executable = '/anaconda3/bin/python'
>   sys.prefix =
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol'
>   sys.exec_prefix =
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol'
>   sys.path = [
>
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol/lib/
> *python38.zip*',
>
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol/lib/
> *python3.8*',
>
>
> '/opt/concourse/worker/volumes/live/47090e1c-dc82-4ab8-7ab0-318c9b6b0520/volume/python_1593706456233/_h_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehol/lib/lib-dynload',
>   ]
> Fatal Python error: init_fs_encoding: failed to get the Python codec of the
> filesystem encoding
> Python runtime state: core initialized
> ModuleNotFoundError: No module named 'encodings'
>
> Current thread 0x7fff72602000 (most recent call first):
> 
> (r-reticulate) MacBook-Pro-de-Diego:~ Diego$
>
> Obviamente no está bien creado el entorno. Además he destacado del
> resultado las referencias python3.8 que he visto, y que no pueden hallarse
> porque yo tengo instalado python3.6.5, aunque sabía que estábamos ya en la
> versión 3.8, pero no me había animado a instalarla porque no sé cómo
> hacerla y tenía muchos otros asuntos entre manos, entre ellos un estudio de
> NLP.
>
> Ahora, visto esto, no sé qué hacer. Creo que voy a actualizar la versión de
> python y también la de anaconda-navigator, que acabo de ver que está
> desactualizada y, luego, voy a empezar de nuevo. Permitidme que os cuente
> luego.
>
>Muchísimas gracias. Atentamente.
>
>Diego Martín Oliva.
>
> El lun., 3 ago. 2020 a las 16:08, Carlos Ortega ( >)
> escribió:
>
> > Hola Diego,
> >
> > Prueba a hacer otra cosa.
> >
> >- Abre una consola y activa ese environment que has creado
> >(r-reticulate)
> >- Y una vez activado escribe "python". Entrarás a  la consola de
> >"python".
> >- Ahí, escribe "import transformers"
> >- Si no te devuelve error, es que en el entorno está bien instalado
> >esa librería y por tanto el problema es de acceso desde "R".
> >
> > Gracias,
> > Carlos Ortega
> > www.qualityexcellence.es
> >
> > El lun., 3 ago. 2020 a las 13:53, Diego Martín ( >)
> > escribió:
> >
> >> Estimados compañeros:
> >>
> >>   Muchísimas gracias a todos por
> >> responderme. Especialmente a Carlos Ortega que, como siempre, me ha
> vuelto
> >> a enseñar el camino.
> >>
> >>   El asunto era exactamente como
> >> decía Carlos; habiendo yo aplicado la siguiente solución:
> >>
> >> library(reticulate)
> >>
> >> conda_create(envname = "r-reticulate",
> >>  

Re: [R-es] Mejor paquete machine learning

[Jorge I Velez]

Hola Jesús,
No sé si mejor o peor, pero caret [1,2] funciona bastante bien.
Espero sea de utilidad.
Saludos,
Jorge.-

[1] https://topepo.github.io/caret/
[2] https://www.machinelearningplus.com/machine-learning/caret-package/



On Fri, Jul 31, 2020 at 2:27 AM Jesús Para Fernández <
j.para.fernan...@hotmail.com> wrote:

> Buenas
>
> He visto el paquete mlr3, pero se que hay algun otro del estilo.
>
> Cuál usais vosotros y por qué?
>
> Un saludo
>
> Get Outlook for Android
>
> [[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

Re: [R-es] Actualizar la versión de glue

[Jorge I Velez]

Lissette,

La forma más fácil es verificar qué versión tienes es haciendo

R>  packageVersion("glue")
## [1] '1.4.1'

Para actualizarla a la _última_ versión, debes hacer

R> install.packages('glue')

De esta forma siempre vas a instalar la versión más reciente disponible en
https://cran.r-project.org/web/packages/glue/index.html

Saludos,
Jorge.-

On Tue, Jul 21, 2020 at 6:15 PM Lissette Fuentes Lorca <
lissette...@gmail.com> wrote:

> Hola! si, pensé que era algo referido a la versión de glue que tengo.
> Cómo puedo saber qué versión tengo instalada?
> cómo puedo actualizar la misma?
>
> Tal como señalan, al parecer tengo la  1.3.1 y necesito la 1.3.2 .
> Muchas 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
>

[[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 al instalar glue

[Jorge I Velez]

Hola Lissette,
Que versión de glue tienes instalada?
Saludos,
Jorge.-

El El mar, 21 de jul. de 2020 a la(s) 3:07 p. m., Lissette Fuentes Lorca <
lissette...@gmail.com> escribió:

> Hola,
> Estoy haciendo un análisis de datos de Twitter y cuando intento correr esta
> línea de código:
>
> timelines %>%
>   dplyr::filter(created_at > "2020-01-01") %>%
>   dplyr::group_by(screen_name) %>%
>   ts_plot("days", trim = 1L) +
>   ggplot2::geom_point() +
>   ggplot2::labs(
> title = "Tuits publicados por cada cuenta"
>   )
> Me da el siguiente error
>
> Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck
> = vI[[j]]) :
>   namespace ‘glue’ 1.3.1 is already loaded, but >= 1.3.2 is required
>
> Instalé y corrí glue, no sé si tenga que ver con la versión de R que tengo.
>
> Alguien tiene idea cuál es el problema?
> Muchas 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
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Paquetes de R para análisis de encuestas y minería de texto

[Jorge I Velez]

Hola Manuel,

Mira

http://r-survey.r-forge.r-project.org/survey/index.html para la primera
pregunta y

   - https://cran.r-project.org/web/views/NaturalLanguageProcessing.html
   - https://www.tidytextmining.com/
   - https://hub.packtpub.com/9-useful-r-packages-for-nlp-text-mining/

para la segunda.

Saludos,
Jorge.-



On Sat, Jul 18, 2020 at 4:46 PM Manuel Spínola  wrote:

> Estimados miembros del grupo,
>
> Quería preguntarles por sugerencias de paquetes para:
>
> 1. Análisis de encuestas o cuestionarios
> 2. Minería de texto
>
> Muchas gracias y slaudos
>
> Manuel Spínola
>
>
> --
> *Manuel Spínola, Ph.D.*
> Instituto Internacional en Conservación y Manejo de Vida Silvestre
> Universidad Nacional
> Apartado 1350-3000
> Heredia
> COSTA RICA
> mspin...@una.cr 
> mspinol...@gmail.com
> Teléfono: (506) 8706 - 4662
> Personal website: Lobito de río <
> https://sites.google.com/site/lobitoderio/>
> Institutional website: ICOMVIS 
>
> [[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

Re: [R-es] Tomar los elementos de una colMeans sin los nombres

[Jorge I Velez]

Quizás

as.numeric(distmeans)

es lo que buscas?

Saludos,
Jorge.-

El El dom, 17 de may. de 2020 a la(s) 7:35 a. m., Manuel Mendoza <
mmend...@fulbrightmail.org> escribió:

> Muy buenas tardes, llevo un rato con una chorrada que casi me da vergüenza
> preguntar, pero no doy con la tecla. Creo un vector con las medias de
> algunas columnas de una df:
>
> distmeans<-colMeans(df[,1:6])
>
> Son 6 elementos, pero siguen manteniendo el nombre de las variables
> originales, incluso cuando creo otros vectores a partir de ellos:
>
> rcp26<-c(distmeans[1], distmeans[4])
> rcp45<-c(distmeans[2], distmeans[5])
> rcp85<-c(distmeans[3], distmeans[6])
>
> Hago un cbind con ellos pero sigue incluyendo los nombres originales de las
> variables y no sale lo que busco:
>
> distdf<-data.frame(cbind(year,rcp45,rcp45,rcp85))
>
> ¡No consigo librarme de ellos!
>
> Estoy seguro de que se puede hacer todo esto de forma sencilla, pero no me
> sale :(
>
> Gracias,
> Manuel
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Mantener nombres de variables

[Jorge I Velez]

Hola Samura,
Realmente no. R no está cambiando tus variables.
El nombre “LeagueN” corresponde al nivel N de la variable League.
Seguramente debes tener otro nivel — como no aparece, éste es el nivel de
referencia.  Lo mismo aplica para las demás variables dentro del modelo de
RLM. Eso si, tus variables son las mismas, pero para efectos de
interpretación, debes tener cuidado con qué nivel tienes como referencia y
para cuál estás estimando el efecto.
Por favor no recodifiques con dummies — R lo hace por tí y es más
eficiente.
Saludos,
Jorge.-

El El dom, 26 de abr. de 2020 a la(s) 3:17 p. m., Samura . <
tontit...@hotmail.com> escribió:

> Hola
> tengo la siguiente duda,
> ¿Cómo hago para que R no me cambie el nombre de las variables en este
> ejemplo?
>
>
> library(ISLR)
> data(Hitters)
> Hitters<-na.omit(Hitters)
>
> modelo=lm(Salary~.,data=Hitters)
> modelo
>
> las variables
> League
> Division
> NewLeague
>
> ahora se llaman
> LeagueN
> DivisionW
> NewLeagueN
>
> quiero pasar esas variables  a 0-1,
> he probado con recode o con fastDummies , pero al final acaba cambiándome
> las variables.
>
>
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Conversión de variable categóricas a Númericas.

[Jorge I Velez]

Edison,
Necesitas
as.numeric(as.character(tusnumeros))
Saludos,
Jorge.-

El El dom, 19 de abr. de 2020 a la(s) 8:53 p. m., Edison Andrés Rincon
Cardoso  escribió:

> Cordial saludo.
>
> Estoy teniendo un problema al convertir una columna categórica a numérica.
> Espero que me puedan colaborar.
> Esta variable está en números escritos, es decir, "uno", "dos" etc. Ya los
> cambie por números (que siguen siendo factores) mediante fct_recode y quise
> de ahí convertirla mediante as.numeric pensando que R iba a identificar los
> números y no iba a alterar el orden, pero no fue así. El siete, que era la
> posición del 1 (en orden), apareció emplazándolo en el resultado.
> También quisiera saber si hay alguna manera de que se empiece desde cero y
> no desde 1.
>
> Espero haber sido claro en la explicación
>
> Agradezco mucho su colaboración.
>
>
> --
>
>
>*   Edison Rincón *
>Tel: 3154574672
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Gráficos "handy-style"

[Jorge I Velez]

Hola Héctor,
Mira https://cran.r-project.org/web/packages/xkcd/index.html
Jorge.-


On Wed, Feb 26, 2020 at 11:53 AM "Hector Gómez Fuerte" 
wrote:

> Buenas tardes,
> ¿conoce alguien algún paquete del R que permita crear gráficos
> estadísticos con un aspecto como si se hubieran dibujado a mano? Busco algo
> similar al paquete de Python  matplotlib. Algunos lo llaman "xkcd comic
> style" (y hay un paquete en R para ello, el xkcd, pero no es exactamente
> igual ni tiene las posibilidades del matplotlib).
> Muchas gracias y saludos.
> ___
> 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

Re: [R-es] Creación de vector con una secuencia de números enteros.

[Jorge I Velez]

Basta con

mi_vector <- 11:30

Si requieres algo mas general, entonces

?seq

es lo que buscas.

Saludos,
Jorge.-

El El dom, 28 de jul. de 2019 a las 11:17 a. m., Jesús Para Fernández <
j.para.fernan...@hotmail.com> escribió:

> Mas fácil que todo eso
>
> mi_vector <- c(11:30)
>
> c es lo que crea el vector y para crear una secuencia usas :
>
> Obtener Outlook para Android
>
> 
> From: R-help-es  on behalf of Daniel
> García Abiétar 
> Sent: Sunday, July 28, 2019 6:13:18 PM
> To: r-help-es@r-project.org 
> Subject: [R-es] Creación de vector con una secuencia de números enteros.
>
> Buenas.
>
> Soy un chaval novato en el uso de R. Ando haciendo un curso online, y he
> llegado a una parte en la que me he atascado. La instrucción es la
> siguiente:
>
> crea una variable llamada
> | 'mi_vector' que contenga un vector con los números enteros del 11 al
> 30. Recuerda
> | que puedes usar el operador secuencia ':'
>
> A lo cual yo he introducido el siguiente comando:
>
> 'mi_vector' <- vector("numeric", length = 11:30)
>
> Pero me da el siguiente error:
>
> Error in vector("numeric", 11:30) : invalid 'length' argument
>
> ¿Alguien podría echarme una mano para resolver por qué me sucede esto?
> ¿Dónde está mi error?
>
> Muchas gracias, 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
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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 un rango de varias variables explicativas a un modelo

[Jorge I Velez]

Hola Rolando,

Quizás subset() sea tu amigo en este caso.  Intenta:

## selecciona la variable dependiente y las variables independientes
d0 <- subset(wageszm14, select = c(lwage, SEXO:marr, pot23:pot317)

## ajusta el modelo
fit <- qregspiv(lwage ~ ., shpfile = zm15, tau = 0.5, nboot = 70, data = d0)

## resultados
summary(fit)

Saludos,
Jorge.-



On Sun, Jun 2, 2019 at 6:00 PM Rolando Valdez  wrote:

> Hola,
>
> Quiero especificar una ecuación con varias variables explicativas de una
> manera eficiente sin necesidad de escribir todas y cada una. Tengo un
> conjunto de variables (junto con otras) dentro de una base de datos que se
> llaman pot23 pot311 pot312 pot 316 pot317... pot80. No
> necesariamente están secuenciadas. Quisiera saber cómo indicar que incluya
> todas las variables de pot23 a pot80 en una ecuación.
>
> He intentado lo siguiente, pero no funciona:
>
> > pots <- paste("pot",23:321, sep="")
> > eqreg2 <- lwage~SEXO+EDAD+HLENGUA+ESCOACUM+marr+wageszm14[,pots]
> > fit <- qregspiv(eqreg2, shpfile = zm15, tau = 0.5, nboot = 70, data =
> wageszm14)
> Error: Can't find columns `pot24`, `pot25`, `pot26`, `pot27`, `pot28`, ...
> (and 273 more) in `.data`.
>
> De igual forma, después estaría interesado en obtener el logaritmo de todas
> esas variables pot~
>
> Gracias de antemano por cualquier tipo de ayuda.
>
> --
> Rol~
>
> [[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

Re: [R-es] Crear nuevo dataframe y eliminar duplicados

[Jorge I Velez]

Hola Jesús,
Dale una mirada a los ejemplos en ?duplicated.
Saludos,
Jorge.-



On Fri, Apr 26, 2019 at 12:21 PM xyg...@gmail.com  wrote:

> Buenas tardes. Tengo un dataframe (df) con varias columnas (A, B, C, D, E,
> F, G, H) de las que quiero quedarme solo con algunas (digamos C, D, E). Lo
> consigo mediante select(df, ˝C˝, ˝D˝, ˝E˝). Hasta ahí todo correcto.
>
> Ahora no consigo ver como convertir el resultado en un nuevo dataframe con
> los datos que me interesan y como eliminar sus datos duplicados.
>
> Espero haberme explicado bien.
>
> Un saludo y gracias anticipadas
> Jesús
> ___
> 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

Re: [R-es] PROBLEMAS NOMBRES DE COLUMNAS CON ESPACIOS CONEXION R-SQL

[Jorge I Velez]

Hola Ana,
Lo que ocurre es que “Corteoptimo” es diferente a “Corte optimo”. Lo mismo
ocurre con “Cortediairo”.  Por eso no reconoce que existen.
Saludos,
Jorge.-

El El mar, 16 de abr. de 2019 a las 8:56 a. m., Ana Jimenez Rebollo <
anaji...@gmail.com> escribió:

> Buenas tardes,
>
> Estoy tratando de realizar un update en SQL desde R:
> sqlUpdate(conexion1, data.frame(AUXILIAR), tablename = "AUXILIAR")
>
> y me devuelve el siguiente error:
> *Error in sqlUpdate(conexion1, data.frame(AUXILIAR), tablename =
> "AUXILIAR",  : *
> *  data frame column(s) Corteoptimo Cortediario not in database table*
>
> El problema es que sí existen esas columnas en SQL pero parece que R está
> almacenando internamente los nombres sin espacios ya que si ejecuto la
> función names en R para el data.frame AUXILIAR (que es con el que quiero
> actualizar en R):
> names(AUXILIAR)
>  [1] "Date""Time""Seasson"
>"TimeH"   "Day" "Holiday"
>
>  [7] "Corte optimo""Corte diario"
>
> Y si obtengo los nombres de la tabla SQL:
> sqlColumns(conexion1, "AUXILIAR")
> TABLE_CAT TABLE_SCHEM TABLE_NAME COLUMN_NAME
> OS MKOnline dbo AUXILIAR Date
> OS MKOnline dbo AUXILIAR Time
> OS MKOnline dbo AUXILIAR Seasson
> OS MKOnline dbo AUXILIAR TimeH
> OS MKOnline dbo AUXILIAR Day
> OS MKOnline dbo AUXILIAR Holiday
> OS MKOnline dbo AUXILIAR Corte optimo
> OS MKOnline dbo AUXILIAR Corte diario
>
>
> Los nombres son idénticos por lo que intuyo que el problema viene en el
> nombre almacenado internamente. He probado a renombrar las columnas del
> data.frame sustituyendo los espacios por puntos y por guiones bajos pero me
> devuelve el mismo error.
> PD: El data.frame con el que quiero actualizar en SQL y la tabla llevan el
> mismo nombre, siento si eso lleva a algún tipo de confusión.
> ¿Alguna sugerencia?
>
> Muchas gracias de antemano.
>
> Saludos,
> AJR.
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Consulta eliminar columna autómatica

[Jorge I Velez]

Quizás

write.table(.., row.names = FALSE)

?

Saludos,
Jorge.-

El El dom, 3 de mar. de 2019 a las 11:18 a. m., Andrés Hirigoyen <
andreshirigo...@gmail.com> escribió:

> Buenas tardes, estoy intentando guardar una salida como archivo .xlsx, pero
> necesito que la primer columna generada (que indica el número de registros)
> en todo data.frame no se guarde. En el ejemplo es tengo 8 registros
>
>y grupo edad
> 1 28 1   30
> 2 24 1   45
> 3 30 1   44
> 4 24 1   41
> 5 22 2   42
> 6 31 2   43
> 7 24 2   42
> 8 28 2   44
>
> ¿Ideas?
> Muchas 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
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] crear un vector con las categorías

[Jorge I Velez]

Estimado Manuel,

Debes definir ecsta como factor usando, por ejemplo,

factor(ecsta, levels = ...)

antes de ajustar el modelo.

Dale una mirada a

?factor

para má detalles.

Saludos,
Jorge.-

El El lun, 18 de feb. de 2019 a las 1:03 p. m., Manuel Mendoza <
mmend...@mncn.csic.es> escribió:

>
> Buenas tardes, tengo un loop que hace un árbol de clasificación cada
> vez y va creando un vector con una predicción que hace. Son 4
> categorías (pongamos a, b, c y d), pero en vez de ir añadiendo la
> categoría predicha me añade al vector el nº (del 1 al 4) al que
> corresponde esa categoría. Supongo que se puede hacer que añada la
> categoría (la letra), pero no sé cómo.
>
> preds <- {}
>
> for (i in 1:nrow(data)) {
>
>training <- data[-i, ]
>
>fitrp <- rpart(ecsta ~ .,data=training, cp=0)
>
>Pred <- predict(fitrp,data[i,], type="class")
>
>preds[i] <- Pred
>
> }
>
> Gracias como siempre,
> Manuel
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> .
>
> --
> Dr Manuel Mendoza
> Department of Biogeography and Global Change
> National Museum of Natural Science (MNCN)
> Spanish Scientific Council (CSIC)
> C/ Serrano 115bis, 28006 MADRID
> Spain
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Donde obtener tablas de datos ejemplo

[Jorge I Velez]

?cor

El El mié, 12 de dic. de 2018 a las 3:51 p. m., jorge della gaspera <
jara...@gmail.com> escribió:

> Buenas tardes, donde puedo encontrar tablas de datos para trabajar con R
> las funciones de correlación de variables.
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Variables: non-numeric argument to 'pairs'

[Jorge I Velez]

> : num  1 1 1 2 3 1 1 1 1 2 ...
>  $ Situacion actual
>: chr  "Contestare basandome en una relacion
> anterior." "Contestare basandome en una relacion anterior." "Contestare
> basandome en una relacion anterior." "Con pareja." ...
>  $ De mi pareja recibo todo el apoyo moral que necesito,
> : num  4 3 5 1 5 4 4 1 4 2 ...
>  $ Mi pareja me trata con respeto,
> : num  5 4 5 5 5 4 5 3 4 5 ...
>  $ Si surge algun conflicto se que podemos solucionarlo mi pareja y yo,
>: num  4 3 3 5 5 4 4 3 3 5 ...
>  $ Por pensar en mi pareja me distraigo de mis actividades cotidianas,
> : num  1 3 3 5 2 1 5 5 2 1 ...
>  $ Puedo hablar libremente de lo que sea con mi pareja,
>: num  3 3 5 5 5 3 5 5 4 5 ...
>  $ Desde que estoy con mi pareja me he conocido mas como persona,
>: num  3 5 5 5 5 4 5 5 3 2 ...
>  $ Mi pareja me aprueba,
> : num  5 3 5 3 5 4 5 5 4 3 ...
>  $ Mi pareja me impulsa para cumplir mis metas,
>: num  5 4 5 5 5 5 4 4 3 2 ...
>  $ Las expectativas que ambos tenemos sobre una relacion son
> complementarias, : num  5 4 5 1 5 4 3 3 2 3 ...
>  $ La relacion con mi pareja me ha transformado como persona,
>: num  4 4 5 5 5 3 5 5 2 1 ...
>  $ Sin mi pareja me sentiria perdido,
>: num  2 2 4 1 2 1 1 2 1 1 ...
>  $ Si mi pareja no esta cerca de mi fisicamente me siento intranquilo/a,
> : num  1 2 3 1 1 2 1 1 1 1 ...
>  $ Mi pareja y yo realizamos varias actividades en conjunto,
> : num  3 4 5 1 5 4 2 5 2 2 ...
>  $ La relacion con mi pareja se ha vuelto rutinaria,
> : num  3 3 1 1 5 2 1 1 3 4 ...
>  $ He dejado de hacer cosas por estar con mi pareja,
> : num  2 2 4 5 1 1 4 4 2 1 ...
>  $ Me siento vulnerable ante mi pareja,
>: num  1 3 1 1 1 1 3 1 1 1 ...
>  $ Mi pareja es una de las personas mas importantes en mi vida,
>: num  5 3 5 5 5 4 5 4 4 3 ...
>   [list output truncated]
>  - attr(*, "na.action")= 'omit' Named int  44 45 47 49 58 62 66 69 72 82
> ...
>   ..- attr(*, "names")= chr  "44" "45" "47" "49" ...
> >
>
> El dom., 25 nov. 2018 a las 16:37, Jorge I Velez (<
> jorgeivanve...@gmail.com>) escribió:
>
>> Kelly,
>> Cual es el resultado de
>> str(Base)
>> ?
>> Saludos,
>> Jorge.-
>>
>> El El dom, 25 de nov. de 2018 a las 4:00 a. m., Kelly Johanna Monroy
>> Malagon  escribió:
>>
>>> Hola a todos
>>>
>>> Estoy cargando una base de datos a R llamada "Base", desde excel. A la
>>> hora
>>> de utilizar funciones como:
>>>
>>> >pairs(Base)
>>> Error in pairs.default(Base) : non-numeric argument to 'pairs'
>>>
>>> Una forma de arreglar este problema es utilizando la funcion
>>> as.numeric(),
>>> pero me toca hacerlo variable por variable:
>>>
>>> >Base$Variable1<-as.numeric( Base$Variable1)
>>> >Base$Variable2<-as.numeric( Base$Variable2)
>>>
>>> etcétera, el problema es que tengo muchas variables y no puedo hacer esto
>>> una por una.
>>>
>>> Por esto vengo aquí a buscar una idea o asesoramiento.
>>>
>>> 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
>>>
>> --
>> Sent from my phone. Please excuse my brevity and misspelling.
>>
> --
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Variables: non-numeric argument to 'pairs'

[Jorge I Velez]

Kelly,
Cual es el resultado de
str(Base)
?
Saludos,
Jorge.-

El El dom, 25 de nov. de 2018 a las 4:00 a. m., Kelly Johanna Monroy
Malagon  escribió:

> Hola a todos
>
> Estoy cargando una base de datos a R llamada "Base", desde excel. A la hora
> de utilizar funciones como:
>
> >pairs(Base)
> Error in pairs.default(Base) : non-numeric argument to 'pairs'
>
> Una forma de arreglar este problema es utilizando la funcion as.numeric(),
> pero me toca hacerlo variable por variable:
>
> >Base$Variable1<-as.numeric( Base$Variable1)
> >Base$Variable2<-as.numeric( Base$Variable2)
>
> etcétera, el problema es que tengo muchas variables y no puedo hacer esto
> una por una.
>
> Por esto vengo aquí a buscar una idea o asesoramiento.
>
> 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
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Visualizar grandes volumenes de datos

[Jorge I Velez]

Jesús,
Intenta con hexbin (
https://cran.r-project.org/web/packages/hexbin/index.html)
Saludos,
Jorge.-


On Thu, Nov 8, 2018 at 7:13 AM Jesús Para Fernández <
j.para.fernan...@hotmail.com> wrote:

> Buenas,
>
> He probado a intentar graficar mas de 700.000 puntos y es una locura. No
> doy con una libreria que consiga graficar grandes volumenes de datos.
> ¿Existe alguna en R?
>
> Gracias
> Jesús
>
> [[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

Re: [R-es] Fwd: Paletas de colores para datos cualitativos

[Jorge I Velez]

Dale una mirada a TraMineR    —JIV

On Wed, Oct 31, 2018 at 7:20 AM Carla Melloni  wrote:

>
>
>
>
>
> Hola, mi nombre es Carla y trabajo en Estadística. Soy nueva en R y estoy
> haciendo un curso para principiantes. Tengo unos datos de trayectorias de
> pacientes de un hospital -categorías cualitativas- y quisiera graficar cada
> trayectoria empleando paletas de colores, para prescindir del texto.
> Adjunto imagen de trabajo realizado en R para un estudio longitudinal, que
> se acerca a lo que quisiera hacer. Gracias!
> ___
> 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

Re: [R-es] Librería para simulación de datos GWAS en estudios de casos y controles

[Jorge I Velez]

Fernando,

Efectivamente hay varias opciones.  Dos alternativas son

https://rdrr.io/bioc/hierGWAS/man/simGWAS.html
https://arxiv.org/pdf/1710.01236.pdf
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5605948/
https://cnsgenomics.com/software/gcta/#Overview (no propiamente en R, pero
tiene interfaz)

Espero sea de utilidad.

Saludos,
Jorge Velez.-


On Mon, Oct 29, 2018 at 2:19 PM Fernando Sanchez via R-help-es <
r-help-es@r-project.org> wrote:

> Hola a todos,
> Estoy buscando alguna librería de R que sea capaz de generar datos
> sintéticos para estudios GWAS (Genome-Wide association studies) de casos y
> controles. Sé que hay unas cuantas, pero me gustaría tener información de
> primera mano de alguien que haya usado alguna y me la pueda recomendar.
> un saludo y muchas gracias,
> Fernando
>
>
> [[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

Re: [R-es] Instalar paquetes antiguos R en Windows 10

[Jorge I Velez]

Claro.  Lo que ocurre es que

C:\Users\Victoria MartÌn

es diferente a

C:/Users/Victoria Martmn/

No sé qué debes hacer para corregir el problema, pero no debe ser difícil.

Saludos,
Jorge Velez.-



El El dom, 21 de oct. de 2018 a las 10:55 a. m., JOSE MARTIN AREVALO <
martin_j...@gva.es> escribió:

> Hola a todos y gracias por responder.
>
> He seguido la sugerencia de ver si era un problema de permisos, y creo que
> no es el caso. Me da la impresión que R busca archivos en un directorio que
> no existe y, además, no creo haberle configurado. Os pongo la salida que me
> da cuando bajo este archivo del repositorio (es la misma que cuando lo hago
> desde el directorio):
>
> > install.packages("devtools")
> Installing package into ëC:/Users/Victoria
> MartÌn/Documents/R/win-library/3.5í
> (as ëlibí is unspecified)
>
>   There is a binary version available but the source version is later:
>  binary source needs_compilation
> devtools 1.13.6  2.0.0 FALSE
>
> installing the source package ëdevtoolsí
>
> probando la URL '
> https://cran.cnr.berkeley.edu/src/contrib/devtools_2.0.0.tar.gz'
> Content type 'application/x-gzip' length 388238 bytes (379 KB)
> downloaded 379 KB
>
> * installing *source* package 'devtools' ...
> ** package 'devtools' successfully unpacked and MD5 sums checked
> Warning in file(file, if (append) "a" else "w") :
>   cannot open file 'C:/Users/Victoria
> Martmn/Documents/R/win-library/3.5/devtools/DESCRIPTION': No such file or
> directory
> Error in file(file, if (append) "a" else "w") :
>   no se puede abrir la conexiÛn
> ERROR: installing package DESCRIPTION failed for package 'devtools'
> * removing 'C:/Users/Victoria MartÌn/Documents/R/win-library/3.5/devtools'
> In R CMD INSTALL
>
> The downloaded source packages are in
> ëC:\Users\Victoria
> MartÌn\AppData\Local\Temp\Rtmp8SKA07\downloaded_packagesí
> Warning message:
> In install.packages("devtools") :
>   installation of package ëdevtoolsí had non-zero exit status
>
>
> Esto me desconcierta bastante: "Warning in file(file, if (append) "a" else
> "w") :  cannot open file 'C:/Users/Victoria
> Martmn/Documents/R/win-library/3.5/devtools/DESCRIPTION': No such file or
> directory” porque ese directorio no existe y no sé de dónde lo ha sacado R.
> Como veréis, he llegado incluso a instalar la versión 3.5.1, y me sigue
> pasando lo mismo. Cuando he intentado cambiar el directorio, no lo hace.
> Al aplicar getwd(), la salida que da es: C:/Users/Victoria
> MartÌn/Documents/R/win-library/3.5, pero creo que no me deja cambiar el
> directorio con “setwd()”.
>
> Espero que con esto os aporte algo más de información, porque yo llevo una
> semana dándole al tema y no hay manera.
>
> Muchas gracias.
>
> Pepe Martín
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Construcción de archivo de texto

[Jorge I Velez]

Gracias, Marcelino.  Alguna idea sobre cómo almacenarlo en un archivo de
texto?
Saludos,
Jorge.-



On Tue, Jul 10, 2018 at 5:13 AM Marcelino de la Cruz Rot <
marcelino.delac...@urjc.es> wrote:

> Hola. A ver esto, qué tal:
>
>
> for (i in 1:length(d)){
>   cat(paste("Pedigree: ",unique(d[[i]]$ped),"\n"))
>   print(d[[i]],row.names=F)
>
> }
>
> Saludos,
>
> Marcelino
>
>
>
> El 10/07/2018 a las 11:31, Jorge I Velez escribió:
> > Hola a todos,
> >
> > A partir de los siguientes datos:
> >
> > d <- list(`1` = structure(list(ped = c(1L, 1L, 1L, 1L, 1L, 1L, 1L),
> >  id = 1:7, father = c(2L, 0L, 0L, 2L, 2L, 2L, 2L), mother = c(3L,
> >  0L, 0L, 3L, 3L, 3L, 3L), sex = c(2L, 1L, 2L, 2L, 2L, 1L,
> >  2L), affected = c(1L, 2L, 1L, 1L, 2L, 2L, 2L)), row.names = c("1",
> > "2", "3", "4", "5", "6", "7"), class = "data.frame"), `2` =
> structure(list(
> >  ped = c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> >  2L), id = 201:214, father = c(0L, 0L, 0L, 201L, 201L, 201L,
> >  201L, 201L, 0L, 203L, 203L, 209L, 209L, 209L), mother = c(0L,
> >  0L, 0L, 202L, 202L, 202L, 202L, 202L, 0L, 204L, 204L, 208L,
> >  208L, 208L), sex = c(1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L,
> >  1L, 1L, 2L, 1L, 1L), affected = c(1L, NA, 1L, 0L, NA, 1L,
> >  1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L)), row.names = c("42", "43",
> > "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54",
> > "55"), class = "data.frame"))
> > d
> >
> > estoy interesado en obtener un archivo de texto "out.txt" con la
> siguiente
> > estructura:
> >
> >
> > pedigree: 1
> > ped id father mother sex affected
> >1  1  2  3   21
> >1  2  0  0   12
> >1  3  0  0   21
> >1  4  2  3   21
> >1  5  2  3   22
> >1  6  2  3   12
> >1  7  2  3   22
> > pedigree: 2
> > ped  id father mother sex affected
> >2 201  0  0   11
> >2 202  0  0   2   NA
> >2 203  0  0   11
> >2 204201202   20
> >2 205201202   1   NA
> >2 206201202   21
> >2 207201202   21
> >2 208201202   20
> >2 209  0  0   10
> >2 210203204   10
> >2 211203204   10
> >2 212209208   20
> >2 213209208   10
> >2 214209208   11
> >
> >
> > La idea es relativamente simple:  guardar cada data.frame de "d" sin row
> > labels y con el encabezado "pedigree: " seguido del número en la columna
> > "ped".
> >
> > Alguna idea?  Por supuesto esto es una versión muy reducida del problema
> > real.
> >
> > Muchísimas gracias por su ayuda.
> >
> > Saludos cordiales,
> > Jorge.-
> >
> >   [[alternative HTML version deleted]]
> >
> > ___
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> > .
> >
>
> --
> Marcelino de la Cruz Rot
> Depto. de Biología y Geología
> Física y Química Inorgánica
> Universidad Rey Juan Carlos
> Móstoles España
>
>

[[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] Construcción de archivo de texto

[Jorge I Velez]

Hola a todos,

A partir de los siguientes datos:

d <- list(`1` = structure(list(ped = c(1L, 1L, 1L, 1L, 1L, 1L, 1L),
id = 1:7, father = c(2L, 0L, 0L, 2L, 2L, 2L, 2L), mother = c(3L,
0L, 0L, 3L, 3L, 3L, 3L), sex = c(2L, 1L, 2L, 2L, 2L, 1L,
2L), affected = c(1L, 2L, 1L, 1L, 2L, 2L, 2L)), row.names = c("1",
"2", "3", "4", "5", "6", "7"), class = "data.frame"), `2` = structure(list(
ped = c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L), id = 201:214, father = c(0L, 0L, 0L, 201L, 201L, 201L,
201L, 201L, 0L, 203L, 203L, 209L, 209L, 209L), mother = c(0L,
0L, 0L, 202L, 202L, 202L, 202L, 202L, 0L, 204L, 204L, 208L,
208L, 208L), sex = c(1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L,
1L, 1L, 2L, 1L, 1L), affected = c(1L, NA, 1L, 0L, NA, 1L,
1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L)), row.names = c("42", "43",
"44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54",
"55"), class = "data.frame"))
d

estoy interesado en obtener un archivo de texto "out.txt" con la siguiente
estructura:


pedigree: 1
ped id father mother sex affected
  1  1  2  3   21
  1  2  0  0   12
  1  3  0  0   21
  1  4  2  3   21
  1  5  2  3   22
  1  6  2  3   12
  1  7  2  3   22
pedigree: 2
ped  id father mother sex affected
  2 201  0  0   11
  2 202  0  0   2   NA
  2 203  0  0   11
  2 204201202   20
  2 205201202   1   NA
  2 206201202   21
  2 207201202   21
  2 208201202   20
  2 209  0  0   10
  2 210203204   10
  2 211203204   10
  2 212209208   20
  2 213209208   10
  2 214209208   11


La idea es relativamente simple:  guardar cada data.frame de "d" sin row
labels y con el encabezado "pedigree: " seguido del número en la columna
"ped".

Alguna idea?  Por supuesto esto es una versión muy reducida del problema
real.

Muchísimas gracias por su ayuda.

Saludos cordiales,
Jorge.-

[[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 un cmeans

[Jorge I Velez]

Manuel,
Mira
https://www.rdocumentation.org/packages/e1071/versions/1.6-8/topics/cmeans
No debes incluir “i” en cmeans().
Saludos,
Jorge.-

Sent from my phone. Please excuse my brevity and misspeling.


El El mié, 27 de jun. de 2018 a las 9:33 a. m., Manuel Mendoza <
mmend...@mncn.csic.es> escribió:

>
> Pues no lo sé. Es algo del código de RandomForest, supongo.
>
>
>
> Quoting Jesús Para Fernández :
>
> > U es un dataframe?
> >
> > Obtener Outlook para Android
> >
> > 
> > From: R-help-es  on behalf of
> > Manuel Mendoza 
> > Sent: Wednesday, June 27, 2018 12:25:10 PM
> > To: r-help-es@r-project.org
> > Subject: [R-es] error en un cmeans
> >
> >
> > Tuve que dejar, de momento, lo del "loop con matriz que cambia de
> > nombre" por algo más urgente, y me salió otro problema.
> >
> > Hago un cmeans con el paquete e1071;
> >
> >cl<-cmeans(Data,i,20,verbose=F,method="cmeans",m=2)
> >
> >   y me da este error: Error in apply(u, 1, which.max) : dim(X) must
> > have a positive length.
> >
> > Estoy harto de usar esa línea de código, con dfs iguales a la de
> > ahora, por lo que no entiendo por qué falla.
> >
> > Gracias,
> > Manuel
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > .
> >
> >
> > --
> > Dr Manuel Mendoza
> > Department of Biogeography and Global Change
> > National Museum of Natural History (MNCN)
> > Spanish Scientific Council (CSIC)
> > C/ Serrano 115bis, 28006 MADRID
> > Spain
> >
> > ___
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
> --
> Dr Manuel Mendoza
> Department of Biogeography and Global Change
> National Museum of Natural History (MNCN)
> Spanish Scientific Council (CSIC)
> C/ Serrano 115bis, 28006 MADRID
> Spain
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Pasar palabras de una lista a una variable del dataframe

[Jorge I Velez]

Una forma es con %in%.  —JIV

El El mar, 22 de may. de 2018 a las 11:46 a. m., Miriam Alzate <
miriam.alz...@unavarra.es> escribió:

> Buenas tardes,
>
> Tengo una lista de 600 palabras. Quiero saber cuántas de esas palabras
> aparecen en cada observación de mi variable "texto". La variable "texto"
> es de tipo caracter. ¿Cómo lo haríais?
>
> Muchas gracias.
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Crear ROC desde la matriz de confusion

[Jorge I Velez]

Jesús,

La función roc.from.table en el paquete epicalc <
https://cran.r-project.org/src/contrib/Archive/epicalc/> hace lo que
necesitas.

El paquete fue retirado de CRAN, pero puedes acceder al *source* en este
link.  El siguiente paso es descargarlo en tu equipo e instalarlo.

Espero sea de utilidad.

Saludos,
Jorge.-


2017-10-27 4:59 GMT-05:00 Jesús Para Fernández :

> Buenas,
>
> ¿Es posible calcular la ROC desde la matriz de confusion? no entiendo muy
> bien como se crea la matriz de confusion y como pueden salir varios puntos
> desde solo la matriz...
>
> Gracias
> Jesús
>
> [[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

Re: [R-es] Lrtest interpretation y error en BIC

[Jorge I Velez]

Buenas tardes, Miriam

Un valor p significativo, digamos al 5%, indica que el modelo "actual" es
mejor que el modelo "anterior". En tu caso, el modelo 2 parece ser mejor
que el modelo 1.

Por otro lado, desafortunadamente no encuentro razón aparente para que
BIC(...) sea NA a partir de lo que nos envías.  Necesitaríamos mas
información para poder ayudarte.

Espero sea de utilidad.

Saludos,
Jorge.-


2017-10-06 14:56 GMT-05:00 :

> Hola a todos,
>
> Estoy trabajando con un modelo negativo binomial de ceros inflados (ZINB).
>
> ¿Cómo podría interpretar el resultado del lrtest del paquete lmtest?
>
>   #Df  LogLik  Df  Chisq Pr(>Chisq)
> 1  24 -4006.7
> 2   3 -5726.4 -21 3439.4  < 2.2e-16 ***
>
> Es el resultado del codigo:
> >lrtest(m10)
>
> m10 es mi modelo ZINB.
>
> Por otro lado, al calcular el test BIC (Bayesian information criterion) el
> resultado que me da del siguiente comando, es NA. ¿Qué puede dar error?
>
> >BIC(m10)
>
> Muchas gracias,
>
> Miriam
>
> ___
> 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

Re: [R-es] gráfico

[Jorge I Velez]

Jose,

Usa

plot(z, col = as.numeric(pred.class))

Saludos,
Jorge.-


2017-09-03 8:29 GMT-05:00 :

> Estimados
> Quisiéramos que la salida final plot (z) fuera de tres colores, uno para
> cada especie, no encuentro el comando para ello. Apreciaría su ayuda
> Saludos
> Jose
>
>
> rm(list = ls())
> library(MASS)
> data(iris)
> datos <- data.frame(iris[,1:4],clase=as.vector(iris[,5]))
> head(datos)
> especies <- iris[,5] # La columna cinco se refiere a las especies(clase)
> z <- lda(clase~.,datos)
> predict(z,iris[1:4, ])$class #Para clasificar las observaciones
> Tabla <- table(iris[,5],predict(z,iris[,1:4])$class)
> Tabla
> (pred.class<-predict(z)$class) # Para clasificar las observaciones
> plot(z)
>
>
>
> --
> Este mensaje le ha llegado mediante el servicio de correo electronico que
> ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
> Nacional de Salud. La persona que envia este correo asume el compromiso de
> usar el servicio a tales fines y cumplir con las regulaciones establecidas
>
> Infomed: http://www.sld.cu/
>
> ___
> 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

Re: [R-es] Area bajo la curva

[Jorge I Velez]

Buenas noches, Wilmer.

Puedes intentar la aproximación sugerida en
https://stackoverflow.com/questions/4954507/calculate-the-area-under-a-curve


## preparación
x <- c(1,2,3,4,5,6,7,8,9,10)
y <- c(0,2,15,30,50,NaN,NaN,10,2,1)
d <- data.frame(x, y)
d <- d[complete.cases(d),]
id <- order(d$x)

## cálculo
require(zoo)
with(d, sum(diff(x[id])*rollmean(y[id],2)))
## [1] 169.5

Espero sea de utilidad.
Jorge.-


2017-08-31 19:17 GMT-05:00 WILMER CONTRERAS SEPULVEDA 
:

> Buenas noches.
>
> Quisiera saber si existe una libreria o función que me calcule el área bajo
> la curva de dos vectores de números. Ademas esta debe tener en cuenta si
> existen datos faltantes o NaN.
>
> Por ejemplo calcular el área bajo la curva del vector en el eje x:
>
> x<- c(1,2,3,4,5,6,7,8,9,10)
>
> y el vector en el eje y:
>
> y<-c(0,2,15,30,50,NaN,NaN,10,2,1)
>
>
> Muchas gracias.
>
> --
>
> *Wilmer Contreras Sepulveda*
>
> *Grupo de Investigación en Desarrollo de Microelectronica Aplicada*
> *Universidad Francisco de Paula Santander *
>
> [[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

Re: [R-es] Chi cuadrado ajustado

[Jorge I Velez]

La prueba de Wald no es correcta en un caso como estos.

Sugeriría, a pesar de la poca información que proporciona Sebastián,
utilizar Regresión Logística o Regresión Poisson.   En R estos modelos
puede ajustarse con la función glm(...).

Dependiendo de las necesidades, también podrías usar Estadística Bayesiana
a través de una distribución multinomial y una Dirichlet, o los modelos GSK
(también conocidos como loglineales --- ver ?loglin en R).

Saludos,
Jorge.-


2017-08-30 16:14 GMT-05:00 Javier Marcuzzi 
:

> Estimados
>
> No estoy siguiendo el hilo, podría cometer un error, o darle un paro
> cardíaco a un estadístico. Pero podría ser otro test y no siempre el famoso
> chi 2.
>
> https://www.rdocumentation.org/packages/aod/versions/1.3/topics/wald.test
>
> https://es.wikipedia.org/wiki/Prueba_de_Wald
>
>
> Javier Rubén Marcuzzi
>
> De: Sebastian Gadea
> Enviado: miércoles, 30 de agosto de 2017 16:51
> Para: eric
> CC: Lista R
> Asunto: Re: [R-es] Chi cuadrado ajustado
>
> Hola Eric,
>
> 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.
>
> Y no me refiero a si quiero probar la independencia de X1 y X2 y tengo
> además la variable *sexo*, realizar la prueba para entre X1 y X2 para las
> diferentes categorías de la variable* sexo*, sino a alguna forma de
> contemplar la informacion en *sexo* sin dividir la tabla.
>
>  source=link_campaign=sig-email_content=webmail>
> Libre
> de virus. www.avast.com
>  source=link_campaign=sig-email_content=webmail>
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> El 29 de agosto de 2017, 14:21, eric  escribió:
>
> > Hola Sebastian, podrias ser mas especifico ? varias pruebas tienen el
> > nombre de chi-cuadrado,
> >
> > segun wikipedia:
> >
> > En estadística  y
> estadística
> > aplicada  se
> > denomina *prueba χ²* (pronunciado como «ji cuadrado» y a veces como «chi
> > cuadrado») a cualquier prueba
> >  en la que el
> > estadístico 
> > utilizado sigue una distribución χ²
> >  si la
> hipótesis
> > nula  es cierta.
> > Algunos ejemplos de pruebas χ² son:
> >
> >- La prueba χ² de Pearson
> >, la
> >cual tiene numerosas aplicaciones:
> >
> >
> >- La prueba χ² de frecuencias
> >- La prueba χ² de independencia
> >- La prueba χ² de bondad de ajuste
> >
> >
> >- La prueba χ² de Pearson con *corrección por continuidad* o
> *corrección
> >de Yates *
> >- La prueba de Bartlett de homogeneidad de varianzas
> >
> >
> >  ... y el termino "ajustar por otra variable" tampoco es tan claro ... en
> > ciencias sociales dicen algo parecido (controlar por otra variable)
> cuando
> > ajustan un modelo en dos pasos, primero con un subconjunto de las
> variables
> > respuesta y luego con todas ... te refieres a esto ?
> >
> > Aclaranos el punto sebastian,
> >
> > Saludos !!!
> >
> > Eric.
> >
> >
> >
> > On 08/29/2017 01:48 PM, Sebastian Gadea wrote:
> >
> > Una consulta desde un Uruguay muy lluvioso,
> >
> > como se realiza una prueba Chi-cuadrado ajustando por otra variable.
> >
> > Muchas gracias!
> >
> > Saludos!
> > Sebastián.
> >
> >   [[alternative HTML version deleted]]
> >
> > ___
> > R-help-es mailing listR-help-es@r-project.orghttps://stat.ethz.ch/
> mailman/listinfo/r-help-es
> >
> >
> > --
> > 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.
> >
> >
>
>
> --
> Saludos!
> Sebastián.
>
> [[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 

Re: [R-es] Promedio elementos vector

[Jorge I Velez]

Hola Manuel,

Dale una mirada a las funciones

?rowSums
?rowMeans

Saludos,
Jorge.-


2017-08-13 22:13 GMT-05:00 Manuel Máquez :

> Buenas noches estimados colegas:
> Adjunto la matriz,de la que quiero obtener, los promedios de cada elemento
> de los vectores que la componen. Me explico, el primer vector, tiene como
> primer elemento 23 y la suma  de todos los elementos es 134; por lo tanto
> el promedio será 0.1716, el segundo elemento de dicho vector es 20,
> entonces el promedio sería 0.1492, del tercero el promedio es 0.1343 y así
> sucesivamente.
> ¿Me podrían hacerme el favor de indicarme, el o los documentos que debo
> consultar, con objeto de solucionar mi problema?
> Por adelantado doy las gracias más cumplidas.
> Atentamente;
> *MANOLO MÁRQUEZ P.*
>
> ___
> 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

Re: [R-es] Paralelizar el cálculo de distancias

[Jorge I Velez]

Buenas tardes, Jesús.

Inicialmente no me funcionó tu código.  Al hacer algunas correcciones

foreach(j = 1:nrow(B),.combine="cbind")%:%
foreach(i = 1:nrow(A),.combine="c") %dopar%{
sqrt(sum((A[i,]-B[j,])^2))
}

obtuve una matriz de 10x10.  Es eso lo que buscas?

Saludos,
Jorge.-


2017-08-13 15:40 GMT-05:00 Jesús Para Fernández <
j.para.fernan...@hotmail.com>:

> Buenas,
>
>
> Quiero ver si he paralelizado correctamente el proceso. Tengo dos
> dataframes, A y B y quiero calcular la distancia euclídea de todas las
> filas de A sobre todas las filas de B. Para ello he hecho lo siguiente
>
> #cargo las librerias
>
> library(foreach)
>
> library(doParallel)
>
> #establezco el numero de clusters, en mi caso 4, ya que el procesador
> tiene 4 nucleos
>
> cl<-makeCluster(4)
> registerDoParallel(cl)
>
>
> #Creo los dataframes
>
> A<-as.data.frame(matrix(rnorm(50,10,2),ncol=5,nrow=10))
>
> B<-as.data.frame(matrix(rnorm(50,10,2),ncol=5,nrow=10))
>
>
> #calculo las distancias
>
> foreach(j in 1:nrow(B),.combine="cbind") %:%
>
> foreach(i in 1:nrow(A),.combine="c") %dopar% {
>
> sqrt(sum((A[i,]-B[j,])^2))
>
> }
>
>
>
> ¿Cómo lo veis?
>
>
> Un saludo
>
> Jesús
>
> [[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

Re: [R-es] FUNCION "SI" en R

[Jorge I Velez]

Javier,

Otra opción sería

require(car)
?recode

Hay varios ejemplos.

Saludos cordiales,
Jorge.-


2017-04-20 15:48 GMT-05:00 Javier Valdes Cantallopts (DGA) <
javier.val...@mop.gov.cl>:

> Hola  a todos
>
> Primero que todo, gracias por todos los consejos, me han servido ene.
>
>
>
> Mi pregunta es la siguiente:
>
> Como armar una función lógica “SI” PARA CAMBIAR LOS NOMBRES A LOS HEADER
>
> Por ejemplo, si tengo una tabla así;
>
> a
>
> b
>
> c
>
> d
>
> 1
>
> 3
>
> 4
>
> 6
>
> 2
>
> 4
>
> 6
>
> 7
>
>
>
> a= “pluvio mensual”
>
> b= “carga de batería”
>
> c= “temperatura”
>
> d= “radiación”.
>
> *La idea es reemplazar, en el caso que el HEADER sea (a), QUE EL NOMBRE
> DEL HEADER CAMBIE A “PLUVIO MENSUAL”*
>
>
>
> *Se me ocurre DAR UNA INSTRUCCIÓN TIPO COLNAMES<-¿*
>
> *Espero sus respuestas, saludos.*
>
>
>
>
>
>
>
>
>
> [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
>
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Encontrar la primera columna no NA

[Jorge I Velez]

Excelente!

Aqui hay otra aproximación en base:

R> t <- Sys.time()
R> out <- apply(as.matrix(dat), 1, function(x) x[!is.na(x)][1])
R> difftime(Sys.time(), t)
## Time difference of 0.656173 secs

Nada mal :)

Saludos,
Jorge.-


2016-10-27 11:42 GMT-05:00 Carlos J. Gil Bellosta :

> Las operaciones con columnas de data.frames (y sus variantes modernas) son
> muy caras. Así que:
>
> t <- Sys.time()
>
> tmp <- (as.matrix(dat))
> cols <- col(tmp)
> cols[is.na(tmp)] <- Inf
> my.cols <- apply(cols, 1, min)
> my.values <- tmp[cbind(1:nrow(tmp), my.cols)]
>
> difftime(Sys.time(), t)
>
> Y una pregunta: ¿alguien programa en R base todavía?
>
> Un saludo,
>
> Carlos J. Gil Bellosta
> http://www.datanalytics.com
>
> El 27 de octubre de 2016, 18:11, Olivier Nuñez  escribió:
>
> >
> > Por último, utilizando la indexación lineal de matriz que propusó luisfo
> > en su momento:
> >
> > > t <- Sys.time()
> > > M=as.matrix(dat)
> > > index <- which(!is.na(M)) - 1
> > > meses<-colnames(M)
> > > M2<- data.table(columna=index %/% nrow(M) +1L, jugador=index %% nrow(M)
> > +1L , valor=M[index+1L])
> > > setkey(M2,jugador,columna)
> > > M2[,.(First_month=meses[columna[1]],Value_First_month=
> > valor[1]),by=jugador]
> > jugador First_month Value_First_month
> >  1:   1 Uno0.93520715
> >  2:   2 Uno0.85930634
> >  3:   3 dos0.13521503
> >  4:   4 Uno0.86996341
> >  5:   5 dos0.65879889
> > ---
> >  6:   6 Uno0.94728423
> >  7:   7 Uno0.24088571
> >  8:   8 Uno0.07458581
> >  9:   9 Uno0.30535050
> > 10:  10 Uno0.54640585
> > > difftime( Sys.time(), t)
> > Time difference of 0.329 secs
> > >
> > - Mensaje original -
> > De: Javier Villacampa González 
> > Para: Olivier Nuñez 
> > CC: R ayuda 
> > Enviado: Thu, 27 Oct 2016 17:17:12 +0200 (CEST)
> > Asunto: Re: [R-es] Encontrar la primera columna no NA
> >
> > Hemos mejorado bastante desde el inicio. Pero aun andamos lejos. Igual es
> > por el merge que hago. Seguire mirando
> > library(microbenchmark)
> > N <- 1e1
> > tabla <-
> >   microbenchmark(
> > # JVG_dplyr ={
> > #   dat %>%
> > # apply( MARGIN = 1, FUN =
> > #  function(x){
> > #which( !is.na(x)  ) %>%  min( na.rm = TRUE ) %>%
> > return()
> > #}
> > # )
> > #   dat[ , First_month := First_month]
> > #   N_for <- length( unique(First_month ))
> > #   for( j in 1:N_for){
> > # x <- dat[  First_month == j,  j,  with = FALSE]
> > # dat[ First_month == j , Value_First_month := x ]
> > #   }
> > # },
> > JVG ={
> >   dat <-
> > data.table( Uno= sample( c(runif(numero) , rep(NA , numero
> /2e0
> > )) , size = numero ) ,
> > dos= sample( c(runif(numero) , rep(NA , numero
> /1e1
> > )) , size = numero ) ,
> > tres   = sample( c(runif(numero) , rep(NA , numero
> /2e1
> > )) , size = numero ) ,
> > cuatro = sample( c(runif(numero) , rep(NA , numero
> /1e2
> > )) , size = numero ) ,
> > cinco  = sample( c(runif(numero) , rep(NA , numero
> /2e2
> > )) , size = numero ) ,
> > seis   = sample( c(runif(numero) , rep(NA , numero
> /1e3
> > )) , size = numero )
> > )
> >
> >   apply(X = dat,  MARGIN = 1, FUN =
> >   function(x){
> > return(   min(  which( !is.na(x)  ),  na.rm = TRUE ) )
> >   }
> >   )
> >   dat[ , First_month := First_month]
> >   N_for <- length( unique(First_month ))
> >   for( j in 1:N_for){
> > x <- dat[  First_month == j,  j,  with = FALSE]
> > dat[ First_month == j , Value_First_month := x ]
> >   }
> > },
> > Olivier ={
> > dat <-
> >   data.table( Uno= sample( c(runif(numero) , rep(NA , numero /2e0
> > )) , size = numero ) ,
> >   dos= sample( c(runif(numero) , rep(NA , numero /1e1
> > )) , size = numero ) ,
> >   tres   = sample( c(runif(numero) , rep(NA , numero /2e1
> > )) , size = numero ) ,
> >   cuatro = sample( c(runif(numero) , rep(NA , numero /1e2
> > )) , size = numero ) ,
> >   cinco  = sample( c(runif(numero) , rep(NA , numero /2e2
> > )) , size = numero ) ,
> >   seis   = sample( c(runif(numero) , rep(NA , numero /1e3
> > )) , size = numero )
> >   )
> >   dat[,First_month   := apply(X = .SD,MARGIN = 1,FUN =
> function(x)
> > colnames(.SD)[min(which(!is.na(x)))])]
> >   dat[,Value_First_month := apply(X = .SD,MARGIN = 1,FUN =
> function(x)
> > x[min(which(!is.na(x)))])]
> > },
> >

Re: [R-es] Consulta sobre la creación de una columna mediante una condición

[Jorge I Velez]

Buenas noches, Oscar.

No entiendo muy bien que problema tines, pero intenta lo siguiente a
ver si produc lo que buscas:

R> #creo un data frame
R> node = c("TipoDTE","TipoDTE","Other")
R> value = c("33", "35", "")
R> j = data.frame(node, value)
R>
R> j
 node value
1 TipoDTE33
2 TipoDTE35
3   Other
R> with(j, ifelse(node == 'TipoDTE', value, ""))
[1] "2" "3" ""
R> j$TipoDTE.0 <- with(j, ifelse(node == 'TipoDTE', value, "")); j
 node value TipoDTE.0
1 TipoDTE33 2
2 TipoDTE35 3
3   Other
R> j$TipoDTE.0
[1] "2" "3" ""

Saludos,
Jorge.-


2016-10-12 10:09 GMT-05:00 Oscar Benitez :

> Hola
>
> Tengo una consulta simple pero que me está costando resolver
>
> #creo un data frame
> node = c("TipoDTE","TipoDTE","Other")
> value = c("33", "35", "")
> j = data.frame(node, value)
>
> en este punto el dataframe tiene dos campos "node" y "value" que R reconoce
> como factores de la siguiente manera:
> node:  Factor w/ 2 levels "Other","TipoDTE": 2 2 1
> value:   Factor w/ 3 levels "","33", "35" : 2 3 1
>
> Quiero crear un campo nuevo mediante una condición:
>
> j$TipoDTE.0<-ifelse (j$node=="TipoDTE",j$value,"")
>
> Y el resultado es:
>
> node value TipoDTE.0
> TipoDTE 33 2
> TipoDTE 35 3
> Other
>
>
> Lo que veo que está haciendo es colocar los niveles del factor en lugar del
> valor del factor
> Evidentemente estoy haciendo algo mal, pero no encuentro qué es...
>
> Cualquier ayuda será apreciada.
>
> Saludos
>
>
> --
> Oscar Benitez
>
> [[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

Re: [R-es] usar un for para cambiar indice

[Jorge I Velez]

Es más fácil si usas matplot.  Mira ?matplot para más información.  Saludos
cordiales.   --JIV


2016-09-06 13:48 GMT-05:00 Sebastian Kruk :

> Estimados usuarios-de-R:
>
> Tengo un problema.
>
> Si por ejemplo tengo una lista conformada por la serie historia del PBI de
> 10 paìses.
>
> >Tiempo <- seq(1:100)
> >plot(Tiempo, PP$"Serie 1", type="l",col="1")
> >lines(Tiempo, PPcapita$"Serie 2")
> >lines(Tiempo, PPcapita$"Serie 3")
> >lines(Tiempo, PPcapita$"Serie 4")
> >lines(Tiempo, PPcapita$"Serie 5")
> >lines(Tiempo, PPcapita$"Serie 6")
> >lines(Tiempo, PPcapita$"Serie 7")
> >lines(Tiempo, PPcapita$"Serie 8")
> >lines(Tiempo, PPcapita$"Serie 9")
> >lines(Tiempo, PPcapita$"Serie 10")
>
> ¡Y quiero usar un for para lo anterior o vectorizar como lo hago?
>
> Me da error hacer lo siguiente:
>
> Tiempo <- seq(1:100)
> plot(Tiempo, PP$"Serie 1", type="l",col="red")
> for (i in 2:10) {lines(Tiempo, PPcapita$"Serie i",col=i)
>
> ¿Hay alguna forma de resolverlo?
>
> Saludos,
>
> Sebas.
>
> [[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

Re: [R-es] Un pequeño favor....

[Jorge I Velez]

Buenas noches, Heber.

La función ls() no permite hacer lo que buscas:

R> ls()[1:10]
 [1] "cols"  "extractRcode"  "is.prime"  "lun"
 [5] "msp"   "myscatterplot" "pubmed""rae"
 [9] "reducePDF" "spanish"

R> ls()[1:2]
[1] "cols" "extractRcode"

Sin embargo, y creo que es lo que estás buscando

R> head(d, n)

es tu mejor opción.  En la expresión anterior "d" corresponde al objeto
que contiene los datos y "n" el número de filas que quieres observar.

Si por el contrario lo que quieres ver es sólo los encabezados, necesitas

R>  colnames(d)

donde "d" es un objeto R que contiene los datos.

Espero sea de utilidad.

Saludos,
Jorge Velez.-



2016-07-14 17:39 GMT-05:00 heber sarmiento via R-help-es <
r-help-es@r-project.org>:

> Cordial saludo a todos y gracias por la ayuda que me puedan brindar,
> El asunto es que tengo una pequeña base de datos cargada en R y al usar
> ls() para listar salen rápidamente todos los datos y no me permite ver los
> encabezados en cada una de las variables, agradezco si alguine me indica si
> ls() tiene algunas opciones para trabajar por pantallas, es decir que me
> deje ver los encabezados.
> Un abarazo freternal y gracias.
> Heber
> [[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

Re: [R] scatter3D and colours

[Jorge I Velez]

Thank you very much, Sarah, for your help.  The colvar argument was
certainly what I needed.  I will follow your suggestion.
Cheers,
Jorge.-



On Thu, Jul 14, 2016 at 1:29 PM, Sarah Goslee <sarah.gos...@gmail.com>
wrote:

> I assume you want some variant of:
> with(d0, scatter3D(X, Z, Y, bty = "b2", col = cols, colvar = c(1, 2, 3),
>pch = 20, cex = 4, ticktype = "detailed", colkey = FALSE, phi =
> 20, theta = -140, zlim = c(-14, 14), xlim = c(-14, 14), ylim = c(-14, 14),
> xlab = "X", ylab = "Z", zlab = "Y", CI = CI2))
> mtext('Experiment 2', line = -0.5, side = 3, cex = 1.3)
>
> You should reread the help for scatter3D, paying close attention to
> col and colvar arguments. The default value of colvar is z, and that's
> what scatter3D() was using. Setting colvar to your chosen levels gives
> you the result you expect.
>
> Thanks for the complete reproducible example. I wouldn't have even
> looked at your problem without something to test.
>
> Sarah
>
> On Thu, Jul 14, 2016 at 1:55 PM, Jorge I Velez <jorgeivanve...@gmail.com>
> wrote:
> > Dear R-help,
> >
> > I am using the plot3D package to produce 3D spheres along with 95% CIs
> > distinguishing each sphere with a predefined colour (see the reproducible
> > example at the end).
> >
> > I have been successful in producing a similar plot using a different data
> > set (kindly see
> > https://www.dropbox.com/s/snpgiqgqaiqgpiv/exp1combined.pdf?dl=0), but in
> > this case I cannot arrange the colours as desired. Specifically, the
> > colours obtaind do not correspond to those in the "cols" object below
> (kindly
> > see https://www.dropbox.com/s/nhljsuare84g811/test_exp2.pdf?dl=0).
> >
> > Any help/tips would be greatly appreciated.
> >
> > Thank you very much in advance.
> >
> > Best regards,
> > Jorge Velez.-
> >
> >
> > ## package needed
> > if(!require(plot3D)) install.packages("plot3D")
> > require(plot3D)
> >
> > ## data to be plotted
> > d0 <- structure(list(X = c(2.5, -2, 1), Z = c(-3.5, 4, -1), Y = c(8,
> > -8.5, -1)), .Names = c("X", "Z", "Y"), row.names = c("high",
> > "low", "medium"), class = "data.frame")
> > d0
> >
> > ## confidence intervals to be added
> > CI2 <- structure(list(z = structure(c(2, 3, 4, 2, 3, 4), .Dim = c(3L,
> > 2L), .Dimnames = list(c("high", "low", "medium"), c("97.5%",
> > "97.5%"))), y = structure(c(4, 5, 5, 4, 5, 5), .Dim = c(3L, 2L
> > ), .Dimnames = list(c("high", "low", "medium"), c("97.5%", "97.5%"
> > ))), x = structure(c(3, 7, 5, 3, 7, 5), .Dim = c(3L, 2L), .Dimnames =
> list(
> > c("high", "low", "medium"), c("97.5%", "97.5%"))), alen = 0,
> > lwd = 2), .Names = c("z", "y", "x", "alen", "lwd"))
> >
> > ## colours I would like to have
> > cols <- c("#0080ff", "#ff00ff", "darkgreen")
> >
> > # this produces the 3D plot, but the colours are not properly assigned
> > with(d0, scatter3D(X, Z, Y, bty = "b2", col = cols,
> >pch = 20, cex = 4, ticktype = "detailed", colkey = FALSE, phi
> =
> > 20, theta = -140, zlim = c(-14, 14), xlim = c(-14, 14), ylim = c(-14,
> 14),
> > xlab = "X", ylab = "Z", zlab = "Y", CI = CI2))
> > mtext('Experiment 2', line = -0.5, side = 3, cex = 1.3)
> >
>

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] scatter3D and colours

[Jorge I Velez]

Dear R-help,

I am using the plot3D package to produce 3D spheres along with 95% CIs
distinguishing each sphere with a predefined colour (see the reproducible
example at the end).

I have been successful in producing a similar plot using a different data
set (kindly see
https://www.dropbox.com/s/snpgiqgqaiqgpiv/exp1combined.pdf?dl=0), but in
this case I cannot arrange the colours as desired. Specifically, the
colours obtaind do not correspond to those in the "cols" object below (kindly
see https://www.dropbox.com/s/nhljsuare84g811/test_exp2.pdf?dl=0).

Any help/tips would be greatly appreciated.

Thank you very much in advance.

Best regards,
Jorge Velez.-


## package needed
if(!require(plot3D)) install.packages("plot3D")
require(plot3D)

## data to be plotted
d0 <- structure(list(X = c(2.5, -2, 1), Z = c(-3.5, 4, -1), Y = c(8,
-8.5, -1)), .Names = c("X", "Z", "Y"), row.names = c("high",
"low", "medium"), class = "data.frame")
d0

## confidence intervals to be added
CI2 <- structure(list(z = structure(c(2, 3, 4, 2, 3, 4), .Dim = c(3L,
2L), .Dimnames = list(c("high", "low", "medium"), c("97.5%",
"97.5%"))), y = structure(c(4, 5, 5, 4, 5, 5), .Dim = c(3L, 2L
), .Dimnames = list(c("high", "low", "medium"), c("97.5%", "97.5%"
))), x = structure(c(3, 7, 5, 3, 7, 5), .Dim = c(3L, 2L), .Dimnames = list(
c("high", "low", "medium"), c("97.5%", "97.5%"))), alen = 0,
lwd = 2), .Names = c("z", "y", "x", "alen", "lwd"))

## colours I would like to have
cols <- c("#0080ff", "#ff00ff", "darkgreen")

# this produces the 3D plot, but the colours are not properly assigned
with(d0, scatter3D(X, Z, Y, bty = "b2", col = cols,
   pch = 20, cex = 4, ticktype = "detailed", colkey = FALSE, phi =
20, theta = -140, zlim = c(-14, 14), xlim = c(-14, 14), ylim = c(-14, 14),
xlab = "X", ylab = "Z", zlab = "Y", CI = CI2))
mtext('Experiment 2', line = -0.5, side = 3, cex = 1.3)


​## R session details
R> sessionInfo()
R version 3.2.4 Patched (2016-03-28 r70416)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.11.5 (El Capitan)

locale:
[1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8

attached base packages:
[1] stats graphics  grDevices utils datasets  parallel  compiler
[8] methods   base

other attached packages:
[1] plot3D_1.1   readxl_0.1.1

loaded via a namespace (and not attached):
[1] tools_3.2.4  misc3d_0.8-4 Rcpp_0.12.5

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R-es] Reemplazar NA con el último valor no NA de una columna en una data.table

[Jorge I Velez]

Hola Patricio,

Intenta

R> x <- c(1, 2, NA, 4, NA, 5)
R> x
[1]  1  2 NA  4 NA  5
R> require(zoo)
Loading required package: zoo

Attaching package: ‘zoo’

The following objects are masked from ‘package:base’:

as.Date, as.Date.numeric

Warning message:
package ‘zoo’ was built under R version 3.2.5

R> na.locf(x)
[1] 1 2 2 4 4 5

Saludos,
Jorge.-



2016-06-08 12:21 GMT-05:00 Patricio Fuenmayor 
:

> Hola, favor denme una mano.
> Tengo una data.table que contiene columnas con algunos valores NA.
> Necesito reemplazar estos NA con el anterior valor no NA de la columna.
> Ejemplor V1=1,2,NA,4,NA,5. Debo obtener V1=1,2,2,4,4,5
> Y este proceso lo debo realizar a varias columnas de la data.table
>
> Gracias por la ayuda...
>
> [[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

Re: [R-es] Ayuda con R Sweave

[Jorge I Velez]

Hola Jesús,
Revisa las opciones en http://yihui.name/knitr/options/  Necesitas
los argumentos fig.width y fig.height.
Saludos,
Jorge.-



2016-06-02 8:56 GMT-05:00 Jesús Para Fernández :

>
> Gracias por las respuestas
>
> La verdad es qeu estoy trabjando con Rstudio + Knitr, lo unico que en R
> studio le doy a crear nuevo Sweave, pero is me carga la libreria knitr.
>
> Solucioné el problema que os comenté, era por un fallo tonto, a la hora de
> poner datos$A, puse dats$A, es decir, le faltaba la o
>
> Ahora una duda que no se como resolver es la de los gráficos. No consigo
> reducir su tamaño...
>
> 

Re: [R-es] Ayuda con R Sweave

[Jorge I Velez]

Buenos dias Jesús,
 Existe alguna razón especial para no trabajar en RStudio + knitr?  IMO,
es mucho mas fácil y flexible.  Más información en
http://yihui.name/knitr/demo/rstudio/
Saludos,
Jorge.-



2016-06-02 6:10 GMT-05:00 Jesús Para Fernández :

> >Buenas,
>
> Estoy empezando a utilizar latex y hay un par de cosas que no se muy bien
> como hacer ni como buscar en google.
>
> Tengo un data.frame llamado datos con 3 variables, A,B y C
>
> Si pongo
> <<>>
> datos$A
> @
>
> me devuelve error
>
> sin embargo si lo hago con
> <<>>=
> datos[,1]
> @
> me da correctamente. No puedo con Sneave al interior de un data frame
> usando el nombre de las variables con el $???
>
> 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
>

[[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] Última (más reciente) posición de un vector (con "n" desconocido)

[Jorge I Velez]

Hola Gemma,

Intenta

x[length(x)]
tail(x, 1)

Saludos,
Jorge.-



2016-05-28 11:52 GMT-05:00 Gemma Ruiz-Olalla :

> Hola,
>
> Hay alguna forma de acceder a la última posición de un vector si se
> desconoce "n"? (por estar dentro de un bucle "for" que va incrementando su
> contador)
>
> La idea es acceder al valor de vector[n] en el caso de que "n" sea una
> incógnita.
>
> Estoy casi segura de que es posible pero mi conocimiento de R es bastante
> limitado.
>
> Muchas gracias,
>
> Gemma
>
> [[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

[R-es] Aplicar una función repetidamente

[Jorge I Velez]

Hola a todos,

Quisiera aplicar una función f(x) un total de k veces de manera recursiva.
En pseudo código sería algo como

Si k = 1, calcular f(x);
Si k = 2, calcular f(f(x));
Si k = 3, calcular f(f(f(x))).

Al final me gustaria tener una función g cuyos argumentos sean x y el valor
de k. Así,

g(x, k = 2)

daría como resultado f(f(x)).

Cualquier ayuda y/o sugerencia será más que bienvenida.

Muchísimas gracias,
Jorge Velez.-

[[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 Error en modelo lineal

[Jorge I Velez]

Buenos dias Jesus,

Esos valores son _aproximados_.  Las estimaciones, de acuerdo con teoria de
regresion, podrian obtenerse de manera puntual y construir intervalos de
_confianza_ y _prediccion_ alrededor de estos.  Ten en cuenta que el 2do
tipo de intervalos de calcula para observaciones _futuras_.

En R puedes calcularlos de la siguiente manera:

## IC de confianza
## ver ?predict.lm para mas detalles
R> data.frame(y, predict(modelo, interval = "confidence"))
y   fit   lwr  upr
#1  8.35  9.938571  6.580445 13.29670
#2 12.42 11.804286  9.134239 14.47433
#3 18.00 15.535714 13.664949 17.40648
#4 17.58 17.401429 15.396872 19.40599
#5 17.97 19.267143 16.798908 21.73538
#6 20.76 21.132857 18.014915 24.25080

## intervalos de prediccion para x = 25
R> predict(modelo, newdata = data.frame(x = 25), interval = "prediction")
#fit  lwr  upr
#1 23.93143 17.69035 30.17251

Lo anterior significa que E[y|x=25] = 23.93 y el intervalo de prediccion
del 95% es (17.69, 30.17).

Espero sea de utilidad.

Saludos cordiales,
Jorge.-




2016-04-21 8:56 GMT-05:00 Jesús Para Fernández :

> Buenas, una pregunta.
>
> Si yo estoy calculando un modelo lineal, el caso más simple, 1 variable
> respuesta y una variable explicativa y creo un modelo, me da un R2 del 80%
> y quiero ver como es esa relacion entre las variables, para calcular el
> error de predicción del modelo, basta con ver el intervalo de confianza del
> modelo e irme a los extremos?
>
> Por si no me he expresado bien, un ejemplo tonto:
>
> y=c(8.35,12.42,18.00,17.58,17.97,20.76)
> x=c(10,12,16,18,20,22)
>
> modelo<-lm(y~x)
> summary(modelo)
> library(MASS)
> Confint(modelo, level=0.95)
>
> Con esto tengo un modelo :
> > Confint(RegModel.1, level=0.95)
>   Estimate   2.5 %  97.5 %
> (Intercept) 0.610  -6.8296312 8.049631
> x 0.9328571  0.4919135  1.373801
>
> Es decir, el intervalo de confianza de la respuesta y en función de los
> valores x, sería de la forma (ymin,ymax), siendo:
> ymin = -6,82 + 0,49*x
> ymax = 8,05 + 1,37*x
>
> Es esto correcto???
>
> Gracias
> Jesús
>
>
>
>
> [[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

Re: [R-es] Cambiar \u

[Jorge I Velez]

Muchas gracias a Luisfo y Milagros por las respuestas.  Ambas soluciones
funcionaron a la perfección.
Saludos cordiales,
Jorge.-



2016-03-17 2:00 GMT-05:00 Luisfo Chiroque <luisf...@yahoo.com>:

> Buenas Jorge,
>
> Creo que esta linea te podria servir:
> gsub('[^[:alnum:]]', '', words)
>
> Lo que hace basicamente es quedarse con los caracteres alfanumericos. Por
> tanto, esos ‘\u00XX’ desapareceran.
>
> Espero que te sea util.
>
> Un saludo,
> Luisfo
>
> P.D. va todo sin acentos u.u
>
>
> > El 17 mar 2016, a las 5:00, Jorge I Velez <jorgeivanve...@gmail.com>
> escribió:
> >
> > Buenas noches a todos,
> >
> > Estoy trabajando con textos y me gustaría realizar algunos cambios dentro
> > de las palabras.
> >
> > Cómo puedo elimiar "\u" y los números en "\u0093\u0085anunció", de tal
> > forma que el resultado sea "anunció"?   En total tengo un poco menos
> > de 50,000 palabras, así que a continuación proporciono un ejemplo:
> >
> > # palabras
> > words <- c("además", "\u0093\u0085anunció", "acerca", "cierto", "el",
> > "en", "es", "exagera", "frente", "hizo")
> > words
> >
> > Muchas gracias por la ayuda.
> >
> > Saludos,
> > Jorge Velez.-
> >
> >
> > # session
> > R version 3.2.3 (2015-12-10)
> > Platform: x86_64-apple-darwin14.5.0 (64-bit)
> > Running under: OS X 10.11.2 (El Capitan)
> >
> > locale:
> > [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8
> >
> > attached base packages:
> > [1] stats graphics  grDevices utils datasets  parallel  compiler
> > [8] methods   base
> >
> > other attached packages:
> > [1] sentiment_1.0  plyr_1.8.3 rjson_0.2.15
> > [4] RCurl_1.95-4.8 bitops_1.0-6   igraph_1.0.1
> > [7] stringr_1.0.0  tm_0.6-2   NLP_0.1-9
> > [10] wordcloud_2.5  RColorBrewer_1.1-2 readxl_0.1.0
> >
> > loaded via a namespace (and not attached):
> > [1] Rcpp_0.12.3   slam_0.1-32   magrittr_1.5  stringi_1.0-1 tools_3.2.3
> >
> >   [[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

Re: [R] Extract every 2 element for a list which are not equal in length

[Jorge I Velez]

Dear Mohammad,
What's wrong with the result?
Best,
Jorge.-

On Monday, March 7, 2016, Mohammad Tanvir Ahamed via R-help <
r-help@r-project.org> wrote:

> Hi,
>
> a <- c(1:5)b <- c(1:3)
> c <- 1
> d <- 5
> e <- list(a,b,c,d)
>
> # To extract every 1st element
> lapply(e,"[[",1)
>
> ## Out-put
> [[1]]
> [1] 1
>
> [[2]]
> [1] 1
>
> [[3]]
> [1] 1
>
> [[4]]
> [1] 5
>
> #To extract every 2nd element (Need help in this case)
> lapply(e,"[[",2)
>
> ## Expected outcome
> [[1]]
> [1] 2
>
> [[2]]
> [1] 2
>
> [[3]]
> [1] NA
>
> [[4]]
> [1] NA
>
>
>
> Any help will be appreciated . Thanks
>
>
> Tanvir Ahamed
> Göteborg, Sweden  |  mashra...@yahoo.com 
>
> __
> R-help@r-project.org  mailing list -- To UNSUBSCRIBE and
> more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R-es] Consulta r studio no funciona

[Jorge I Velez]

Amalia,
Instalaste R?  Necesitas R para que RStudio funcione.
Saludos,
Jorge.-



2016-03-01 21:45 GMT-05:00 Amalia Carolina Guaymas Canavire <
acarolin...@gmail.com>:

> Hola les escribo porque R studio en windows 8 no me funciona y no sé como
> solucionarlo. El error que arroja es que R no puede iniciar, la sección de
> environment  no se carga, aparece como procesando y cuando deseo probar
> algo en la consola la misma no devuelve resultado.
> Re instale R studio y aún persiste el error.
> Gracias por su ayuda.
>
> [[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

[R-es] Ejecutar una función automáticamente

[Jorge I Velez]

Buenas tardes a todos,

Me gustaría ejecutar una función

f(x, when)

cada cierto tiempo; "x" es un string y "when" es la frecuencia de ejecución
(cada hora, cada media hora).

Alguna sugerencia?

Estoy trabajando en Windows, pero tambien tengo acceso a OS X, asi que
agradecería si pudieran orientarme sobre qué hacer en ambos sistemas
operativos.

Gracias a todos de antemano por la ayuda.

Saludos cordiales,
Jorge Velez.-

[[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] Determinación del punto de corte óptimo

[Jorge I Velez]

Buenas tardes Fernando,

No entiendo muy bien los detalles de tu mensaje (no hay espacios y es
dificil leerlo), pero te sugiero que le des una mirada a
http://biostat.mc.vanderbilt.edu/wiki/Main/CatContinuous antes de proceder.

Saludos cordiales,
Jorge.-


2016-01-26 14:04 GMT-05:00 Fernando Sanchez :

>  Hola a todos, en estos momentos me encuentro inmerso en laresolución del
> siguiente problema. Resulta que dispongo de una variable
> categóricadenominada Severity y que consta de tres categorías (low, medium,
> high). Además,tengo otra variable que se denomina ZX y que puede tomar
> cualquier valorcomprendido entre 0 y 10. Quiero determinar los dos puntos
> de corte óptimos demanera que me dividan a la variable ZX en tres
> categorías (llamadas tambiénlow, medium y high) y que exista el máximo
> nivel de acuerdo entre las variablesZX y Severity. Por acuerdo me refiero a
> que el porcentaje de casos (filas) queen ambas variables caiga en la misma
> categoría sea el máximo posible.Dado que no se me ocurría cómo resolverlo
> de formasistemática, he hecho un pequeño código en R que fija dos puntos de
> cortealeatorios y calculo la tasa de acuerdo. A base de repetirlo 500.000
> veces, mehe hecho una idea de hasta qué nivel de concordancia puedo llegar
> pero buscabauna forma más sistemática. En el caso de un único punto de
> corte, sé queexisten métodos y librerías en R como OptimalCutPoints que
> permitenhacerlo.Debajo os pongo un ejemplo de los datos que estoy
> manejandopara que os hagáis mejor una idea de su
> aspecto.ZX  Severity2.818181818high2.084242424
> high5.32667medium4.758484848low4.795454545high3.367878788
> high5.734848485high3.417575758medium3.16000
> medium4.307272727lowCualquier sugerencia es bienvenida. Si mi
> código le puedeayudar a alguien, estaré encantado en pasárselo.
> [[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

Re: [R-es] Borrar cada fila 400

[Jorge I Velez]

Hola Jesús,

No es necesario un loop para ello.  A continuación una idea utilizando
seq():

datos[-seq(400, NROW(datos), by = 400), ]

Saludos cordiales,
Jorge Velez.-


2015-11-17 9:14 GMT-05:00 Jesús Para Fernández :

> Buenas, tengo un csv [csv final] con 5 filas, que es unión de varios
> csv [csv particular].
>
> Cada csv [csv particular] tiene en la última fila, la 400, una serie de
> valores que quiero eliminar, por lo que del [csv filan] quiero borrar la
> linea 400,800,1200,
>
>
> Lo he intentado con un bucle for:
>
> for(i in 1:5){
> if(i%%400 == 0) {datos[-i,]}
> }
>
>
> Pero no me funciona. Además me han dicho que con apply puede ser mucho más
> eficiente el algoritmo. ¿Alguna idea?
>
> 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
>

[[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] Extract entries from matrix

[Jorge I Velez]

Dear Dr. Steipe,

Thank you for the code.  It does exactly what I needed.

Best regards,
Jorge.-


On Wed, Oct 28, 2015 at 1:18 PM, Boris Steipe <boris.ste...@utoronto.ca>
wrote:

> Your code does not produce the matrix in your image.
> The first three rows contain all-zeros and the last row is missing.
> The following line fixes that:
>
> m <- rbind(m[-(1:3), ], 1:5)
>
> Given that matrix, the following code produces the output
> you have illustrated. It's so trivial however that I suspect
> something must be missing in your problem description.
>
>
>
> v <- numeric(nrow(m))
> j <- 1
> for (i in 1:nrow(m)) {
> if (j > ncol(m) || m[i,j] == 0) {
> j <- 1
> }
> v[i] <- m[i,j]
> j <- j+1
> }
>
> v
>
> Note that this puts 0 in the output if
> there is a zero in your first column and
> the "diagonal". Your example didn't have that.
>
>
>
> B.
>
>
> On Oct 28, 2015, at 1:17 PM, Jorge I Velez <jorgeivanve...@gmail.com>
> wrote:
>
> > Dear all,
> >
> > I thought I would better send an image illustrating that the problem is
> > (hope the file gets through).  In the picture, the matrix "m" is given by
> >
> > ## input
> > m <- structure(c(0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 2, 2,
> > 2, 2, 2, 2, 2, 2, 2, 0, 0, 0, 3, 3, 3, 3, 3, 3, 3, 3, 3, 0, 0,
> > 0, 0, 0, 0, 0, 0, 4, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5,
> > 5), .Dim = c(12L, 5L))
> >
> > We start from the entry [1,1] and select all values in the diagonal until
> > we find the first zero. That is, we select the values in purple.  Because
> > the first zero is found in [4,4], the next matrix begins in the 4th row.
> > The values of interest would be 1, 2, 3, 4, 5 (light blue).  The last
> value
> > in this resulting vector is entry m[8,5], which is located at the maximum
> > number of columns of m. Hence, the next matrix to work with starts at
> row 9
> > and the values of interest are 1, 2 (in orange).
> >
> > The output vector would then be of the same length as the number of rows
> in
> > "m", and would contain the elements previously selected.
> >
> > Any ideas on how to proceed?
> >
> > Thank you very much in advance.
> >
> > Best regards,
> > Jorge Velez.-
> >
> >
> >
> > On Tue, Oct 27, 2015 at 4:38 PM, Jorge I Velez <jorgeivanve...@gmail.com
> >
> > wrote:
> >
> >> Thank you all for your solutions and comments.
> >>
> >> As Dr. Carlson mentioned, we leave rows 1 to 3 out as they are all
> zeroes.
> >> Then, the entries I need to select from m are
> >>
> >> 
> >> entry value
> >> 
> >> 4,1  ---> 1
> >> 5,2  ---> 2
> >> 6,3  ---> 3
> >> 7,1  ---> 1
> >> 8,2  ---> 2
> >> 9,3  ---> 3
> >> 10,4  ---> 4
> >> 11,5  ---> 5
> >> 12,1  ---> 1
> >>
> >> Note that the entry [7,4] is zero, so we start from the first column in
> >> the 7th row and then select entry [7,1] instead.  That's what I meant
> by  "...
> >> the idea is to extract the diagonal elements until a zero is found."  I
> >> should have said *entries* instead of  _diagonal elements_. I am sorry
> Dr.
> >> Turner for the confusion.
> >>
> >> Starting with m
> >>
> >> R> m
> >> #   [,1] [,2] [,3] [,4] [,5]
> >> # [1,]00000
> >> # [2,]00000
> >> # [3,]00000
> >> # [4,]12300
> >> # [5,]12300
> >> # [6,]12300
> >> # [7,]12300
> >> # [8,]12300
> >> # [9,]12340
> >> #[10,]12340
> >> #[11,]12345
> >> #[12,]12345
> >>
> >> the first submatrix to work with is
> >>
> >> # [4,]12300
> >> # [5,]12300
> >> # [6,]12300
> >>
> >> from which the elements of interest are 1, 2, 3.  Note that the 7th row
> of
> >> m is not included here because m[7, 5] = 0.
> >>
> >> Further, the second submatrix is
> >>
> >> # [7,]12300
> >> # [8,]12300
> >> # [9,]12340
> >&g

Re: [R] Extract entries from matrix

[Jorge I Velez]

Dear all,

I thought I would better send an image illustrating that the problem is
(hope the file gets through).  In the picture, the matrix "m" is given by

## input
m <- structure(c(0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 2, 2,
2, 2, 2, 2, 2, 2, 2, 0, 0, 0, 3, 3, 3, 3, 3, 3, 3, 3, 3, 0, 0,
0, 0, 0, 0, 0, 0, 4, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5,
5), .Dim = c(12L, 5L))

We start from the entry [1,1] and select all values in the diagonal until
we find the first zero. That is, we select the values in purple.  Because
the first zero is found in [4,4], the next matrix begins in the 4th row.
The values of interest would be 1, 2, 3, 4, 5 (light blue).  The last value
in this resulting vector is entry m[8,5], which is located at the maximum
number of columns of m. Hence, the next matrix to work with starts at row 9
and the values of interest are 1, 2 (in orange).

The output vector would then be of the same length as the number of rows in
"m", and would contain the elements previously selected.

Any ideas on how to proceed?

Thank you very much in advance.

Best regards,
Jorge Velez.-



On Tue, Oct 27, 2015 at 4:38 PM, Jorge I Velez <jorgeivanve...@gmail.com>
wrote:

> Thank you all for your solutions and comments.
>
> As Dr. Carlson mentioned, we leave rows 1 to 3 out as they are all zeroes.
> Then, the entries I need to select from m are
>
> 
> entry value
> 
> 4,1  ---> 1
> 5,2  ---> 2
> 6,3  ---> 3
> 7,1  ---> 1
> 8,2  ---> 2
> 9,3  ---> 3
> 10,4  ---> 4
> 11,5  ---> 5
> 12,1  ---> 1
>
> Note that the entry [7,4] is zero, so we start from the first column in
> the 7th row and then select entry [7,1] instead.  That's what I meant by  "...
> the idea is to extract the diagonal elements until a zero is found."  I
> should have said *entries* instead of  _diagonal elements_. I am sorry Dr.
> Turner for the confusion.
>
> Starting with m
>
> R> m
> #   [,1] [,2] [,3] [,4] [,5]
> # [1,]00000
> # [2,]00000
> # [3,]00000
> # [4,]12300
> # [5,]12300
> # [6,]12300
> # [7,]12300
> # [8,]12300
> # [9,]12340
> #[10,]12340
> #[11,]12345
> #[12,]12345
>
> the first submatrix to work with is
>
> # [4,]12300
> # [5,]12300
> # [6,]12300
>
> from which the elements of interest are 1, 2, 3.  Note that the 7th row of
> m is not included here because m[7, 5] = 0.
>
> Further, the second submatrix is
>
> # [7,]12300
> # [8,]12300
> # [9,]12340
> #[10,]12340
> #[11,]12345
>
> and the corresponding elements are 1, 2, 3, 4, 5.
>
> And the last matrix is
>
> #[12,]12345
>
> from which the position [12,1] is selected.
>
> So, the resulting entries from this process are 1, 2, 3, 1, 2, 3, 4, 5, 1.
>
> Thank you in advance for any additional insight you may provide.
>
> Regards,
> Jorge Velez.-
>
>
>
> On Tue, Oct 27, 2015 at 4:06 PM, David L Carlson <dcarl...@tamu.edu>
> wrote:
>
>> I don't see how you are getting the result you provide.
>>
>> > m
>>   [,1] [,2] [,3] [,4] [,5]
>>  [1,]00000
>>  [2,]00000
>>  [3,]00000
>>  [4,]12300
>>  [5,]12300
>>  [6,]12300
>>  [7,]12300
>>  [8,]12300
>>  [9,]12340
>> [10,]12340
>> [11,]12345
>> [12,]12345
>> > t(sapply(1:8, function(x) diag(m[x:12, ])))
>>  [,1] [,2] [,3] [,4] [,5]
>> [1,]00000
>> [2,]00300
>> [3,]02300
>> [4,]12300
>> [5,]12300
>> [6,]12340
>> [7,]12345
>> [8,]12345
>>
>> These are all of the diagonals from the 1st through 8th rows. The first 3
>> begin with 0 so we leave them out, but then we have 4th: 1, 2, 3; 5th: 1,
>> 2, 3; 6th: 1, 2, 3, 4, etc so you must have some additional rule in mind to
>> get your answer.
>>
>> -
>> David L Carlson
>> Department of

Re: [R-es] Extraer elementos diagonales de submatrices

[Jorge I Velez]

Estimado Javier,

Gracias por tu mensaje.

No, lo unico que requiero es la lista de números (i1, 2, 3, 1, 2, 3, 4, 5,
1, 2).

Saludos cordiales,
Jorge.-



2015-10-28 14:35 GMT-05:00 Javier Rubén Marcuzzi <
javier.ruben.marcu...@gmail.com>:

> Estimado Jorge I Velez
>
>
>
> No comprendo un punto, dices que deseas construir sub matrices y extraer
> elementos de sub matrices, en el ejemplo en output no hay sub matrices
> (tres matrices como resultado) sino una cadena de números. ¿Cómo necesitas
> el resultado?, ¿Cómo se ve en el ejemplo?, ¿Cómo matrices de matrices?
>
>
>
> Ejemplo
>
> 123
>
> 12345
>
> 12
>
>
>
> O en una forma
>
> 1231234512
>
>
>
> ¿Hay que dejar algo (índice) como para que accedas a algo reconocido o
> especificado, ej, segundo valor de la segunda sub matriz (2,2), ¿o con
> tener la lista de números le es útil?
>
>
>
> Javier Rubén Marcuzzi
> Técnico en Industrias Lácteas
> Veterinario
>
>
>
>
>
>
> *De: *Jorge I Velez
> *Enviado: *miércoles, 28 de octubre de 2015 12:15
> *Para: *R-help-es
> *Asunto: *[R-es] Extraer elementos diagonales de submatrices
>
>
>
>
>
> Buenos dias a todos,
>
>
>
> Quisiera extraer algunas entradas de una matrix "m" teniendo en cuenta
> algunas restricciones.  El siguiente ejemplo ilustra la situacion:
>
>
>
> ## input
>
> m <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2,
>
> 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3,
>
> 3, 0, 0, 0, 0, 0, 4, 4, 4, 4, 4, 0, 0, 0, 0,
>
> 0, 0, 0, 5, 5, 5), .Dim = c(10L, 5L))
>
> m
>
>
>
> ## output
>
> output <- c(1:3, 1:5, 1:2)
>
> output
>
>
>
> Si el numero de filas es r y el numero de columnas k, la idea es construir
> submatrices de dimension k x k y extraer los elementos diagonales NO
> ceros.  Ahora, en caso de encontrar un cero, debe desplazarse a la
> siguiente fila, y construir una nueva matriz k x k.  Graficamente esto
> seria:
>
>
>
> [image: Inline image 2]
>
>
>
> Observe que en este caso, r = 10 y k = 5.  En la primera submatriz, la
> diagonal tiene los valores 1, 2, 3, 0, 0, de los cuales SOLO deben
> seleccionarse 1, 2 y 3 (en lila).  El primer cero se encuentra en la
> posicion [1, 4], asi que la siguiente submatriz debe construirse COMENZANDO
> en la fila 4, columna 1.
>
>
>
> A partir de esta segunda submatriz se obtienen los elements 1, 2, 3, 4 y
> 5, que corresponden a su diagonal (en azul celeste).  La posicion del
> ultimo elemento de esta diagonal es [8, 5], asi que la siguiente submatriz
> debe comenzar en la fila 9.  A partir de esta submatriz se obtienen los
> valores 1 y 2 (en naranja).  El resultado final de todo este proceso es el
> vector
>
>
>
> # [1] 1 2 3 1 2 3 4 5 1 2
>
>
>
> Agradezco a todos el tiempo que les tomo leer este mensaje, y por supuesto
> por cualquier sugerencia que me permita obtener este vector.
>
>
>
> Saludos cordiales,
>
> Jorge Velez.-
>
>
>
>
>
>
>

[[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] Extract entries from matrix

[Jorge I Velez]

Dear R-help,

I am working with a matrix "m" from which I would like to extract some
elements.  An toy example is as follows:

## input matrix
m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 0L,
0L, 0L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L), .Dim = c(22L, 5L))

R> m
# [,1]  [,2] [,3] [,4] [,5]
#  [1,]00000
#  [2,]00000
#  [3,]00000
#  [4,]12300
#  [5,]12300
#  [6,]12300
#  [7,]12300
#  [8,]12300
#  [9,]12340
# [10,]   12340
# [11,]   12345
# [12,]   12345

>From "m", I would like to extract the entries

4, 1
5, 2
6, 3
7, 1
8, 2
9, 3
10, 1
11, 2
12, 3

so at the end of applying a function "f" to "m" I get

1, 2, 3, 1, 2, 3, 4, 1, 2, 3


Basically the idea is to extract the diagonal elements until a zero is
found.

In the real problem the dimensions of "m" are much bigger, but this smaller
version of "m" illustrate what needs to be done.

I would greatly appreciate any ideas on how to do this.

Thanks in advance,
Jorge Velez.-

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Extract entries from matrix

[Jorge I Velez]

Dear Jim,

Thank you very much for your quick reply.

I am sorry for the confusion it may have caused, but I messed up the
indexes in my example.  I would like, from the following matrix "m"

## input
m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,
0L, 0L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 0L, 0L, 0L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 4L,
4L, 4L, 4L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 5L, 5L), .Dim = c(12L,
5L))

to obtain

 1 2 3 1 2 3 4 5 1

Sure using m[idx] will give the desired result.  The problem is that idx is
not known and needs to be determined from "m".  I would like to use
something like

extractDiagonals(m)
## [1]  1 2 3 1 2 3 4 5 1

I look forward to your reply.  Thanks in advance.

Best regards,
Jorge Velez.-



On Tue, Oct 27, 2015 at 2:31 PM, jim holtman <jholt...@gmail.com> wrote:

> If you want to use the numbers you gave a the index into the matrix, then
> you can create a matrix with the values and then index into 'm'.  I don't
> see a '4' in the output example you gave using your index values:
>
> > m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
> +  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 2L, 2L, 2L, 2L,
> +  2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 0L,
> +  0L, 0L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
> +  3L, 3L, 3L, 3L, 3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 4L, 4L, 4L,
> +  4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 0L, 0L, 0L, 0L, 0L,
> +  0L, 0L, 0L, 0L, 0L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
> +  5L), .Dim = c(22L, 5L))
> > # create index matrix
> > indx <- matrix(c(4, 1,
> +  5, 2,
> +  6, 3,
> +  7, 1,
> +  8, 2,
> +  9, 3,
> +  10, 1,
> +  11, 2,
> +  12, 3), ncol = 2, byrow = TRUE)
> >
> >
> > m
>   [,1] [,2] [,3] [,4] [,5]
>  [1,]00000
>  [2,]00000
>  [3,]00000
>  [4,]12300
>  [5,]12300
>  [6,]12300
>  [7,]12300
>  [8,]12300
>  [9,]12340
> [10,]12340
> [11,]12345
> [12,]12345
> [13,]12345
> [14,]12345
> [15,]12345
> [16,]12345
> [17,]12345
> [18,]12345
> [19,]12345
> [20,]12345
> [21,]12345
> [22,]12345
> > indx
>   [,1] [,2]
>  [1,]41
>  [2,]52
>  [3,]63
>  [4,]71
>  [5,]82
>  [6,]93
>  [7,]   101
>  [8,]   112
>  [9,]   12    3
> > m[indx]
> [1] 1 2 3 1 2 3 1 2 3
>
>
> Jim Holtman
> Data Munger Guru
>
> What is the problem that you are trying to solve?
> Tell me what you want to do, not how you want to do it.
>
> On Tue, Oct 27, 2015 at 2:43 PM, Jorge I Velez <jorgeivanve...@gmail.com>
> wrote:
>
>> Dear R-help,
>>
>> I am working with a matrix "m" from which I would like to extract some
>> elements.  An toy example is as follows:
>>
>> ## input matrix
>> m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 2L, 2L, 2L, 2L,
>> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 0L,
>> 0L, 0L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
>> 3L, 3L, 3L, 3L, 3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 4L, 4L, 4L,
>> 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 0L, 0L, 0L, 0L, 0L,
>> 0L, 0L, 0L, 0L, 0L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
>> 5L), .Dim = c(22L, 5L))
>>
>> R> m
>> # [,1]  [,2] [,3] [,4] [,5]
>> #  [1,]00000
>> #  [2,]00000
>> #  [3,]00000
>> #  [4,]12300
>> #  [5,]12300
>> #  [6,]12300
>> #  [7,]12300
>> #  [8,]12300
>> #  [9,]12340
>> # [10,]   12340
>> # [11,]   12345
>> # [12,]   12345
>>
>> >From "m", I would like to extract the entries
>>
>> 4, 1
>> 5, 2
>> 6, 3
>> 7, 1
>> 8, 2
>> 9, 3
>> 10, 1
>> 11, 2
>> 12, 3
>>
>> so at the end of applying a function "f" to "m" I get
>>
>> 1, 2, 3, 1, 2, 3, 4, 1, 2, 3
>>
>&g

Re: [R] Extract entries from matrix

[Jorge I Velez]

Thank you all for your solutions and comments.

As Dr. Carlson mentioned, we leave rows 1 to 3 out as they are all zeroes.
Then, the entries I need to select from m are


entry value

4,1  ---> 1
5,2  ---> 2
6,3  ---> 3
7,1  ---> 1
8,2  ---> 2
9,3  ---> 3
10,4  ---> 4
11,5  ---> 5
12,1  ---> 1

Note that the entry [7,4] is zero, so we start from the first column in the
7th row and then select entry [7,1] instead.  That's what I meant by  "...
the idea is to extract the diagonal elements until a zero is found."  I
should have said *entries* instead of  _diagonal elements_. I am sorry Dr.
Turner for the confusion.

Starting with m

R> m
#   [,1] [,2] [,3] [,4] [,5]
# [1,]00000
# [2,]00000
# [3,]00000
# [4,]12300
# [5,]12300
# [6,]12300
# [7,]12300
# [8,]12300
# [9,]12340
#[10,]12340
#[11,]12345
#[12,]12345

the first submatrix to work with is

# [4,]12300
# [5,]12300
# [6,]12300

from which the elements of interest are 1, 2, 3.  Note that the 7th row of
m is not included here because m[7, 5] = 0.

Further, the second submatrix is

# [7,]12300
# [8,]12300
# [9,]12340
#[10,]12340
#[11,]12345

and the corresponding elements are 1, 2, 3, 4, 5.

And the last matrix is

#[12,]12345

from which the position [12,1] is selected.

So, the resulting entries from this process are 1, 2, 3, 1, 2, 3, 4, 5, 1.

Thank you in advance for any additional insight you may provide.

Regards,
Jorge Velez.-



On Tue, Oct 27, 2015 at 4:06 PM, David L Carlson <dcarl...@tamu.edu> wrote:

> I don't see how you are getting the result you provide.
>
> > m
>   [,1] [,2] [,3] [,4] [,5]
>  [1,]00000
>  [2,]00000
>  [3,]00000
>  [4,]12300
>  [5,]12300
>  [6,]12300
>  [7,]12300
>  [8,]12300
>  [9,]12340
> [10,]12340
> [11,]12345
> [12,]12345
> > t(sapply(1:8, function(x) diag(m[x:12, ])))
>  [,1] [,2] [,3] [,4] [,5]
> [1,]00000
> [2,]00300
> [3,]02300
> [4,]12300
> [5,]12300
> [6,]12340
> [7,]12345
> [8,]12345
>
> These are all of the diagonals from the 1st through 8th rows. The first 3
> begin with 0 so we leave them out, but then we have 4th: 1, 2, 3; 5th: 1,
> 2, 3; 6th: 1, 2, 3, 4, etc so you must have some additional rule in mind to
> get your answer.
>
> -
> David L Carlson
> Department of Anthropology
> Texas A University
> College Station, TX 77840-4352
>
>
>
> -Original Message-
> From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Jorge I
> Velez
> Sent: Tuesday, October 27, 2015 2:44 PM
> To: jim holtman
> Cc: R-help
> Subject: Re: [R] Extract entries from matrix
>
> Dear Jim,
>
> Thank you very much for your quick reply.
>
> I am sorry for the confusion it may have caused, but I messed up the
> indexes in my example.  I would like, from the following matrix "m"
>
> ## input
> m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,
> 0L, 0L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 0L, 0L, 0L, 3L, 3L,
> 3L, 3L, 3L, 3L, 3L, 3L, 3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 4L,
> 4L, 4L, 4L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 5L, 5L), .Dim = c(12L,
> 5L))
>
> to obtain
>
>  1 2 3 1 2 3 4 5 1
>
> Sure using m[idx] will give the desired result.  The problem is that idx is
> not known and needs to be determined from "m".  I would like to use
> something like
>
> extractDiagonals(m)
> ## [1]  1 2 3 1 2 3 4 5 1
>
> I look forward to your reply.  Thanks in advance.
>
> Best regards,
> Jorge Velez.-
>
>
>
> On Tue, Oct 27, 2015 at 2:31 PM, jim holtman <jholt...@gmail.com> wrote:
>
> > If you want to use the numbers you gave a the index into the matrix, then
> > you can create a matrix with the values and then index into 'm'.  I don't
> > see a '4' in the output example you gave using your index values:
> >
> > > m <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L

Re: [R] Number of digits to display?

[Jorge I Velez]

Dear Judson,

Perhaps the following example would give you some ideas:

R> options(digits = 4)
R> rnorm(5)
#[1] -0.57089 -0.14759  0.05717 -0.04935  2.22123
R> options(digits = 3)
R> rnorm(5)
#[1]  0.789  0.616  0.156 -1.315 -1.090
R> options(digits = 2)
R> rnorm(5)
#[1] -1.04 -0.58 -0.50  0.30  0.60

​Best regards,
Jorge.-


On Mon, Oct 26, 2015 at 11:43 AM, Judson  wrote:

> How do I control the number of digits to display,
> say, in a matrix, without rounding or losing accuracy
> in subsequent calculations?
> round() of course reduces accuracy.
>
> Also, if all my results
> are really fractions, is there a way to display
> them as fractions of integers rather than
> decimal expressions?
>
> ... judson blake
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Update dataframe based on some conditions

[Jorge I Velez]

Dear R-help,

I am working on what it seems to be a simple problem, but after several
hours trying to come up with a solution, unfortunately I have not been able
to.

I would like to go from "datain" to "dataout", that is, create the NEWREF
variable according with some restrictions, and update the values for the
remaining variables in the original data set (which is way more bigger than
this example). The problem can be described as having products (coded as
REF) in stock. Here, the total nomber of units in stock are named TOENDREF
and those required for the customer are given by TIMEREF. The idea is to
use as many units of the previous REF as possible before using a new REF.

## input
datain <- structure(list(REF = c("999", "999", "999", "1099", "731", "731",
"731", "731", "1442", "1442", "1442", "1442"), TIMEREF = c(120,
240, 360, 30, 30, 60, 90, 120, 30, 60, 90, 120), TOENDREF = c(390,
270, 150, 480, 480, 450, 420, 390, 480, 450, 420, 390)), .Names = c("REF",
"TIMEREF", "TOENDREF"), row.names = c(NA, 12L), class = "data.frame")
datain

## output
dataout <- structure(list(REF = c(999L, 999L, 999L, 1099L, 731L, 731L,
731L,
731L, 1442L, 1442L, 1442L, 1442L), TIMEREF = c(120L, 240L, 360L,
30L, 30L, 60L, 90L, 120L, 30L, 60L, 90L, 120L), TOENDREF = c(390L,
270L, 150L, 120L, 90L, 30L, 420L, 300L, 270L, 210L, 120L, 0L),
NEWREF = c(999L, 999L, 999L, 999L, 999L, 999L, 731L, 731L,
731L, 731L, 731L, 731L)), .Names = c("REF", "TIMEREF", "TOENDREF",
"NEWREF"), row.names = c(NA, 12L), class = "data.frame")
dataout


I what follows I will try to explain what I want to accomplish:

* Example 1
Take rows 3 and 4 of "datain"

#REF TIMEREF TOENDREF
#3   999 360  150
#4  1099  30  480

As 150 units of REF 999 are available, we could substitute the 30 units of
REF 1099 with them. Hence, the 4th row of the _updated_ "datain" becomes

#REF TIMEREF TOENDREF NEWREF
#3   999 360  150  999
#4  1099  30  120  999

* Example 2
Now, let's take rows 3 to 8 of the _updated_ "datain":

#REF TIMEREF TOENDREF
#3   999 360  150
#4   999  30  120
#5   731  30  480
#6   731  60  450
#7   731  90  420
#8   731 120  390

In row 4, there 120 units available to be used. The number of units
required of REF 731 is 30, which can be easily covered by the remaining 120
units of REF 999. By doing so, the remaining units of REF 999 would then be
90.  Hence, the newly _updated_ "datain" becomes

#REF TIMEREF TOENDREF
#3   999 360  150
#4   999  30  120
#5   999  30   90
#6   999  60   30
#7   731  90  420
#8   731 120  300

Finally, the updated "datain" file after processing the remaining REF would
be

#REF TIMEREF TOENDREF
#9  731  30  270
#10 731  60  210
#11 731  90  120
#12 731 1200

Hopefully I have explained well what I would like to end up with.  If this
is not the case, I will be more than happy to provide more information.

Any help would be very much appreciated.  Thanks in advance.

Best regards,
Jorge Velez.-

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R-es] Crear variable con condiciones

[Jorge I Velez]

Muchas gracias Javier por tu respuesta.

Si. Para obtener "dataout" se utilizan filas anteriores de acuerdo con la
disponibilidad de la variable TOENDREF para cada valor de la variable REF.
Por ejemplo, las filas 3 y 4 de "datain" son

#REF TIMEREF TOENDREF
#3  999 360  150
#4 1099  30  480

En la fila 3, el valor de TOENDREF es 150. Esto indica que hay 150
unidades disponibles de esa referencia. Ahora, en la fila 4, TIMEREF es 30
para REF = 1099.  Como en esta fila TIMEREF es menor que TOENDREF para la
referencia anterior, entonces la nueva variable NEWREF debe ser 999 y no
1099.  El nuevo valor de TOENDREF en esta fila sera 150 - 30 = 120.  Esta
seria la fila 4 de "dataout":

REF TIMEREF TOENDREF NEWREF
#4  1099  30  120999

Para la fila 5 de "dataout", los recursos disponibles corresponden al
_nuevo_ valor de TOENDREF en NEWREF (i.e., 120).  Siguiendo la misma logica
anterior, obtenemos entonces las filas 5 a 12 de "dataout":

REF TIMEREF TOENDREF NEWREF
#5   731  30   90999
#6   731  60   30999
#7   731  90  420731
#8   731 120  300731
#9  1442  30  270731
#10 1442  60  210731
#11 1442  90  120731
#12 1442 1200731

Observa que en la ultima fila se agotaron todos los recursos de TOENDREF
para NEWREF = 731, por lo que no fue necesario utilizar la REF = 1442.

Espero que esta vez las cosas sean un poco mas claras.

Los datos se pueden agrupar por la variable REF, que basicamente se refiere
a la referencia de un producto.  Si aun tengo disponibilidad de ese
producto (variable TOENDREF) entonces lo utilizo y cancelo la referencia
siguiente.  Las unidades que se piden de cada producto corresponden a la
variable TIMEREF.

Gracias a todos de antemano por sus sugerencias.

Saludos,
Jorge
​Velez.-


2015-10-20 22:30 GMT-05:00 Javier Rubén Marcuzzi <
javier.ruben.marcu...@gmail.com>:

> Estimado Jorge I Velez
>
>
>
> Yo hace unos años tuve un problema parecido y creo que había dos posibles
> soluciones y una era aportada por usted.
>
>
>
> No alcanzo a entender correctamente, al procesar el ejemplo hay números
> que no me coinciden con lo que comenta, le sugiero lo siguiente, cree
> nuevamente el ejemplo, con un cambio, un data.frame donde está lo que ahora
> es datain y dataout, que entiendo que es lo que quiere y lo que desea, pero
> al describir el ejemplo ref 1099, indique las coordenadas entre las filas
> (que son 12) y las columnas.
>
>
>
> Yo no alcanzo a comprender correctamente, ¿usa una fila con valores de
> otra fila anterior?, O se me acomodaron mal los datos y me perdí.  El caso
> parecido que yo tuve que resolver, de acuerdo a ciertos valores tomaba los
> de filas anteriores para realizar el cálculo. También necesito saber si hay
> una identificación tipo id de una base de datos o son números que no se
> pueden agrupar. Le pregunto porque en mi caso que puede ser parecido, pude
> realizar agrupación y columnas auxiliares (para procesar las condiciones).
>
>
>
> Javier Rubén Marcuzzi
> Técnico en Industrias Lácteas
> Veterinario
>
>
>
>
>
>
> *De: *Jorge I Velez
> *Enviado: *martes, 20 de octubre de 2015 19:17
> *Para: *R-help-es
> *Asunto: *[R-es] Crear variable con condiciones
>
>
>
>
>
> Buenas tardes a todos,
>
>
>
> Quisiera crear una variable de acuerdo a ciertas condiciones.  Me gustaria
>
> llegar de "datain" a "dataout".
>
>
>
> ## entrada
>
> datain <- structure(list(REF = c("999", "999", "999", "1099", "731", "731",
>
> "731", "731", "1442", "1442", "1442", "1442"), TIMEREF = c(120,
>
> 240, 360, 30, 30, 60, 90, 120, 30, 60, 90, 120), TOENDREF = c(390,
>
> 270, 150, 480, 480, 450, 420, 390, 480, 450, 420, 390)), .Names = c("REF",
>
> "TIMEREF", "TOENDREF"), row.names = c(NA, 12L), class = "data.frame")
>
> datain
>
>
>
> ## salida
>
> dataout <- structure(list(REF = c(999L, 999L, 999L, 1099L, 731L, 731L,
>
> 731L,
>
> 731L, 1442L, 1442L, 1442L, 1442L), TIMEREF = c(120L, 240L, 360L,
>
> 30L, 30L, 60L, 90L, 120L, 30L, 60L, 90L, 120L), TOENDREF = c(390L,
>
> 270L, 150L, 120L, 90L, 30L, 420L, 300L, 270L, 210L, 120L, 0L),
>
> NEWREF = c(999L, 999L, 999L, 999L, 999L, 999L, 731L, 731L,
>
> 731L, 731L, 731L, 731L)), .Names = c("REF", "TIMEREF", "TOENDREF",
>
> "NEWREF"), row.names = c(NA, 12L), class = "data.frame")
>
> dataout
>
>
>
>
>
> A continuacion describo dos casos puntuales para ilu

[R-es] Crear variable con condiciones

[Jorge I Velez]

Buenas tardes a todos,

Quisiera crear una variable de acuerdo a ciertas condiciones.  Me gustaria
llegar de "datain" a "dataout".

## entrada
datain <- structure(list(REF = c("999", "999", "999", "1099", "731", "731",
"731", "731", "1442", "1442", "1442", "1442"), TIMEREF = c(120,
240, 360, 30, 30, 60, 90, 120, 30, 60, 90, 120), TOENDREF = c(390,
270, 150, 480, 480, 450, 420, 390, 480, 450, 420, 390)), .Names = c("REF",
"TIMEREF", "TOENDREF"), row.names = c(NA, 12L), class = "data.frame")
datain

## salida
dataout <- structure(list(REF = c(999L, 999L, 999L, 1099L, 731L, 731L,
731L,
731L, 1442L, 1442L, 1442L, 1442L), TIMEREF = c(120L, 240L, 360L,
30L, 30L, 60L, 90L, 120L, 30L, 60L, 90L, 120L), TOENDREF = c(390L,
270L, 150L, 120L, 90L, 30L, 420L, 300L, 270L, 210L, 120L, 0L),
NEWREF = c(999L, 999L, 999L, 999L, 999L, 999L, 731L, 731L,
731L, 731L, 731L, 731L)), .Names = c("REF", "TIMEREF", "TOENDREF",
"NEWREF"), row.names = c(NA, 12L), class = "data.frame")
dataout


A continuacion describo dos casos puntuales para ilustrar las condiciones
que deben satisfacerse:

* Ejemplo 1
En REF = '1099', el TIMEREF es inferior al valor TOENDREF para REF = 99
(i.e., 30 < 150). Por lo tanto, la nueva variable "NEWREF" debe tomar el
valor de '99'.  Al realizar esta asignacion, se tiene que la fila 4 de
"datain" se convierte en

#REF TIMEREF TOENDREF NEWREF
#4 1099  30  120999

Aqui, 120 es el valor de TOENDREF para REF = '99'.

* Ejemplo 2
El proceso continua reasignando el REF hasta que el valor resultante de
TOENDREF sea inferior (o igual) al valor de TIMEREF del siguiente valor de
REF como ocurre en la fila 6 de dataout

#REF TIMEREF TOENDREF NEWREF
#6 731  60   30999
#7 731  90  420731

Puesto que el valor TOENDREF de 30 en TIMEREF 99 no es suficiente para
cubrir TIMEREF en TIMEREF 731, en la fila 7 NEWREF sera igual a REF.

De antemano muchas gracias por la ayuda.

Saludos cordiales a todos,
Jorge Velez.-

[[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 fread {data.table}

[Jorge I Velez]

Hola Ma Luz,
Es dificil saber que salgo mal sin saber que hiciste.  Podrias por favor
enviarnos el codigo que utilizaste?
Saludos cordiales,
Jorge.-


2015-09-29 9:55 GMT-05:00 MªLuz Morales :

> Buenas tardes,
>
> intento almacenar el contenido de un archivo .docx en un data.table, pero
> solo me devuelve 1855 filas cuando deberñian ser 6821.
>
> Sin embargo, el mismo archivo se me descarga completo usando read_docx
> {qdapTools}, pero este devuelve un vector carácter y no es lo que quiero.
>
> ¿Alguien sabe donde puede estar el problema?
>
> Nota: El docx procede de un archivo .rel que descargué de internet, al cual
> cambié la extensión por .doc y una vez abierto guardé como docx.
>
> Gracias!!
> Un saludo
> MªLuz
>
> [[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

Re: [R-es] Una preguntilla

[Jorge I Velez]

Heber,

Dale una mirada a

?cbind
?file.choose

"cbind" es el acronimo de "column binding" o "pegado por columnas".
Básicamente toma una columna (o varias) de una matrix o data.frame y las
pega con otra(s).

La función fle.choose() abre la ventana que necesitas.  Dentro de la
funcion puedes usar read.table() o similares para leer los datos que el
usuario necesita incluir en el script.

Espero haber podido ayudar.

Saludos,
Jorge.-


2015-09-22 4:12 GMT+10:00 heber sarmiento via R-help-es <
r-help-es@r-project.org>:

> Alguien me puede ilustrar en relación con la manera en como funciona la
> orden cbind y adicional a esto me gustaría saber existe alguna orden que me
> despliegue una ventana para leer datos de un usuario en un scrip.
> De antemano gracias.
> Heber
> [[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

Re: [R-es] intercalar elementos de vectores

[Jorge I Velez]

Hola Carmen,

Si, hay varios paquetes.  De ellos, el mas famoso es sem:
https://cran.r-project.org/web/packages/sem/index.html

Hay varios tutoriales en internet.  Podrias comenzar con
http://socserv.socsci.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf
https://www.google.com/url?sa=trct=jq=esrc=ssource=webcd=1cad=rjauact=8ved=0CCsQFjAAahUKEwizq-rixLrHAhXKlR4KHczkA-Murl=http%3A%2F%2Fsocserv.socsci.mcmaster.ca%2Fjfox%2FMisc%2Fsem%2FSEM-paper.pdfei=9knXVbPYMsqreszJj5gOusg=AFQjCNFd1Wjy6WKJPEybeLNnOq5Mz-VEEwsig2=-BHZgieK0dD9zFuP5xrM0g
y
sus versiones mas recientes.

Saludos cordiales,

Jorge I. Vélez
JCSMR, Canberra
​


2015-08-22 1:54 GMT+10:00 carmen cecilia otero villadiego 
ccovilladi...@gmail.com:

 Hola buenos días,

 Les agradecería mucho si me pueden ayudar con la sgte consulta:

 R tiene paquetes donde podamos hacer análisis de ecuaciones estructurales o
 análisis path??

 si es así me pueden dar una guía???

 De antemano mil gracias,

 Saludos.

 El 24 de febrero de 2015, 7:49, Fernando Macedo ferm...@gmail.com
 escribió:

   Buenas a todos.
  Relato el problema:
 
  - tengo un archivo de 316 columnas por 562000 filas (aprox.).
  - esas 316 columnas representan 158 sujetos, o sea dos columnas por cada
  individuo conteniendo información que debe ser condensada en una sola.
 
  Lo que necesito es ir tomando las dos columnas de cada individuo e
  intercalar los elementos de los vectores formando uno solo.
 
  Ejemplificando sería algo así:
 
   a
  [1] a a a a a
   b
  [1] b b b b b
   c
   [1] a b a b a b a b a b
 
 
  Estoy haciendo con un loop for pero es realmente muy lento. He buscado
 por
  algún paquete que ya lo haga directamente pero no he tenido mucho éxito.
 Me
  imagino que con sapply o apply pueda ser mucho más efectivo pero me ha
  resultado complicado para entender la sintaxis de estas funciones cuando
  involucra más de un objeto (vector, matriz, etc...).
 
  Desde ya agradezco las sugerencias que puedan verter sobre este problema.
 
  --
  Fernando Macedo
 
  [[alternative HTML version deleted]]
 
  ___
  R-help-es mailing list
  R-help-es@r-project.org
  https://stat.ethz.ch/mailman/listinfo/r-help-es
 



 --
 *Cecilia Otero Villadiego*

 *Estadístico - Esp. mercadeo Gerencial*
 *Cel:300 661 92 16 - 301 280 74 16*

 [[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

Re: [R-es] función cantidad mayor de valores

[Jorge I Velez]

Hola Javier,

Intenta los siguientes cambios

g4 - function(codigo_llega, n_caracteres){
  codigo_llega - as.character(codigo_llega)
  if(n_caracteres == 6) res - substr(codigo_llega, start=0, stop=4)
  else res - codigo_llega
  res
}
g4 - Vectorize(g4)
x - c('Jorge Velez','Javier Marcuzzi','Daniel Merino','Susana deus
Alvarez', 'Carlos Ortega')
g4(x, 6)
g4(x, 4)

Saludos cordiales,
Jorge.-



2015-07-28 13:08 GMT+10:00 javier.ruben.marcu...@gmail.com:

 Señores



 Tengo un problema, donde use distintas alternativas y el informe de error
 es el mismo.



 Mi ultima alternativa es una función que copio y pego junto con un ejemplo
 de uso.

  g3 - function(n_caracteres, codigo_llega){
 +   if (n_caracteres == 6)
 +   {
 + res - substr(codigo_llega, start=0, stop=4)
 +   }
 +   else
 +   {
 + res - codigo_llega
 +   }
 +   res
 + }
  g3(6,Javier Marcuzzi)
 [1] Javi

 Pero cuándo uso los datos reales, el mensaje de error es el siguiente
 (copio y pego lo último que se ve junto con el mensaje)


 [9989] QS01EE01 QS01EE02 QS01EE03 QS01EE04 QS01EE05 QS   QS01
  QS01EQS01EX   QS01EX01 QS01EX02
 [1] QS
  [ reached getOption(max.print) -- omitted 846 entries ]
 7148 Levels: QA QA01 QA01A QA01AA QA01AA01 QA01AA02 QA01AA03 QA01AA04
 QA01AA30 QA01AA51 QA01AB ... QV10XX03
 Warning message:
 In if (n_caracteres == 4) { :
   the condition has length  1 and only the first element will be used


 ¿Alguna idea sobre ese mensaje, aparentemente hasta el elemento 10.000
 funciona, luego tengo un problema (con la función que envío como con otras
 alternativas)








 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


[[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] procesamiento de textos con R

[Jorge I Velez]

Hola Ma Luz,

He trabajado con las siguientes opciones:

stylo:  https://sites.google.com/site/computationalstylistics/stylo
tm:  http://cran.r-project.org/web/packages/tm/index.html

Espero te sirvan.

Saludos cordiales,

Jorge
​ Velez

JCSMR, Canberra


2015-07-07 19:14 GMT+10:00 MªLuz Morales mlzm...@gmail.com:

 Buenos días,

 quisiera saber si existe algún paquete en R para procesamiento de texto,
 búsqueda de similitudes y ese tipo de cosas. He estado buscando pero no he
 encontrado nada al respecto.

 Gracias
 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

Re: [R-es] pr

[Jorge I Velez]

​Hola a todos,

Tomando el ejemplo de Carlos, podrías construir una funcion que realice lo
que necesitas:

vacunar - function(pop.size, incid, durEst, ...){
resEst - est.R0.AR(pop.size=pop.size, incid=incid,...)
ro - resEst$R
vac - 1 - (1/ro)
vac
}

vacunar(pop.size = 500, incid = c(1,5,12,13,15,19,22,34,41,53,70), durEst =
7)
## [1] 0.3246206

Saludos cordiales,
Jorge.-


2015-06-28 2:30 GMT+10:00 Carlos Ortega c...@qualityexcellence.es:

 Hola,

 Puedes hacerlo de esta forma:

 #-
 library(R0)

 durEst -7 #7 dias estimacion duracion
 resEst - est.R0.AR(pop.size=500, incid=c(1,5,12,13,15,19,22,34,41,53,70))
 #el resultado es  1.48

 ro - resEst$R

 vac -1-(1/ro) #  (1-1/Ro)
 vac

 IT - ro/durEst   #  Ro/7
 IT
 #-

 est.R0.AR devuelve un objeto en el que uno de los elementos es R que
 almacena la estimación.

 Saludos,
 Carlos Ortega
 www.qualityexcellence.es


 El 27 de junio de 2015, 16:01, jbetancourt jbetanco...@iscmc.cmw.sld.cu
 escribió:

  Estimados
  No tengo experiencia en programación y  quisiera preparar el  siguiente
  script para los estudiantes
  Una supuesta enfermedad, al parecer de transmisión aérea y que suele
 durar
  7 días, entra en una localidad de 500 habitantes, la  incidencia diaria
 es
  c(1,5,12,13,15,19,22,34,41,53,70)
  Al calcular Ro de la manera siguiente:
 
  library(R0)
  est.R0.AR(pop.size=500, incid=c(1,5,12,13,15,19,22,34,41,53,70))  #el
  resultado es  1.48
 
  yo quisiera programar el rango de personas a vacunar (1-1/Ro) sin tener
  que escribir 1.48 sino que el script lo tome del resultado de la
  estimación  realizada en est.R0.AR, lo mismo para  calcular las personas
  a vacunar Ro/7
  vac -1-(1/1.48) #  (1-1/Ro)
  vac
  IT-1.48/7   #  Ro/7
  IT
  Mis saludos cordiales
  José
 
  [[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

 [[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

Re: [R-es] extraño error con seq

[Jorge I Velez]

Francisco,
Tiene que ver con aritmetica de punto flotante y como los numeros se
representan en R. En la FAQ 7.31 esta explicado.
Espero sea de utilidad.
Jorge.-

On Tuesday, June 2, 2015, Francisco Rodriguez Sanchez 
f.rodriguez.s...@gmail.com wrote:

 Estimados colisteros,

 Acabo de encontrar un extraño error usando la función seq:

 seq(from = 0.6 + 0.1, to = 0.7, by = 0.1)
 [1] 0.7
 # todo bien

 seq(from = 2.2 + 0.1, to = 2.3, by = 0.1)
 Error in seq.default(from = 2.2 + 0.1, to = 2.3, by = 0.1) :
wrong sign in 'by' argument


 La ayuda de seq establece que Specifying to - from and by of opposite
 signs is an error. Y en efecto:

 2.3 - (2.2 + 0.1)
 [1] -4.440892e-16

 arroja un resultado muy pequeño pero negativo, mientras que 'by' es
 positivo (0.1).

 He probado con otros muchos números y esto solo me ocurre con el caso de
 2.2 y 2.3. Por ejemplo:

 0.7 - (0.6 + 0.1)
 [1] 0


 A qué se debe esto? O más importante, cómo puedo evitar este error?

 Muchas gracias de antemano

 Paco

 --
 Dr Francisco Rodriguez-Sanchez
 Integrative Ecology Group
 Estacion Biologica de Doñana - CSIC
 Avda. Americo Vespucio s/n
 41092 Sevilla (Spain)
 http://bit.ly/frod_san


 [[alternative HTML version deleted]]

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



-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Count number in r

[Jorge I Velez]

Dear Arshad,

Here is a possibility using tapply():

with(d, tapply(Rain, list(Month, Year), function(x) sum(x  .1)))
##1971
#1   12
#20


where d is your data.frame().   See also ?aggregate and ?ave.

Best,
Jorge.-


On Tue, May 19, 2015 at 8:10 PM, Hafizuddin Arshad 
hafizuddinarsha...@gmail.com wrote:

 Dear R users,

 Could someone help me on this? I have this kind of data set:

 structure(list(Year = c(1971L, 1971L, 1971L, 1971L, 1971L, 1971L,
 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L,
 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L,
 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L,
 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L, 1971L,
 1971L, 1971L), Month = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
 2L, 2L, 2L), Rain = c(58.9, 74.6, 17.7, 7.8, 1.2, 1, 5.3, 0.7,
 1.2, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
 0, 10.4, 17.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)), .Names = c(Year,
 Month, Rain), class = data.frame, row.names = c(NA, -44L
 ))

 I want to count data in Rain that is greater and equal to 0.1 mm
 according to their Month and Year. I have used this code, but it seems
 so wrong.

 raindat - read.csv('my data set',header=TRUE)
 yearcorr-min(raindat$Year)-1
 years-unique(raindat$Year)
 rainmonth-as.data.frame(matrix(0,nrow=2,ncol=12))
 for(year in years) {
   for(month in 1:12) {
 if(any(raindat$Year==yearraindat$Month==month))
   rainmonth[year-yearcorr,month]-
   length((which(raindat$Rain =
 0.1))[raindat$Year==yearraindat$Month==month])
   }
 }
 rownames(rainmonth)-years
 names(rainmonth)-month.abb
 rainmonth

 Thank you so much.


 Arshad

 [[alternative HTML version deleted]]

 __
 R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide
 http://www.R-project.org/posting-guide.html
 and provide commented, minimal, self-contained, reproducible code.


[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Help with pedigree() function in kinship2

[Jorge I Velez]

Dear R-help,

I am interested in plotting some pedigrees and came across the kinship2
package. What follows is an example of the pedigrees I am working with.

Now, when running

## check package availability
if(!require(kinship2)) install.packages('kinship2')
require(kinship2)


## data to plot
d - structure(list(FamilyID = c(1, 1, 1, 1, 1, 1, 1,
1, 1), PatientID = structure(c(2L, 3L, 5L, 11L, 12L, 15L,
16L, 17L, 6L), .Label = c( 1,  2,  3,  4,  5,  6,
 7,  9, 10, 11, 13, 14, 18, 20, 23, 24, 25,
27, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39,
41, 43, 45, 50, 62, 63, 64, 65, 66, 67, 85,
88), class = factor), FatherID = structure(c(1L, 1L, 6L,
1L, 5L, 6L, 1L, 7L, 6L), .Label = c(0, 1, 10, 11, 13,
2, 23, 27, 28, 3, 33, 34, 35, 38, 5, 62,
64, 66, 9), class = factor), MotherID = structure(c(1L,
1L, 7L, 1L, 14L, 7L, 1L, 5L, 7L), .Label = c(0, 10, 18,
2, 24, 29, 3, 30, 33, 34, 39, 4, 43, 5,
6, 63, 65, 9), class = factor), Sex = structure(c(2L,
1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L), .Label = c(Female, Male), class =
factor),
AffectionStatus = structure(c(1L, 1L, 2L, 1L, 2L, 2L, 1L,
2L, 2L), .Label = c(1, 2), class = factor)), .Names =
c(FamilyID,
PatientID, FatherID, MotherID, Sex, AffectionStatus
), row.names = c(NA, 9L), class = data.frame)

## plotting
ped - with(d, pedigree(PatientID, FatherID, MotherID, Sex,  affected =
AffectionStatus, famid = FamilyID))

## Error in pedigree(PatientID, FatherID, MotherID, Sex, affected =
AffectionStatus,  :
##Value of 'dadid' not found in the id list 1/0 1/0 1/2
1/0 1/2

I get an error.  My sessionInfo() is at the end.  I was wondering if
someone could help me to dissect what the cause of this error is and how it
can be fixed.

Thank you very much for your help.

Best regards,
Jorge Velez.-


R sessionInfo()
R version 3.2.0 Patched (2015-05-04 r68320)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.10.3 (Yosemite)

locale:
[1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8

attached base packages:
[1] stats graphics  grDevices utils datasets  parallel  compiler
 methods
[9] base

other attached packages:
[1] kinship2_1.6.0 quadprog_1.5-5 Matrix_1.2-0   readxl_0.1.0

loaded via a namespace (and not attached):
[1] Rcpp_0.11.6 grid_3.2.0  lattice_0.20-31

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R-es] ubicación de paquete

[Jorge I Velez]

Estimado Ricardo,
Si entiendo bien, necesitas ubicar el paquete 'stats'.  Este paquete hace
parte de la instalacion regular de R.  No uso Windows regularmente, pero
creo que el paquete esta en C:\Program Files\R\R-3.1.3patched\library
(cambia R-3.1.3patched por tu version actual).
Espero te sirva.
Jorge.-

2015-04-30 10:07 GMT+10:00 Alva Valiente, Ricardo (RIAV) 
r...@cajatrujillo.com.pe:

 Estimados una �ltima consulta como para cerrar el d�a. Estoy buscando el
 paquete Stats el cual esta en la versi�n 2.15.3, pero en la nueva versi�n
 de R no lo encuentrodonde lo puedo ubicar?, o como puedo saber el nuevo
 nombre del paquete.

 Atte.
 Ricardo Alva Valiente

 Aviso Legal: La informaci�n de este correo electr�nico, as� como de sus
 archivos adjuntos, es confidencial y est� dirigida exclusivamente a �l o
 los destinatarios. Si Usted ha recibido este correo por error, por favor
 av�senos inmediatamente por este medio y elim�nelo de su sistema. Se
 encuentra prohibido cualquier uso, reproducci�n, divulgaci�n o distribuci�n
 por otras personas distintas de �l o los destinatarios. Cualquier opini�n
 emitida en este correo electr�nico es propia del autor o remitente y no
 representa necesariamente la opini�n de la Caja Trujillo. A pesar de
 esfuerzos razonables en el control de virus y programas maliciosos, la Caja
 Trujillo no puede asegurar que �stos no se encuentren en este correo por
 causas ajenas a su control, por lo que usted debe analizar este correo y
 sus archivos adjuntos antes de abrirlos. Caja Municipal de Ahorro y Cr�dito
 de Trujillo www.cajatrujillo.com.pe 

 [[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

Re: [R-es] Loop sobre muchos data frames

[Jorge I Velez]

Gracias, Oscar.

Al parecer corresponde a un vector de caracteres.  Puesto que lo
que necesitamos es la informacion asociado a cada uno de esos nombres,
creo que necesitas _leer_ la informacion que esta en cada uno de ellos.
Por ejemplo,

y - scan(bNames[1], what = 'character')
y

tendra la informacion contenida en el archivo qB001.txt.  Una vez con esa
informacion puedes construir el corpus como

TermDocumentMatrix(y, control = list(tokenize = nGram1Tok))

y continuar el proceso que tienes en mente.

Ten en cuenta que no puedo probar nada de lo que sugiero, asi que es solo
una sugerencia.

Saludos cordiales,
Jorge.-


2015-04-13 10:24 GMT+10:00 Oscar Benitez oscar.benitez1...@gmail.com:

 Jorge, amigos de R-help

 obtengo lo siguiente:


  str(bNames)
  chr [1:150] qB001 qB002 qB003 qB004 qB005 qB006 qB007 ...


 Gracias y saludos
 Oscar

 El 12 de abril de 2015, 20:15, Jorge I Velez jorgeivanve...@gmail.com
 escribi�:

 Hola Oscar,

 A lo mejor este equivocado, pero creo que necesitas

 read.table(bNames[[i]])

 en lugar de

 bNames[[i]]

 Como esta, bNames[[i]] corresponde al nombre del corupus, no al corpus
 como tal.

 Nos ayudaria muchisimo el resultado de

 str(bNames)

 Saludos cordiales,
 Jorge.-


 2015-04-13 9:25 GMT+10:00 Oscar Benitez oscar.benitez1...@gmail.com:

 Jorge, estimados colaboradores de R-help

 Estuve tratando de utilizar un script para uno de los pasos en mi
 an�lisis, que es transformar cada uno de los corpus en mi espacio de
 trabajo en un objeto TermDocumentMatrix

 Tengo un vector llamado bNames que lista todos los corpus que quiero
 pasar a TDM, y constru� los siguientes comandos:

 tdm.n1 - vector('list', length = length(bNames))

 for(i in seq_along(tdm.n1)){
   tdm.n1.[[i]] -
 TermDocumentMatrix(bNames[[i]],control=list(tokenize=nGram1Tok))
 }

 pero obtengo el siguiente error:

 Error in UseMethod(TermDocumentMatrix, x) :
   no applicable method for 'TermDocumentMatrix' applied to an object of
 class character

 entonces revis� qu� es lo que bNames[[i]] le est� pasando al for

 bNames[[1]]
 [1] qB001

 Una cadena de texto, no el corpus llamado qB001!!, por eso obtengo el
 error.
 Alguien se le ocurre c�mo hacer para que lea cada uno de los objetos en
 la lista y se los pase de a uno a la funci�n?
 Muchas gracias por adelantado
 Oscar






 El 10 de abril de 2015, 11:28, Jorge I Velez jorgeivanve...@gmail.com
 escribi�:

 Oscar,

 Cambia

 (l[i])

 por

 read.table(l[i])

 Olvide leer cada archivo en el codigo anterior.

 Saludos,
 Jorge.-

 2015-04-11 0:02 GMT+10:00 Oscar Benitez oscar.benitez1...@gmail.com:

 Jorge
 Gracias por el consejo. Aparentemente no lo estoy aplicando bien, pues
 el objeto que obtengo no contiene lo que quiero.
 Me explico, al ejecutar

 txt - vector('list', length = length(names)) #names el el vector
 donde ya ten�a almacenada la lista de txt's
 for(i in seq_along(txt)){
   txt[[i]] - Corpus(VectorSource(names[i]))
 }

 obtengo el objeto txt:
  class(txt)
 [1] list

 si extraigo solamente el primer objeto de esa lista:
  txt[[1]]
 VCorpus (documents: 1, metadata (corpus/indexed): 0/0)

 si quiero ver el contenido del primer copus

  inspect(txt[[1]])
 VCorpus (documents: 1, metadata (corpus/indexed): 0/0)

 [[1]]
 PlainTextDocument (metadata: 7)
 qB001.txt

 me informa cosas sobre el objeto, pero los datos no est�n all�...
 deber�a mostrarme algo as� como:

 inspect(cbD02[1:1]) #inspecciono el corpus cbD120, creado a mano por
 la sentencia cbD120-Corpus(VectorSource(qT120))

 #..contenido del corpus..
 I went to go see some gypsy to tell me my future, but she asked for my
 photograph instead.
 itz me the one who was talkin to u last time !!!
 starts at 4pm. come get sunny munny :)
 kind of scary to imagine what needs military wiping
 #.m�s.contenido del corpus..


 si quiero aplicarle un funci�n propia de limpieza de corpus, por
 ejemplo eliminar los n�meros presentes en el corpus
  tm_map(txt[[1]], removeNumbers)
 VCorpus (documents: 1, metadata (corpus/indexed): 0/0)

 no hace nada de nada...

 Saludos

 Oscar





 El 10 de abril de 2015, 1:15, Jorge I Velez jorgeivanve...@gmail.com
 escribi�:

 Oscar,

 Una forma de trabajar con este tipo de archivos es utilizando listas:

 # directorio del proyecto
 setwd('~/proyecto')

 # archivos de texto
 l - list.files(pattern = '.txt')

 # procesamiento
 txt - vector('list', length = length(l))
 for(i in seq_along(txt)){
 txt[[i]] - Corpus(VectorSource(l[i]))
 }

 # para acceder a la informacion del primero archivo, solo debes
 escribir
 txt[[1]]


 Espero sea de utilidad.

 Saludos,
 Jorge.-


 2015-04-10 14:14 GMT+10:00 Oscar Benitez oscar.benitez1...@gmail.com
 :

 Hola a todos!
 Estoy en un proyecto de text mining y por razones de los recursos
 con que
 cuento tuve que separar los archivos de texto de input del proyecto
 en
 muchos archivos peque�os.
 Luego de transformar cada uno de estos archivos en un corpus
 separado,
 puedo aplicar limpieza sobre cada corpus, buscar n

Re: [R-es] DUDA_R code_Modelos

[Jorge I Velez]

Lorena,
Solo por curiosidad, cual es la frecuencia y el porcentaje de ceros en tus
datos?
Saludos,
Jorge.-


2015-04-10 0:04 GMT+10:00 Lorena Tudela Marco lorenatudelama...@gmail.com:

 Muchas gracias por las r�pidas respuestas.

 Sin duda, es un gustazo contar con un foro de debate y apoyo como este! ;)

 Si, hab�a actualizaciones! as� que ya esta solucionado!

 Un abrazo y gracias! Que teng�is buen d�a,

 Lorena

 El 9 de abril de 2015, 14:20, Ruben Tobalina Ramirez 
 lagrimaescr...@gmail.com escribi�:

  Hola, a lo mejor es una tonter�a, pero el paquete ha tenido
  actualizaciones?
 
  No conozco que paquetes usas pero recuerdo un script de text minig con
  el paquete tm que me daba error tras unos meses guardado por una
  actualizaci�n.
 
  un saludo
 
  El d�a 9 de abril de 2015, 13:26, Marcelino de la Cruz
  marcelino.delac...@upm.es escribi�:
   puedes comprobar esto:
  
   range(as.numeric(Nijt))
  
  
   El valor m�nimo deber�a ser 0.
  
  
  
   El 09/04/2015 a las 13:09, Lorena Tudela Marco escribi�:
  
   Hola a todxs,
  
   Hace un tiempo con una base de datos count data y guarde el c�digo R
   utilizado.
   Ahora intento sacar de nuevo las estimaciones de los modelos, y me
   aparece el siguiente error:
  
   /Error in zeroinfl(as.numeric(Nijt) ~ X8 + X20 + SPAIN +
 United.Kingdom
   +  : /
   /  invalid dependent variable, minimum count is not zero/
   /
   /
   No consigo averiguar que esta fallando y poder reestimar los modelos.
  
   Adjunto el c�digo R al email, por si sirviese de ayuda.
  
   *Cualquier comentario o sugerencia es bienvenido!!!*
  
   Un abrazo grande y buen d�a,
  
   Lorena
  
  
   ___
   R-help-es mailing list
   R-help-es@r-project.org
   https://stat.ethz.ch/mailman/listinfo/r-help-es
  
  
   ___
   R-help-es mailing list
   R-help-es@r-project.org
   https://stat.ethz.ch/mailman/listinfo/r-help-es
 
  ___
  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

Re: [R-es] Dimensiones arreglo.-

[Jorge I Velez]

Freddy,
El comportamiento es de esperar y corresponde a un tema de dise�o.  La
explicacion, creo, esta en ?[
Saludos,
Jorge.-


2015-04-10 0:28 GMT+10:00 Freddy Omar L�pez Quintero 
freddy.vat...@gmail.com:

 Hola muchachos,

 Espero que est�n muy bien.

 Estoy trabajando con unos arreglos y encontr� un comportamiento que no
 conoc�a. Me pregunto si es algo de esperar.

 Es lo siguiente. Supongamos el siguiente arreglo:

 arreglo - array(runif(10*2*2), dim=c(10, 2, 2)) # dim: 10x2x2

 y que sobre �l selecciono las primeras 5 filas (o las primero 5 l�neas de
 la primera dimensi�n, no s�):

 arreglo[1:5, , ] # dim: 5x2x2

 Pero si ahora hago lo propio con un arreglo como el que sigue:

 arreglo - array(runif(10*1*2), dim=c(10, 1, 2))

 me encuentro con estas dimensiones:

 arreglo[1:5, ,]  # dim: 5x2

 Naturalmente, esto produc�a inconsistencias en mi c�digo. Ide� un
 condicional para el caso en el que la segunda dimensi�n es 1 y otro para
 dimensiones mayores.

 �Es esto de esperar?�Es una inconsistencia?�Hay una funci�n que me estoy
 perdiendo en estas situaciones?

 Gracias y disculpen las molestias.

 --
 �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


[[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-foro: Re: Bases de datos R

[Jorge I Velez]

Estimada Veronica,

El comando

R  ls()

deberia darte todos los objetos que tienes en el espacio de trabajo.  Si
quieres borrar un objeto en particular, por ejemplo el llamado datos,
solo debes escribir

R rm(datos)

donde rm() es la funcion remove.  Para mas ayuda consulta ?rm

Ahora, si quieres borrar TODOS los objetos de la sesion de R sin reiniciar
R, debes escribir

R  rm(list = ls())

Espero sea un poco mas claro.

Saludos cordiales,
Jorge.-



2015-04-08 15:58 GMT+10:00 pepeceb pepe...@yahoo.es:

 Hola Ver�nicasi le das la orden ls() te aparecer�n todos los objetos que
 has creadoLuego para eliminar el que quieras, basta con
 rm (objeto a eliminar)
 De todas formas, si vuelves a cargar una base de datos con el mismo
 nombre, elimina a la anterior y ya puedes trabajar con la nueva.


   De: Veronica Escorcia veronic...@hotmail.com
  Para: jose cebrian pepe...@yahoo.es
  Enviado: Mi�rcoles 8 de abril de 2015 4:30
  Asunto: R-foro: Re: Bases de datos R

 R-foro � Foros � Foro de debate � Bases de datos R
 |  | Re: Bases de datos Rde Veronica Escorcia - mi�rcoles, 8 de abril de
 2015, 03:56 |
 |   | Hola Manuel.Ya te expliqu� lo que m�s pude, no s� qu� m�s decirte
 :/ Me podr�as indicar por favor c�mo borrar todo para iniciar un nuevo
 proyecto?
 Es con Rcmdr, s�lo necesito saber c�mo eliminar las bases de datos que ya
 tengo cargadas y as� inicio de nuevo todo el proceso. Muchas gracias.
 Mostrar mensaje anterior | ResponderVer el mensaje en su contexto |

  Darse de baja de este foro  Dar de baja de todos los foros


 [[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

Re: [R-es] Hoja de Referencia (CheatSheet) para gglplot2 traducida al español...

[Jorge I Velez]

Feliciraciones a ambos, muy buen trabajo!  --JIV


2015-04-06 23:55 GMT+10:00 Carlos Ortega c...@qualityexcellence.es:

 Hola,

 Durante esta pasada Semana Santa, Santiago Mota y yo mismo hemos traducido
 al espa�ol la cheatsheet de ggplot2 que recientemente public� RStudio.

 Nos acaban de comunicar que la hoja de referencia ya est� disponible:

 http://www.rstudio.com/resources/cheatsheets/

 --
 Saludos,
 Carlos Ortega
 www.qualityexcellence.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

Re: [R-es] Informes periódicos con R

[Jorge I Velez]

Hola Jesus,

La idea es que tengas un documento markdown en RStudio y que el encabezado
sea similar a

---
title: Aqu� va tu t�tulo
author: Aqu� quien escribi� el documento
date: Y la fecha...
output: word_document
---

Despues de los segundos --- escribes el texto que necesitas siguiendo la
sintaxis clasica de markdown.

Una vez compilado el documento desde RStudio (tambien puede hacerse
directamente desde R con la funcion source) se generara el documento .docx
que contiene toda la informacion que necesitas.

Espero sea un poco mas claro.

Saludos,
Jorge.-



2015-04-07 21:50 GMT+10:00 Jesus Herranz jesus.herr...@imdea.org:

 Hola, Jorge

 Lo he mirado y parece que es una opci�n bastante buena, pero no logro
 saber c�mo se salva en Word, solo lo logro en html

 Muchas gracias

 Jes�s







 *De:* Jorge I Velez [mailto:jorgeivanve...@gmail.com]
 *Enviado el:* martes, 07 de abril de 2015 11:40
 *Para:* Jesus Herranz
 *CC:* R-help-es
 *Asunto:* Re: [R-es] Informes peri�dicos con R



 Hola Jesus,

 Una forma es RStudio + Markdown.  Hay infinidad de tutoriales en internet,
 pero la referencia basica es http://rmarkdown.rstudio.com/   Puedes
 organizar las salidas para que sean en Word, PDF o HTML.  La escongencia
 depende de lo que quieras hacer posteriormente con el informe generado.

 Espero sea de utiilidad.

 Saludos cordiales,

 Jorge.-





 2015-04-07 19:37 GMT+10:00 Jesus Herranz jesus.herr...@imdea.org:

 Hola



 Necesito elaborar un informe que se ejecutar� mensualmente, ya que los
 datos
 en los que se basa ir�n cambiando. El informe estar� en Word, y contiene
 texto fijo, tablas, gr�ficos y resultados obtenidos con R.

 El informe ser� ejecutado despu�s por una persona que no sabe nada de R, al
 que le suministrar� los scripts de R, pero deber�a ser bastante trasparente
 para �l.



 �Qu� paquetes me aconsej�is? y �alguna web con un tutorial?



 Prefiero algo sencillo, ya que no tengo mucho tiempo para dedicar a esto



 Gracias



 Jes�s






 [[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

Re: [R-es] BUCLE

[Jorge I Velez]

David:Por que no pruebas generando las 30 muestras _unicas_ de antemano
y sobre estas calcular lo que necesitas? --JIV


2015-03-25 16:22 GMT+11:00 David Contreras davidcontrera...@gmail.com:

 Buena noche a todos,

 Nuevamente requiero de su ayuda con algo puntual:

 1. Tengo en el vector Muestras:

   rm(list = ls())
   set.seed(2085089)# Fija Datos

 data-matrix(-1,30,ncol=10)

  MuestraS
  [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
 [1,]100011011 1

 Ahora, dentro de un bucle for voy a generar muestras aleatorias, pero
 requiero que en ese mismo proceso, apenas se tenga una muestra duplicada
 esta se sobrescriba, no he logrado que el bucle se quede en una posici�n y
 de esta manera generar una nueva muestra que no este duplicada.
 Para hace lo descrito en el parrafo anterior uso el siguiente c�digo:

   a-1
   for (a in 1:30) {
 data [a,]-sample(MuestraS,replace=T)
 if (sum(duplicated(data[1:a, ]))!=0){a-a}
 else {a-a+1}
   }

 No esta haciendo lo que requiero, se detecta el duplicado pero la muestra
 se generar y continua hasta llegar a 30 dejando el duplicado.
 No se si deba usar el bucle for o deberia cambiar por un while, pero no lo
 he podido lograr de ninguna de las dos formas.

 Un saludo, agradezco su colaboraci�n con el tema.


 David C.

 [[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

Re: [R-es] Combinatoria

[Jorge I Velez]

Hola Miguel,

Para que crear la funcion mCRn que propones si con

R  choose(4, 2)
[1] 6

obtienes el mismo resultado?

Aunque reconozco que son importantes, no se si en realidad las funciones
que mencionas son imprescindibles.

En cuanto al sistema de ecuaciones no lineales que planteas, dale una
mirada al paquete nleqslv.

Saludos cordiales,
Jorge.-


2015-03-23 17:30 GMT+11:00 Our Utopy ourut...@gmail.com:

 Hola Jorge, gracias por ayudarme antes de empezar. Ayer me fui a cama con
 tan solo escribir el t�tulo Ecuaciones no lineales

 B�sicamente lo que voy a buscar es lo cl�sico, es decir, resolver cualquier
 ecuaci�n no lineal, ya sea polin�mica o trascendente, con una o varias
 inc�gnitas.

 Ejemplo t�pico:

 x^2 + y^2 = 1
 y = sin x

 Si me van surgiendo dudas, aqu� estar�, no lo dud�is.

 Por cierto Carlos, veo que las funciones *combinations* y *permutations*
 del paquete *gtools*  hacen lo que yo buscaba, pero si estuviese en medio
 de un mar de c�lculos, lo que no me hace es calcular cuantas son, me veo
 obligado a buscar mis propias variables, por ejemplo el numero de
 combinaciones con repetici�n de tres banderas tomadas de dos en dos ser�a:

  mCRn - nrow(combinations(n=3, r=2, v=x, repeats=TRUE))

 Sin embargo, el paquete* combinat* s� que posee la funci�n *mCn*(4,2) que
 nos da 6

 No s�, lo que intento deciros es que me resulta muy chocante que un
 programa tan potente como R haya construido funciones variadas para otros
 temas y para combinatoria lo haya obviado, cuando las veo imprescindibles
 para probabilidad.

 De todas maneras me hab�is ayudado much�simo. Gracias.

 Miguel

 El 23 de marzo de 2015, 4:18, Jorge I Velez jorgeivanve...@gmail.com
 escribi�:

  Hola Miguel,
  Exactamente que necesitas hacer con ecuaciones no lineales?  Encontrar la
  raiz?  Si este es el caso, no necesitas, en principio, cargar ningun
  paquete.  Con la funcion ?optim en base podrias hacer lo basico.
  Saludos cordiales,
  Jorge.-
 
 
  2015-03-23 10:22 GMT+11:00 Our Utopy ourut...@gmail.com:
 
  Ya me extra�aba a mi!
  Ma�ana a primera hora miro todo lo que me comentas.
  Ahora estoy con ecuaciones no lineales y ya veo que hay que cargar
  librerias
  Muchas gracias. Un saludo
 
  El 23 de marzo de 2015, 0:08, Carlos Ortega c...@qualityexcellence.es
  escribi�:
 
   Hola,
  
   Por precisar un par de detalles:
  
  
  - Con el paquete gtools se pueden generar:
  - las variaciones, permutaciones, combinaciones, variaciones con
 repetici�n y combinaciones con repetici�n (mira el ejemplo
  adjunto con las
 combinaciones con repetici�n).
 - Quedan sin cubrir las permutaciones con repetici�n.
  
  
  - Y dentro del paquete base, s� que hay formas de calcular
 algunas
 
  cosas de combinatoria. Mira la funci�n choose() para el c�lculo
 de
  las
  combinaciones. En el ejemplo de la funci�n incluso aparece una
 forma
  de
  generar el tri�ngulo de Pascal.
  
   #- Ejemplo de Combinaciones con repetici�n (gtools) -
  
x - c('rojo', 'azul', 'verde')
# Combinaciones sin repetici�n
combinations(n=3, r=2, v=x, repeats=FALSE)
[,1]   [,2]
   [1,] azul rojo
   [2,] azul verde
   [3,] rojo verde
# Combinaciones con repetici�n
combinations(n=3, r=2, v=x, repeats=TRUE)
[,1][,2]
   [1,] azul  azul
   [2,] azul  rojo
   [3,] azul  verde
   [4,] rojo  rojo
   [5,] rojo  verde
   [6,] verde verde
  
  
 #---
  
   Saludos,
   Carlos Ortega
   www.qualityexcellence.es
  
   El 22 de marzo de 2015, 22:50, Our Utopy ourut...@gmail.com
 escribi�:
  
   Hola amigos, muchas gracias por vuestra ayuda.
  
   Entonces veo que mi sorpresa era leg�tima. Por todos vuestros mails
 la
   conclusi�n es que:
  
  - En el m�dulo base de R no incluye combinatoria elemental, ni
  siquiera
  el n�mero combinatorio Cm,n hay que cargar el paquete *combinat*
  - Y para las variaciones con repetici�n el paquete* gtools*
  - Y a�n as� no tenemos ni las combinaciones ni las permutaciones,
  ambas
  con repetici�n.
  
  
   Pues tendremos que ponernos a ello �no cre�is? y hacer un paquete que
  lo
   resuelva y que genere incluso un tri�ngulo de Pascal y otras
 cuestiones
   del
   binomio de Newton.
   Ya s� que es matem�tica elemental pero repito que estoy sorprendido
 por
   este hecho y m�s porque estas *T�cnicas de conteo* son
 imprescindibles
   para
   el c�lculo b�sico de probabilidades.
  
   Seguimos en la lucha. Un saludo
  
  
   El 22 de marzo de 2015, 19:32, Francisco Rodr�guez 
 fjr...@hotmail.com
  
   escribi�:
  
Carlos, creo que el ejemplo que mandas se refiere a las Variaciones
  con
Repetici�n de 3 elementos tomados de 3 en 3 y cuyo n�mero asciende
 a
  27
   
En las Permutaciones con repetici�n con las observaciones se forman
   grupos
disjuntos cuya suma de cardinales es el n�mero total, as� pues
 puede
hablarse de

Re: [R-es] Combinatoria

[Jorge I Velez]

Hola Miguel,
Exactamente que necesitas hacer con ecuaciones no lineales?  Encontrar la
raiz?  Si este es el caso, no necesitas, en principio, cargar ningun
paquete.  Con la funcion ?optim en base podrias hacer lo basico.
Saludos cordiales,
Jorge.-


2015-03-23 10:22 GMT+11:00 Our Utopy ourut...@gmail.com:

 Ya me extra�aba a mi!
 Ma�ana a primera hora miro todo lo que me comentas.
 Ahora estoy con ecuaciones no lineales y ya veo que hay que cargar
 librerias
 Muchas gracias. Un saludo

 El 23 de marzo de 2015, 0:08, Carlos Ortega c...@qualityexcellence.es
 escribi�:

  Hola,
 
  Por precisar un par de detalles:
 
 
 - Con el paquete gtools se pueden generar:
 - las variaciones, permutaciones, combinaciones, variaciones con
repetici�n y combinaciones con repetici�n (mira el ejemplo adjunto
 con las
combinaciones con repetici�n).
- Quedan sin cubrir las permutaciones con repetici�n.
 
 
 - Y dentro del paquete base, s� que hay formas de calcular algunas
 cosas de combinatoria. Mira la funci�n choose() para el c�lculo de
 las
 combinaciones. En el ejemplo de la funci�n incluso aparece una forma
 de
 generar el tri�ngulo de Pascal.
 
  #- Ejemplo de Combinaciones con repetici�n (gtools) -
 
   x - c('rojo', 'azul', 'verde')
   # Combinaciones sin repetici�n
   combinations(n=3, r=2, v=x, repeats=FALSE)
   [,1]   [,2]
  [1,] azul rojo
  [2,] azul verde
  [3,] rojo verde
   # Combinaciones con repetici�n
   combinations(n=3, r=2, v=x, repeats=TRUE)
   [,1][,2]
  [1,] azul  azul
  [2,] azul  rojo
  [3,] azul  verde
  [4,] rojo  rojo
  [5,] rojo  verde
  [6,] verde verde
 
  #---
 
  Saludos,
  Carlos Ortega
  www.qualityexcellence.es
 
  El 22 de marzo de 2015, 22:50, Our Utopy ourut...@gmail.com escribi�:
 
  Hola amigos, muchas gracias por vuestra ayuda.
 
  Entonces veo que mi sorpresa era leg�tima. Por todos vuestros mails la
  conclusi�n es que:
 
 - En el m�dulo base de R no incluye combinatoria elemental, ni
 siquiera
 el n�mero combinatorio Cm,n hay que cargar el paquete *combinat*
 - Y para las variaciones con repetici�n el paquete* gtools*
 - Y a�n as� no tenemos ni las combinaciones ni las permutaciones,
 ambas
 con repetici�n.
 
 
  Pues tendremos que ponernos a ello �no cre�is? y hacer un paquete que lo
  resuelva y que genere incluso un tri�ngulo de Pascal y otras cuestiones
  del
  binomio de Newton.
  Ya s� que es matem�tica elemental pero repito que estoy sorprendido por
  este hecho y m�s porque estas *T�cnicas de conteo* son imprescindibles
  para
  el c�lculo b�sico de probabilidades.
 
  Seguimos en la lucha. Un saludo
 
 
  El 22 de marzo de 2015, 19:32, Francisco Rodr�guez fjr...@hotmail.com
  escribi�:
 
   Carlos, creo que el ejemplo que mandas se refiere a las Variaciones
 con
   Repetici�n de 3 elementos tomados de 3 en 3 y cuyo n�mero asciende a
 27
  
   En las Permutaciones con repetici�n con las observaciones se forman
  grupos
   disjuntos cuya suma de cardinales es el n�mero total, as� pues puede
   hablarse de Permutaciones con Repetici�n de 10 elementos tomados de 5
  en 5,
   2 en 2 y de 3 en 3, en este caso, el n�mero de combinaciones que sale
   seria: PR5;3,2,3 = 10!/(5!�2!�3!)
  
   Este caso es el que digo que seria mas dificil de preparar y quizas si
   requiere algo de programacion, porque el numero de grupos es variable,
  pero
   a lo mejor en lo que mandas se puede introducir algun tipo de vector,
 en
   todo caso lo miro ma�ana
  
  
   Un saludo
  
   --
   Date: Sun, 22 Mar 2015 13:22:36 +0100
   Subject: Re: [R-es] Combinatoria
   From: c...@qualityexcellence.es
   To: fjr...@hotmail.com
   CC: ourut...@gmail.com; r-help-es@r-project.org
  
  
   S�, tambi�n...
   Para las permutaciones, n=r.
   Y con el par�metro repeats.allowed controlas si son con o sin
  repetci�n:
  
   #--
#Permutaciones *con repetici�n*
 
permutations(n=3, r=3, v=x, repeats.allowed=TRUE)
 [,1][,2][,3]
[1,] azul  azul  azul
[2,] azul  azul  rojo
[3,] azul  azul  verde
[4,] azul  rojo  azul
[5,] azul  rojo  rojo
[6,] azul  rojo  verde
[7,] azul  verde azul
[8,] azul  verde rojo
[9,] azul  verde verde
   [10,] rojo  azul  azul
   [11,] rojo  azul  rojo
   [12,] rojo  azul  verde
   [13,] rojo  rojo  azul
   [14,] rojo  rojo  rojo
   [15,] rojo  rojo  verde
   [16,] rojo  verde azul
   [17,] rojo  verde rojo
   [18,] rojo  verde verde
   [19,] verde azul  azul
   [20,] verde azul  rojo
   [21,] verde azul  verde
   [22,] verde rojo  azul
   [23,] verde rojo  rojo
   [24,] verde rojo  verde
   [25,] verde verde azul
   [26,] verde verde rojo
   [27,] verde verde verde
#Permutaciones *sin repetici�n*
 
permutations(n=3, r=3, v=x, repeats.allowed=FALSE)
[,1][,2][,3]
   [1,] azul  rojo  verde
   [2,] azul  

Re: [R-es] Familia *pply

[Jorge I Velez]

Hola Fernando,

No puedo ver las negritas, pero intuyo que lo que quieres es calcular la
media por columnas?  Si es asi, hay dos maneras:

1.  Usa colMeans(x), donde x es tu matriz de datos
2.  Usa apply(x, 2, mean) donde x es tu matriz de datos

Existe una tercera pero menos conocida posibilidad que es usando el paquete
matrixStats.  Esta implementado en C en su mayoria y, de acuerdo con el
autor, es mucho mas rapido que la familia *apply.  En
http://cran.r-project.org/web/packages/matrixStats/vignettes/matrixStats-methods.html
puedes encontrar mas informacion.

Saludos cordiales,
Jorge.-



2015-03-19 11:01 GMT+11:00 Fernando Macedo ferm...@gmail.com:

 Buenas a todos. Desde hace un tiempo estoy tratando de aplicar las
 funciones de la familia *pply en todo lo que puedo, pero todav�a no es
 algo que me surja tan r�pidamente o naturalmente al momento de los loops
 como usar for().
 Conozco las ventajas de usar estas funciones y por eso mi intento de
 hacerme de ellas.

 Por ejemplo en este problema:

 data=matrix(rnorm(100*20),20,100)
 col=sample(1:100,100)

 t1=Sys.time()

 medias=replicate(1000,{
sel=sample(1:20,10)
pareja=sample(sel,100,replace = T)
ta=Sys.time()
 *recep=NULL**
 **  for(i in 1:100){**
 **n=col[i]**
 **m=pareja[i]**
 **c=data[m,n]**
 **recep=c(recep,c)**
 **  }**
 *  tb=Sys.time()
media=mean(recep)
tt=tb-ta
c(media,tt)
 })

 t2=Sys.time()

 diftime=(t2-t1)[[1]]

 sum(medias[2,])/diftime


 la parte que est� en negrita (si us� bien los Sys.time()) me representa
 (hice varias pruebas) aprox un 60% del tiempo total empleado.

 Mi pregunta es, para este ejemplo �c�mo plantear�an una soluci�n usando
 funciones *pply?
 Y luego ver cuanto aumenta en el rendimiento del uso del tiempo.

 De paso, la salida que obtengo es una matriz de 2 por 1000, cuando ser�a
 m�s interesante una matriz de 1000 por 2. Si se usa simplify = F como
 argumento de replicate() resulta en una lista. �Existe alg�n argumento
 que directamente obtenga una matriz de 1000 por 2? (Esto �ltimo pensando
 en de repente 10 o 100 de repeticiones y salidas m�s complejas).


 Saludos!

 --
 Fernando Macedo


 [[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

Re: [R-es] ELIMINAR FILAS DUPLICADAS DENTRO DE UNA MATRIZ

[Jorge I Velez]

?duplicated   --JIV

On Wednesday, March 18, 2015, David Contreras davidcontrera...@gmail.com
wrote:

 Buen d�a,

 Necesito quitar filas duplicadas dentro de una matriz, estaba usando la
 funci�n unique pero �sta no lleva a cabo el proceso que requiero, para
 ilustrar un poco m�s lo que requiero, por ejemplo:

 Tengo la matriz A, que contiene 3 filas y 8 columnas:

 0 1 0 1 0 1 0 1
 1 1 1 1 1 1 0 0
 0 1 0 1 0 1 0 1

 Necesitar�a que la misma matriz dejando solo por ejemplo las dos primeras
 filas ya que la fila 1 y 3 est�n duplicadas.

 Agradezco su colaboraci�n.

 Saludos,

 David C0

 [[alternative HTML version deleted]]

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



-- 
Sent from my phone. Please excuse my brevity and misspelling.

[[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] Calculos matrices diferentes tamaños

[Jorge I Velez]

Hola David,

Gracias por la aclaracion.  Confio en que esta vez si se obtiene lo que
buscas:

R set.seed(123)
R A - matrix(rpois(20, 4), ncol = 4)
R B - matrix(rpois(274*3, 4), ncol = 3)
R A
# [,1] [,2] [,3] [,4]
#[1,]3187
#[2,]6443
#[3,]3751
#[4,]6443
#[5,]7428

R # el resultado es una matriz de orden 5x3
R sapply(1:ncol(B), function(.col) sapply(A[, .col + 1], function(x)
sum(B[, .col]  x)))
#[,1] [,2] [,3]
#[1,]  2569   13
#[2,]   98  102  157
#[3,]   14   59  244
#[4,]   98  102  157
#[5,]   98  2066

Ahora hagamos los calculos para los elementos de A:

R sum(B[,1]  A[1,2])
# [1] 256
R sum(B[,1]  A[2,2])
# [1] 98
R sum(B[,1]  A[3,2])
# [1] 14
R sum(B[,1]  A[4,2])
# [1] 98
R sum(B[,1]  A[5,2])
# [1] 98
R sum(B[,2]  A[1,3])
# [1] 9
R sum(B[,2]  A[2,3])
# [1] 102
R sum(B[,2]  A[3,3])
# [1] 59
R sum(B[,2]  A[4,3])
# [1] 102
R sum(B[,2]  A[5,3])
# [1] 206
R sum(B[,3]  A[1,4])
# [1] 13
R sum(B[,3]  A[2,4])
# [1] 157
R sum(B[,3]  A[3,4])
# [1] 244
R sum(B[,3]  A[4,4])
# [1] 157
R sum(B[,3]  A[5,4])
# [1] 6

Saludos,
Jorge.-



2015-03-10 1:40 GMT+11:00 David Contreras davidcontrera...@gmail.com:

 Buen D�a Jorge,

 Es algo parecido a lo que necesito, pero no debo tener en cuenta la suma
 de la columna sino por ejemplo:

 Tengo el elemento [1,2] en la matriz A que contiene el n�mero 5 y por otra
 parte tengo la columna [,2] de la Matriz B con valores por ejemplo
 (1,2,3,4,5,6,7,8,9,10,11,12). El resultado deber�a ser un vector que en la
 posici�n uno debe indicar que existen 7 valores mayores a 5 en la columna
 [,2] de la matriz B en comparaci�n con el elemento [1,2] e la Matriz A.

 Espero ser muy claro ya que a veces se torna un poco complicado la
 explicaci�n de lo que se quiere programar.

 Saludos,

 David Contreras.

 El 9 de marzo de 2015, 1:12, Jorge I Velez jorgeivanve...@gmail.com
 escribi�:

 Hola David,

 Si entiendo correctamente, lo siguiente deberia hacer lo que describes:

 set.seed(123)
 A - matrix(rpois(20, 4), ncol = 4)
 B - matrix(rpois(274*3, 4), ncol = 3)

 # el resultado es una matriz de orden 5x3
 sapply(1:ncol(B), function(.col) sapply(A[, .col + 1], function(x) sum(x
  B[, .col])))

 # [,1] [,2] [,3]
 #[1,]5  253  243
 #[2,]  112  119   69
 #[3,]  243  1643
 #[4,]  112  119   69
 #[5,]  112   18  260

 Este resultado indica que el elemento A[1, 2] es 5 veces mayor que todos
 los elementos de B[,1]; A[2, 2] es 112 veces mayor; A[3,2] es 243
 veces mayor; A[4,2] es 112 veces mayor y finalmente que A[5,2] es 112
 veces mayor que B[, 1].  Similarmente, A[5, 4] es 260 veces mayor que
 los elementos de B[,3].

 Un por uno seria algo como

 R sapply(A[,2], function(x) sum(x  B[,1]))
 [1]   5 112 243 112 112
 R sapply(A[,3], function(x) sum(x  B[,2]))
 [1] 253 119 164 119  18
 R sapply(A[,4], function(x) sum(x  B[,3]))
 [1] 243  69   3  69 260

 Saludos cordiales,
 Jorge.-


 2015-03-09 16:04 GMT+11:00 David Contreras davidcontrera...@gmail.com:

 Buen d�a amigos,

 Antes que nada, reitero mis agradecimientos por sus respuestas tan
 precisas
 y oportunas, sigo en mi trabajo de grado y tengo dificultad con un
 proceso:

 Tengo dos matrices A y B, A es de dimensi�n 5x4 y B es de dimensi�n
 274x3.
 Requiero de un proceso en el que tome el primer elemento de la Matriz A
 columna 2 y lo compare con todos los elementos de la Matriz B columna 1 y
 me diga en un vector de dimensi�n (1x5) cuantas veces fue el elemento de
 A
 mayor que cada uno de los 274 elementos de B y asi hasta llegar a los 5
 elementos de A, es decir el resultado de este proceso sera un vector con
 5
 elementos y en cada uno me dice cuentas veces fue mayor el elemento [i,j]
 mayor que los 274 elementos de B.

 Agradezco su colaboraci�n con el tema que como siempre es muy valiosa.


 Saludos,

 David Contreras.

 [[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

Re: [R] How to access https page

[Jorge I Velez]

Hi Hui,

I have used the source_url function in the devtools package with good
results.  Give it a shot!

Best,
Jorge.-


On Tue, Mar 10, 2015 at 9:39 AM, Hui Du hui...@savvyrookies.com wrote:

 Hi All,

 I am trying to parse some information from website, say, a linkedin page.
 The linkedin url was

 url = http://www.linkedin.com/in/huidu;

 I had no problem to use readLines and XML package to collect the
 information I need. However, that url became 
 https://www.linkedin.com/in/huidu; now.

 url = https://www.linkedin.com/in/huidu;

 It failed readLines function.

  readLines(url)
 Error in file(con, r) : cannot open the connection
 In addition: Warning message:
 In file(con, r) : unsupported URL scheme


 Do you know any way to read-in web information if the url is https? Thanks
 a lot.

 Hui

 [[alternative HTML version deleted]]

 __
 R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide
 http://www.R-project.org/posting-guide.html
 and provide commented, minimal, self-contained, reproducible code.


[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R-es] nueva en R

[Jorge I Velez]

Hola, Gloria.  Bienvenida!

Hay algunos manuales en la parte inferior de
http://cran.r-project.org/other-docs.html

Saludos cordiales,
Jorge.-


2015-02-26 22:36 GMT+11:00 Gloria Perez Fuentes fuper...@gmail.com:

 hola!soy nueva aqui, estoy haciendo un trabajo sobre R y toda la
 informacion que encuentro esta en ingles. me gustaria saber algun libro
 sobre R o donde podria conseguir buena informacion sobre el programa en
 espa�ol, ya sea de su funcionamiento, finalidad del programa,
 caracteristicas etc. Nunca he utilizado este programa. 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


[[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] CÀLCULOS DE PRUEBAS ESTADÍSTICAS

[Jorge I Velez]

David,

En el primer caso puedes consultar

?wilcox.test

Para la prueba de rachas, dale una mirada al ejemplo 6.2 en
http://knuth.uca.es/repos/ebrcmdr/pdf/actual/ebrcmdr.c6.pdf

Saludos cordiales,
Jorge.-


2015-02-26 16:47 GMT+11:00 David Contreras davidcontrera...@gmail.com:

 Buen d�a comunidad,

 Sigo aqu� en mi trabajo de grado en horarios extras, me pueden ayudar
 indicando como se calcula el test de los signos y el de rachas en R.

 Agradezco nuevamente su colaboraci�n que ha sido muy valiosa.

 Saludos,

 DC.

 [[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

Re: [R-es] como eliminar los nombres de las series en el grafico

[Jorge I Velez]

Eric,

Creo que necesitas Hmisc:::xYplot(..., label.curves = FALSE)

Saludos,
Jorge.-


2015-02-26 11:15 GMT+11:00 eric ericconchamu...@gmail.com:

 Estimada comunidad, estoy haciendo un grafico con la funcion xYplot
 (paquete Hmisc) para agregar barras de error a los puntos del grafico.
 Pero pasa algo que no me gusta: aparece el nombre de las series en el
 grafico.

 Alguien sabe como puedo eliminarlas ? no quiero usar ggplot porque ya
 tengo hechos todos los graficos con lattice.

 Aqui el codigo que uso y los datos (estan separados por tab, no por ,)
 y el grafico los envio adjuntos.

 xYplot(Cbind(ave,ul,ll) ~ con
 , groups=sol,data=ca.med.sincon ,xlab=Solvent concentration (mM)
 , ylab=Contact Angle (�), method=lower bars,  cex=1.1, type=b
 , scales=list(x=list(log=10), equispaced.log=FALSE)
 , abline=list(h=ca.med[ca.med$sol==con,3])
 )

 Espero que alguien sepa porque ya he gastado mucho tiempo buscando y
 probando cosas.

 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



[[alternative HTML version deleted]]

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