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] poner un código en Internet

[Jorge I Velez]

Dropbox?—JIV

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

> Buenas tardes (de nuevo), ¿sabéis si hay alguna forma sencilla de poner un
> código en Internet, al que pueda acceder cualquiera desde un link?
> 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] 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",
> >>  pack

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]

ro hacer,
> : 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&utm_campaign=sig-email&utm_content=webmail>
> Libre
> de virus. www.avast.com
>  source=link&utm_campaign=sig-email&utm_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 list
R-help-es@r-project.org
https://sta

Re: [R-es] Chi cuadrado ajustado

[Jorge I Velez]

Sebastián,

Una forma es vía Regresión Logística.  En https://stats.idre.ucla.
edu/r/dae/logit-regression/ podrías encontrar algunas ideas.

Saludos,
Jorge.-


2017-08-29 11:48 GMT-05:00 Sebastian Gadea :

> 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 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] 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] ezANOVA

[Jorge I Velez]

Buenas tardes, Miguel.

Basta con que instales el paquete "quantreg" y luego el paquete que quieres
usar.  Seria algo como

install.package('quantreg')
require(quantreg)
library(devtools)
dev_mode()
install_github('mike-lawrence/ez')
library(ez)

Espero sea de utilidad.

Saludos,
Jorge.-


2017-01-29 13:55 GMT-05:00 Miguel Lázaro via R-help-es <
r-help-es@r-project.org>:

> Hola a todos/as
> después de algún tiempo si trabajar con R, estoy con la idea de probar
> ezANOVA pero no soy capaz de descargarlo. He mirado y he seguido distintos
> métodos según he leído en diversos foros pero nada me funciona. Por
> ejemplo
> library(devtools)
> dev_mode()
> install_github('mike-lawrence/ez')
> library(ez)
>
> pero, cuando después de todo hago
>
> library(ez)
>
> la respuesta es
>
> Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck
> = vI[[j]]) :
>   there is no package called ‘quantreg’
> Error: package or namespace load failed for ‘ez’
>
> Imagino que será cosa de compatibilidades, versiones...
> He probado con R 3.3.2, que es el último. Agradecería alguna orientación
> al respecto.
>
> 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] 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)))])]
> > },
> > Olivier2={
> >   dat <-
> > data.table( Uno= sample( c(runif(numero) , rep(NA , numero
> /2e0
> > )) 

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] Chaid

[Jorge I Velez]

Do
factor(Survived)
when calling chaid(...)
Best,
Jorge.-

On Sunday, August 21, 2016, MIKE DE LA HOZ  wrote:

> Hi,
>
>
> I am running a chaid tree using titanic dataset (see attachment)
>
>
>
> setwd("C:/Users/miguel")
>
> titanic <- read.csv("train.csv")
> titanic.s <- subset( titanic, select = -c(PassengerId, Name ) )
>
> ctrl <- chaid_control(minsplit = 20, minbucket = 5, minprob = 0)
> chaidTitanic <- chaid(Survived ~ ., data = titanic, control = ctrl)
>
>
>
> It looks like I get the following error
>
> Error: is.factor(x) is not TRUE
>
>
> can you please help me here? I am not able to follow this type of error.
>
>
> Thanks
>
>
>

-- 
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 nombres a una matriz

[Jorge I Velez]

Luis,
Dale una mirada a ?within y ?transform
Saludos,
Jorge.-


2016-08-03 12:13 GMT-05:00 LAE :

>
> Hola a todos,
>
>
> Estoy teniendo problemas para cambiar el nombre de los factores.
>
> Mi matriz original es de la siguiente forma:
> > alfa
>
>local mes amb Fase sp1   sp2 sp3 sp4
> L1-P1-C Ago.10 Lentico LFP 22.111664 0 0 0
> L3-P3-M Ago.10 Lentico LFP 22.111664 5.527916 5.527916 0
> P-1-C Ago.10 Lotico LFP 11.055832 0 0 0
> P-3-M Ago.10 Lotico LFP 55.27916 5.527916 11.055832 0
> L1-P1-C Dic.10 Lentico LFP 5.527916 0 0 0
> L3-P3-M Dic.10 Lentico LFP 0 0 0 0
> P-1-C Dic.10 Lotico LFP 0 16.583748 11.055832 0
> P-3-M Dic.10 Lotico LFP 5.527916 0 0 0
> L1-P1-C Mar.12 Lentico LFP 55.27916 0 0 0
> L3-P3-M Mar.12 Lentico LFP 71.862908 0 0 0
> P-1-C Mar.12 Lotico LFP 5.527916 0 11.055832 0
> P-3-M Mar.12 Lotico LFP 5.527916 0 5.527916 0
> L1-P1-C Dic.12 Lentico LFP 60.807076 5.527916 16.583748 0
> L3-P3-M Dic.12 Lentico LFP 0 0 0 0
> P-1-C Dic.12 Lotico LFP 0 0 0 0
> P-3-M Dic.12 Lotico LFP 11.055832 0 5.527916 0
> L1-P1-C Dic.13 Lentico LFP 93.974572 0 0 0
> L3-P3-M Dic.13 Lentico LFP 55.27916 0 16.583748 0
> P-1-C Dic.13 Lotico LFP 0 5.527916 5.527916 0
> P-3-M Dic.13 Lotico LFP 0 0 11.055832 0
> L1-P1-C Abr.15 Lentico LFP 127.142067 0 0 0
> L3-P3-M Abr.15 Lentico LFP 49.751244 5.527916 5.527916 0
> P-1-C Abr.15 Lotico LFP 22.111664 0 22.111664 0
> P-3-M Abr.15 Lotico LFP 0 0 16.583748 0
> L1-P1-C Ago.11 Lentico HLP 0 0 0 0
> L3-P3-M Ago.11 Lentico HLP 0 0 0 0
> P-1-C Ago.11 Lotico HLP 0 0 0 5.527916
> P-3-M Ago.11 Lotico HLP 0 0 0 0
> L1-P1-C Ago.12 Lentico HLP 22.111664 0 0 0
>
>
> A esta matriz la tengo que cortar y luego ordenar
>
>  NH3.001<-alfa[5:16,]
> NH3.001
>
> NH3.002<-alfa[1:4,]
> NH3.002
>
>
> Para armar la matriz final tengo que reemplazar el nombre de los factores
> de la variable categórica FASE.
>
> NH3.0002<-replace(NH3.002,NH3.002=="LFP","HFP")
> NH3.0002
>
> Pero no lo estoy consiguiendo. Ya entiende pasarlo como as.character o
> as.factor y no me funciona.
>
> Alguien me podría ayudar resolver este problema.
>
> Desde ya muchas gracias atentamente, Luis
>
>
> Dr. Luis A. Espínola.
> Biólogo - Ecologia de Sistemas Acuáticos
> Instituto Nacional de Limnología (INALI)
> Ciudad Universitaria - Paraje "El Pozo"
> CP: 3000, Santa Fe. Argentina
> TE: + 54 342 4511645
>
>
>
> El 3 ago 2016, a las 1:32 p.m., Carlos Ortega 
> escribió:
>
> Hola,
>
> Lo que estás comentando es la aplicación de un t-test para las dos muestras
> (tus dos series temporales) y con este test determinar si son las medias
> iguales o no.
>
> Pero con las series temporales, no puedes aplicar un t-test porque en las
> series temporales los valores están auto-correlados (hay dependencias entre
> unos y otros). Por eso tienes que confirmar tu hipótesis de que son o no la
> misma serie temporal de otra forma.
>
> Mira primero esto, que te dice cómo hacerlo:
> https://stat.ethz.ch/pipermail/r-help/2006-September/112363.html
>
> Y como supongo que aparecerán conceptos nuevos, mira estas dos referencias
> más:
>
> http://perfdynamics.blogspot.com.es/2014/04/melbournes-weather-and-cross.html
> http://ellisp.github.io/blog/2015/09/19/timeseries-same-acf
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
>
> El 3 de agosto de 2016, 4:27, Javier Gómez Gonzalez 
> escribió:
>
> Hola a todos;
>
> Quisiera saber si la frunción timeVariation cuando los intervalos de
> confianza se superponen aunque no sea en todo el perido temporal; es
> correcto decir que: “No podemos afirmar que los valores medios horarios de
> concentración de NO de la serie temporal jun-dic 10-13 son distintos de los
> de la serie temporal jun-dic 14 debido a que los intervalos de confianza al
> 95% de ambas series temporales se superponen el x% de las horas del día y
> que el intervalo de confianza al 95% de la diferencia de las medias incluye
> al cero un y% de las veces”. También tengo otra interpretación que los
> valores del periodo jun-dic 14 son la mayoría de las horas del día
> superiores a los del periodo jun-dic 10-13 pero como los intervalos de
> confianza se superponen las mediciones pueden que no sean
> significativamente distintas. ¿Son acertadas algunas de las dos
> interpretaciones? ¿Sí los intervalos de confianza se superponen en una sola
> hora del día es suficiente para afirmar que las mediciones no son
> significativamente distintas?
>
> En el caso de las medias semanales y mensuales cuando se superponen los
> intervalos de confianza para un día o un mes solamente; ¿es correcto decir
> que para ese día los valores de contaminación no son significativamente
> distintos y para el resto 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
>
>
>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> [[alternative HTML version deleted]]
>
> ___

Re: [R-es] Duda sobre construccion de vectores

[Jorge I Velez]

Buenas tardes, Francisco.

Quizás no sea la manera mas eficiente, pero proporciona el resultado que
buscas:

R> M <- matrix(c(2,3,3,4,5,6,7,7,8,11,10,5), nrow = 3, ncol = 4)
R> index <- c(2,3,3)
R> sapply(1:NROW(M), function(.row) M[.row,][index[.row]])
[1] 4 7 8

Saludos,
Jorge.-



2016-07-25 11:54 GMT-05:00 Francisco Rodríguez :

> Hola buenas tardes comunidad, seguro que hay algo fácil, pero se me escapa
> y es que estoy hoy algo torpe seguramente, a ver me gustaría hacer lo
> siguiente:
> Dada una matriz M (que puede ser muy grande), me gustaria crear un vector,
> sin usar bucles (ya que cuando es muy grande, todo puede ser muy lento) del
> siguiente modo.
> Para cada una de sus filas, me gustaría seleccionar un único elemento de
> una columna dada por otro vector, así por ejemplo:
> -Si la matriz de partida es:
> M <- matrix(c(2,3,3,4,5,6,7,7,8,11,10,5), nrow = 3, ncol = 4)
> -Así pues se tiene que:
> M
>  [,1] [,2] [,3] [,4][1,]247   11[2,]357
>  10[3,]3685
> -Y el vector indicador es:
> Indicador <- c(2,3,3)
> -Me gustaría tener una nuevo vector, cuyos componentes fueran (si el
> vector se llamase Resultado) los siguientes:
> Resultado{1} 4 7 8
> -Es decir sería la segunda componente del primer registro, la tercera
> componente del segundo registro y la tercera componente del tercer registro
> -No recuerdo si se digo hace tiempo algo de con la librería data.table, me
> da igual tener que usar cualquier librería, pero el tema es tratar de hacer
> que sea una operación vectorial que por ejemplo evite un bucle tal como:
> Resultado <- vector("numeric", 3)
> for (i in (1:3)){
> Resultado[i] = M[i, Indicador[i]]}
> -El cuál, como he comentado sería lento si el número de filas de matriz
> fuese muy elevado
> Un saludo y 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] Calcular media cada x tramos

[Jorge I Velez]

Buenos dias, Jesús,

Hay varias formas, pero te sugeriría trabajar con la función rollapply
del paquete zoo.  En

R> require(zoo)
R> ?rollapply

podrás encontrar varios ejemplos.

Espero sea de utilidad.

Saludos cordiales,
Jorge Velez.-


2016-07-15 7:31 GMT-05:00 Jesús Para Fernández :

> Buenas,
>
> Estoy intentando caracterizar una curva de 180 puntos y sacar su vector de
> caracteristica:
>
> Punto de inicio,
> Máximo
> Minimo
> 
>
> Para ello tengo un dataframe con las distintas curvas, de la forma
> siguiente:
>
>
> segundos   valor curva
> 1 1701
> 2 175 1
> ..
> 180  125 1
> 1  190 2
> 2   1932
> .
>
>
> Como puedo hacer para sacar el promedio de cada 10 segundos de la curva???
>
> 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] 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-es] par() y ggplot2

[Jorge I Velez]

Buenas tardes, Mauricio,

El sistema ggplot2 no funciona como el base, y por eso se requiere utilizar
una aproximacion diferente.  En [1] y [2] hay dos posibles soluciones;
generalmente utilizo la segunda.

Espero sea de utilidad.

Saludos cordiales,
Jorge Velez.-

[1]  http://www.cookbook-r.com/Graphs/Multiple_graphs_on_one_page_(ggplot2)/
[2]
http://rstudio-pubs-static.s3.amazonaws.com/2852_379274d7c5734f979e106dcf019ec46c.html




2016-07-13 12:28 GMT-05:00 Mauricio Monsalvo :

> Hola.
> Tengo 4 gráficos:
> a <-   qqnorm(total$ImpTotal)#con lattice
> y b, c y d son variaciones de este tipo​, creados con ggplot2:
> ​b <- g <- ggplot(data = total, aes(x=ImpTotal, fill=Convenio))
>   g + geom_histogram(binwidth = 90) ​
> Los quiero representar de la forma:
> par(mfrow=c(2,2))
> a
> b
> c
> d
> Pero los que crea con ggplot2 aparecen en par(mfrow=c(1,1)), es decir:
> ignora la instrucción.
> ¿No puedo forzar a ggplot2 a que me haga caso? ¿El problema radica en
> utilizar dos paquetes gráficos distintos? ¿Cuál sería la función
> equivalente a par(mfrow=c(2,2)) en ggplot2? Hasta donde puedo entender, los
> parámetros de faceting (http://ropensci.github.io/plotly/ggplot2/)
> funcionan para subconjuntos de datos, pero no para indicar una función
> equivalente a par().-
> Desde ya, muchas gracias!!
> --
> Mauricio
>
> [[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] 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...
>
> <>=
> plot(datos)
> @
>
> si cambio el valor del width sigue saliendo siempre el mismo tamaño
>
>
> --
> From: jorgeivanve...@gmail.com
> Date: Thu, 2 Jun 2016 08:52:14 -0500
> Subject: Re: [R-es] Ayuda con R Sweave
> To: j.para.fernan...@hotmail.com
> CC: r-help-es@r-project.org
>
>
> 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 <
> j.para.fernan...@hotmail.com>:
>
> >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] 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]

Exacto.  Y por eso escribi: "Sin embargo, y por tu ejemplo, este parece no
ser tu caso."   --JIV



2016-04-22 4:02 GMT-05:00 Jesús Para Fernández :

> Es imposible la multicolinealidad severa cuando tienes una variablre
> respuesta y una variable explicativa nada mas no???
>
> --
> From: jorgeivanve...@gmail.com
> Date: Thu, 21 Apr 2016 11:35:43 -0500
>
> Subject: Re: [R-es] Calcular Error en modelo lineal
> To: j.para.fernan...@hotmail.com
> CC: r-help-es@r-project.org
>
> Hola Jesús,
> Lo que mencionas es posible en presencia de multicolinealidad severa. Solo
> en ese caso los intervalos de confianza y prediccion son equivalentes.  Sin
> embargo, y por tu ejemplo, este parece no ser tu caso.
> Saludos,
> Jorge.-
>
>
> 2016-04-21 9:49 GMT-05:00 Jesús Para Fernández <
> j.para.fernan...@hotmail.com>:
>
> Enun ejemplo real estoy viendo como el intervalo de confianza usando lo
> que me comentas me ha salido mucho más pequeño de lo que la realidad luego
> refleja. ¿Cómo es esto posible??
>
> Es decir, veo que para valores de 2,70 obtengo una respuesta de entre 2,69
> y 2,90 y sin embargo luego en la realidad tengo valores entre 2,20 y 3
>
>
>
> Gracias
> Jesús
>
> --
> From: jorgeivanve...@gmail.com
> Date: Thu, 21 Apr 2016 09:09:03 -0500
> Subject: Re: [R-es] Calcular Error en modelo lineal
> To: j.para.fernan...@hotmail.com
> CC: r-help-es@r-project.org
>
>
> 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 <
> j.para.fernan...@hotmail.com>:
>
> 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] Una inquietud

[Jorge I Velez]

Heber,

Intenta con

R> qf(0.025,4,3, lower.tail = FALSE)

Saludos,
Jorge.-



2016-04-21 11:58 GMT-05:00 heber sarmiento via R-help-es <
r-help-es@r-project.org>:

> Estoy calculando algunos valores con qf(0.025,4,3) en R y no me coinciden
> los valores con los de la tabla que trae un texto, alguien me puede
> colaborar explicándome?
> 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] Calcular Error en modelo lineal

[Jorge I Velez]

Hola Jesús,
Lo que mencionas es posible en presencia de multicolinealidad severa. Solo
en ese caso los intervalos de confianza y prediccion son equivalentes.  Sin
embargo, y por tu ejemplo, este parece no ser tu caso.
Saludos,
Jorge.-


2016-04-21 9:49 GMT-05:00 Jesús Para Fernández :

> Enun ejemplo real estoy viendo como el intervalo de confianza usando lo
> que me comentas me ha salido mucho más pequeño de lo que la realidad luego
> refleja. ¿Cómo es esto posible??
>
> Es decir, veo que para valores de 2,70 obtengo una respuesta de entre 2,69
> y 2,90 y sin embargo luego en la realidad tengo valores entre 2,20 y 3
>
>
>
> Gracias
> Jesús
>
> --
> From: jorgeivanve...@gmail.com
> Date: Thu, 21 Apr 2016 09:09:03 -0500
> Subject: Re: [R-es] Calcular Error en modelo lineal
> To: j.para.fernan...@hotmail.com
> CC: r-help-es@r-project.org
>
>
> 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 <
> j.para.fernan...@hotmail.com>:
>
> 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] 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

[R-es] Cambiar \u

[Jorge I Velez]

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

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 :

> 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 
> 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-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-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-es] Extraer elementos diagonales de submatrices

[Jorge I Velez]

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.-
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

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

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

> 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] Pasar un listado de variables como argumento de una función

[Jorge I Velez]

Llego un poco tarde a la discusión, pero...

Intenta, a todo costo, evitar subset().  Utilizar indexacion es muchisimo
mas eficiente (aunque no igual de claro) que subset() y los resultados son
los mismos.  Si la base de datos es "pequeña", el tiempo de ejecucion es
similar, pero no ocurre lo mismo cuando la base de datos es "grande".

Saludos cordiales,

Jorge
​Velez

JCSMR, Canberra


2015-07-07 21:52 GMT+10:00 Carlos Ortega :

> Hola,
>
> Pero con lo que quieres hacer, estás definiendo el objetivo de la función
> "subset()".
> Otra cosa es que no te funcione alguno de los casos.
>
> A mí, el primer caso de "Excluir", me funciona:
>
> > DATOS <- data.frame(SE= c("M", "H", "M", "M", "H"),
> + EDAD  = c(50,  60,  20,  18,  30),
> + GRP_S = c("B", "0", "B", "A", "B"),
> + HTA   = c("N", "S", "N", "N", "N"))
> > # Excluir
> > DATOS_S <- subset(DATOS, select = -c(EDAD, GRP_S))
> > names(DATOS_S)
> [1] "SE"  "HTA"
>
> Para el caso de "Incluir", lo que definies como "c(SE, GRP_S:HTA)", ¿qué
> sentido tiene?.
> Estás definiendo una secuencia "GRP_S:HTA"  cuando ":" se utiliza en el
> contexto que quieres como secuencia entre dos números. Sí, los ":" también
> se utilizan como factores, pero para indicar una interacción de un modelo,
> que no es tu caso.
>
> Cometiendo un "pecadillo" con la función "attach()" , porque de otra forma
> SE, EDAD y GRP_S no se entienden como variales en el entorno, mira el
> resultado de lo que obtienes:
>
> > attach(DATOS)
> > c(SE, EDAD:GRP_S)
>  [1]  2  1  2  2  1 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33
> 32 31 30 29 28 27 26 25 24 23 22 21 20
> [37] 19 18 17 16 15 14 13 12 11 10  9  8  7  6  5  4  3
> Warning messages:
> 1: In EDAD:GRP_S :
>   numerical expression has 5 elements: only the first used
> 2: In EDAD:GRP_S :
>   numerical expression has 5 elements: only the first used
>
> Creo que para la elección de columnas, sigue siendo más recomendable el uso
> de "[" que de "subset()", en esa dirección va lo que te proponía. Y el caso
> de "Incluir", es el complementario de "Excluir", por lo que con esa
> solución puedes funcionar para cualquier caso.
>
> Mira estas referencias adicionales:
>
> http://stackoverflow.com/questions/9860090/in-r-why-is-better-than-subset
> http://adv-r.had.co.nz/Subsetting.html#subsetting-operators
>
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El 7 de julio de 2015, 11:16, Griera  escribió:
>
> > Hola:
> >
> > Gracias de nuevo por la ayuda!
> >
> > La solución, como no, funciona. Pero yo quería alguna cosa más flexible y
> > universal que le pudiera pasar como parámetro diferentes opciones de
> > incluir i excluir variables. Si estuviera fuera de la función seria:
> >
> > =
> > DATOS <- data.frame(SE= c("M", "H", "M", "M", "H"),
> > EDAD  = c(50,  60,  20,  18,  30),
> > GRP_S = c("B", "0", "B", "A", "B"),
> > HTA   = c("N", "S", "N", "N", "N"))
> > # Excluir
> > DATOS_S <- subset(DATOS, select = -c(EDAD, GRP_S))
> > names(DATOS_S)
> > DATOS_S <- subset(DATOS, select = -c(SE, EDAD:GRP_S))
> > names(DATOS_S)
> > # Incluir
> > DATOS_S <- subset(DATOS, select = c(SE, GRP_S:HTA))
> > names(DATOS_S)
> > =
> >
> > Pero cuando le paso a la función -c(EDAD, GRP_S)) o c(SE, GRP_S:HTA)), no
> > funciona. Existe alguna solución que pueda incorporar esta flexibilidad?
> >
> > Gracias Carlos y saludos.
> >
> >
> > On Tue, 7 Jul 2015 00:26:54 +0200
> > Carlos Ortega  wrote:
> >
> > > Hola,
> > >
> > > Puedes hacerlo de esta otra forma:
> > >
> > > #-
> > > DES = function(XDATOS, XDROP) {
> > >   #print(names(XDATOS))
> > >   #print(XDROP)
> > >   DATOS_S <- XDATOS[, setdiff(names(XDATOS), XDROP) ]
> > >   return(DATOS_S)
> > > }
> > >
> > > DES(DATOS, c("EDAD", "GRP_S"))
> > > #-
> > >
> > >
> > > Saludos,
> > > Carlos Ortega
> > > www.qualityexcellence.es
> > >
> > >
> > > El 6 de julio de 2015, 21:59, Griera  escribió:
> > >
> > > > Hola:
> > > >
> > > > Quiero que una función realice una serie de cálculos pero eliminando
> > las
> > > > variables que no interesan (diferentes según e fichero a analizar).
> > Intento
> > > > pasar esta lista como argumento con un c("VAR1", "VAR2", etc), pero
> no
> > lo
> > > > consigo. Un ejemplo seria:
> > > >
> > > > DATOS <- data.frame(SE=c("M", "H", "M", "M", "H"),
> > > > EDAD=c(50, 60, 20, 18, 30),
> > > > GRP_S=c("B", "0", "B", "A", "B"))
> > > > DES = function(XDATOS, XDROP=F)
> > > > {
> > > > print(names(XDATOS))
> > > > DATOS_S <- subset(XDATOS, select = -c(XDROP))
> > > > }
> > > >
> > > > Sin "" da el error:
> > > >
> > > > > DES(DATOS, c(EDAD, GRP_S))
> > > > [1] "SE""EDAD"  "GRP_S"
> > > > Error in print(XDROP) : object 'EDAD' not found
> > > >
> > > > Con "" da el error:
> > > >
> > > > > DES(DATOS, c("EDAD", "GR

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 :

> 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] Problema con loop

[Jorge I Velez]

De nada, Fernando.

Otra alternativa es usando mclapply() en el paquete "parallel". Estoy en
Mac/Linux y funciona muy bien. Seguramente hay algunos equivalentes en
Windows.

Saludos,

Jorge Velez
JCSMR, Canberra



2015-07-04 12:17 GMT+10:00 Fernando Macedo :

>  Muchas gracias Jorge!
>
> Si la opción del loop con for ya la tenía implementada y me demora
> bastante por eso quería probar con apply o en este caso sapply, muchas
> gracias!
>
> Saludos!!
>
> F Macedo
>
> El 03/07/15 a las 23:09, Jorge I Velez escribió:
>
>  Hola Fernando,
>
>  Podrias considerar las siguientes opciones:
>
>  R2 <- vector("list", x)
> for(i in 1:x){
> modelo <- lm(y ~ efectos[, 1:i])
> R2[[i]] <- summary(modelo)$r.squared
> }
> R2
>
>  ​ opt2 <- sapply(1:x, ​function(i){
> modelo <- lm(y ~ efectos[, 1:i])
>summary(modelo)$r.squared
> })
> opt2
>
>  Para mas información revisa ?sapply, ?lapply y ?lm.  En caso de que
> necesites otros parámetros del modelo de regresión, te sugiero revisar el
> resultado de
>
>  modelo <- lm(y ~ efectos[, 2])
>  names(summary(modelo))
>
>  Saludos cordiales,
>
>  Jorge Velez
> JCSMR, Canberra
>
>
> 2015-07-04 12:01 GMT+10:00 Fernando Macedo :
>
>> Buenas a todos, acá estoy yo de nuevo con problemas de loops.
>>
>> Tengo el siguiente problema: un vector de datos (y) y una serie de
>> efectos. El loop lo que intenta es evaluar el R² de modelos incrementando
>> por vuelta una variable efecto.
>>
>> Seria algo así:
>>
>> for(i in 1:x) {
>> modelo=lm(y~efectos[,1:i])
>> ... codigo para guardar R² y otros por cada vuelta...
>> }
>>
>> apply() puede ir columna a columna pero no sé como hacer para que me vaya
>> "acumulando columnas" por así decir. De echo no se si puede hacer.
>>
>> Espero haberme explicado y como siempre agradezco de antemano la atención.
>>
>> Un abrazo!
>>
>> --
>> Fernando Macedo
>>
>> ___
>> 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] Problema con loop

[Jorge I Velez]

Hola Fernando,

Podrias considerar las siguientes opciones:

R2 <- vector("list", x)
for(i in 1:x){
modelo <- lm(y ~ efectos[, 1:i])
R2[[i]] <- summary(modelo)$r.squared
}
R2

​opt2 <- sapply(1:x, ​function(i){
   modelo <- lm(y ~ efectos[, 1:i])
   summary(modelo)$r.squared
})
opt2

Para mas información revisa ?sapply, ?lapply y ?lm.  En caso de que
necesites otros parámetros del modelo de regresión, te sugiero revisar el
resultado de

modelo <- lm(y ~ efectos[, 2])
names(summary(modelo))

Saludos cordiales,

Jorge Velez
JCSMR, Canberra


2015-07-04 12:01 GMT+10:00 Fernando Macedo :

> Buenas a todos, acá estoy yo de nuevo con problemas de loops.
>
> Tengo el siguiente problema: un vector de datos (y) y una serie de
> efectos. El loop lo que intenta es evaluar el R² de modelos incrementando
> por vuelta una variable efecto.
>
> Seria algo así:
>
> for(i in 1:x) {
> modelo=lm(y~efectos[,1:i])
> ... codigo para guardar R² y otros por cada vuelta...
> }
>
> apply() puede ir columna a columna pero no sé como hacer para que me vaya
> "acumulando columnas" por así decir. De echo no se si puede hacer.
>
> Espero haberme explicado y como siempre agradezco de antemano la atención.
>
> Un abrazo!
>
> --
> Fernando Macedo
>
> ___
> 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 :

> 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 
> 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 
> 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] 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 :

> 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 
> 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 :
>>
>>> 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 
>>> 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 :
>>>>
>>>>> 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]]
>>>>> <>
>>>>>
>>>>> si quiero ver el contenido del primer copus
>>>>>
>>>>> > inspect(txt[[1]])
>>>>> <>
>>>>>
>>>>> [[1]]
>>>>> <>
>>>>> qB001.txt
>>>>>
>>>>> me informa cosas sobre el objeto, pero los datos no est�n all�...
>>>>

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

[Jorge I Velez]

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 :

> 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 
> 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 :
>>
>>> 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]]
>>> <>
>>>
>>> si quiero ver el contenido del primer copus
>>>
>>> > inspect(txt[[1]])
>>> <>
>>>
>>> [[1]]
>>> <>
>>> 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)
>>> <>
>>>
>>> no hace nada de nada...
>>>
>>> Saludos
>>>
>>> Oscar
>>>
>>>
>>>
>>>
>>>
>>> El 10 de abril de 2015, 1:15, Jorge I Velez 
>>> 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.-
>>>

Re: [R-es] 'Nancycats' en R

[Jorge I Velez]

Javier,
Como mencione, es solo una codificacion y no un "peso" como mencionas.
Entiendo perfectamente a lo que te refieres con los heterocigotos y hay
que tener cuidado en la interpretacion.  Sin embargo, "0.5" solo significa
que es heterocigoto en ese marcador, nada mas.
Saludos,
Jorge.-

2015-04-12 1:11 GMT+10:00 Javier Marcuzzi :

> Estimado Gemma Ruiz Olalla
>
> De curioso mir� una b�squeda en internet, y se la comparto, de esta
> observe la p�gina 8
>
> http://adegenet.r-forge.r-project.org/files/tutorial-basics.pdf
>
> Estoy de acuerdo con Jorge y Carlos, pero por las dudas, es solo una
> codificaci�n, yo sin mirar m�s no interpretar�a que el heterocigota tiene
> un "peso" de 0.5, porque si hay sobredominancia, dominancia, etc ...,
> aunque desconozco el caso en particular.
>
> Javier Rub�n Marcuzzi
>
> El 11 de abril de 2015, 11:52, Jorge I Velez 
> escribi�:
>
>> Gemma,
>>
>> Tiene que ver con cuantos alelos hay en cada locus.  Generalmente esto
>> se codifica como 0, 1 y 2, pero los autores optaron por una codificacion
>> 0,
>> 0.5 y 1 (que corresponde a tener 0, 1 o 2 alelos, respectivamente, en
>> el locus de interes).  Asumiendo un microsatelite con alelos "A" y "a",
>> por
>> ejemplo, 0 corresponderia al genotipo "aa", 0.5 a "Aa" y 1 a "AA".
>>
>> Saludos cordiales,
>> Jorge.-
>>
>>
>>
>> 2015-04-12 0:37 GMT+10:00 Gemma Ruiz-Olalla :
>>
>> > Buenas tardes,
>> >
>> > Estamos intentando hacer un estudio sobre diversidad en gatos callejeros
>> > ('stray cats') de la base de datos 'nancycats' (librer�a 'adegenet') de
>> R.
>> > Sin embargo, no entendemos muy bien los datos. �Alguien puede
>> explicarnos a
>> > qu� corresponde la variable respuesta?
>> >
>> > Entendemos que las filas son los 237 gatos (observaciones), y las
>> columnas
>> > se corresponden con los 9 microsat�lites (�dentro de los que hay las 17
>> > colonias de gatos?)
>> > Pero a qu� se corresponden las cifras 0.0, 0.1, 0.5, etc.?
>> >
>> > Muchas gracias de antemano.
>> >
>> > --
>> > Gemma
>> > gemma.ruizola...@gmail.com
>> >
>> > [[alternative HTML version deleted]]
>> >
>> > ___
>> > R-help-es mailing list
>> > R-help-es@r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>> [[alternative HTML version deleted]]
>>
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>>
>

[[alternative HTML version deleted]]

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

Re: [R-es] 'Nancycats' en R

[Jorge I Velez]

Gemma,

Tiene que ver con cuantos alelos hay en cada locus.  Generalmente esto
se codifica como 0, 1 y 2, pero los autores optaron por una codificacion 0,
0.5 y 1 (que corresponde a tener 0, 1 o 2 alelos, respectivamente, en
el locus de interes).  Asumiendo un microsatelite con alelos "A" y "a", por
ejemplo, 0 corresponderia al genotipo "aa", 0.5 a "Aa" y 1 a "AA".

Saludos cordiales,
Jorge.-



2015-04-12 0:37 GMT+10:00 Gemma Ruiz-Olalla :

> Buenas tardes,
>
> Estamos intentando hacer un estudio sobre diversidad en gatos callejeros
> ('stray cats') de la base de datos 'nancycats' (librer�a 'adegenet') de R.
> Sin embargo, no entendemos muy bien los datos. �Alguien puede explicarnos a
> qu� corresponde la variable respuesta?
>
> Entendemos que las filas son los 237 gatos (observaciones), y las columnas
> se corresponden con los 9 microsat�lites (�dentro de los que hay las 17
> colonias de gatos?)
> Pero a qu� se corresponden las cifras 0.0, 0.1, 0.5, etc.?
>
> Muchas gracias de antemano.
>
> --
> Gemma
> gemma.ruizola...@gmail.com
>
> [[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] 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 :

> 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
> >  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] 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 :

> 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 
>  Para: jose cebrian 
>  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] 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 :

> 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 :
>
> 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] Informes periódicos con R

[Jorge I Velez]

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 :

> 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] 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 :

> 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] 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 :

> 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] Aleatoriedad

[Jorge I Velez]

Hola Miguel,

Quizas lo siguiente ayude un poco a entender que ocurre:

# para reproducir los resultados
set.seed(123)

# proporcion de 'todos cero'
all0 <- mean(rbinom(1000, 8, .5) == 0)
all0

# proporcion de 'todos uno'
all1 <- mean(rbinom(1000, 8, .5) == 0)
all1

# distribution de 'todos cero'
hist(replicate(1000, mean(rbinom(1000, 8, .5) == 0)), las = 1, main =
'Todos cero', xlab = 'Proporci�n')
abline(v = all0, col = 2, lwd = 2)

# distribution de 'todos uno'
hist(replicate(1000, mean(rbinom(1000, 8, .5) == 8)), las = 1, main =
'Todos uno', xlab = 'Proporci�n')
abline(v = all1, col = 2, lwd = 2)

Como veras, la proporcion de 'todos cero' o 'todos uno' varia de manera
considerable, pero en promedio (linea vertical roja) es cercano al valor
que mencionas.

Mas informacion en ?set.seed,  ?rbinom, ?hist t ?replicate.

Saludos cordiales,
Jorge.-



2015-03-25 6:46 GMT+11:00 Our Utopy :

> Hola de nuevo, ya empiezo a ser pesado �no? bueno, no importa porque
> aprendemos todos. Eso, al menos, me parece.
>
> Hoy estuve estudiando en R el tema de la aleatoriedad. Veo que hay
> m�ltiples posibilidades pero me est�n chocando mucho. Encuentro que el
> generador de n�meros pseudo aleatorios  es m�s pseudo de lo que deber�a.
>
> Me explico, quiero generar 0 y 1 aleatorios. Estoy trabajando con una
> martingala.
>
> Uso x <- runif(1,0,1) y despu�s seg�n cual sea la probabilidad de ganar
> elijo el punto de divisi�n del intervalo apropiado.
>
> Bien, los modelos que extraigo no me parecen coherentes. Me gustaban mucho
> m�s los obtenidos con el aleatorio() de Microsoft Excel. O eso me parece
>
> �Por qu�? pues porque la probabilidad de sacar 8 "ceros" u 8 "unos"
> seguidas es 0.003 y a mi me est� pareciendo que lo hace mucho m�s a menudo
> de lo debido.
>
> Y ora cosa, �sab�is si cuando apago el ordenador me aparecen al iniciarlo
> de nuevo las mismas series aleatorias, en el mismo orden que en el trabajo
> previo?  y de ser as� �existe alguna funci�n en R que las haga iniciar por
> otro lugar de la serie interna predefinida?
>
> �Sab�is si hay test de aleatoriedad para series predefinidos en R? Es que
> el buscador de funciones ?? me busca en los paquetes incluidos pero no
> instalados pero ... �hay algo para buscar en el repositorio?
>
> Gracias ya anticipadas. Un saludo
>
>
>
>
>
> --
> Our Utopy
> http://utopicmaths.blogspot.com.es/
> http://financialmeth.blogspot.com.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