similar to: Private data libraries (idea & trick)

Dear Rusers, My question is about "recode variables". First, i'd like to say something about the idea of recoding: My dataset have three variables:type,soiltem and airtem,which means grass type, soil temperature and air temperature. As we all known, the change of air temperature is greater than soil temperature,so the values in those two different temperaturemay represent different

Greetings - I am working on some data that contain chemical names with air concentrations, and I am creating a data frame with date/time and each chemical having its own column. However, these are organic chemicals (e.g. 1-butene, 2,3,4-trimethylbenzene etc). The package I am going to be using the data with is openair, and many of the great functions require you to specify a column name which

I am recoding some data. Many values that should be 1.5 are recorded as 1-2. Some example data and my solution is below. I am curious about better approaches or any other suggestions. Thanks! # example input data myData <- read.table(textConnection("id, v1, v2, v3 a,1,2,3 b,1-2,,3-4 c,,3,4"),header=TRUE,sep=",") closeAllConnections() # the first column is IDs so remove

Kjetil Kjernsmo <kjetil.kjernsmo@astro.uio.no> writes: > On 8 Nov 2000, Peter Dalgaard BSA wrote: > > >> > range(sapply(1:2000, function(n) mean(rpois(10000, c(15,15+1e-8))))) > >> [1] 14.8692 15.1200 > > > >AHA! Spotted, I think. > > Wow! Great, that was fast! > > >It is possible to return from rpois in step N, > >in which case

Just to get these issues recorded in the bug repository: a) We don't handle multiple window closes well My "favourite" version of this is when I ^Z out of an R session with an open graphics window to do something on the command line, then kill the graphics window (which does nothing), kill it again, realise that it's because the R session was stopped, and: [pd@blueberry

I've just gotten through translating my notes for the course in basic statistics for health researchers. I've placed a copy on ftp://blueberry.kubism.ku.dk/priv/R-intro.tgz This is essentially a translation of the Danish notes that I wrote during the Autumn of 1997 for use with Rsept, i.e. roughly R v.50. They lean heavily on the textbook, Altman: Practical Statistics for Medical

Dear R-helpers, I am helping a SAS user run some analyses in R that she cannot do in SAS and she is complaining about R's peculiar (to her!) way of recoding variables. In particular, she is wondering if there is an R package that allows this kind of SAS recoding: IF TYPE='TRUCK' and count=12 THEN VEHICLES=TRUCK+((CAR+BIKE)/2.2); Thanks for any help or suggestions you might be able

ExpeRts, I'm trying to compare three survival curves using the function "survdiff" in the survival package. Following is my code and corresponding error message. > survdiff(Surv(st_months, status) ~ strata(BOR), data=mydata) Error in survdiff(Surv(st_months, status) ~ strata(BOR), data = mydata) : No groups to test When I check the "strata" of the variable. I get .

I've rolled up R-0.65.1.tgz a moment ago. As usual, you shouldn't get it from Auckland, but wait for it to be mirrored at a CRAN site near you within a day or two. For those who are too curious to wait, I'll put a copy in ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.65.1.tgz There's also a version split in two for floppies and a patch file if you prefer that. (If you use the

I've rolled up R-0.65.1.tgz a moment ago. As usual, you shouldn't get it from Auckland, but wait for it to be mirrored at a CRAN site near you within a day or two. For those who are too curious to wait, I'll put a copy in ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.65.1.tgz There's also a version split in two for floppies and a patch file if you prefer that. (If you use the

Dear R-forumites, I want to apply a Wilcoxon test on subsets of the data frame mydata, splitted using the myindice variable. When I send : wilcoxtest <- by(mydata, mydata$myindice, function(x) {wilcox_test(x$value~x$fact)}) I get : Error in eval(expr, envir, enclos) : object "x" not found whereas wilcoxtest <- by(mydata, mydata$myindice, function(x)

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130130/7bc4cd37/attachment.pl>

Dear Developers, I've been struggling through writing R extension in C. I've been using an external pointer to store my data (please see sample below). I encountered a very weird erroneous behaviour: when I tried to use my external pointer to a structure holding several types of data, including SEXPs, I discovered that SEXPs change their types between returning from initialization

That's not an R question but a stats question, but I wouldn't do it. For one thing: The variance of binary data is a function of the mean, so the research question is dubious in the first place. Secondly, the test is based on ranking and comparing absolute differences from the group median, which for binary data is generally 0 or 1, so all absolute differences will be 1.... Put

I am hoping for some advice regarding the difficulties I have been having recoding variables which are contained in a csv file. Table 1 (below) shows there are two types of blanks - as reported in the first two columns. I am using windows XP & the latets version of R. When blanks cells are replaced with a value of n using syntax: > affect [affect==""] <- "n" there

Hi, I have used the mgcv library to generate a simple additive model. I want to know how to plot the function on the raw data with confidence intervals whan I have TWO variables in the model. I get it to work with one variable but not with two. I am on the limit for what I understand in R, so be gentle. I have read the help file on predict.gam, but did not get any help out of it. #My model:

> For those who *are* desperate, I've left a copy in > ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.64.0.tgz > (Be gentle, that's my desktop PC!) Well, *now* it's there anyway... -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918

> For those who *are* desperate, I've left a copy in > ftp://blueberry.kubism.ku.dk/pub/R-devel/R-0.64.0.tgz > (Be gentle, that's my desktop PC!) Well, *now* it's there anyway... -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918

Dear All, I am trying to reproduce the example that I found online here http://bit.ly/11VG4ha However, when I run my script (pasted at the end of the email), I notice that there is a factor 2 between the values for the coefficients for the categorical variable female calculated by my script and in the online example. Any idea about where this difference comes from? Besides, how can I

Dear R-helpers, I know of two ways to reading data within an R program, using textConnection and stdin (demo program below). I've Googled about and looked in several books for comparisons of the two approaches but haven't found anything. Are there any particular advantages or disadvantages to these two approaches? If you were teaching R beginners, which would you present? Thanks, Bob