similar to: Hmisc summary.formula.reverse

Hello, I was using snow to parallel-process some code in R-2.10 (32-bit windows. ). The code is as follows: require(foreach) require(doSNOW) cl <- makeCluster(6, type='SOCK') registerDoSNOW(cl) bl2 <- foreach(i=icount(length(unqmrno))) %dopar% { (some code here) } stopCluster(cl) When I run the same code in Windows R-2.11 (either 32-bit or 64-bit), R hangs at

Hi, I think there is a minor bug in odfWeave. In the function odfStyleGen, the following line has an extra "=": if(length(grep("italic", thisStyle$fontType))) fontText <- c(fontText, tagattr("fo:font-style=", "italic")) This is creating an error if some text needs to be formatted as italic, since the corresponding entry in

Hi, I'm trying to customize some formating for a table using odfWeave. I'm having some trouble getting the formatting to work. I've defined the following format defns: origStyleDefs <- getStyleDefs() StyleDefs=origStyleDefs StyleDefs$ArialLeftBold$parentStyleName='' StyleDefs$ArialLeftBold$type='Paragraph' StyleDefs$ArialLeftBold$textAlign='left'

Hi, Is there a function similar to the layout function in base graphics in either lattice or ggplot? I'm hoping someone has written a function wrapper to the appropriate commands in grid that would make this easier :) Abhijit [[alternative HTML version deleted]]

Hi, I would like to use odfWeave to output some contingency tables (the output of "table") into OOo. I know I can do this in LaTex (using "latex" in the Hmisc package), but I was wondering if it is possible in OdfWeave. My documentation to odfTable says inputs can only be vector, matrix or data.frame, and I'm having a hard time converting my table into one of these

Hi All, I have registered "abhijit" for SIP in asterisk Server. I am able to register my softphone (SJPhone) to the server using the name "abhijit". But whenever I try to make any calls I am gettinh the following error message:- *CLI> -- Registered SIP 'abhijit' at 172.20.28.85 port 5060 expires 120 May 26 07:34:52 NOTICE[2761]: pbx.c:1738 pbx_extension_helper:

Dear useRs, When I'm using Hmisc's latex() function with ctable=TRUE on a summary.formula with method="reverse" object and saving in a .tex file, the latter contains an unwanted empty line which makes compilation fail. Here is a brief example : library(Hmisc) test <- data.frame(a=sample(1:30,10),b=sample(c("good","bad","ugly"),10,replace=T))

> From: Abhijit Dasgupta <adasgupta@araastat.com> > Date: October 31, 2010 1:30:02 AM EDT > To: Matt Shotwell <shotwelm@musc.edu> > Subject: Re: [R] ForestPlot or similar > > I just did something very similar using ggplot's pointrange geom. In the following, I'm plotting hazard ratios, for which the nominal value is 1 and not 0. x has 5 columns: drug, hr,

Hi Gurus: If I have a large dataset of the form of: > x <- data.frame(V1 = runif(10), V2 = sample(c('A','B','C'),10,T)) > x V1 V21 0.2691580 A2 0.8711267 B3 0.2674728 C4 0.3278876 A5 0.1809152 A6 0.2499651 C7 0.9155174 A8 0.8004974 B9 0.7885516 A10 0.9301630 A And I want a V3 that =V2 if V2=A, and =D if V2=B or C. In other words I

Hi Abhijit I can't answer your questions but I would be very interested in a fully working binary of Ices2 for Win32 once it's ready. Thanks for the effort. Best regards, Ross Levis ==================================== StationPlaylist.com http://www.stationplaylist.com Low-cost music scheduling, live assist & automation software. ==================================== Abhijit

Dear R-users: I can not reproduce figure 4 of *Statistical Tables and Plots using S and LATEX* by Prof. Frank Harrell with the following code: rm(list=ls()) library(Hmisc) getHdata(pbc) attach(pbc) age.groups <- cut2(age, c(45,60)) g <- function(y) apply(y, 2, quantile, c(.25,.5,.75)) y <- with(pbc, cbind(Chol=chol,Bili=bili)) # You can give new column names that are not legal S names

I think qqline does not do exactly what it is advertised to do ("`qqline' adds a line to a normal quantile-quantile plot which passes through the first and third quartiles."). Consider the graph: tmp <- qnorm(ppoints(10)) qqnorm(tmp) qqline(tmp) The line (which I expected go through all the points), has a slightly shallower slope than does the points plotted by qqnorm. I think

Hi, I have a simple question about quartiles in R, especially how they are calculated using the boxplot. Quartiles (.25 and .75) in boxplot are different from the summary function and also don't match with the 9 types in the quantile function. See attachment for details. Can you give me the details on how the boxplot function does calculate these values? Cheers, Rene Brinkhuis

Why am I getting a wrong result for quartiles? here is my code: > cbiomass = c(910, 1058, 929, 1103, 1056, 1022, 1255, 1121, 1111, 1192, > 1074, 1415) > summary(cbiomass) > IQR(cbiomass) The result R gives me is: For the summary > Min. 1st Qu. Median Mean 3rd Qu. Max. 910 1048 1088 1104 1139 1415 For IQR > 91.25 ********* The true Q1 is 1039

Hello, I'm currently using a simple plot to visualize some mean values. I'm having ~200 datapoints on the x-axis, each has 10 records. I'm currently plotting only the mean value of each of the datapoints. What I need is a way to visualize the quartiles/error/whatever of these points. I thought about boxplots, but I have to many points on the xaxis - it would be impossible to see

I'm trying to devide x into tertiles, but ends up with integer limits even x holds one decimal. The analysis is extremely sensitive to the limits and I like to keep them right. How can that be done? quartiles <- quantcut( x[x >= 0], q=seq(0,1, by=(1/3)) > table(quartiles) quartiles [180,344] (344,448] (448,644] 16467 16476 16452 [[alternative HTML version deleted]]

Hi - What is the limit (rows and columns) on the size of a data set that R will process? Thanks. Abhijit Dr. Abhijit Roy Citi - Global Consumer Group - Business Analytics and Methods O: 91 80 4041 6398 Fax: 91 80 2211 0827 [[alternative HTML version deleted]]

Dear all, The incredibly useful Hmisc package provides a method to generate summary tables that can be typeset in latex. The Alzola and Harrell book "An introduction to S and the Hmisc and Design libraries" provides an example that generates mean and quartiles for continuous variables, and numbers and percentages for count variables: summary() with method = 'reverse'. I

Hello, I am modeling some survival data wih cph (Design). I have modeled a predictor which showed non linear effect with restricted cubic splines. I would like to retrieve the se(coef) for other, linear, predictors. This is just to make nice LateX tables automatically. I have the coefficients with coef(). How do I do that? Thanks, David Biau. [[alternative HTML version deleted]]

When data contains both factor and numeric variables, how to get quartiles for all numeric variables? n <- 100 x1 <- runif(n) x2 <- runif(n) x3 <- x1 + x2 + runif(n)/10 x4 <- x1 + x2 + x3 + runif(n)/10 x5 <- factor(sample(c('a','b','c'),n,replace=TRUE)) x6 <- factor(1*(x5=='a' | x5=='c')) data1 <- cbind(x1,x2,x3,x4,x5,x6) data