similar to: R alternative to SAS PROC REPORT

I am trying to simulate a series of ones and zeros (1 or 0) and I am using "rbinom" but realizing that the number of successes expected is not accurate. Any advice out there. This is the example: N<-500 status<-rbinom(N, 1, prob = 0.15) count<-sum(status) 15 percent of 500 should be 75 but what I obtain from the "count" variable is 77 that gives the probability of

Hi friends, I need suggestions/directions on how to producing a waterfall plot for present extend of change in tumour size for a set of respondents in a study.  Example of use of waterfall plot is in the following slides presented at ASCO 2007 by Axel Grothey. Link is

I am new to more radical programming in R. I am trying to write a nested 'for' loop to produce output that takes subscripts like: for i taking values 1,2,3,4,5 and j taking values 1,2,3 I want to output for a computation using the combination values of i and j a value x like this; i j x 1 1 x11 1 2 x12 1 3 x13 2 1 x21 2 2 x22 2 3 x23 3 1 x31 3 2 x32

I am trying to define groupings from levels of factor variables and this the warning message that R give "& not meaningful for factors". The nature of my task is this. I have a variable stage which has the levels (1B, 2A, 2B) - these are the AJCC TNM stages of cancer, and another variable diameter with factor levels ("=< 4", "4 - 6.5, > 6.5; limit values are

Hello, This may sound crazy, but I have a large number of Multiplan data files I'm in the process of recuperating, and I'm hoping to avoid having to open them one by one to convert them into a modern, directly usable format. So I was wondering if someone somewhere had encountered this and maybe had a way to read Multiplan files in R, or, at least, a way to read SYLK files directly in R.

Dear All, A friend of mine seeks help in converting the attached ascii to netcdf. Any assistance will be sincerely appreciated. Regards, Marshall -------------- next part -------------- An embedded message was scrubbed... From: Sepo Hachigonta <shachigo at csag.uct.ac.za> Subject: asscii2netcdf Date: Fri, 05 Sep 2008 19:15:49 +0200 Size: 1292765 URL:

Hi, Sorry for the simple question. Is there a way to change the name of only one column of an existing data frame? I know colnames allows you to set the name of all the columns, but only one column in the middle of my data frame needs a new name. Thanks, -Nina

Hi, Given a vector, say x=sample(0:9,10) x [1] 0 6 3 5 1 9 7 4 8 2 I can find the location of an element by which(x==2) [1] 10 but what if I want to find the location of more than one number? I could do c(which(x==2),which(x==3)) but isn't there something more streamlined? My first guess was y=c(2,3) which(x==y) integer(0) which doesn't work. I haven't found any clue in the R

Dear latex/R-Sweavers, Using the codel below, I can color text in individual cells for latex output. Is there a similar way to get a background shading? My attempts failed because I did not get the closing brace at the right place with Hmisc/latex. library(Hmisc) x <- as.data.frame(diag(rnorm(3),nrow=3)) cellTex <- matrix(rep("", NROW(x) * NCOL(x)), nrow=NROW(x)) cellTex[2,2]

HI there, i know this is a basic question, though i need some help because this is somewhat away from my current issue, but nevertheless interesting to me... Lets assume i have some estimated probabilities, say estimated by a logit model. i know i can also state them as an odds ratio. Now i?d like to state these odds ratios as a reference to a specific outcome of my investigated

Dear All: I am working on writing some R functions to make statistical reports automatically. Dr. Harrell's Hmisc has all the wonderful stuff. But sometimes I need change some formats, so I want to read through it and make some modifications to fit my project. Ideally, I want proceed as following: 1. change some source of Hmisc 2. compile and install the modified Hmisc 3. debug my

Dear R-users, Attach with is my data..what i want to do is finding a suitable distribution for my data..I want to run a few test like the poisson and the exponential distribution. Please help me on how to find the p-value for poisson as well as the exponential distribution without knowing the parameter. Is it possible?? Thanks in advance. love, Anisah -------------- next

Hi, I'll apreciate your guidance on how can I re-create the output from SAS PROC VARCLUS in R. I've found the varclus function in Hmisc. However, is it possible to use that function to compute for each variable the 1-R**2 ratio (this is the ratio of 1 minus the R-squared with Own Cluster to one minus the R-squared in the Next Closest cluster)? Thanks in advance for any help, Lars.

To anyone that may be able to help, First I applogize if this message ends up in your inbox twice. The first one seemed to get stuck in a pending status so I deleted and started over. I am relatively new to R, however I know the basics through some classes I have taken. Unfortunately the classes left off the more complicated analysis that I need to learn and use. I am attempting to replicate

SAS proc rank has ties options of high and low that would allow producing ranks of the type found in the sports pages, e.g., rank (c(1,1,2,2,2,2,3)) == 1 1 3 3 3 3 7 Could R support these ties.methods?

I am running a 3-parameter nonlinear fit using the default Gauss-Newton method of nls. initialValues.L = list(b=4,d=0.04,t=180); fit.nls.L = nls( myModel.nlm , fData.L, start = initialValues.L, control = nls.control(warnOnly = TRUE), trace=T ); summary.nls.L = summary(fit.nls.L); I run the same analysis in SAS proc NLIN. proc nlin data=apples outest=a; parms b=4 d=.04 t=180; model Y =

I am trying to use lme to fit a mixed effects model to get the same results as when using the following SAS code: proc mixed; class refseqid probeid probeno end; model expression=end logpgc / ddfm=satterth; random probeno probeid / subject=refseqid type=cs; lsmeans end / diff cl; run; There are 3 genes (refseqid) which is the large grouping factor, with 2 probeids nested within each refseqid,

Example: I have the following model > model <- lmer(response ~ time * trt * bio + (time|id), data = dat) where time = time of observation trt = treatment group (0-no treatment / 1-treated) bio = biological factor (0-absent / 1-present) and I would like to obtain an estimate (with standard error) of the change in response over time for individuals in the

PROC MIXED is used to fit mixed effects model for correlated data. Usually we can use either a REPEATED statment or a RANDOM statement. The random statement is corresponding to lme function in R -- specifying a random effect term. The repeated statement actually directly specifies the covariance structure -- is there a similar function in R to do this? I currently want to specify a unstructured

Dear R friends A while a go I sent an email to the epi-list and later to the help-list and no answer could fully illuminate my question. So Im trying again with a more specific matter. Im trying to work on a script (function) to analyse data from a diagnostic test meta-analysis with random effects. This was first described by an author using SAS witn PROC NLMIXED. Im not an expert in R and much