similar to: Representar datos longitudinales mediante splines

Buenos días, Soy Cristian Rodelo y quería haceros una pregunta. Tengo un conjunto de datos en los que a un grupo de sujetos se les realizaron medidas de tensión arterial a lo largo de 24 horas. En promedio 4 determinaciones horarias. Es decir que cada sujeto se le realizaron 4 determinaciones horarias durante 24 horas. Bien, quiero determinar la trayectoria longitudinal de las mediciones de

Buenas noches a todos las personas del foro, Me dirijo a vosotros con la siguiente cuestión: ¿Cómo es posible obtener, para cada individuo (identificado con la variable ID) de un data frame, la CANTIDAD ACUMULADA hasta una determinada fecha (día, mes y año), que es diferente para cada individuo y que se denota por la variable ENTRADA? Ejemplo: require(data.table) datos2 <-

Hola a todos: quiero consultarles para estar seguro de que estoy entendiendo bien el funcionamiento de la función smoothSpline() del paquete 'timeSeries'. Tengo una serie temporal con datos mensuales a la cual quiero suavizar usando splines para, por ejemplo, comparar con otras series temporales. Por lo que estuve viendo, me conviene usar la función smoothSpline() que se basa en

Buenas tardes a todo el mundo, El motivo de que me dirija al foro es que no sé cómo añadir unas determinadas filas en un data table de una forma eficaz y según el identificador de sujeto. Además, quisiera hacerlo sin tener que recurrir a ningún paquete de R adicional del tipo plyr, dplyr, etc. (no sé si esto es posible). EJEMPLO: require(data.table) DT <- data.table( id =

Estimados, Necesito estudiar los efectos de dos factores sobre una variable de conteo con un rango estrecho de variación (es decir, probablemente sin distribución Normal). El experimento tiene más de 10 (fechas de muestreo) observaciones repetidas en el tiempo. En primer lugar, considero que debería aplicar un modelo lineal generalizado (¿Mixto?) usando Poisson como distribución de referencia.

Saludo cordial para cada uno. Les pido ayuda para generar números aleatorios de una mixtura k-puntos. Sabemos que la función de distribución F es una mixtura k-puntos si es de la forma F(x) = p_1 F_1(x) + p_2 F_2(x) + … + p_k F_k(x), donde F_j es una función de distribución de probabilidad, p_j > 0 y suma(p_j) = 1, para j = 1, 2, …, k. En mi caso particular F es la suavización de la

Hola. Tengo que representar en un mapa, valores de proporciones de extranjeros con tarjetas de residencia en vigor, según nacionalidad y país de procedencia. Los datos son de este estilo Rumania Andalucía 0,200365823 Aragón 0,377192073 Asturias 0,205353012 I. Baleares 0,06592893 Canarias 0,03050312 Cantabria 0,190397393 Castilla-León 0,211670421 Castilla-La Mancha

Estimados, tengo una consulta acerca de gráficas con xyplot. Tengo una base de datos con dos perfiles de concentracion por sujeto, ya que a cada sujeto se le administraron dos formulaciones. Con la función xyplot he podido graficar cada perfil separando sujetos y formulaciones: pl <- xyplot (DV~TIME | factor(ID)+factor(FORM), data=tab, type = 'l',

Hola, compañeros. Pido ayuda con algo que sé que tiene que ser simple, pero la presión de tener que sacarmelo de encima me simplifica a mí demasiado y no me doy cuenta. Tengo una matriz de datos en la que tengo características tipo factor, necesito trasponer esa información a una matriz de datos binarios, en función de algunas combinaciones de esas variables. Unas combinaciones tienen que pasar a

Hola chic en s, querría construir mi primera función, y tengo una duda respecto al comando par( mfrow =c(3,3)). Primero de todo, tengo una tabla con 10 variables, para cada variable, unos 145 datos. Quiero representar para cada variable su gráfica de dispersión respecto a las demás. Es decir, coger la primera variable y la segunda, y hacer gráfica, coger la primera variable y la tercera, y hacer

Hola:   Tengo una lista de 2 elementos, cada uno de los cuales contiene información relativa a un sujeto (Pablo y Carlos).   lSujetos <- list() lSujetos[[1]] <- list(nomfich="Pablo", colTime=5, colVars=c(6,7,8)) lSujetos[[2]] <- list(nomfich="Carlos", colTime=5, colVars=c(6,7,8)) A continuación, leo las series temporales correspondientes a cada individuo. En este caso

Hi, all, I want to use lowess(x, y) where x and y are vectors of length of 4000+. In fact, x and y are log of some vectors. So, some of the elements are NaN. lowess() can not take away those elements then do the fitting. It will give the error message and do nothing. 1. Can anybody tell me how to get rid of those NaN's and use lowess()? 2. How to get the LOWESS fitting values for any

Colleagues ---------------------------------- System info: R version rw1022 on NT ESS v. 5.1.18 using emacs ver. 20.4 ---------------------------------- I would like to access the residuals from a lowess fit (function lowess() ) as an aid to assessing whether I am over-smoothing or under-smoothing a dataset (as in Cleveland's book "Visualizing data", section 3.3) Unless I

The lowess() help page refers to documentation in "src/appl/lowess.doc". This was moved to "src/library/stats/src/lowess.doc" in 2007. This patch fixes it: Index: src/library/stats/man/lowess.Rd =================================================================== --- src/library/stats/man/lowess.Rd (revision 79137) +++ src/library/stats/man/lowess.Rd (working copy) @@

Hi, A lowess method is defined in our package for one of our S4 classes. To explicitely define the generic this is being used if (!isGeneric("lowess")) setGeneric("lowess", useAsDefault = lowess) This works fine for many other methods inherited from various R packages. In this case a warning is issued in R CMD check: * checking R code for possible problems ... NOTE

I'm using this code to plot a smoothed line. These two columns of data really represent 4 groups and I'd like to plot a separate line for each group but have them all in the same plot. The R-Docs for lowess do not seem to indicate some type of "GROUPS=var_name" option. What would be the syntax for this? plot(AWGT ~ lipid ) lines(lowess(lipid , AWGT, f=.8)) -- Dean

Dear Colleagues: I need to make a graphic that uses the Nimbus rather than Helvetica font family so that the font can be embedded in the encapsulated postscript file. This is to satisfy a requirement from a journal for electronic submission of figures. I do the following: postscript('lowess.eps', onefile=FALSE, pointsize=18, horizontal=FALSE, family='NimbusSan')

Hi there, I'm plotting some glass RI values just by plotting plot(x) then I put on my lowess smoother lines(lowess(x)) now I want to put on some 95% Confidence Interval bands of the lowess smoother, but don't know how?? Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E gareth.campbell@esr.cri.nz gcam032@gmail.com [[alternative

x <- c(0,7,8,14,15,120,242) y <- c(122,128,130,158,110,110,92) lowess(x,y) $x [1] 0 7 8 14 15 120 242 $y [1] 122.0000 128.0000 132.2857 158.0000 110.0000 -4930.0000 110.0000 R version 2.2.1, 2005-12-20, i486-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils"

Hello Masters, I run the loess() function to obtain local weighted regressions, given lowess() can't handle NAs, but I don't improve significantly my situation......, actually loess() performance leave me much puzzled.... I attach my easy experiment below #------SCRIPT---------------------------------------------- #I explore the functionalities of lowess() & loess() #because I have