Re: [R-es] Representar datos longitudinales mediante splines

[Carlos Ortega]

Hola,

El suavizado con "lowess()" te da error:

> for(i in unique(df$id)) {
+   lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id))  }
Error in lowess(df$x[df$id == i], df$y[df$id == i], col = df$id) :
  unused argument (col = df$id)

Y el mensaje de error es bastante explícito. Incluyes el parámetro "col"
dentro de lowess que no lo considera.

Es un simple problema de dónde has colocado un paréntesis de cierre.

Así sí que funciona (destaco en rojo la posición del paréntesis de cierre):

#-
for(i in unique(df$id)) {
  lines(lowess(df$x[df$id == i], df$y[df$id == i]*)*, col = df$id)  }
  points(x = df$x, y = df$y, pch=16, col = df$id)
# Los datos asociados a cada id con un color de referencia
title("Scatter Plot of value vs. time")
#-

[image: Imágenes integradas 1]


Saludos,
Carlos Ortega
www.qualityexcellence.es



El 3 de marzo de 2016, 21:35, Francisco Javier 
escribió:

> Buenas noches a todos,
>
>
> Me dirijo a vosotros porque estoy trabajando con medidas repetidas sobre
> un grupo de sujetos y quisiera graficar la evolución temporal de cada uno
> de ellos pero suavizada mediante un spline (o en su defecto algún tipo de
> línea suavizadora que no necesariamente tenga que pasar por los puntos de
> cada sujeto). Además, quisiera hacerlo en R básico, esto es, sin recurrir
> al uso de paquetes adicionales tipo lattice,  ggplot2.
>
>
> A modo de ejemplo, tengo el siguiente data frame:
>
>
> df <- data.frame(
>
> id = factor(rep(c(1,2,3), c(3,4,3))),
> x = c( 1, 14, 22,  2, 9,  20, 25,  4,  15,  22),
> y = c(35, 28, 52, 79, 64, 70, 95, 45, 102, 105))
>
>
> E intento (sin conseguir resultado):
>
>
> plot(df$x, df$y, xlab = "time", ylab="value", type="n", xlim = c(0,25),
> ylim=c(20,120))
> for(i in unique(df$id)) {
>  lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id))  }
> points(x = df$x, y = df$y, pch=16, col = df$id)  # Los datos
> asociados a cada id con un color de referencia
> title("Scatter Plot of value vs. time")
>
>
> ¡Muchas gracias por cualquier ayuda al respecto!
>
>
> Francisco
>
> [[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
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Fwd: confirm 3d318a03e0b489e0782feadf5364bf2186e9d7b1

[teresa valero]

confirm 3d318a03e0b489e0782feadf5364bf2186e9d7b1
[[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] Script con error SOLUCIONADO

[Manuel Máquez]

Carlos, gracias por tu rápida respuesta.
Mira prefiero que aparezcan *todos* los ceros porque de esa manera, me doy
cuenta de cuando ocurren valores fuera del rango. Ahora que si se eliminan
todos me sirve pero menos.
Te repito mi agradecimiento.
*MANOLO MÁRQUEZ P.*

[[alternative HTML version deleted]]

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

Re: [R-es] nueva distribución de R y problema solucionado

[Javier Marcuzzi]

Estimados

Dos partes, la primera y sin importancia, “por donde van los tiros” de este 
lado del océano se comprende con una probabilidad del 99,85 % y un error del 
0.0032, con R, otro software da error de presupuesto.

La parte importante: del comentario Otra de las mejoras que han introducido y 
está disponible como paquete para todos, es el uso de un repositorio fiable

¿Cómo se usa esa posibilidad?

Javier Rubén Marcuzzi

De: miguel.angel.rodriguez.mui...@sergas.es
Enviado: miércoles, 2 de marzo de 2016 11:33
Para: c...@datanalytics.com
CC: javier.ruben.marcu...@gmail.com; R-help-es@r-project.org
Asunto: Re: [R-es] nueva distribución de R y problema solucionado

Hola Carlos.

Veo que has pillado perfectamente por dónde iban los tiros
:-)

(otro localismo, sorry... pero creo que este se entiende perfectamente)

Un saludo,
Miguel.


El 02/03/2016 a las 15:30, Carlos J. Gil Bellosta escribió:
Hola, ¿qué tal? 

Sobre

El 2 de marzo de 2016, 11:06,  
escribió: 
Que Microsoft tenga su propia versión de R y (si es el caso) su propia
versión de los paquetes... con lo dados que han sido en el pasado a
"tirar por su cuenta".. no sé yo...

Opiniones?

creo que, en primer lugar, deberíamos felicitarnos con que el superpoderoso 
Microsoft se haya fijado en R.  ¡Quién nos lo iba a decir hace 10 o 12 años!

Tiene otra ventaja. Ayer mismo hablaba con una gente de un banco y me contaban 
las pegas (reales o potenciales) que ponía el Banco de España a modelos de 
riesgos desarrollados en R. Ahora se le puede decir al regulador que se está 
usando código bajado del mismo sitio que sus actualizaciones del sistema 
operativo. 

Finalmente, creo que mientras no se pruebe lo contrario, Microsoft está 
haciendo un uso tal vez inhabitual pero permitido de la GPL y licencias de uso. 
La GPL en concreto no prohíbe el uso comercial del código. NI siquiera venderlo 
por tanto dinero como esté dispuesto a pagar un tercero. Solo tienes que 
redistribuir el código fuente. Si Microsoft consigue hacer dinero vendiendo R, 
mejor para él. Si alguien quiere darle dinero a Microsoft por algo que puede 
conseguir gratis en otra parte, peor para él. 

El único riesgo sería no tanto un "fork" --que estaría cubierto por la GPL y 
tendría que ser libre-- como una reimplementación de cero parecida pero no 100% 
compatible con R. Sería un caso parecido al de las diversas máquinas virtuales 
de Java o compiladores de C. Si esa reimplementación no fuese libre pero 
funcionase mejor que la habitual, tendríamos un problema.

Mientras tanto, creo que no hay motivo para preocuparse.

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com 



Nota: A información contida nesta mensaxe e os seus posibles documentos 
adxuntos é privada e confidencial e está dirixida únicamente ó seu 
destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe, por 
favor elimínea. A distribución ou copia desta mensaxe non está autorizada.

Nota: La información contenida en este mensaje y sus posibles documentos 
adjuntos es privada y confidencial y está dirigida únicamente a su 
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje, 
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.

See more languages: http://www.sergas.es/aviso-confidencialidad

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

Re: [R-es] Representar datos longitudinales mediante splines

[Javier Marcuzzi]

Estimado Francisco Javier

Puede realizar eso con xy, o la parte gráfica que prefiera (hay que intentar y 
se llegaría), pero en su código no veo alguna función de ajuste de curva, 
spline permite ajustar curvas y fijar algunos puntos, hay varias posibilidades 
al respecto, debería seleccionar una librería con la palabra clave spline que 
cree apropiada. 
 Yo cuándo uso splines al data.frame le agrego las referencias al spline (una 
cantidad de columnas de acuerdo al modelo), y al realizar la gráfica la 
predicción correspondiente al ajunte de curvas. 

Javier Rubén Marcuzzi

De: Francisco Javier
Enviado: jueves, 3 de marzo de 2016 17:36
Para: r-help-es@r-project.org
Asunto: [R-es] Representar datos longitudinales mediante splines

Buenas noches a todos,


Me dirijo a vosotros porque estoy trabajando con medidas repetidas sobre un 
grupo de sujetos y quisiera graficar la evoluci�n temporal de cada uno de ellos 
pero suavizada mediante un spline (o en su defecto alg�n tipo de l�nea 
suavizadora que no necesariamente tenga que pasar por los puntos de cada 
sujeto). Adem�s, quisiera hacerlo en R b�sico, esto es, sin recurrir al uso de 
paquetes adicionales tipo lattice,  ggplot2.


A modo de ejemplo, tengo el siguiente data frame:


df <- data.frame(

id = factor(rep(c(1,2,3), c(3,4,3))),
x = c( 1, 14, 22,  2, 9,  20, 25,  4,  15,  22),
y = c(35, 28, 52, 79, 64, 70, 95, 45, 102, 105))


E intento (sin conseguir resultado):


plot(df$x, df$y, xlab = "time", ylab="value", type="n", xlim = c(0,25), 
ylim=c(20,120))
for(i in unique(df$id)) {
 lines(lowess(df$x[df$id == i], df$y[df$id == i], col = df$id))  }
points(x = df$x, y = df$y, pch=16, col = df$id)  # Los datos 
asociados a cada id con un color de referencia
title("Scatter Plot of value vs. time")


�Muchas gracias por cualquier ayuda al respecto!


Francisco

[[alternative HTML version deleted]]



[[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 al importar una BD que esta en formato de SPSS

[eric]

Perfecto, muchas gracias Javier.

Saludos, Eric.






On 03/03/2016 02:36 AM, Javier Marcuzzi wrote:

Estimado Eric

Recién veo este sitio, nunca importe SPSS, no puedo compartir
experiencia al respecto más que enviarle el link por si lo cree útil.

https://github.com/hadley/haven

Javier Rubén Marcuzzi

*De: *eric 
*Enviado: *lunes, 29 de febrero de 2016 14:51
*Para: *Lista R 
*Asunto: *[R-es] problema al importar una BD que esta en formato de SPSS

Estimados, tengo que hacer un calculo muy simple, pero con una BD mas o

menos grande (250mil filas x 500 columnas) ... esta BD esta en formato

de SPSS y la importo asi:

library(foreign)

bdr <- read.spss("CASEN_2013_MN_B_Principal.sav",

use.value.labels=FALSE, to.data.frame=TRUE)

luego, quiero transformar el DF en un data.tale pues el calculo requiere

obtener promedios de acuerdo a ciertos criterios, lo que es muy facil

especificar con DT y el calculo es bastante rapido tambien ...

pero al tratar de transformar bdr a data.table

bdr <- as.data.table(bdr)

me transforma las columnas con los datos numericos que debo usar, en

characteres, y al tratar de volverlas a datos numericos con as.numeric()

me reemplaza todos los datos con NA

que estoy haciendo mal ?

adjunto algunas filas del archivo, muchas gracias,

Eric.

--

Forest Engineer

Master in Environmental and Natural Resource Economics

Ph.D. student in Sciences of Natural Resources at La Frontera University

Member in AguaDeTemu2030, citizen movement for Temuco with green city

standards for living

Nota: Las tildes se han omitido para asegurar compatibilidad con algunos

lectores de correo.



--
Forest Engineer
Master in Environmental and Natural Resource Economics
Ph.D. student in Sciences of Natural Resources at La Frontera University
Member in AguaDeTemu2030, citizen movement for Temuco with green city 
standards for living


Nota: Las tildes se han omitido para asegurar compatibilidad con algunos 
lectores de correo.


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

[R-es] Suscripcion a ayudas en español

[Tony Castillo via R-help-es]

Hola me interesa suscribirme a las ayudas en español, estoy siguiendo un
curso de coursera y quiero acceder a las ayudas en mi idioma.
Gracias

-- 
Tony Castillo Calzadilla, Universidad de Deusto, Bilbao

[[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] Script con error SOLUCIONADO

[Carlos Ortega]

Hola,

Pero, ¿qué es lo que quieres hacer?

   - ¿Que por un lado vaya contando las veces que el valor no es mayor que
   cinco y lo vaya guardando en sma?
   - ¿Y cuando el valor sea mayor, guarde un cero?..

Gracias,
Carlos.

El 3 de marzo de 2016, 7:47, Manuel Máquez  escribió:

> Gracias por la ayuda que me proporcionaron Luisfo Llador, Karel López
> Quintero y Carlos Ortega ya está caminando el script; me queda sin embargo
> una duda en su funcionamiento, que quedó de esta manera:
> j <- 1
> sma <- 0
> dat <- read.csv('1A.csv')
> for(i in 1:length(dat$d)){
>  if(dat$d[i] > 5){
> j <- j + 1
> sma[j] <- 0}
>  else{
> sma[j] <- sma[j] + 1}
> }
> La duda consiste en que dat$d[i] tiene así los primeros 11 datos: 5, 3, 5,
> 7, 1, 1, 4, 2, 14, 10, 3 y sma me da 3, 4, 0, 1, 0, 0, 0, 6, 3, 3,1; pero
> entonces entre 3 y 4 debería estar un 0; y después entre 4 y 1 debería
> haber 0, 0 y no sólo 0.
> Entonces pienso que probablemente 'metí la pata', pero no encuentro donde.
> ¿Me podrían hacer el favor de volverme a ayudar?
> Nuevamente anticipo las gracias más cumplidas.
>
> *MANOLO MÁRQUEZ P.*
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



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