similar to: appending similar data frames?

Dear R Gang, I'm interested in using R and the nls package for fitting kinetic models. I'm having some difficulty getting a model specified for nls though. The math for the model that I want to fit is dg(t)/dt = K1 f(t) - k2 g(t) where g(t) and f(t) are measured data at a sequence of times t. K1 and k2 are the parameters of the model. If I solve this, the solution is g(t) = K1

I'd like to use stripplot for some plots because I want to use the jitter parameter. On the other hand, I'd like to use dotplot because I'd like to have the horizontal lines that it includes. dotplot doesn't have a jitter option and I'm not having any success with getting panel.grid(h=-1) with stripplot. Can anyone show me how to make dotplot-like lines on a stripplot? Or

Dear R-gang, I have a question about handling underscores in names in R 1.8.1 and 1.9.0. I recently installed 1.9.0 on a machine and found that many codes no longer work as a result of the changed behavior in make.names. I have numerous data files that have dashes, periods and underscores in the header row. I've got numerous R codes that read those files with read.table and read.csv and

Can anyone show me how to connect a scrollbar to a text widget using the tcltk library? I'm confused about how to make sure the two widgets know about each other and how to set up the appropriate call-backs for scrolling. Mike -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

I'd like to capture the results of a print() command so that I can put the results in a tktext widget. I know I can redirect the output to a file with sink. I'm looking for something like that, but where I can redirect the output into a variable. Is there a way to do that? Mike -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'd like to make a plot showing frequency of an event. The data is in a data from that includes Year, Month and Day (of month) fields, so I created a Date with ISOdate(Year, Month, Day, tz=''). I can plot frequencies for the year 2002 with > thisyear <- Date[Year==2002] > hist( thisyear, xaxt='n' ) > axis.POSIXct( 1, at=seq(min(thisyear), max(thisyear),

I have some Bernoulli data something like this: x<-sort(runif(100,1,20)) p<-pnorm(x,10,3) y<-as.numeric(runif(x)<p) plot(x,y) lines(x,p) This plot is not very satisfactory because the ogive does not visually fit the (0,1) points very well, and also because the points tend to fall on top of one another. The second problem can be eliminated by adding vertical jitter. However I was

I'm experiencing some confusion with the ellipsis argument (...). In my .Rprofile, I have the following functions: stderr <- function(x, ...) { sqrt( var(x, ...) / length(x) ) } se <- stderr I can use tapply to calculate some means: > tapply( Diameter, factor(Region), mean, na.rm=TRUE ) I II III IV V 0.02896429

Hi I'm fairly new to R and the list, so please take what I say accordingly! Far as I can see, strptime gives you a string in some specified format. In order to do any kind of numerically-based modelling with that, you need to obtain a number to work with. One way to do this is by getting the time with Sys.time() instead and coercing it to a number using as.integer(): >

Hallo! GENERAL QUESTION: I'm trying to change the tick marks of the x-axis in a grouped data plot (nlme). CONCRETE EXAMPLE: In the example (see below) I want the x-axis to have tick marks at 0, 6, 12, 18, 24. How can I do this? WHAT I TRIED I tried "normal" methods like axis(...) but this does not work with this plot. And I also tried xlim=c(0,24) but the ticks are unchanged and

Hi there folks I have R 1.3.0 running on Red Hat Linux 7.1 with GhostView 3.5.8. When I produce a file like: > bitmap("foo.pdf", type="pdfwrite") > plot(foo) > dev.off() and then try to read it with gv, I get an unrecoverable error and nothing displayed. Any help as to whether this ***should*** work would be appreciated. cheers, John -- John Williams

Dear R-help, I have a data set consisting of measurements made on multiple subjects. Measurement sessions are repeated for each subject on multiple dates. Not all subjects have the same number of sessions. To create a factor that represents the session, I do the following: data <- read.csv('test-data.csv') # data appended below data$date <- as.Date(data$date,

R help, How can I plot the below data table by Data ~ Site and group by Grinder and Equip Id so I get a chart like this Excel version? I have tried coplot with little success and lattice makes a pretty good chart like I want, but I am using the DCOM so it does not display correctly. I would like to make this using the base library if possible. <<...OLE_Obj...>> DATA TABLE

Hello, I have multiple data files. Each file contains a single column and 1.5 million rows. I want to create normalized pdfs (area under curve is 1) and histograms to compare with one another. Could anybody suggest if there exists an easy way or built in function in R. At present I am using Origin and Excel together to do this. A single file needs 10 minutes and I have a total of 929 files! So

I would like to draw smooth curves instead of polygons. I could not find any spline function to do that : given an x and a y vectors, they all take the x in increasing order. Is there a function to draw a smooth curve through a set of points in any order ? with many thanks in advance Alain Guerreau CNRS Paris [[alternate HTML version deleted]]

I have a question about using cumsum on subsets of a data frame. Suppose I have a frame that looks something like this > tmp f x y 1 left 1 0 2 left 2 0 3 left 3 9 4 left 4 10 5 left 5 23 6 left 6 45 7 left 7 13 8 left 8 2 9 left 9 6 10 right 1 10 11 right 2 26 12 right 3 9 13 right 4 50 14 right 5 78 15 right 6 20 16 right 7 7 17 right 8 20 18 right 9 19

The FAQ Section 7 is a very useful place for new users to find out any number of R idiosycracies. However there is no numbering on the FAQ Table of Content or on the Sections Tables of Contents. An R-help list reply of "Read FAQ 7.10" in response to a question about converting a factor to numeric is a bit cryptic. The only time 7.10 appears is after the searcher has found the entry.

Hi I have the following formula (I hope it is clear - if no, I can try to do better the next time) h(x, a, b) = integral(0 to pi/2) ( ( integral(D/sin(alpha) to Inf) ( ( f(x, a, b) ) dx ) dalpha ) and I want to do an mle with it. I know how to use mle() and I also know about integrate(). My problem is to give the parameter values a and b to the

Full_Name: Mike Miller Version: 1.9.0 OS: win2000 Submission from: (NULL) (134.68.121.109) The docs for make.names imply that the result of make.names(c("a and b", "a_and_b"), unique=TRUE) should be "a.and.b" "a.and.b.1" when it is actually "a.and.b" "a_and_b" The docs do not reflect this very major change in

This message is about a curious difference in timing between two ways of computing the same function. One uses expm, so is expected to be a bit slower, but "a bit" turned out to be a factor of >1000. The code is below. We would be grateful if anyone can point out any egregious bad practice in our code, or enlighten us on why one approach is so much slower than the other. The problem