similar to: plotting time series for particular months!

Hello all, I am trying to use "clim.pact" package for my work, but since this is the beginning for me to use gridded datasets in "R", I am having some trouble. I want to do seasonal analyses like trends, anomalies, variograms, EOF and probably kriging too to downscale my 1 degree gridded data to 0.5.  So, as a first step, I compiled my entire dataset (with 25

Hello All. I would like to repeat a multi-screen, multivariate time series plot across multiple units. Consider the following data: <r-help@stat.math.ethz.ch> UNIT TRIAL Y1 Y2 Y3 Y4 Y5 Y6 1 109 1 NA NA NA 13718 3655 15240 2 109 2 15000 6000 30000 13057 6198 21962 3 109 3 10000 10000 20000 8356 12929 27992 4 109

Dear useR's, I have a problem to read data (daily precipitation). The data are formated in months by columns wich leads to empty items in months with less then 31 days. What i did was reading the file and try to reshape the object a<- read.table("nds.txt", header=TRUE, sep=";", dec = ".", na.string="NA") ar<-reshape(a,

Dear Users, I have a monthly data for a number of years(1960-2007) for a number of stations and i wish to extract sesonal time-series for the months of March-May and October-November for very station. I have read this data with read.table in R with stations as columns and time (months) as rows. My attempt to aggregate with the zoo package using the function as.yearqtr failed since this sums

Dear all, I'm wanting to do a series of comparisons among 4 categorical variables: a <- aggregate(y, list(var1, var2, var3, var4), sum) This gets me a very nice 2-dimensional data frame with one column per variable, BUT, as help for aggregate says, <<empty subsets are removed>>. I don't see in help(aggregate) how I can change this. In contrast, a <- tapply(y,

Hi All. I'd like to change the default plotting colors used to construct a scatterplot with regression line in the CAR package. So, scatterplot(y~pred,smooth=FALSE, xlab="X", ylab="Y", lwd=2) If I change the palette (e.g., palette(ranbow(6)), I can change the color of the lines and points. However, the axes and labels remain in black (i.e., the first color in the

Greetings, I'm trying to find a function to generate a random landscape consisting of x by y cells, allows a specification of number of classes and their proportional abundance in the landscape. Does R support such a function ? Thanks in advance Steve Steve Friedman, Assistant Professor Center for Global Change and Earth Observations Michigan State University East Lansing, Michigan

Hi All. I need a bit of help extracting the residual error variance from the VarCorr structure from lmer. #Here's a 2-way random effects model lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat) #Get the structure vc.fit <- VarCorr(lmer.1) #results in..... $person 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.7755392 $rater 1 x 1 Matrix

Hello R users- I am new to R, and tried searching the archives and literature for an answer to this - please be patient if I missed something obvious. I am fitting a logistic regression model, and would like to include variance functions (specifically the varIdent function). I cannot figure out how to do this either in glmmPQL (or something similar) for the model with random effects, or in glm

Hi R helpers, I was wondering if anybody knows if is possible to generate bootstrap values for a cluster analysis in R. What I am trying to do is obtain some confidence on the clusters formed by resampling the data set. A similar type of analysis is used in molecular taxonomy and the confidence values of each cluster are placed in the nodes of the dendogram. Any ideas on how to do this

Hi, I was trying to do a multi-level modeling of the data using the nlme package. I encounter the following message: > addRandomSlopeRange<- lme(FAN.range ~ CHN.range + Age + Block.T, data=data.block, random = ~CHN.range|CHI, method = "ML", na.action=na.omit) Error in lme.formula(FAN.range ~ CHN.range + Age + Block.T, data = data.block, : nlminb problem, convergence error code

I understand that Rcmdr needs the tcltk package in order to run. However, when I attempt to using the library (Rcmdr) command, I get the message below. I have attempted to start X in R as well as using starting it from the Utilities folder on my Mac to no avail. I can't figure out if there is an issue with X11 or with R itself (I just updated to 1Is there a terminal command I can use

Hi all, I want to let linux server join ad by using a trust ad's child domain user, but failed with error. below is my env and what I have try I have 3 domain controller: test.com,demo.com and chn.demo.com test.com with demo.com is two way trust. and chn.demo.com is the child domain of demo.com demo at demo.com chn at chn.demo.com can join ad member to test.com I have tested demo at

On Sat, 01 Oct 2016 23:20:01 +0000 Fay zhang via samba <samba at lists.samba.org> wrote: > Hi all, > I want to let linux server join ad by using a trust ad's child domain > user, but failed with error. > below is my env and what I have try > > I have 3 domain controller: test.com,demo.com and chn.demo.com > > test.com with demo.com is two way trust. and

Dear all, I am working with a data file which is the record of precipitation measurement normaly done every 10 minutes. I would like to check if there are missing times in my data file. Is there a function existing able to check for that in R ? Thanks by advance, Jessica

Hello, I've noticed that with g++ 4.3 I get the following "warning: deprecated conversion from string constant to ?char*?" when I call e.g., Rprintf("Hello world"); or when I try to define a const char e.g., char lower = 'L'; I am wrapping the c++ code in extern "C" {}. This seems to be something new with g++ 4.3, because the same code compiled

Hi, I'd like to unsubscribe from the list. Thanks Neyra ________________________________ De: "r-help-request@r-project.org" <r-help-request@r-project.org> Para: r-help@r-project.org Enviado: sáb, octubre 30, 2010 5:30:07 AM Asunto: R-help Digest, Vol 92, Issue 31 Send R-help mailing list submissions to r-help@r-project.org To subscribe or unsubscribe via the

Hello, all; I get an unexpected result when trying to plot a probability histogram with R1.9.1 on windows xp: #with the following code: > x <- runif(100,0,1) > hist(x) > hist(x, freq=F) > h <- hist(x, freq=F) > summary(h) # Length Class Mode #breaks 11 -none- numeric #counts 10 -none- numeric #intensities 10 -none- numeric #density 10

Dear R helpers I have a daily rainfall dataset of 5 stations, but the time-series are of different lengths. Three stations measured precipitation from 1917 to 2004, one from 1930 to 1979, and other from 1975 to 1998. I'd like to know how to analyse if rainfall is similiar at all stations. Is it possible to compare time-series like these? Best Regards!

Btw, all the pbinom() examples posted to S-news today return the proper result in R Martin >>>>> "JS" == James Stapleton <stapleton at stt.msu.edu> writes to S-news : JS> More interesting behavior of pbinom: >> pbinom(3,2,.8) JS> [1] 1 >> pbinom(30,2,.8) JS> [1] 0.9967653 >> pbinom(40,2,.8) JS> [1]