similar to: pairs and panel.smooth for two groups

Hi, First of all, thanks for the efforts of all the developers and contributors - I'm very new to R and at the moment am just using it to create some graphics, but it seems to be quite powerful. I've googled the archives and wasn't able to find anyhting that dealt with this problem, so would appreciate any suggestions/tips. In a pairs plot I'd like to plot a linear regression line

Hello, I'm trying to pull off a certain graph using splom, and can't quite get my panel functions right. Basically, the equivalent using pairs would be something like this (using iris data set as an example): panel.corval <- function(x, y, digits=2, prefix="", cex.cor,col,pch) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1))

Hello, Below is some example code that should reproduce an error I'm encountering while trying to create a tiff plot with two subplots. If I run just the following bit of code through the R GUI the result is what I'd like to have appear in the saved tiff image: x<-seq(0:20) y<-c(1,1,2,2,3,4,5,4,3,6,7,1,1,2,2,3,4,5,4,3,6) plot(x,y,type="l",las=1,ylim=c(0,12))

Hello, I have a customized pairs () fonction as follows that displays correctely my data. ------------------------------------------------------------------------ panel.cor1 <- function (x, y, digits=2, prefix="") { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- cor(x, y,use="pairwise.complete.obs",

Hello list, I am looking to create a figure for my dataset using splom, where there is a splom subplot for each level of a factor within, for example, a 2x2 layout. For each subplot, I wish to put the r-value between each variable pair in the lower panel. The code I have thus far is below, using the iris dataset so that it is reproducible. This almost does what I want, but for some reason

Paul, I was trying to make a minimal self contained example, but I guess I went too far on the minimizing. The original problem came from code more like: library(TeachingDemos) hist(rexp(1000), main='') abline( v=1, col='red') sp.par <- subplot(hist(rnorm(100), main=''), x='topright') op <- par(sp.par[c('usr', 'plt')]) abline(v=0,

Hi, I am trying to use function pairs to plot the scatterplot, but I only want to keep the upper triangle, what's the argument to do this? Thanks, -Jack [[alternative HTML version deleted]]

Hi, I making a figure with six sub-plots using par(mfcol=c(2,3)). In the last sub-plot I want to have two graphs instead of one. I have tried using par(fig=x,y,z,v) but this par seems to overwrite the first par. Is there a simple solution? Thanks! Anna -- View this message in context: http://r.789695.n4.nabble.com/Multiple-plots-in-one-subplot-tp4203525p4203525.html Sent from the R help

Hello, For an experiment, I selected plots of land within a forest either with honeysuckle or without honeysuckle. Thus, my main factor is fixed, with 2 levels: "honeysuckle present"(n=11) and "honeysuckle absent"(n=8). Within each plot of land, I have a "trenched" subplot and an "untrenched" subplot. Within each subplot of every plot, I measured soil

Greetings R users, I am an R graphics newbie trying to produce a custom trellis plot using pairs() with R 2.7.2. I have spatial data on which I run a geographically weighted regression (gwr, using the -spgwr- package). I want to check the gwr coefficients for multicollinearity and spatial association, following Wheeler and Tiefelsdorf (2005), and I would like to summarize the results of this

I ran into this issue when trying to modify a subplot (subplot function in the TeachingDemos package), but here is a more minimal example of the issue (I don't know that it is serious enough to call a bug): This code works how I expect: dev.new() hist(rexp(100)) par(plt=c(0.5,0.9,0.5,0.77), usr=c(0,1,0,1)) box() # show new plt region points(c(0,1), c(0,1), pch=16, col='red', cex=3)

In lattice #toy data library(ggplot2) library(lattice) x <- rnorm(100) y <- rnorm(100) k <- sample(c("Weak","Strong"),100,replace=T) j <- sample(c("Tall","Short"),100,replace=T) w <- data.frame(x,y,j,k) xyplot(y~x|j+k,scales=list(y=list(relation="free"))) will give you a scale in each subplot with a range equal to the range of y

Hi, I'am trying to make a multiple bar plot over a map and I'm having difficulties with the distance between axes labels and the axis. Trying to control this with mgp does not help because it controls both axes simultaneously. For example, with default values (mgp = c(3, 1, 0)) y-axis labels are ok, but x-axis labels are not. Setting mgp = c(3, 0, 0) gives good x-axis labels but the

Hello everyone, I hope this is a fair enough question, but I don’t have access to a copy of Bates and Pinheiro. It is probably quite obvious but the answer might be of general interest. If I fit a fixed effect with an added quadratic term and then do it as an orthogonal polynomial using maximum likelihood I get the expected result- they have the same logLik.

I'm new to R - and trying to create a plot similiar to the spider plot at http://www.statgraphics.com/eda.htm#radar . I can't figure out several things... most of which I would think would be straightforward.... How can I change the lines for each series plotted instead of creating a filled area? How can I get the labels for each of the radial axes at the outside of the plot

Hi, is it possible to reverse the order in which panel.lmline() or panel.smooth() operation in xyplot() ? This type of situation might occur when plotting some variable with depth, but the relation we want to describe is variable ~ depth, and not depth ~ variable, as the plotting formula would suggest. # an example: d <- 1:100 v <- d * rnorm(100) xyplot(d ~ v, ylim=c(100,0),

Hello, can anyome tell me how to access the full script of the panel.smooth function so that I can change the thickness of the smoothing line or its colour ? All I could access is : > panel.smooth function (x, y, col = par("col"), bg = NA, pch = par("pch"), cex = 1, col.smooth = "red", span = 2/3, iter = 3, ...) { points(x, y, pch = pch, col = col, bg =

Hello all, I would like to recreate the plot shown here (from a useR 2009 presentation): http://www.agrocampus-ouest.fr/math/useR-2009/abstracts/pdf/Hocking.pdf I downloaded the code for that image, and discovered that it relies on external web services, and also having PERL installed on the computer. I believe this could be done "locally" using the "seqLogo"

I'm just learning to work with R, and am having some difficulty controlling the output of my dotplot. The default order of the function dotchart seems to be the largest value on the bottom and the smallest on the top. I wanted it in the opposite order (which is how my data frame is currently sorted), so I used the following approach, which works orderidx <- order(b2000$RR)

Hello, I have come across some unexpected behaviour of the function heatmap in the stats package. This looks like a bug to me, but I might have misunderstood something. When calling the function in symmetric mode, the ColSideColors are plotted correctly, but RowSideColors appear in reverse order. This code (modified from the example on the help page) demonstrates the problem: cU <-