similar to: shading between to curves that intersect repeatedly

Hello, I cannot figure out how to shade between two lines in a plot. For example, if I am trying to plot a confidence envelope and I would like to shade the interior of the envelope grey. Is this possible in R? Thanks in advance. Heather Lynch

Hello all, could anyone explain the exact meaning of parameter nproc? Why different values of nproc give so different critical values, i.e. meanL2BB$computeCritval(0.05,nproc=3) [1] 0.9984853 meanL2BB$computeCritval(0.05,nproc=1) [1] 0.4594827 The strucchange-package description gives "integer specifying for which number of processes Brownian motions should be simulated" - do I need

To R programming experts, I am a undergraduate student, and now doing research personally. I apply diagonal bivariate poisson (R package "bivpois") with stochatics weighted function (refer to dixoncoles97 section 4.5 to 4.7). However I dont know how to fit this stochatical weighted function to the completed bivariate poisson model. I know that some other references for dynamic soccer

Dear all, I know that my question is somewhat special but I tried several times to solve the problems on my own but I am unfortunately not able to compute the following test statistic using the quantreg package. Well, here we go, I appreciate every little comment or help as I really do not know how to tell R what I want it to do^^ My situation is as follows: I have a data set containing a

I am trying to learn how to use the graphics from the lattice package ( and am very new to R). I am trying to replicate the example plot referenced below, by using the lattice xyplot & lpolygon to create panels. I get what appears to be the correct shape of the filled region, but cannot get the position to overlay properly. I have attempted with various settings of position. ( i.e.

Hello again With the poylgon function it's possible to plot shaded areas under a curve. But somehow it connects the start and the endpoint of a line and fills whats between them. I would actually like to set the boarders of the shading by two min and max curves, but I failed. i tried stl like: polygon(min, max, border=0), where min and max are my curves plotet by plot(max) lines(min)

I would like to plot multiple random walks onto the same graph. My p variable dictates how may random walks there will be. par(mfrow=c(1,1)) p <- 100 N <- 1000 S0 <- 10 mu <- 0.03 sigma <- 0.2 nu <- mu-sigma^2/2 x <- matrix(rep(0,(N+1)*p),nrow=(N+1)) y <- matrix(rep(0,(N+1)*p),nrow=(N+1)) t<- (c(0:N))/N for (j in 1:p) { z <- rnorm(N,0,1) x[1,j] <- 0 y[1,j]

Dear Listers, I am drawing a plot of two density curves, for male and female incomes. I would like to shade/hatch/color (whatever) the areas under the curves which are distinctive for each gender. This is the code I have tried so far: m <- density(topmal.d$y, bw = "sj") f <- density(topfem.d$y, bw = "sj") par(mfrow = c(1,1)) plot(x = c(0,400), y = c(0,0.02), type =

Hello, I have a question regarding shading regions under curves to display 95% confidence intervals. I generated bootstrap results for the slope and intercept of a simple linear regression model using the following code (borrowed from JJ Faraway 2005): > attach(allposs.nine.d) > x<-model.matrix(~log(d.dist,10))[,-1] > bcoef<-matrix(0,1000,2) > for(i in 1:1000){ +

Hi, I have been using maxLik to do some MLE of Geometric Brownian Motion Process and everything has been going fine, but know I have tried to do it with jumps. I have create a vector of jumps and then added this into my log-likelihood equation, know I am getting a message: NA in the initial gradient My codes is hear # n<-length(combinedlr) j<-c(1,2,3,4,5,6,7,8,9,10)

I have a set of 52 weekly values, one is the desired high lake stage for a week, and the other is the desired low lake stage for each week. It looks like this: week High Low 1 16 14.5 2 16 14.5 3 15.95 14.45 4 15.84 14.34 5 15.73 14.23 6 15.61 14.11 7 15.5 14 8 15.38 13.88 9 15.25 13.75 10 15.13 13.63 .... 52 15 13.5 I can plot out two nice curves

Is there a package that will allow me to fit Brownian motion and Ornstein?Uhlenbeck models of evolution for discrete traits? I know that geiger and ouch have commands for fitting these models for continuous traits, but these aren't suitable for discrete trait evolution, correct? -- View this message in context:

Hi, we are running Samba 2.2.2 on Solaris8; one client (MS 2000) is not able to map a drive: The tcpdump shows the correct user name <user>@<dom>. smbd -d 3 shows <user>_<dom> and thus the authentication on the PDC fails. Does anybody know why Samba translates the @ into an _? Thanks a lot, Roman __________________________________________________________________ The NEW

Dear all, I'm using the curve() function to plot discharge Q against water depth a. However, I would like to have a graph of water depth a plotted against discharge Q. How can this be done? Minimal working example: S0 = 0.004 n = 0.04 tanalpha = 1.4/1.5 par(mar = c(5,5,1,1)) # b, l, t, r

Dear R Helpers, I have searched for any R package or code for simulating multivariate fractional Brownian motion (mFBM) or multivariate fractional Gaussian noise (mFGN) when a covariance matrix are given. Unfortunately, I could not find such a package or code. Can you suggest any solution for multivariate FBM and FGN simulation? Thank you for your help. Best Regards, Ryan ----- Wonsang You

Forum, I have two time series that I plot up using the data/code below. The x-axis shows weekly time step starting on 6/3/1999 (SP=10) and ending on 10/25/2007 (SP=448). I'm trying to shade between the two lines in a similar fashion as the graph found here: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=122 Can the polygon function be used to do this? There are lots of

Hi: I was wondering whether I can find some help for computer simulation of stochastic processes (e.g. Brownian motion), for pedagogicl/instructional purposes. Any help would be appreciated. thanks, Ravi.

Hi, I'm writing a short course tutorial to Browninan Motion/ Langevin Equation. At the end of the theory section I wanted to add a short GNU R example, so the students can play a little around. I already looked in the MASS book (by Venables and Ripley) but I couldn't find any Brownian Motion/ Langevin Equation package. Are there any good packages or tutorials available which cover R and

i have the following constructed and running very well,, thanks to Gabor Grothendieck for his help. >data.info <- c("station.id", "year", "date", "max.discharge") > > for(i in 1:num.files) { + station.id <- substring(data[i], 1,8) + DF <- read.table(data[i], sep=",", blank.lines.skip = TRUE) + z <- zoo(DF[,4],

note: whole function is below- I am sure I am doing something silly. when I use it like USGS(input="precipitation") it is choking on the precip.1 <- subset(DF, precipitation!="NA") b <- ddply(precip.1$precipitation, .(precip.1$gauge_name), cumsum) DF.precip <- precip.1 DF.precip$precipitation <- b$.data part, but runs fine outside of the function: days=7