Displaying 20 results from an estimated 9000 matches similar to: "shade between two curves"
2006 Sep 29
2
GLM information matrix
Is there a function that provides the Fisher information matrix for a
generalized linear model? I do not see how to access the off-diagonal
matrix elements of the value returned by glm. (I'm particularly
interested in logistic regression.)
If not, what is a good way to use R to compute Hessians or other partial
derivatives of log likelihoods?
I would appreciate any guidance.
David
2007 Aug 05
1
Selectively shading areas under two density curves
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 =
2009 Aug 13
3
shading between two smoothed curves
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
2007 Oct 03
2
Shading area under density curves
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){
+
2004 Sep 03
2
debugging an S4 method
Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify.
Thanks,
David
_____________________________
David Bickel http://davidbickel.com
Research Scientist
Pioneer Hi-Bred International
Bioinformatics & Exploratory Research
7250 NW 62nd Ave., PO Box 552
Johnston,
2010 Jul 14
3
Mathematica and R
What are some effective ways to leverage the strengths of R and
Mathematica for the analysis of a single data set?
More specifically, are there any functions that can assist with any of
the following?
1. Calling an R function from Mathematica.
2. Calling a Mathematica function from R.
3. Using XML or another reliable data format to pass vectors, matrices,
and/or lists from one environment to
2005 Feb 22
2
estimate the parameter of exponential distribution, etc.
Given a numeric vector of observations, does R have any generic way to estimate the parameters of commonly used distributions (exponential, gamma, etc.) without numerically optimizing the likelihood function?
Thanks,
David
_______________________________________
David R. Bickel http://davidbickel.com
Research Scientist
Pioneer Hi-Bred International
Bioinformatics & Exploratory Research
7250
2004 Aug 31
2
enter browser on error
Is there a way I can get R to automatically enter the browser inside a user-defined function on the generation of an error? Specifically, I'm trying to debug this:
Error in as.double.default(sapply(lis, FUN)) :
(list) object cannot be coerced to double
In addition: There were 38 warnings (use warnings() to see them)
> traceback()
8: as.double.default(sapply(lis, FUN))
7:
2005 Dec 23
2
convolution of the double exponential distribution
Is there any R function that computes the convolution of the double
exponential distribution?
If not, is there a good way to integrate ((q+x)^n)*exp(-2x) over x from
0 to Inf for any value of q and for any positive integer n? I need to
perform the integration within a function with q and n as arguments. The
function integrate() is giving me this message:
"evaluation of function gave a
2004 May 05
1
Installation on Windows XP
R never gave me any serious problems until I tried to use it on Windows XP. I had 1.8.1 installed and working successfully for a few weeks. This morning when I launched it, I received this error message: "R for Windows GUI front-end has encountered a problem and needs to close. We are sorry for the inconvenience."
I downloaded 1.9.0 and installed it, but when I launch it, I get this
2008 Oct 24
2
Shading underneath a line plot.
Hi,
I'm sure this should be fairly simple to do but I haven't found an example.
I have 10 data points which are plotted as a line, I want to shade under
this line with a colour. Other examples have shading under curves where
there are two sets of data use polygons.
mydata <- c(268,251,254,250,244,246,247,243,241,243)
plot(mydata, type="o")
2003 Jul 22
3
curves with shaded areas?
Hi,
I want to make a plot with abline where the area below or above the
curve is shaded. I can't find any documentation on that. Can anybody
help me with that?
thanks, Juffi
2004 Jun 14
4
interrupt in Linux
Does anyone know how to interrupt R in RedHat? Neither control-c nor Esc is working. What I don't want to do is close the window or kill the entire R process.
Thanks,
David
This communication is for use by the intended recipient and ...{{dropped}}
2007 Oct 29
1
shading of curves with polygon
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)
2010 Mar 09
3
Shade area under curve
I want to shade the area under the curve of the standard normal density.
Specifically color to the left of -2 and on. How might i go about doing
this?
Thanks
--
View this message in context: http://n4.nabble.com/Shade-area-under-curve-tp1586439p1586439.html
Sent from the R help mailing list archive at Nabble.com.
2002 Jan 30
2
Shade area under curve?
Hi all,
I've got this graphics question which really should be easy. I want to shade
an area between bounds under a curve. A suitable beginning seems to be the
following:
> plot(dnorm,-4,4)
> segments(-4,0,4,0)
> segments(-2,0,-2,dnorm(-2))
> segments(2,0,2,dnorm(2))
It is the area between -2 and 2 which I want to shade (or something
similar). Hints anyone?
Robert
2009 Dec 15
1
shading between to curves that intersect repeatedly
Dear R-Community,
I have the following problem and was wondering whether someone could help
me out. I have a number of time series of discharge (up to 3000 data
points in some cases) which I would like to compare on the same line plot
by shading the area between some of them. Can anyone tell me how to do
that? I have tried to fiddle around with the polygon function, but the
curves interact a
2012 Aug 15
1
shade overlapping portions of circles (or other shapes)
I'm making some illustrations and it would be convenient to
automatically shade the overlapping portions of circles. These
illustrations are for Social Choice theory, a field in political
science and economics. I've wrestled together some examples so you
can see what I mean, but have not mastered the "color overlapping
sections" problem (as you will see):
2004 Aug 02
2
random seed puzzle
After reading the help page on set.seed, I am unsure about how to solve the following problem. I need to call function f a thousand times. The list of values returned by f, should be as random as possible. f calls g twice:
f <- function(){g1 <- g(1); g2 <- g(2); c(g1; g2)}
The function g in turn calls sample and returns a number, but also depends on its argument, so, starting from the
2002 May 19
3
How to shade part of a density plot
I'm trying to shade part of a density plot. The code I'm trying (using
the Old Faithful data as an example) is something like this:
# The Old Faithful geyser data
data(faithful)
d <- density(faithful$eruptions, bw = "sj")
plot(d)
polygon(d[d$x>4], col = "wheat")
I expected that the part of the curve to the right of 4 on the x axis
should be shaded, but nothing