similar to: Help on multi-line plot

Dear R users, I am currently working in subsetting a zooreg() object using either window or subset. I have a solution but it may be a bit cumbersome when I start working with actual data. Your inputs would be greatly appreciated. Example: I have a zooreg() object that starts in 1997 and ends in 2001. This object contains daily data for the 4 years

Hi there I'm working with ncdf data. I have different dataset for 4 runs, 4 seasons and 3 timeslices (48 datasets in total). The datasets have the following dimensions: 96 longitudes, 48 latitudes and 30 time steps. To read all of them in, I wrote the following loop: runs <- c("03","04","05","06") years <- c(1851,1961,2061) seasons <-

Mathew Brown Institute of Bioclimatology University of G?ttingen B?sgenweg 2 37077 G?ttingen, Germany t: +49 551 39 9359 mathew.brown at forst.uni-goettingen.de On 9/24/2011 6:00 PM, r-help-request at r-project.org wrote: > Send R-help mailing list submissions to > r-help at r-project.org > > To subscribe or unsubscribe via the World Wide Web, visit >

I think this is an interesting discussion -- I've learned from both Steve's and Brian's comments, and I'm broadening it to R-help since I think others will be interested as well. The problem up for comment is: result <- apply(array.3D, 1:2, sum) Where array.3D is 3000 by 300 by 3. The original poster already had a perfectly good replacement for this problem that was

I am trying to extract only the winters (defined to be 01-Dec through 28-Feb) of daily data from 1948-2002. There are 90 days in each winter season. I wrote the following code to gather the winter dates into a single vector: DJF <- NULL for(year in 1949:1999) { temp.begin <- strptime(paste("01/12", year-1, sep="/"), "%d/%m/%Y") temp.end <-

Dear all Im new in R I have a necdf data set that I want to plot : this is my data set [1] "file Tmax.DJF.daily.1981_1999.echama2.nc has 4 dimensions:" [1] "longitude Size: 127" [1] "latitude Size: 110" [1] "ht Size: 1" [1] "t Size: 1680" [1] "------------------------" [1] "file Tmax.DJF.daily.1981_1999.echama2.nc has 1

Dear all, I would like to know how can I concatenate 2 data.frames into a single one. Both data frames have the same number of columns and the same class type in each correspondent column. So what I want is to have a new data.frame where I have first the values from one data.frame and then the values from a second data.frame would came after in this new data.frame. Thanks in advance. Med

When I try to convert the zoo object to a timeSeries object, which would allow me to utilize Rmetrics packages, I get an error message. > Data<-read.zoo("c:\\DOWUBSPRICING.txt,na.strings="NA",sep="\t",header=T) > is(Data) "zoo" > as.timeSeries.zoo(Data) Error in .local (.Object, . ) Is this happening because I am using daily data?

Hi all, I've got a list of files from 1914 to 2000. For each file, I can call the read.table function as follows. > file <- read.table("1914.txt") But I want to make a script so that at the end of the loop, the script selects the file 1915.txt, 1916.txt and so on. I've used seq() to create the increment seq <- paste(seq(1914, 2006, by=1), "*.txt.h",

Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp, zlim=c(1000,1020),color.palette = colorRampPalette(c("blue", "white",

Hi List, I have (otherwise successfull) migrated a server from cyrus to dovecot (Mac OS X Server 10.5.8 to 10.6.4 [dovecotd version to be 1.1.20apple0.5]). There are ca. 300+ users of which a very small handfull (ca. 1 %) are having reported a usage size that is very different to what 'du' reports. How can I rebuild the maildirsize and why might it be wrong in the first place(?) Will the

Dear all, Can I do a permutation test if the number of individuals in one group is much bigger than in the other group? I searched the literature but I didin´t find any assumption that refers to this subject for permutation tests. Best regards João Fadista Ph.d. student UNIVERSITY OF AARHUS Faculty of Agricultural Sciences Dept. of Genetics and Biotechnology Blichers Allé 20, P.O.

Hi all, How can I overcome the error "object of type 'closure' not subsettable" I ran the following script seq <- paste(seq(1914, 1916, by=1), "*.y", sep=".") # make sequence c <- 3 # total number of files d2 <- file # creates dummy file # Input sequence in loop for (i in 1:3){ list <- list.files("~/ukcp09/txt/x.djf", seq[[i]]) file

Hi all, I have made a plot with ssplot, using a SpatialPointsDataFrame. The content is quite simple, as I have 9 grid points with lon/lat coordinates and 9 values attached to these coordinates. They are in a square area of 3 by 3 gridboxes. I would like to lay a map from maps() over these values, but when I try this, the grids of the maps (smaller) do not overlap with the grids I have

Dear all, I would like to know if I can compare by a significance test 2 models with different kind of parameters. Perhaps I am wrong but I think that we can only compare 2 models if one is a sub model of the other. Med venlig hilsen / Regards João Fadista Ph.d. studerende / Ph.d. student AARHUS UNIVERSITET / UNIVERSITY OF AARHUS Det Jordbrugsvidenskabelige Fakultet / Faculty of

If I have a contourplot (in the lattice package) and I want to add straight lines to it, how do I do this? I see that there are llines() and lsegement() functions for lattice plots, but they don't seem to do anything in this case: library(lattice) library(KernSmooth) x=rnorm(10000) y=x+rnorm(x,0,.5) a=bkde2D(cbind(x,y),.7) z=as.vector(a$fhat) grid=expand.grid(x=a$x1,y=a$x2) grid$z=z

Good day everyone, Can anyone suggest an effective method to solve the following problem: I have 2 dataframes D1 and D2 as follows: D1: dates ws wc pwc 2005-10-19:12:00 10.8 80 81 2005-10-20:12:00 12.3 5 15 2005-10-21:15:00 12.3 3 15 2005-10-22:15:00 11.3 13 95 2005-10-23:12:00 12.3 13 2 2005-10-24:15:00 10.3 2 95 2005-10-25:15:00 10.3 2 2 D2:

Dear List-Members, since 3 weeks I have been heavily working on reproducing the results of an economic paper. The method there uses the numerical solution of an integral within nonlinear least squares. Within the integrand there is also some parameter to estimate. Is that in the end possible to implement in R [Originally it was done in GAUSS]? I'm nearly into giving up. I constucted an

Hello, I am estimating the following model: so2.lme<-lme(so2~1+I(alcadakm^2)+dia,data=subjectes2,na.action=na.omit) And when I try to plot the random effects of the quadratic term with respect to a covariate (mam) I get an error: > so2.lmeRE<-ranef(so2.lme,augFrame=T) > plot(so2.lmeRE,form=I(alcadakm^2)~mam) Error in plot.ranef.lme(so2.lmeRE, form = I(alcadakm^2) ~ mam ) : Only

Is there a form of with (or an equivalent function) that can be applied to objects? I'd like to be able to write something like (drawing from bioconductor + trellis as an example) xyplot(maA ~ maM | maPrintTip, object = swirl[,1]) which would be interpreted as xyplot(swirl[,1]@maA ~ swirl[,1]@maM | swirl[,1]@maPrintTip) (or even better, match functions not slots eg. as