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

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)

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){ +

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

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

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

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.

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]

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)

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 R users, I still feel new to R so please apologize if I am doing something stupid here. My use of the polygon() function produces a result that I cannot comprehend: In a plot, I would like to shade the area below a normal distribution. However, I do not want the entire area to be shaded, but just the area on the right side of a vertical line that I draw through the distribution (in

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

Is there a function that can shade the region between two arbitrary curves? For example, how can I fill in the area between these two plotted curves? > x <- 1:10 > y <- x + rnorm(10) > z <- x + 10 + rnorm(10) > plot(x, y, type = "l", ylim = c(0,20)) > lines(x, z) I would appreciate any help. David

I have 2 datafiles 'target' and 'observed' as shown below (I will gladly email these 2 small files to whomever). X25. And X75. Indicate the value of 25th and 75th-percentile of the target ('what should be') and the observed ('what is'). The i.value is simply the month. > target X i.value X25. X75. 1 one.month 1 10.845225 17.87237 2

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