search for: djf

Hello R-Users I am struggling with this line plot, it might be simple but I am missing something here. First of all I want to make multiple line plots across seasons (DJF,MAM,JJA,SON) for 12 variables (here, called nodes) and fill them with the node. So that season=x-axis, node=line col and freq=y-axis. My plot currently links all DJFs, MAMs, JJAs and SONs, however I will like them to be linked across seasons, so that for node 1, I should have a line plot of DJF,...

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 <- strptime(paste("28/02", year, sep="/"), "%d/%m/%Y") ## Gives wrong answer: temp.seq <- seq(temp.begin, tem...

...: I have a zooreg() object that starts in 1997 and ends in 2001. This object contains daily data for the 4 years aa<-zooreg(1:1825,start=as.Date("1997-01-01")) My aim is to subset the data according to seasons (Southern Hemisphere) for continuous years: December - January - February (DJF-Summer: 1997-2001), March - April - May (MAM-Autumn: 1997-2001), June - July - August (JJA-Winter: 1997-2001), September - October - November (SON-Spring: 1997-2001) thereby analysing the seasons data only for all years. The example below is only for DJF but I would like to replicate the analysis f...

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 variables:&quot...

...s and Biotechnology Blichers Allé 20, P.O. BOX 50 DK-8830 Tjele Tel: +45 8999 1900 Direct: +45 8999 1900 Mobile: +45 E-mail: Joao.Fadista@agrsci.dk <mailto:Joao.Fadista@agrsci.dk> Web: www.agrsci.dk <http://www.agrsci.dk/> ________________________________ Tilmeld dig DJF's nyhedsbrev / Subscribe Faculty of Agricultural Sciences Newsletter <http://www.agrsci.dk/user/register?lan=dan-DK> . Denne email kan indeholde fortrolig information. Enhver brug eller offentliggørelse af denne email uden skriftlig tilladelse fra DJF er ikke tilladt. Hvis De ikke er de...

...rt of my code: minlon<-4.5 maxlon<-12 minlat<-46.5 maxlat<-52.5 gt=GridTopology(cellcentre.offset=c(5.75,47.5),cellsize=c(2.5,2),cells.dim=c(3,3)) grd=SpatialGrid(gt,proj4string = CRS(as.character(NA))) gridparameters(grd) gd<-as.data.frame(grd,"SpatialGrid") h<-season_djf h<-as.data.frame(h) djf.att=SpatialPointsDataFrame(gd,h) gridded(djf.att)=TRUE spplot(djf.att, col.regions = colorRampPalette(c("red","orange", "yellow", "lightblue", "blue", "purple")), xlab="Longitude(?)", ylab="Lati...

...s and Biotechnology Blichers Allé 20, P.O. BOX 50 DK-8830 Tjele Tel: +45 8999 1900 Direct: +45 8999 1900 Mobile: +45 E-mail: Joao.Fadista@agrsci.dk <mailto:Joao.Fadista@agrsci.dk> Web: www.agrsci.dk <http://www.agrsci.dk/> ________________________________ Tilmeld dig DJF's nyhedsbrev / Subscribe Faculty of Agricultural Sciences Newsletter <http://www.agrsci.dk/user/register?lan=dan-DK> . Denne email kan indeholde fortrolig information. Enhver brug eller offentliggørelse af denne email uden skriftlig tilladelse fra DJF er ikke tilladt. Hvis De ikke er de...

...sed seq() to create the increment seq <- paste(seq(1914, 2006, by=1), "*.txt.h", sep=".") But it didn't work. Below is the script; seq <- paste(seq(1914, 2006, by=1), "*.txt", sep=".") for (i in 1:93){ list <- list.files("~/ukcp09/txt/x.djf", seq) file <- lapply(list, read.table) mean <- (Reduce("+", file))/3 } When I call > objects() > list [1] "2006.01.txt.h" "2006.02.txt.h" "2006.12.txt.h" The listed files are the last files in the sequence. Am I doing the loop correct...

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

...ences Dept. of Genetics and Biotechnology Blichers Allé 20, P.O. BOX 50 DK-8830 Tjele Phone: +45 8999 1900 Direct: +45 8999 1900 E-mail: Joao.Fadista@agrsci.dk <mailto:Joao.Fadista@agrsci.dk> Web: www.agrsci.org <http://www.agrsci.org/> ________________________________ DJF now offers new degree programmes <http://www.agrsci.org/content/view/full/34133> . News and news media <http://www.agrsci.org/navigation/nyheder_og_presse> . This email may contain information that is confidential. Any use or publication of this email without written permission from...

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

...lices (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 <- c("DJF","MAM","JJA","SON") for (i in runs) { for (j in years) { for (k in seasons) { data <- open.ncdf(paste("/data09/cstocker/data/atm/tr_1500_", i,"/PRECT/PREC_tr1500_", i,"_", j,"_", k,".nc", sep="")) p...

...#39; 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 <- lapply(list, read.table) # Calculations mean <- (Reduce("+", file))/c d2[[i]] <- file[[i]] - mean Apparently, the following command is the source of the error. d2[[i]] <- file[[i]] - mean It works OK when I typed into the terminal the following aft...

...ends in 2001. This object contains daily data for the 4 years >>> >>> aa<-zooreg(1:1825,start=as.Date("1997-01-01")) >>> >>> My aim is to subset the data according to seasons (Southern Hemisphere) for continuous years: December - January - February (DJF-Summer: 1997-2001), March - April - May (MAM-Autumn: 1997-2001), June - July - August (JJA-Winter: 1997-2001), September - October - November (SON-Spring: 1997-2001) thereby analysing the seasons data only for all years. The example below is only for DJF but I would like to replicate the analysis f...

Hello there folks, I have been using the latest alpha version, I believe it is 5 , which proves to be a bit better than the previous work, but I do face some serious problems. This is my first time installing samba on a SunOS 4.1.4 and I do have some concern about samba works. Assuming the existance of 2 Unix systems that are exporting filesystems between them, and 1 PC that requests home

...46 +0100 From: S?ren H?jsgaard <Soren.Hojsgaard at agrsci.dk> Subject: Re: [R] Summary information by groups programming assitance To: "Ranney, Steven" <steven.ranney at montana.edu>, <r-help at r-project.org> Message-ID: <C83C5E3DEEE97E498B74729A33F6EAEC05630089 at DJFPOST01.djf.agrsci.dk> Content-Type: text/plain; charset="iso-8859-1" Maybe summaryBy (or lapplyBy/splitBy) in the doBy package might help you. Regards S?ren ________________________________ Fra: r-help-bounces at r-project.org p? vegne af Ranney, Steven Sendt: ma 22-12-2008 22:51 Til...