similar to: strsplit with dataframes

I have data from several sensors that recorded data at hourly intervals during seven months. I want to separate daily variation from the trend, and also be able to zoom in on only one month of data. I have not been able what functions to use, I can not figure out from the help for 'ts' how to use hourly data. I guess this is routine-work for a lot of people so I hope someone can

I want to create a character vector based on the table (shortened for display) below: Where there are ones in the matrix I want the column name to appear and where there are zeros nothing, which would make the vector in this shortened case: combinations <- ("A B","A C","A E H","A F G","B C D","E G H",A C D E H","A C D F

I have been using R for a few months to plot my data, and fit statistical models and so on. It is stated that R is not only a package for statistcs and graphics but also a programming language. Currently I am working with Fortran 90 to do numerical simulations (1D reactive -transport), and subsequently visualize it in some external software (spreadsheet). As I now started using R I will plot my

I really like to use R for all my graphs, and as I work with stable isotopes I want to have a proper chemical notation in my plots I have looked at ?plotmath, but didn't find the answer and also searched the R website. ------------------------------------------------------------------------ -- plot(1:10,xlab=expression(^{14}*C)) # I want to have a superscript with nothing in front, but it

Dear R-gurus, I would like to zoom in a plot, e.g. I select a region on the x-axis and then I would like the ranges on the y-axis to change accordingly. Is it possible to do this with existing functions, or do I have to invent some data selection before plotting? See below a short example, where I select ylim with trial and error, which I want to avoid. Cheers, Henrik Andersson

Just a small suggestion since Windows have a file system not designed for command line use... Would it be possible to add the possibility of automatically adding/modifying the path to the R executables in the windows installer program? --------------------------------------------- Henrik Andersson Netherlands Institute of Ecology - Centre for Estuarine and Marine Ecology P.O. Box 140 4400 AC

Dear helpers, I wrote a code which estimates a multi-equation model with generalized least squares (GLS). I can use GLS because I know the covariance matrix of the residuals a priori. However, it is a bit slow and I wonder if anybody would be able to point out a way to make it faster (it is part of a bigger code and needs to run several times). Any suggestion would be greatly appreciated. Carlo

Do I have the right impression that it's currently not possible to produce confidence intervals for the nls predictions using R? I had a course were we used SAS PROC nlin and there you could get intervals for the parameters and the prediction but I do not have access to SAS. Would it be difficult to implement, I tried to dig into the help pages of nls, vcov and nlsModel but I could not

I would like to display output to a tcltk window from e.g. a call to summary(). I tried to get something else than oneliners into a text window of the kind found at: http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/TextWindows.html But without success. Henrik ------------------------------------------------------------- Henrik Andersson Netherlands Institute of Ecology - Centre for

If I want to fiddle with the axes, I normally plot things like this for instance plot(1:5,axes=F) axis(2,labels=letters[1:5]) BUT, if I want to do remove the axes from a scatter plot matrix using i.e. pairs(cbind(runif(10),runif(10)),axes=F) I get Error in plot.default(x[, j], x[, i], xlab = "", ylab = "", axes = FALSE, : formal argument "axes" matched by

I would like to have arrow heads at the end of my axes, since I am plotting variable where the absolute amount is irrelevant, there is not supposed to be numbers on the axes. An imperfect example: plot(rnorm(10),bty='l',xaxt='n',yaxt='n',ylab='',xlab='',type='l') abline(h=0) Like this but without, the xaxis and with arrrowheads More like this in

Hello R-people, I have a function which plots two figures next to each other using if(dev.cur()==1) x11(width=14,height=7) to create a new window and set the size of the window automatically if no device is open. This works fine in interactive mode but I don't know how Sweave gets the parameters for the figures, and right now the default behavior causes this figure to be square and the

I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3) 17.01 18.15 I want 17.0 18.2 Not 17 18.2 Why is the last digit stripped off in the case when it is zero! Is this a "feature" of R or did I miss something? --------------------------------------------- Henrik Andersson Netherlands Institute of

I have used Sweave a lot the latest year, but never really used any long function calls. If I have code which look like this ------------------------------------------------------------- gof <- benthic.flux(ID="Gulf of Finland", meas.conc=conc, bw.conc=bw.conc, time=times,

Here's another way: > dat <- read.table("clipboard", header=TRUE)[,-(1:3)] > dat A B C D E F G H 1 1 1 0 0 0 0 0 0 2 1 0 1 0 0 0 0 0 3 1 0 0 0 1 0 0 1 4 1 0 0 0 0 1 1 0 5 0 1 1 1 0 0 0 0 6 0 0 1 0 1 0 1 1 7 0 0 0 1 1 1 0 1 8 1 0 1 1 1 0 0 1 9 1 0 1 1 0 1 1 0 > apply(dat > 0, 1, function(idx) paste(names(dat)[idx], collapse=" ")) 1 2

I am about to write general functions to read the output of simulations models. These model generate output files with different sections which I want to analyze plot etc. Since this will be used many people at the department I wanted to make sure that will do this in the best way. For instance I want to read a snippets of data from a text that look like this.

Dear R listers, I know I'm breaking the rules by asking a "homework" related question-- I hope you'll forgive me. I am a social psychology graduate student, and the only one in my department who uses R. I successfully completed my multiple regression and structural equation modeling courses using R (John Fox's car and sem packages were a big help, as was his book).

I'm hoping that this is a relatively easy question for someone familiar with the lme4 package. I'm accustomed to using HLM software and writing a simple 2 level [null] equation like this: L1 - Yij = b0 + e L2 - b0 = B00 + u0 The following command in R provides results that are identical to the HLM program. results <- lmer( Y ~ 1 |id , PanelData4) I can't

R, I want to get means and confidence limits on the original scale for the treatment effect after running a mixed model. The data are: response<-c(16,4,5,8,41,45,10,15,11,3,1,64,41,23,18,16,10,22,2,3)

Hi everyone, I have posted a similar question to this list, but I don't get a reply. I really want to solve this problem, so I post it again... I am trying to use R to fit some mixed-effects models for a nested data. The data is a simulated data with 111 subjects. Each subject has 6 waves' data. Below are the first two subjects' data : > simu1[1:12,] Grouped Data: gf ~ age |