similar to: Adding a Rotated Density Plot to an Existing Plot

hi R-users: I want to draw three plot into one figure by layout and the script has been shown below. But I find R does not show the xlab and ylab information completely as shown the figure attached. How can I midify the script.? thank you . xxlab<-paste(cpmd," (",ro,"%)",sep=" ") yylab<-paste(rfmd," (",co,"%)",sep=" ")

The code I'm using is shown below. I would like to have a larger median at the top of the plot so that I can show the entirity of "title_text". Several times I tried messing with "par(mar", but that seemed to make matters worse. By any chance can anyone provide any insight as to the best way to increase the top/upper/northern margin, so the entirity of the title is

Hi Does anyone know how I might pick out diagonal elements of a matrix using a vector? If I create a matrix a: a <- matrix(c(1:16), 4, byrow=TRUE) and I want to pick out the elements (1,1),(2,2),(3,3), or another arbitrary diagonal (upper or lower), is there any way I can use a vector to do this? So if I want a diagonal of size 3, I could create a vector like x <- c(0:2) and then pick

Hello, I'm doing some graphics for a paper and a need customize such with filled region above the density curve. My attempts I get something very near what I need, but I don't solve the problem of use semitransparent filled. Below a minimal reproducible code. Someone has any idea? require(lattice) # toy data... dt <- expand.grid(A=1:2, B=1:3, y=1:50) dt$y <- rnorm(nrow(dt), dt$B,

Dear R users, Using the library(mgcv) and running R under MacOSX, I have fitted a generalised additive model with binomial errors in order to check the linearity of two continuous variables ap2mm and diffdaysm in a glm: > mymodel.gam <- gam(diedhos~ s(ap2mm) + Dweekm + s(diffdaysm) + Dweekm:diffdaysm + ap2mm:Dweekm, binomial) I would like postscript gam plots for the two smoothed

Hi all If I have a string, say "ABCDA", and I want to convert this to the sum of the letter values, e.g. A -> 1 B -> 2 etc, so "ABCDA" = 1+2+3+4+1 = 11 Is there an elegant way to do this? Trying something like which(LETTERS %in% unlist(strsplit("ABCDA", ""))) is not quite correct, as it does not count repeated characters. I guess what I need is

Hi all Does anyone know if it is possible when plotting a line or scatter plot, to selectively color the data points based on the data value? i.e. if plotting say the percentage change in stock price movements, to color +ve points in green and -ve points in red? And extending this to a user-defined range of colors based on the quartile of the data points? Thanks Rory Rory Winston RBS Global

Hi (If you're wondering, this is a Project Euler question :)) If I wanted to calculate the sum of the digits in the decimal representation of 2^1000, what would be a good way to go about that? I've tried the following methods: # Calculate the sum of digits in the decimal representation of 2^n # Only works for smaller values of n bsum <- function(n) { s <- 0 e <-

Hi all (I'm sure this question has been asked before, but I cant find it). If I have two character vectors: > x <- c("aaa","bbb","ccc") > y <- c("1","2","3") How can I get the cartesian product of the string values? > expand.grid(x,y) Gives me a data frame with separate columns...however, I cant seem to get *apply

Hi Firstly , happy Christmas to R-Help! Secondly, I wonder if anyone can help me with the following query: I am trying to reproduce some explicit probability calculations performed in APPL (a Maple extension for computational probability). For instance, in APPL, to compute the probability that the sum of 10 iid uniform variables [0,1] will be between 4 and 6, (i..e Pr( 4 < \sum_{i=1}^{10}X_i

Hi all On page 39 of this paper [1] by Andrew Lo there is a very interesting correlation network diagram (sorry I dont have a more convenient link to the type of diagram I'm talking about). does anyone know of any package in R that can generate these types of diagrams? Cheers -- Rory [1] http://web.mit.edu/alo/www/Papers/august07.pdf [[alternative HTML version deleted]]

Hi R, Can we download Reuters (3000 Xtra) data from R? Does ODBC package help me in this? Or otherwise, is there a way to extract daily closing prices data of Reuters from R? Thank you very much, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}}

Hi I am having a bit of difficulty with changing the canvas size on a trellis/lattice plot. I am plotting two "cubes" of 3-dimensional random numbers, as follows: library(gsl) library(lattice) q <- qrng_alloc(type="sobol", 3) npoints <- 200 rs <- qrng_get(q,npoints) # Plot the normal variates in a 3-dim cube p1 <- cloud(rnorm(npoints) ~ rnorm(npoints) +

Hi all I have a question about correct usage of persp(). I have a simple neural net-based XOR example, as follows: library(nnet) xor.data <- data.frame(cbind(expand.grid(c(0,1),c(0,1)), c(0,1,1,0))) names(xor.data) <- c("x","y","o") xor.nn <- nnet(o ~ x + y, data=xor.data, linout=FALSE, size=1) # Create an (x.y) surface and predict over all points d <-

Hi I am currently (for pedagogical purposes) writing a simple numerical analysis library in R. I have come unstuck when writing a simple Newton-Raphson implementation, that looks like this: f <- function(x) { 2*cos(x)^2 + 3*sin(x) + 0.5 } root <- newton(f, tol=0.0001, N=20, a=1) My issue is calculating the symbolic derivative of f() inside the newton() function. I cant seem to get R to

Hi I have a dataset that looks like this (dput'd below): > head(x, 20) time status 1 2009-07-02 10:32:37 1 2 2009-07-02 10:32:43 0 3 2009-07-02 10:32:43 1 4 2009-07-02 10:32:44 0 5 2009-07-02 10:32:44 1 6 2009-07-02 10:32:48 0 7 2009-07-02 10:32:48 1 8 2009-07-02 10:32:54 0 9 2009-07-02 10:33:04 1 10 2009-07-02 10:33:04 0 11 2009-07-02 10:33:05 1 12 2009-07-02 10:33:07 0 13 2009-07-02

Hi all Here is a rather naive implementation of the SEPA algorithm for generating lexical permutations: lexperm3 <- function(x, n=length(x)) { perms <- list() k <- 1 perms[[k]] <- x k <- k + 1 for (y in 1:(factorial(n)-1)) { i <- n-1 while (x[i] > x[i+1] && i > 0) { i <- i - 1 } # i is largest index st x[i] > x[i+1] j <- n #

Hallo, I would like to add the marginal distributions along the X and the Y axis to a scatter plot. Can anybody help me, please? Thank you, Paolo -- Paolo Bulla Istituto di Metodi Quantitativi Universit?? "L. Bocconi" viale Isonzo 25 20136 Milano paolo.bulla at unibocconi.it

Sender: r-help-bounces at r-project.org On-Behalf-Of: comtech.usa at gmail.com Subject: Re: [R] [R-SIG-Finance] forecasting earnings, sales and gross margin of a company... Message-Id: <b1f16d9d0810231239k506d582i7ecb908b84bc1642 at mail.gmail.com> Recipient: ngottlieb at marinercapital.com -------------------------------------------------------- This information is being sent at the

Hi I have two irregular time series, which are of different lengths and being and end at different times. For the common subset of time that they both span, they should have the same values, but the values may occur at slightly different time intervals. I am trying to "line up" the identical values and reconcile them. I have merged the two series into a zoo object which looks