Search Archives

Solved my own problem by using: odfTable.matrix( as.matrix ( with (mydata, table (site_id, reaction)) ) ) ******************************************************************************************************************** This message may contain confidential information. If yo...{{dropped:21}}

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130606/3af2455c/attachment.pl>

So the issue is something to do with the [['xxx']] construction of your data. Can you explain what thats' all about - as it errors all over the shop when I try using that... You've set me on a mission...

Duh! did it again! the variables need str in front of them don't they!! sub = sprintf('Seasonal station with natural streamflow - Lat: %s Lon: %s Gross Area %s km\UB2 - Effective Area %s km\UB2 ', round( str[['metadata']][['latitude...

>> col=c("blue","red")mydfr$[treatment] > > Yes, but I would like to use the function for lots of other dataframes > as well, so embedding 'mydfr' in the function is not the ideal > solution... In that case I'd try something...

> > # I tried defining a function like this > myplot <- function(...)plot(..., pch=19, col=c("blue","red")[treatment]) > > # So i can call it like this: > with(mydfr, myplot(Xmeas, Ymeas)) > > # but: > Error in plot.xy(xy, type, ...) : object 'treatment' not found...

I have some data which is censored and I want to determine the median. Its actually cost data for a cohort of patients, many of whom are still on treatment and so are censored. I can do the same sort...

> 2. The answer will be wrong. The reason is that the censoring occurs on a time scale, not a $ scale: you don't stop observing someone because > total cost hits a threshold, but because calendar time does. The KM routines...

I'm anything but an expert in R however if I'm labeling a graph axis with a superscript I have tended to use: > plot (x , y , xlab = expression ("label"^2)) But when you try to have more than one...

> 2. The answer will be wrong. The reason is that the censoring occurs on a time scale, not a $ scale: you don't stop observing someone because > total cost hits a threshold, but because calendar time does. The KM routines...

> I run the script and it exports a PDF called "version 1". > I want it to check if "version 1" already exists. If so, > then I want the new graphs to be exported as > "version 2", and so on. > > Is...

I might be missing something really obvious, but is there an easy way to locate all non-unique values in a data frame? Example mydata <- numeric() mydata$id <- 0:8 mydata$unique <- c(1:5, 1:4) mydata$result <- c...

This may be a really obvious question but I just can't figure out how to do it. I have a small dataset that I am trying to compare to some controls. It is essential that the controls are matched...

>Ah, I think I see what you want. Try this on each pair of exclusive sets: > n_total<-dim(mydataset)[1] under65<-mydataset$age <= 65 n_under65<-sum(under65) under65row<-c(sum(mydataset$dose[under65] == "FD"), sum(mydataset$dose[under65...

>Ah, I think I see what you want. Try this on each pair of exclusive sets: <snip> >Then under65row and over65row should be the first two rows of your result. >Can't test this at the moment, but I don...

I've been tweaking code for several days on and off in R, cut and pasting in from a text editor (I just leave them open all the time). I think I got something that was usable but then a...

Hi, I'm sure this should be simple but I can't figure it out! I want to get the median survival calculated by the survfit function and use the value rather than just be able to print it. Something...

I am using seq with the expression seq(1.4, 2.1, by=0.001) to create a sequence of references from 1.4 to 2.1 in 0.001 increments. They appear to be created correctly. They have a...

I may be doing this wrong! but I have a function which I have simplified a lot below. I want to pass some 'field names' of a data-frame to the function for it to then do some manipulation of...

Slightly confused because if I try: > newdata.yaxis = c(2.0000473, 3.123456, 3.23456, 2.67890, 1.56789) > newdata.yaxis_4 = round (newdata.yaxis, digits = 4) > newdata.yaxis [1] 2.000047 3.123456 3.234560 2.678900 1.567890...