search for: updata

Dear listserve members, especially Prof. Harrell: I am trying to create a factor variable that has fewer levels than the original. I have a factor: >rosa$risk1 [1] 2 2 5 1 ... [1799] 3 3 1 3 1 6 3 3 1 5 3 5 3 3 3 0 3 3 3 1 1 3 Levels: 0 1 2 3 4 5 6 8 But when I do this: rosa2 <- upData(rosa, levels=list(risk1=list(1='1',2='2',3='3',0=c('0','4','6','8')))) I get: Error: syntax error But I am able to reproduce the dat2 <- upData(dat, etc...) example in Prof. Harrell's Hmisc notes. Thanks! Tanya Murphy

How do I add data to a .rdata file? In my case, I have a time series that needs to get updated every day. Thanks, Jason [[alternative HTML version deleted]]

...############################### # BELOW THIS LINE, NOTHING MORE SHOULD NEED TO BE CONFIGURED. our @ifconf = (); our @iprout = `/sbin/ip route`; our $p1; our $p2; our $ppp0=0; our $ppp1=0; our $FAILED='FALSE'; our $FAIL='TRUE'; our $attempt=0; our $date=''; our @data=(); our @updata=(); our @log=(); our @total_log=(); our $cur_day=0; our $cur_month=0; our $cur_year=0; our $late_day=0; our $late_month=0; our $late_year=0; our %months=(Jan,1,Feb,2,Mar,3,Apr,4,May,5,Jun,6,Jul,7,Aug,8,Sep,9,Oct,10,Nov,11,Dec,12); our $start_time=`/bin/date`; our $stop_time; #####################...

Something strange is happening with update.packages on my machine (Win2000, R-1.5.0 installed from binaries). Some of the packages (see below) keep getting updated every time. When I issue update.packages() I get this every time for the last few days (I think since I installed 1.5 and the packages but I am not sure). This happens with CRAN shown below and also with CRAN set to

...y summaryM * show.html, print.html: removed; conflicted with htmltools/rmarkdown * label: required exact match on attribute name "label" to not retrieve the "labels" attribute (thanks: Ivan Puzek) * htmlSN: new function to convert floating point to scientific format usint html * upData, cleanup.import: fix NA becoming new factor levels (thanks: Beau Bruce) * htmlTranslate: new function * plotlySave: new function * histSpikep: new function * upData: new argument html * plot.describe: added ggplot2 graphics, improved continuous display using ggplotly * labelPlotmath: cleaned up by...

...y summaryM * show.html, print.html: removed; conflicted with htmltools/rmarkdown * label: required exact match on attribute name "label" to not retrieve the "labels" attribute (thanks: Ivan Puzek) * htmlSN: new function to convert floating point to scientific format usint html * upData, cleanup.import: fix NA becoming new factor levels (thanks: Beau Bruce) * htmlTranslate: new function * plotlySave: new function * histSpikep: new function * upData: new argument html * plot.describe: added ggplot2 graphics, improved continuous display using ggplotly * labelPlotmath: cleaned up by...

When updataing packages from CRAN, I got the following error message: updating HTML package descriptions Warning message: DLL attempted to change FPU control word from 8001f to 9001f Everything seems fine, the packages updated, just wondering what this was. R 1.6.2 on Windows 98. Details on the packages...

...code can work with xyplot from the lattice package (does not have error bar though). Can somebody kindly offer some help? Thanks, Wen dfr <- expand.grid(month=1:12, continent=c('Europe','USA'), sex=c('female','male')) set.seed(1) dfr <- upData(dfr, y=month/10 + 1*(sex=='female') + 2*(continent=='Europe') + runif(48,-.15,.15), lower=y - runif(48,.05,.15), upper=y + runif(48,.05,.15)) myPanel <- function(x, y, subscripts, ...) { panel.xYplot(x, y, ...)...

Dear list, I've tried to search the archives but found nothing, although I may use the wrong wording in my searches. I've also double-checked the upData function in Hmisc, but it does something else. I'm wondering if one can update a dataframe by "forcing into" it a shorter dataframe containing the corrections, like the "update" provided in SAS data steps. In this simple example: a <- data.frame(id=c(1:5),x=rnorm(5)...

...quot;%d-%B-%y %H:%M") whereas I can get what I want if I do it in 2 steps and rbinding ? Thanks and best regards, Jean-Louis This is the R console output > bsamp<-read.table("bsampl2.csv",header=T,sep=";") > names(bsamp)<-tolower(names(bsamp)) > bsamp<-upData(bsamp,drop=c("study")) Input object size: 23896 bytes; 15 variables Dropped variable study New object size: 23016 bytes; 14 variables > bsamp$visitdat<-as.character(bsamp$visitdat) > bsamp$samtime<-as.character(bsamp$samtime) > bsamp$admtime<-...

...rary(Hmisc) getHdata(prostate) # Variables in prostate had units in ( ) inside variable labels. Move # these units of measurements to separate units attributes # wt is an exception. It has ( ) in its label but this does not denote units # Also make hg have a legal R plotmath expression prostate<-upData(prostate, moveUnits=TRUE,units=c(wt="", hg="g/100*ml"),labels=c(wt="Weight Index = wt(kg)-ht(cm)+200")) attach(prostate) stage<- factor(stage, 3:4, c("Stage 3","Stage 4")) s6<-summary(stage~rx+age+wt+pf+hx+sbp+dbp+ekg+hg+sz+sg+ap+bm,method=&q...

...directory ERROR: compilation failed for package 'RBGL' ** Removing '/usr/local/lib/R/library/RBGL' %%%%%%%%%%%%%%%%%%%%%%%%%%%%% Is it maybe because I use redhat9 with libxml2-2.5.4-1, which doesn't neet the not meets the requirment of RBGL_1.12.0 ? But I do not know how to updata them, even after downloading the file libxml2-2.5.4-3.rh9.i386.rpm, I failed to install it: [root at yxnlinux software]# rpm -i libxml2-2.5.4-3.rh9.i386.rpm file /usr/bin/xmlcatalog from install of libxml2-2.5.4-3.rh9 conflicts with file from package libxml2-2.5.4-1 file /usr/bi...

...trouble writing my own simple latex code and converting it to pdf using pdftex in miktex). Thanks, Mike Babyak Department of Psychiatry and Behavioral Science Duke University Medical Center ************************** Here's the code: library(Hmisc) library(survival) getHdata(pbc) pbc<-upData(pbc, moveUnits=TRUE, labels=c(stage='Histologic Stage\nLudwig Criteria')) kmsurv <- function(S, times) { f <- survfit.km(factor(rep(1,nrow(S))), S) tt <- c(0, f$time) ss <- c(1, f$surv) # add first point to survival curve approx(tt, ss, xout=times, method='constant', f=...

Hi all, I am trying to do this type of graph (code taken from the example of xYplot help). dfr <- expand.grid(month=1:12, continent=c('Europe','USA'), sex=c('female','male')) set.seed(1) dfr <- upData(dfr, y=month/10 + 1*(sex=='female') + 2*(continent=='Europe') + runif(48,-.15,.15), lower=y - runif(48,.05,.15), upper=y + runif(48,.05,.15)) xYplot(Cbind(y,lower,upper) ~ month,groups=sex, subset=continent=='Europe', data=dfr...

Hello, The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0 and beta1 via GLM) 2) with lambda, estimate alpha via ML estimation 3) with updataed alpha, replicate 1) and 2) until alpha is converged to a value I coded 1) and 2) (it works), but faced some problems with the loop. It produced some errors. Is there someone to help me to figure this problem for the loop. Here are two codes: [I] one is for 1) and 2) (working), and [II] another o...

I am using xYplot to plot data and CIs. How do I add several lines to the figure? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011 +1-434-982-4729 Lab: B019

Hi, I want to add error bars to a plot generated with xyplot. I've tried both errbar() and plotCI(), but in both cases the points are not in the same place. It's as if the two functions are using a different frame of reference for the plotting area. for example: means <- c(92.5, 92.25, 90.9, 91.0, 94.15, 90.05) #means time <- c(1,1,2,2,3,3) #occasion variable group <-

Dear R Gurus Just started on R ! Using xYplot from Hmisc (R 1.9, W2K) I get a grey/blue background that I would like to change to white (ie no background) or may be to another color. Tried to do that with par(bg) but only changed the color of the trellis heading. What's the right command to do that ? Kind regards, JL PS if anyone has nice default settings for win device please let me

Hi R users, I am new to using DBMS with R for large datasets. Thanks to all who responded with useful suggestion to my earlier postings about using large datasets and DBMS with R. I am writing to get some help about how to design good tables in DBMS to take full advantage of the wonderful built-in facilities in R, like labels. I am using RMySQL client. Because R makes good use of variable

...na(c$X4)] <- c$X4.y[is.na(c$X4)] c <- c[,1:11] One quickly becomes tired with this solution with my real-life dataframes where different columns would require "repair" from one case to the next. I think I still prefer making the narrower dataframe like the wider one: b2 <- upData(b, X1=NA, X2=NA, X3=NA, X5=NA, X6=NA, X7=NA, X8=NA, X9=NA, X10=NA) b2 <- b2[,c(1, 3:5, 2, 6:11)] d <- rbind(a, b2) But again this requires quite a bit of fine-tuning from one case to the next in my real-life dataframes. I suspect R has a neat way to do this and I just cannot come up wi...